Faster short-read mapping with strobemer seeds constructed from syncmers Read alignment to genomes is a fundamental computational step used in many bioinformatic analyses, and often, it is the computational bottleneck. Therefore, it is desirable to perform the alignment step as fast as possible without compromising accuracy. Most alignment algorithms consider a seed-and-extend approach, where the time-consuming seeding step identifies and decides on candidate mapping locations. Recently, several advances have been made on seeding methods for fast sequence comparison.

We combine two such methods, syncmers and strobemers, in a novel seeding approach for constructing dynamic-sized fuzzy seeds and implement the method in a short-read aligner, strobealign. We benchmark strobealign to other traditional and recently proposed aligners on simulated and biological data and show that strobealign is an order of magnitude faster than traditional aligners such as BWA and Bowtie2 while being both faster and more accurate than more recently proposed aligners. Our aligner can free up substantial time and computing resources needed for read alignment in many pipelines.

Availabilityhttps://github.com/ksahlin/strobealign.
GEM-DeCan: Improved tumor immune microenvironment profiling through novel gene expression and DNA methylation signatures predicts immunotherapy response Quantifying the proportion of the different cell types present in tumor biopsies remains a priority in cancer research. So far, a number of deconvolution methods have emerged for estimating cell composition using reference signatures, either based on gene expression or on DNA methylation from purified cells. These two deconvolution approaches could be complementary to each other, leading to even more performant signatures, in cases where both data types are available. However, the potential relationship between signatures based on gene expression and those based on DNA methylation remains underexplored.

Here we present five new deconvolution signature matrices, based on DNA methylation or RNAseq data, which can estimate the proportion of immune cells and cancer cells in a tumour sample. We test these signature matrices on available datasets for in-silico and in-vitro mixtures, peripheral blood, cancer samples from TCGA, bone marrow from multiple myeloma patients and a single-cell melanoma dataset. Cell proportions estimates based on deconvolution performed using our signature matrices, implemented within the EpiDISH framework, show comparable or better correlation with FACS measurements of immune cell-type abundance and with various estimates of cancer sample purity and composition than existing methods.

Using publicly available data of 3D chromatin structure in haematopoietic cells, we expanded the list of genes to be included in the RNAseq signature matrices by considering the presence of methylated CpGs in gene promoters or in genomic regions which are in 3D contact with these promoters. Our expanded signature matrices have improved performance compared to our initial RNAseq signature matrix. Finally, we show the value of our signatures in predicting patient response to immune checkpoint inhibitors in three melanoma and one bladder cancer cohorts, based on bulk tumour sample gene expression.

We also provide GEM-DeCan: a snakemake pipeline, able to run an analysis from raw sequencing data to deconvolution based on various gene expression signature matrices, both for bulk RNASeq and DNA methylation data. The code for producing the signature matrices and reproducing all the figures of this paper is available on the GEM-DeCan repository.
The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. SufB intein splicing in Mycobacterium tuberculosis is influenced by two remote conserved N-extein Histidines Inteins are auto-processing domains that implement a multi-step biochemical reaction termed protein splicing, marked by cleavage and formation of peptide bonds. They excise from a precursor protein, generating a functional protein via covalent bonding of flanking exteins. We report the kinetic study of splicing and cleavage reaction in [Fe-S] cluster assembly protein SufB from Mycobacterium tuberculosis. Although it follows a canonical intein splicing pathway, distinct features are added by extein residues present in the active site. Sequence analysis identified two conserved histidines in the N-extein region; His-5 and His-38. Kinetic analyses of His-5Ala and His-38Ala SufB mutants exhibited significant reductions in splicing and cleavage rates relative to the SufB wild-type precursor protein. Structural analysis and molecular dynamics simulations suggested that Mtu SufB displays a unique mechanism where two remote histidines work concurrently to facilitate N-terminal cleavage reaction. His-38 is stabilized by the solvent-exposed His-5, and can impact N-S acyl shift by direct interaction with the catalytic Cys1. Development of inteins as biotechnological tools or as pathogen specific novel antimicrobial targets requires a more complete understanding of such unexpected roles of conserved extein residues in protein splicing. Low photorespiratory capacity is sufficient to block the induction of photosynthetic capacity in high light The induction of high photosynthetic capacity in high light (HL) is a common response among many herbaceous dicot plants, however, the signals that control this response remain largely unknown. Here, multiple independent lines of evidence are presented in support of the conclusion that low photorespiratory capacity acts a negative signal to limit photosynthetic capacity acclimation in HL in Arabidopsis thaliana. Using a panel of natural accessions, primary nitrogen (N) assimilation and photorespiration rates early after a shift to growth in HL, as well as activities for key enzymes in these pathways, were shown to positively correlate with the magnitude of the subsequent induction of photosynthetic capacity, which occurred several days later. Time-resolved metabolomic data during acclimation to HL were collected using a strongly acclimating ecotype and a weakly acclimating ecotype, revealing in greater detail the differences in N assimilation, photorespiration, and triose-phosphate utilization pathways underlying efficient photosynthetic capacity acclimation. When shifted into HL growth conditions under non-photorespiratory conditions, weakly acclimating ecotypes and even photorespiratory mutants gained the ability to strongly induce high photosynthetic capacity in HL. Thus, a negative, photorespiration-dependent signal early in the HL shift appears to block photosynthetic capacity acclimation in accessions with low photorespiratory capacity, whereas accessions with high photorespiratory capacity are licensed to increase photosynthetic capacity. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Priority conservation areas and a global population estimate for the Critically Endangered Philippine Eagle derived from modelled range metrics using remote sensing habitat characteristics Many range-restricted taxa are currently experiencing population declines yet lack fundamental information regarding distribution and population size. Establishing baseline estimates for both these key biological parameters is however critical for directing conservation planning for at-risk range-restricted species. The International Union for the Conservation of Nature (IUCN) Red List uses three range metrics that define species distributions and inform extinction risk assessments: extent of occurrence (EOO), area of occupancy (AOO) and area of habitat (AOH). However, calculating all three metrics using standard IUCN approaches relies on a geographically representative sample of locations, which for rare species is often spatially biased. Here, we apply model-based interpolation using Species Distribution Models (SDMs), correlating occurrences with remote-sensing covariates, to calculate IUCN range metrics, protected area coverage and a global population estimate for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi). Our final range wide continuous SDM had high predictive accuracy (Continuous Boyce Index = 0.927) and when converted to a binary model estimated an AOH = 23,185 km2, a maximum EOO = 605,759 km2, a minimum EOO = 272,272 km2, with an AOO = 53,867 km2. Based on inferred habitat from the AOH metric, we estimate a global population of 318 breeding pairs (range: 258-362 pairs), or 636 mature individuals, across the Philippine Eagle global range. Protected areas covered 34 % of AOH, 15 % less than the target representation, with the continuous model identifying key habitat as priority conservation areas. We demonstrate that even when occurrences are geographically biased, robust habitat models can be built that enable quantification of baseline IUCN range metrics, protected area coverage, and a population size estimate. In the absence of adequate location data for many rare and threatened taxa, our method is a promising spatial modelling tool with widespread applications, in particular for island endemics facing high extinction risk. Admp regulates tail bending by controlling ventral epidermal cell polarity via phosphorylated myosin translocation. The tails transient but pronounced ventral bending is found in many chordate embryos and provides an interesting model of how tissue interactions contribute to the overall shape (Lu et al., 2020). Here, we identify one of the upstream regulators of tail bending. Admp knock-down embryo (Admp MO) completely inhibited ventral tail bending. Admp MO showed specific inhibition of the cell-cell intercalation ventral epidermis. In the early stages of the intercalation, each ventral epidermal cell exhibited a boat-shaped morphology (boat cell) with a narrow apical surface in which phosphorylated myosin (pMLC) accumulated. In the Admp MO, the number of the boat cells reduced and pMLC was observed at the basal side in the Admp MO embryo. Finally, laser cutter experiments showed a pMLC-dependent relaxation of the ventral tail bending, which was observed when each ventral midline epidermal cell was cut along the apico-basal axis. Based on these results, we propose a novel function for Admp, in which Admp controls the pMLC translocation to the apical side of the ventral epidermal cells and leads the tail bending by ventral resistance to the notochord extension.

Summary StatementAdmp is an upstream regulator of the tail bending in the protochordate Ciona tailbud embryo, regulating tissue polarity of the ventral midline epidermis by the translocation of the phosphorylated myosin. Contribution of arm movements to balance recovery after tripping in older adults. Falls are common in daily life, often caused by trips and slips and, particularly in older adults, with serious consequences. Although arm movements play an important role in balance control, there is limited research into the role of arm movements during balance recovery after tripping in older adults. We investigated how older adults use their arms to recover from a trip and the difference in the effects of arm movements between fallers (n=5) and non-fallers (n=11).

Sixteen older males and females (69.7{+/-}2.3 years) walked along a walkway and were occasionally tripped over suddenly appearing obstacles. We analysed the first trip using a biomechanical model based on full-body kinematics and force-plate data to calculate whole body orientation during the trip and recovery phase. With this model, we simulated the effects of arm movements at foot-obstacle impact and during trip recovery on body orientation.

Apart from an increase in sagittal plane forward body rotation at touchdown in fallers, we found no significant differences between fallers and non-fallers in the effects of arm movements on trip recovery. Like earlier studies in young adults, we found that arm movements during the recovery phase had most favourable effects in the transverse plane: by delaying the transfer of angular momentum of the arms to the body, older adults rotated the tripped side more forward thereby allowing for a larger recovery step. Older adults that are prone to falling might improve their balance recovery after tripping by learning to prolong ongoing arm movements. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. SARS-CoV-2 spike-glycoprotein processing at S1/S2 and S2' and shedding of the ACE2 viral receptor: roles of Furin and TMPRSS2 and implications for viral infectivity and cell-to-cell fusion The spike (S)-protein of SARS-CoV-2 binds ACE2 and requires proteolytic "priming" at PRRAR685{downarrow} into S1 and S2 (cleavage at S1/S2), and "fusion-activation" at a S2 site for viral entry. In vitro, Furin cleaved peptides mimicking the S1/S2 cleavage site more efficiently than at the putative S2, whereas TMPRSS2 inefficiently cleaved both sites. In HeLa cells Furin-like enzymes mainly cleaved at S1/S2 during intracellular protein trafficking, and S2 processing by Furin at KPSKR815{downarrow} was strongly enhanced by ACE2, but not for the optimized S2 KRRKR815{downarrow} mutant (S2), whereas individual/double KR815AA mutants were retained in the endoplasmic reticulum. Pharmacological Furin-inhibitors (Boston Pharmaceuticals, BOS-inhibitors) effectively blocked endogenous S-protein processing in HeLa cells. Furthermore, we show using pseudotyped viruses that while entry by a "pH-dependent" endocytosis pathway in HEK293 cells did not require Furin processing at S1/S2, a "pH-independent" viral entry in lung-derived Calu-3 cells was sensitive to inhibitors of Furin (BOS) and TMPRSS2 (Camostat). Consistently, these inhibitors potently reduce infectious viral titer and cytopathic effects, an outcome enhanced when both compounds were combined. Quantitative analyses of cell-to-cell fusion and spike processing revealed the key importance of the Furin sites for syncytia formation. Our assays showed that TMPRSS2 enhances fusion and proteolysis at S2 in the absence of cleavage at S1/S2, an effect that is linked to ACE2 shedding by TMPRSS2. Overall, our results indicate that Furin and TMPRSS2 play synergistic roles in generating fusion-competent S-protein, and in promoting viral entry, supporting the combination of Furin and TMPRSS2 inhibitors as potent antivirals against SARS-CoV-2.

IMPORTANCESARS-CoV-2 is the etiological agent of COVID-19 that resulted in >5 million deaths. The spike protein (S) of the virus directs infection of the lungs and other tissues by binding the angiotensin-converting enzyme 2 (ACE2) receptor. For effective infection, the S-protein is cleaved at two sites: S1/S2 and S2. Cleavage at S1/S2, induces a conformational change favoring the recognition of ACE2. The S2 cleavage is critical for cell-to-cell fusion and virus entry into host cells. Our study contributes to a better understanding of the dynamics of interaction between Furin and TMPRSS2 during SARS-CoV-2 entry and suggests that the combination of a non-toxic Furin inhibitor with a TMPRSS2 inhibitor could significantly reduce viral entry in lung cells, as evidenced by an average synergistic [~]95% reduction of viral infection. This represents a powerful novel antiviral approach to reduce viral spread in individuals infected by SARS-CoV-2 or future related coronaviruses. BRCA2 BRC missense variants disrupt RAD51-dependent DNA repair Pathogenic mutations in the BRCA2 tumor suppressor gene predispose to breast, ovarian, pancreatic, prostate, and other cancers. BRCA2 maintains genome stability through homology-directed repair (HDR) of DNA double-strand breaks (DSBs) and replication fork protection. Nonsense or frameshift mutations leading to truncation of the BRCA2 protein are typically considered pathogenic, however, missense mutations resulting in single amino acid substitutions can be challenging to functionally interpret. The majority of missense mutations in BRCA2 have been classified as Variants of Uncertain Significance (VUS) with unknown functional consequences. In this study, we identified three BRCA2 VUS located within the BRC repeat region to determine their impact on canonical HDR and fork protection functions. We provide evidence that S1221P and T1980I, which map to conserved residues in the BRC2 and BRC7 repeats, compromise the cellular response to chemotherapeutics and ionizing radiation, and display deficits in fork protection. We further demonstrate biochemically that S1221P and T1980I disrupt RAD51 binding and diminish the ability of BRCA2 to stabilize RAD51-ssDNA complexes. The third variant, T1346I, located within the spacer region between BRC2 and BRC3 repeats, is fully functional. We conclude that T1346I is a benign allele whereas S1221P and T1980I are hypomorphic disrupting the ability of BRCA2 to fully engage and stabilize RAD51 nucleoprotein filaments. Our results underscore the importance of correctly classifying BRCA2 VUS as pathogenic variants can impact both future cancer risk and guide therapy selection during cancer treatment.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=191 SRC="FIGDIR/small/461749v2_ufig1.gif" ALT="Figure 1">
View larger version (53K):
org.highwire.dtl.DTLVardef@1e5e3c6org.highwire.dtl.DTLVardef@1cd0e60org.highwire.dtl.DTLVardef@16adf4corg.highwire.dtl.DTLVardef@16a0089_HPS_FORMAT_FIGEXP  M_FIG C_FIG Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. 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. Protection induced by a human monoclonal antibody recognizing two different epitopes in a conserved region of streptococcal M proteins Group A streptococci have evolved multiple strategies to evade human antibodies, making it challenging to create effective vaccines or antibody treatments. Here, we have generated antibodies derived from the memory B cells of an individual who had successfully cleared a group A streptococcal infection. The antibodies bind with high affinity in the central region of the surface-bound M protein. Such antibodies are typically non-opsonic. However, one antibody could effectively promote vital immune functions, including phagocytosis and in vivo protection. Remarkably, this antibody primarily interacts through a bivalent dual-Fab cis mode, where the Fabs bind to two distinct epitopes in the M protein. The dual-Fab cis binding phenomenon is conserved across different groups of M types. In contrast, other antibodies binding with normal single-Fab mode to the same region can not bypass the M proteins virulent effects. A broadly binding, protective monoclonal antibody could be a candidate for anti-streptococcal therapy. Our findings highlight the concept of dual-Fab cis binding as a means to access conserved, and normally non-opsonic regions, for protective antibody targeting. Divergent Acyl Carrier Protein Decouples Mitochondrial Fe-S Cluster Biogenesis from Fatty Acid Synthesis in Malaria Parasites Plasmodium falciparum malaria parasites are early-diverging eukaryotes with many unusual metabolic adaptations. Understanding these adaptations will give insight into parasite evolution and unveil new parasite-specific drug targets. Most eukaryotic cells retain a mitochondrial fatty acid synthesis (FASII) pathway whose acyl carrier protein (mACP) and 4-phosphopantetheine (Ppant) prosthetic group provide a soluble scaffold for acyl chain synthesis. In yeast and humans, mACP also functions to biochemically couple FASII activity to electron transport chain (ETC) assembly and Fe-S cluster biogenesis. In contrast to most eukaryotes, the Plasmodium mitochondrion lacks FASII enzymes yet curiously retains a divergent mACP lacking a Ppant group. We report that ligand-dependent knockdown of mACP is lethal to parasites, indicating an essential FASII-independent function. Decyl-ubiquinone rescues parasites temporarily from death, suggesting a dominant dysfunction of the mitochondrial ETC followed by broader cellular defects. Biochemical studies reveal that Plasmodium mACP binds and stabilizes the Isd11-Nfs1 complex required for Fe-S cluster biosynthesis, despite lacking the Ppant group required for this association in other eukaryotes, and knockdown of parasite mACP causes loss of both Nfs1 and the Rieske Fe-S protein in ETC Complex III. This work reveals that Plasmodium parasites have evolved to decouple mitochondrial Fe-S cluster biogenesis from FASII activity, and this adaptation is a shared metabolic feature of other Apicomplexan pathogens, including Toxoplasma and Babesia. This discovery also highlights the ancient, fundamental role of ACP in mitochondrial Fe-S cluster biogenesis and unveils an evolutionary driving force to retain this interaction with ACP independent of its eponymous function in FASII.

Significance StatementPlasmodium malaria parasites are single-celled eukaryotes that evolved unusual metabolic adaptations. Parasites require a mitochondrion for blood-stage viability, but essential functions beyond the electron transport chain are sparsely understood. Unlike yeast and human cells, the Plasmodium mitochondrion lacks fatty acid synthesis enzymes but retains a divergent acyl carrier protein (mACP) incapable of tethering acyl groups. Nevertheless, mACP is essential for parasite viability by binding and stabilizing the core mitochondrial Fe-S cluster biogenesis complex via a divergent molecular interface lacking an acyl-pantetheine group that contrasts with other eukaryotes. This discovery unveils an essential metabolic adaptation in Plasmodium and other human parasites that decouples mitochondrial Fe-S cluster biogenesis from fatty acid synthesis and evolved at or near the emergence of Apicomplexan parasitism. Targeted Gene Editing in Plants using CRISPR-Cas9, Single-Stranded DNA Oligonucleotides, and Protoplast Regeneration Versatile genome editing can be facilitated by the insertion of DNA sequences into specific locations. Current protocols involving clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins rely on low efficiency homology-directed repair or non-homologous end joining with modified double-stranded DNA oligonucleotides as donors. Our simple protocol eliminates the need for expensive equipment, chemical and enzymatic donor DNA modification or plasmid construction by using polyethylene glycol-calcium to deliver non-modified single-stranded DNA oligonucleotides and CRISPR-Cas9 ribonucleoprotein into protoplasts. Plants regenerated via edited protoplasts achieved targeted insertion frequencies of up to 50.0% in Nicotiana benthamiana and 13.6% in rapid cycling Brassica oleracea without antibiotic selection. Using a 60-nt donor containing 27 nt in each homologous arm, 6 of 22 regenerated N. benthamiana plants showed targeted insertions, and one contained a precise insertion of a 6-bp HindIII site. The inserted sequences were transmitted to the next generation and invite the possibility of future exploration of versatile genome editing by targeted DNA insertion in plants. Essential Regression - a generalizable framework for inferring causal latent factors from multi-omic human datasets High-dimensional cellular and molecular profiling of human samples highlights the need for analytical approaches that can integrate multi-omic datasets to generate predictive biomarkers and prioritized causal inferences. Current methods are limited by high dimensionality of the combined datasets, the differences in their data distributions and their integration to infer causal relationships. Here we present Essential Regression (ER), an interpretable machine learning approach for high-dimensional multi-omic datasets, that addresses these problems by identifying latent factors and their likely cause-effect relationships with the system-wide outcome/properties of interest. ER is a novel data-distribution-free latent-factor regression model that integrates multi-omic datasets and identifies latent factors significantly associated with an outcome. ER outperforms a range of state-of-the-art methods in terms of prediction performance on simulated datasets. ER can be coupled with probabilistic graphical modeling thereby strengthening the causal inferences. ER generates novel cellular and molecular predictions, using multi-omic human systems immunology datasets, pertaining to immunosenescence and immune dysregulation. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Accounting for stimulations that do not elicit motor-evoked potentials when mapping cortical representations of multiple muscles The representation of muscles in the cortex can be mapped using navigated transcranial magnetic stimulation. The mapping can be quantified via measures like the centre of gravity or the size of the region of excitability. Determining these measures typically relies only on stimulation points that yield motor-evoked potentials (MEPs); stimulations that do not elicit an MEP, i.e., non-MEP points, are ignored entirely. In this study, we show how incorporating non-MEP points may affect the estimates of the size and centroid of the excitable area in eight hand and forearm muscles after mono-phasic single-pulse TMS. We performed test-retest assessments in twenty participants and estimated the reliability of centroids and sizes of the corresponding areas using inter-class correlation coefficients. For most muscles, the reliability turned out good, if not excellent. As expected, removing the non-MEP points significantly decreased area sizes and area weights, suggesting that conventional approaches that do not account for non-MEP points are likely to overestimate the regions of excitability. Disc and Actin-Associated Protein 1 Influence Attachment in the Intestinal Parasite Giardia lamblia The deep-branching eukaryote Giardia lamblia is an extracellular parasite that attaches to the host intestine via a microtubule-based structure called the ventral disc. Control of attachment is mediated in part by the movement of two regions of the ventral disc that either permit or exclude the passage of fluid under the disc. Several known disc-associated proteins (DAPs) contribute to disc structure and function, but no force-generating protein has been identified among them. We recently identified several Giardia actin (GlActin) interacting proteins at the ventral disc, which could potentially employ actin polymerization for force generation and disc conformational changes. One of these proteins, Disc and Actin Associated Protein 1 (DAAP1), is highly enriched at the two regions of the disc previously shown to be important for fluid flow during attachment. In this study, we investigate the role of both GlActin and DAAP1 in ventral disc morphology and function. We confirmed interaction between GlActin and DAAP1 through coimmunoprecipitation, and used immunofluorescence to localize both proteins throughout the cell cycle and during trophozoite attachment. Similar to other DAPs, the association of DAAP1 with the disc is stable, except during cell division when the disc disassembles. Depletion of GlActin by translation-blocking antisense morpholinos resulted in both impaired attachment and defects in the ventral disc, indicating that GlActin contributes to disc-mediated attachment. Depletion of DAAP1 through CRISPR interference resulted in intact discs but impaired attachment, gating, and flow under the disc. As attachment is essential for infection, elucidation of these and other molecular mediators is a promising area for development of new therapeutics against a ubiquitous parasite.

Author SummaryGiardia lamblia is a single-celled organism and one of the most common gastrointestinal parasites worldwide. In developing countries, recurrent Giardia infections are common, due to lack of access to clean water. Giardia infections can lead to diarrhea, vomiting, dehydration, disruption of the intestinal microbiome, and chronic infections can lead to irritable bowel syndrome. Because existing drug treatments have side effects and Giardias resistance to drugs is increasing, new treatment strategies are needed. The parasites attachment to the hosts intestine is mediated by a Giardia-specific structure that resembles a suction cup and is called the ventral adhesive disc. We previously identified DAAP1, a protein which interacts with Giardia actin and localizes to the ventral disc. Here, we explore the relationship between these two proteins and investigate their role in disc-based attachment. Most disc proteins, including DAAP1, are unrelated to any human proteins, making them appealing drug targets to inhibit parasite attachment and infection. Phenotypic plasticity is aligned with phenological adaptation on both micro- and macroevolutionary timescales In seasonally-variable environments, phenotypic plasticity in phenology may be critical for adaptation to fluctuating environmental (temperature) conditions. Using an 18-generation longitudinal dataset from natural insect (damselfly) populations, we show that phenology has strongly advanced. Individual fitness data suggest this is likely an adaptive response towards a thermally-dependent fitness optimum. A laboratory experiment revealed that developmental plasticity qualitatively matches environmental-dependence of selection, partially explaining observed phenological advance. Expanding our analysis to the macroevolutionary level, we use a database of over 1-million occurrence records and spatiotemporally-matched temperature data from 49 Swedish Odonate species to infer macroevolutionary dynamics of phenology. Phenological plasticity was more closely aligned with adaptation for species that have recently colonized northern latitudes, with more phenological mismatch at lower latitudes. Our results show that phenological plasticity plays a key role in microevolutionary dynamics within a single species, and such plasticity may have facilitated post-Pleistocene range expansion in this insect clade.