Genome-wide local ancestry and direct evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus) The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of a human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here we assess whether ancestry at mitochondrially-targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially-targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

GRAPHICAL SUMMARY

O_FIG O_LINKSMALLFIG WIDTH=170 HEIGHT=200 SRC="FIGDIR/small/457829v2_ufig1.gif" ALT="Figure 1">
View larger version (33K):
org.highwire.dtl.DTLVardef@1031fb0org.highwire.dtl.DTLVardef@1575faorg.highwire.dtl.DTLVardef@f4795corg.highwire.dtl.DTLVardef@4ba74e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Highlights* Using high-density genome-wide SNP data, we present evidence for mitonuclear coevolution in hybrid African cattle.
* We observe a significant increase of taurine ancestry across multiple hybrid populations at mitochondrially-targeted nuclear genes.
* Our results provide support for the hypothesis of mitonuclear incompatibility between the zebu nuclear genome and the taurine mitochondrial genome.
Analysis of the leaf metabolome in Arabidopsis thaliana mutation accumulation lines reveals association of pleiotropy and fitness consequences Understanding the mechanisms by which mutations affect fitness and the distribution of mutational effects are central goals in evolutionary biology. Mutation accumulation (MA) lines have long been an important tool for understanding the effect of new mutations on fitness, phenotypic variation, and mutational parameters. However, there is a clear gap in predicting the effect of specific new mutations to their effects on fitness. Here, we complete gene ontology analysis and metabolomics experiments on Arabidopsis thaliana MA lines to determine how spontaneous mutations directly affect global metabolic output in lines that have measured fitness consequences. For these analyses, we compared three lines with relative fitness consistently higher than the unmutated progenitor and three lines with lower relative fitness as measured in four different field trials. In a gene ontology analysis, we find that the high fitness lines were significantly enriched in mutations in or near genes with transcription regulator activity. We also find that although they do not have an average difference in the number of mutations, low fitness lines have significantly more metabolic subpathways disrupted than high fitness lines. Taken together, these results suggest that the effect of a new mutation on fitness depends less on the specific metabolic pathways disrupted and more on the pleiotropic effects of those mutations, and that organisms can explore a considerable amount of physiological space with only a few mutations.

Significance StatementAs the source of all genetic variation, new mutations are crucial for understanding how organisms adapt to their environment. However, the ways in which new mutations affect the range of metabolic reactions that occur in the cell is unknown. With a combination of gene functional analyses and measurement of the small molecules that drive cellular function and physiology, we find that mutations associated with high fitness are disproportionately found in or near proteins-coding genes that regulate the specific timing and location of gene expression and are less disruptive of cellular physiology.
Alkaline phosphatase activity of serum affects osteogenic differentiation cultures Fetal bovine serum (FBS) is a widely used supplement in cell culture medium, despite its known variability in composition which greatly affects cellular function and consequently the outcome of studies. In bone tissue engineering, the deposited mineralized matrix is one of the main outcome parameters, but using different brands of FBS can result in large variations. Alkaline phosphatase (ALP) is present in FBS. Not only is ALP used to judge the osteogenic differentiation of bone cells, it may affect deposition of mineralized matrix. The present study focused on the enzymatic activity of ALP in FBS of different suppliers and its contribution to mineralization in osteogenic differentiation cultures. It was hypothesized that culturing cells in a medium with high intrinsic ALP activity of FBS will lead to higher mineral deposition compared to media with lower ALP activity. The used FBS types were shown to have significant differences in enzymatic ALP activity. Our results indicate that the ALP activity of the medium not only affected the deposited mineralized matrix but also the osteogenic differentiation of cells as measured by a changed cellular ALP activity of human bone marrow derived mesenchymal stromal cells (hBMSC). In media with low inherent ALP activity, the cellular ALP activity was increased and played the major role in the mineralization process; while, in media with high intrinsic ALP activity contribution from the serum, less cellular ALP activity was measured and the ALP activity of the medium also contributed to mineral formation substantially. Our results highlight the diverse effects of ALP activity intrinsic to FBS on osteogenic differentiation and matrix mineralization and how FBS can determine the experimental outcomes, in particular for studies investigating matrix mineralization. Once again, the need to replace FBS with more controlled and known additives is highlighted.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=55 SRC="FIGDIR/small/462693v2_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@902078org.highwire.dtl.DTLVardef@43b6fborg.highwire.dtl.DTLVardef@1a244fborg.highwire.dtl.DTLVardef@5ef819_HPS_FORMAT_FIGEXP  M_FIG C_FIG A bioactive phlebovirus-like envelope protein in a hookworm endogenous retrovirus Endogenous retroviruses (ERVs), accounting for 14% of our genome, serve as a genetic reservoir from which new genes can emerge. Nematode ERVs are particularly diverse and informative of retrovirus evolution. We identify Atlas virus - a novel, intact ERV in the human hookworm Ancylostoma ceylanicum with an envelope protein genetically related to GN-GC glycoproteins from phleboviruses. A cryo-EM structure of Atlas GC reveals a class II viral membrane fusion protein fold not previously seen in retroviruses. Atlas GC has the structural hallmarks of an active fusogen. Atlas GC trimers insert into membranes with endosomal lipid compositions and low pH. When expressed on the plasma membrane, Atlas GC has cell-cell fusion activity. Atlas virus is transcriptionally active at specific stages of hookworm development. The preserved activities and expression pattern of Atlas GC suggest it has acquired a cellular function. Our work reveals an unexpected degree of structural and genetic plasticity in retroviruses. Reconstructing SNP Allele and Genotype Frequencies from GWAS Summary Statistics The emergence of genomewide association studies (GWAS) has led to the creation of large repositories of human genetic variation, creating enormous opportunities for genetic research and worldwide collaboration. Methods that are based on GWAS summary statistics seek to leverage such records, overcoming barriers that often exist in individual-level data access while also offering significant computational savings. Here, we propose a novel framework that can reconstruct allelic and genotypic counts/frequencies for each SNP from case-control GWAS summary statistics. Our framework is simple and efficient without the need of any complicated underlying assumptions. Illustrating the great potential of this framework we also propose three summary-statistics-based applications implemented in a new software package (ReACt): GWAS meta-analysis (with and without sample overlap), case-case GWAS, and, for the first time, group polygenic risk score (PRS) estimation. We evaluate our methods against the current state-of-the-art on both synthetic data and real genotype data and show high performance in power and error control. Our novel group PRS method based on summary statistics could not be achieved prior to our proposed framework. We demonstrate here the potential applications and advantages of this approach. Our work further highlights the great potential of summary-statistics-based methodologies towards elucidating the genetic background of complex disease and opens up new avenues for research. A theoretical model of mitochondrial ATP synthase deficiencies. The role of mitochondrial carriers. The m.8993T>G mutation of the mitochondrial MT-ATP6 gene is associated with NARP syndrome (Neuropathy, Ataxia and Retinitis Pigmentosa). The equivalent point mutation introduced in yeast Saccharomyces cerevisiae mitochondrial DNA considerably reduced the activity of ATP synthase and of cytochrome-c-oxidase preventing yeast growth on oxidative substrates. The over-expression of the mitochondrial oxodicarboxylate carrier (Odc1p) is able to rescue the growth on oxidative substrate by increasing the substrate-level phosphorylation of ADP coupled to conversion of -ketoglutarate (AKG) into succinate with an increase in Complex IV activity. Previous studies showed that equivalent point mutations in ATP synthase behave similarly and can be rescued by Odc1p overexpression and/or uncoupler of OXPHOS from ATP synthesis. In order to better understand the mechanism of ATP synthase mutation bypass, we developed a core model of mitochondrial metabolism based on AKG as respiratory substrate. We describe the different possible metabolite output and the ATP/O ratio values as a function of ATP synthase inhibition. Secretory cells are the primary source of pIgR in small airways BackgroundLoss of secretory immunoglobulin A (SIgA) is common in COPD small airways and likely contributes to disease progression. We hypothesized loss of SIgA results from reduced expression of pIgR, a chaperone protein needed for SIgA transcytosis, in the COPD small airway epithelium.

MethodspIgR-expressing cells were defined and quantified at single-cell resolution in human airways using RNA in-situ hybridization, immunostaining, and single-cell RNA sequencing. Complementary studies in mice utilized immunostaining, primary murine tracheal epithelial cell (MTEC) culture, and transgenic mice with secretory or ciliated cell-specific knockout of pIgR. SIgA degradation by human neutrophil elastase or secreted bacterial proteases from non-typeable Haemophilus influenzae (NTHi) was evaluated in vitro.

ResultsWe found that secretory cells are the predominant cell type responsible for pIgR expression in human and murine airways. Loss of SIgA in small airways was not associated with a reduction in secretory cells but rather a reduction in pIgR protein expression despite intact PIGR mRNA expression. Neutrophil elastase and NTHi-secreted proteases are both capable of degrading SIgA in vitro and may also contribute to a deficient SIgA immunobarrier in COPD.

InterpretationLoss of the SIgA immunobarrier in small airways of patients severe COPD is complex and likely results from both pIgR-dependent defects in IgA transcytosis and SIgA degradation.

Key MessagesWhat is the key question? Localized SIgA deficiency in small airways is an established driver of COPD pathogenesis, but the mechanism of loss remains unclear. We hypothesized loss of SIgA is due to reduced numbers of pIgR-expressing cells in SIgA-deficient small airways.

What is the bottom line? pIgR is primarily expressed by secretory cells in human and murine airways. Although numbers of secretory cells are similar between SIgA-deficient and SIgA-replete airways in COPD, there is reduced expression of pIgR protein, but not mRNA, in SIgA-deficient airways. Additionally, host and bacterial proteases degrade SIgA in vitro, suggesting loss of SIgA may relate to both impaired transcytosis and increased degradation.

Why read on? This study highlights the complexity of SIgA immunobarrier maintenance and suggests that strategies aimed at restoring the SIgA immunobarrier will need to account for both impaired transcytosis and degradation by host and/or bacterial proteases. Dramatic consequences of reducing erythrocyte membrane cholesterol on Plasmodium falciparum Cholesterol is the most abundant lipid in the erythrocyte. During its blood stage development, the malaria parasite establishes an active cholesterol gradient across the various membrane systems within the infected erythrocyte. Interestingly, some antimalarial compounds have recently been shown to disrupt cholesterol homeostasis in intraerythrocytic stages of Plasmodium falciparum. These studies point to the importance of cholesterol for parasite growth. Previously, reduction of cholesterol from the erythrocyte membrane by treatment with methyl-{beta}-cyclodextrin (M{beta}CD) was shown to inhibit parasite invasion and growth. In addition, M{beta}CD treatment of trophozoite stage P. falciparum was shown to result in parasite expulsion from the host cell. We have revisited these phenomena by using live video microscopy, ultrastructural analysis, and response to antimalarial compounds. By using time-lapse video microscopy of fluorescently tagged parasites, we show that M{beta}CD treatment for just 30 min causes dramatic expulsion of the trophozoite stage parasites. This forceful expulsion occurs within 10 sec. Remarkably, the plasma membrane of the host cell from which the parasite has been expelled does not appear to be compromised. The parasitophorous vacuolar membrane (PVM) continued to surround the extruded parasite, but the PVM appeared damaged. Treatment with antimalarial compounds targeting PfATP4 or PfNCR1 prevented M{beta}CD-mediated extrusion of the parasites, pointing to a potential role of cholesterol dynamics underlying the expulsion phenomena. We also confirmed the essential role of erythrocyte plasma membrane cholesterol for invasion and growth of P. falciparum. This defect can be partially complemented by cholesterol and desmosterol but not with epicholesterol, revealing stereospecificity underlying cholesterol function. Overall, our studies advance previous observations and reveal unusual cell biological features underlying cholesterol depletion of infected erythrocyte plasma membrane.

ImportanceMalaria remains a major challenge in much of the world. Symptoms of malaria are caused by the growth of parasites belonging to Plasmodium spp. inside the red blood cells (RBC), leading to their destruction. The parasite depends upon its host for much of its nutritional need. Cholesterol is a major lipid in the RBC plasma membrane, which is the only source of this lipid for malarial parasites. We have previously shown that certain new antimalarial compounds disrupt cholesterol homeostasis in P. falciparum. Here we use live time-lapse video microscopy to show dramatic expulsion of the parasite from the host RBC when the cholesterol content of the RBC is reduced. Remarkably, this expulsion is inhibited by the antimalarials that disrupt lipid homeostasis. We also show stereospecificity of cholesterol in supporting parasite growth inside RBC. Overall, these results point to a critical role of cholesterol in physiology of malaria parasites. Neprosin belongs to a new family of glutamic peptidase family based on in silico evidence Neprosin was first discovered in the insectivorous tropical pitcher plants of Nepenthes species as a novel protease with prolyl endopeptidase (PEP) activity. Neprosin has two uncharacterized domains of neprosin activation peptide and neprosin. A previous study has shown neprosin activity in hydrolyzing proline-rich gliadin, a gluten component that triggers celiac disease. In this study, we performed in silico structure-function analysis to investigate the catalytic mechanism of neprosin. Neprosin sequences lack the catalytic triad and motifs of PEP family S9. Protein structures of neprosins from Nepenthes x ventrata (NvNpr) and N. rafflesiana (NrNpr1) were generated by ab initio methods and comparatively assessed to obtain high-quality models. Structural alignment of models to experimental structures in the Protein Data Bank (PDB) found a high structural similarity to glutamic peptidases. Further investigations reveal other resemblances to the glutamic peptidases with low optimum pH that activates the enzyme via autoproteolysis for maturation. Two highly conserved glutamic acid residues, which are stable according to the molecular dynamics simulation, can be found at the active site of the substrate cleft. Protein docking demonstrated that mature neprosins bind well with potent antigen I-gliadin at the putative active site. Taken together, neprosins represent a new glutamic peptidase family, with a putative catalytic dyad of two glutamic acids. This study illustrates a hypothetical enzymatic mechanism of the neprosin family and demonstrates the useful application of an accurate ab initio protein structure prediction in the structure-function study of a novel protein family. 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. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. A family of unusual immunoglobulin superfamily genes in an invertebrate histocompatibility complex Most colonial marine invertebrates are capable of allorecognition, the ability to distinguish between themselves and conspecifics. One long-standing question is whether invertebrate allorecognition genes are homologous to vertebrate histocompatibility genes. In the cnidarian Hydractinia symbiolongicarpus, allorecognition is controlled by at least two genes, Allorecognition 1 (Alr1) and Allorecognition 2 (Alr2), which encode highly polymorphic cell surface proteins that serve as markers of self. Here, we show that Alr1 and Alr2 are part of a family of 41 Alr genes, all of which reside a single genomic interval called the Allorecognition Complex (ARC). Using sensitive homology searches and highly accurate structural predictions, we demonstrate that the Alr proteins are members of the immunoglobulin superfamily (IgSF) with V-set and I-set Ig domains unlike any previously identified in animals. Specifically, their primary amino acid sequences lack many of the motifs considered diagnostic for V-set and I-set domains, yet they adopt secondary and tertiary structures nearly identical to canonical Ig domains. Thus, the V-set domain, which played a central role in the evolution of vertebrate adaptive immunity, was present in the last common ancestor of cnidarians and bilaterians. Unexpectedly, several Alr proteins also have immunoreceptor tyrosine-based activation motifs (ITAMs) and immunoreceptor tyrosine-based inhibitory motifs (ITIMs) in their cytoplasmic tails, suggesting they could participate in pathways homologous to those that regulate immunity in humans and flies. This work expands our definition of the IgSF with the addition of a family of unusual members, several of which play a role in invertebrate histocompatibility.

Significance StatementThe immunoglobulin superfamily (IgSF) is one of the largest and most functionally versatile domain families in animal genomes. Although their amino acid sequences can vary considerably, IgSF domains have been traditionally defined by conserved residues at several key positions in their fold. Here, we sequenced an invertebrate histocompatibility complex and discovered a family of IgSF genes with amino acid sequences that lack most of these residues yet are predicted to adopt folds virtually identical to canonical V-set and I-set IgSF domains. This work broadens the definition of the IgSF and shows that the V-set domain was present earlier in animal evolution than previously appreciated. Game of theories: an agent-based model of how variation in contest cost-benefits determines the evolution of contest resolution rules Game-theory based models are used to understand rules that animals use to settle contests over indivisible resources. However, the empirical literature of contests indicates controversial support to models, with some species supporting different models and other species showing no support to any model. Since strategies used to resolve contests may have different associated costs, it is possible that different conditions have determined the evolution of distinct assessment strategies used by animals. We used an agent-based model to explore the importance of the following conditions: resource availability, probability of reproduction with resource, and damage costs on evolution of assessment strategies. We used self- and mutual-assessment models as a heurist framework to build agents with different assessment strategies. In our model, agents competed for resources in scenarios with different combinations of resource availability, probability of reproduction with resource, and damage costs. We found that agents following the self-assessment with damage strategy were prevalent in scenarios with no probability of reproduction without the resource, independently of other variables. We also found that agents following the non-aggressive strategy occurred in all scenarios. However, agents using the non-aggressive strategy were prevalent only in scenarios with probability of reproduction with the resource. Finally, we observed that agents using mutual-assessment occurred only in a scenario with high risk of damage, low availability of resources, and with probability of reproduction without the resource. These results indicate that agents following the self-assessment with damage and non-aggressive strategies may be able to stay at most scenarios. Inhibitory proteins block substrate access to the membrane-embedded active site of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma. Cell type matching in single-cell RNA-sequencing data using FR-Match Reference cell type atlases powered by single cell transcriptomic profiling technologies have become available to study cellular diversity at a granular level. We present FR-Match for matching query datasets to reference atlases with robust and accurate performance for identifying common and novel cell types and suboptimally clustered cell types in the query data. FR-Match shows excellent performance for cross-platform, cross-sample type, and cross-tissue region cell type matching. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

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

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

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

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. State-dependent network interactions differentially gate sensory input at the motor and command neuron level in Caenorhabditis elegans Neural responses are influenced by both external stimuli and internal network states. While network states have been linked to behavioral and stimulus states, little is known about how sensory inputs are filtered by whole-brain activity to affect motor and command neurons. Here, we recorded whole-brain activity of Caenorhabditis elegans experiencing bacterial food stimuli, and modeled how sensory inputs affect motor and command neurons in a network state-dependent manner. First, we classified active neurons into six functional clusters: two sensory neuron clusters (ON, OFF), and four motor/command neuron clusters (AVA, RME, SMDD, SMDV). Using encoding models, we found that ON and OFF sensory neurons that respond to onset and removal of bacteria, respectively, employ different adaptation strategies. Next, we used decoding models to show that bacterial onset and removal differentially drive AVA and RME cluster activity. To explore state-dependent effects on AVA and RME clusters, we developed a model that identified network states and fitted submodels for each state to predict how each of the six functional clusters drive AVA and RME cluster activity. We also identified network states in which AVA and RME clusters were either largely unperturbed by or receptive to bacterial sensory input. Furthermore, this model allowed us to disentangle the state-dependent contributions of stimulus timescales and bacterial content to neural activity. Collectively, we present an interpretable approach for modeling network dynamics that goes beyond implication of neurons in particular states, and moves toward explicitly dissecting how neural populations work together to produce state dependence.

Significance StatementA major function of the brain is to transform sensory information into behavior. As the first receiver of sensory input, sensory neuron activity is often most correlated with stimulus features. However, this high-fidelity representation of sensory input becomes diluted as it travels to downstream neurons, where sensory information is integrated with network activity. By the time sensory information reaches motor neurons, it is often difficult to dissociate the influence of sensory input from the influence of network activity. Here, we describe a method that is fully interpretable such that we can show how neural populations on a whole-brain scale interact to produce network states. From there, we can attribute motor neuron activity to network history and sensory input. 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. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. Force profile of the two-handed hardstyle kettlebell swing in novice older adults: an exploratory profile BackgroundUnderstanding the force profile of an exercise increases clinical confidence when assessing the benefits and potential risks of a prescribed exercise. This exploratory study presents the force profile of the hardstyle kettlebell swing in novice older adults and compares peak force with kettlebell deadlifts. These data will help inform healthcare providers and coaches who are considering prescribing kettlebell exercises for older adults.

MethodsThirty-five community-dwelling males and females (59-79 years) were recruited, from applicants to participate in the BELL trial. Two-handed hardstyle swings were performed with 8-16 kg. Deadlifts were performed with 8-24 kg and 8-32 kg for females and males, respectively. Ground reaction force was obtained from a floor-mounted force platform. Pairwise comparisons of peak force, forward force, rate of force development, swing cadence, sex, and kettlebell mass, were investigated for the kettlebell swing, with representative force-time curves described. Pairwise comparisons of peak force, sex and kettlebell mass were investigated for the deadlift, with comparisons of peak force, kettlebell mass, and sex, between swings and deadlifts.

ResultsFor kettlebells up to 16 kg, paired samples T-tests show a large exercise effect ({delta} > 1.4) with peak force higher for swings than deadlifts. Data shows: (i) higher peak force during swings than deadlifts ({delta} = 1.77), reaching 4.5 (1.0) N.kg-1, (ii) peak force during an 8 kg swing was greater than a 32kg deadlift, (iii) negligible difference in normalised peak force between males and females performing kettlebell swings, but a moderately large effect size during deadlifts (males > females, {delta} = 0.69), (iv) mean rate of force development of 19.9 (4.7) N.s-1.kg-1 with a very weak, positive correlation with kettlebell mass (y = 14.4 + 0.32x), and trivial effect of sex, (v) mean forward force equal to 5.5% of vertical force during swings, increasing from 3.8 (1.6) % with 8 kg to 7.1 (2.6) % with 16 kg.

ConclusionDuring kettlebell swings, there is negligible difference in normalised net peak force between novice males and females using the same absolute loads. Where ground reaction force is a therapeutic target, kettlebell swings with an 8 kg kettlebell could have similar effects to much heavier deadlifts (>24 kg). Kettlebell swings performed with lighter loads, could provider similar therapeutic value to much heavier deadlifts, and may be a more appealing, affordable, and convenient option for older adults. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Identification of CD4+ Sub-population of Resident Cardiac Fibroblasts Linked to Myocardial Fibrosis Infiltration with inflammatory T-cells and accumulation of cardiac myofibroblasts are hallmarks of cardiac fibrosis and maladaptive remodeling. The origin, identity, and functions of the resident cardiac cells involved in this process are, however, unclear. To determine the identity of cells contained in regions exhibiting fibrosis, mass cytometry profiling was performed using resident human ventricular cardiac fibroblasts and right ventricle autopsy tissues from individuals diagnosed with pulmonary hypertension and SUGEN/hypoxia rats. Results showed that a subpopulation of resident myocardial fibroblasts expresses increased levels of CD4+, a helper T-cell surface marker, in addition to mesenchymal markers in humans and rats. Characterization of the resident cardiac fibroblast subpopulation, both structurally and functionally, using transcriptome and secretome analysis of the secreted cytokines, chemokines, proteins, and metabolites, evidenced that IL-1{beta} induces a phenotypic switch of human cardiac fibroblasts from mesenchymal to CD4+ lymphoidal lineage in vitro. RNA sequencing (RNA-seq) analysis of FACS-sorted CD4-expressing cardiac fibroblasts further revealed that the transcriptome of such IL-1{beta}-induced CD4+ fibroblast population exhibited classical lymphoidal and stem cell-like signatures. Lastly, reversal of cell clustering, phosphorylation of MAPK p38 and NF-{kappa}B p65, and phenotypic switching was achieved with the administration of an IL-1R antagonist. In conclusion, we have identified a subpopulation of cardiac fibroblasts which exhibits structural and functional attributes of both mesenchymal and lymphoid cells which is induced by IL-1{beta}-IL-1R-NFkB pathway for differentiation of cardiac fibroblast cells. These data suggest that cardiac fibroblast transdifferentiation during inflammation may form the basis for maladaptive remodeling during myocardial fibrosis.