High-quality genome assembly of Cinnamomum burmannii (chvar. Borneol) provides insights into the natural borneol biosynthesis Cinnamomum burmannii (chvar. Borneol) is a well-known medicinal and industrial plant cultivated in the Lingnan region of China. It is the key source from organism of natural borneol (D-borneol), one of the precious and widely used Chinese herbal medicines with a variety of medicinal effects. Here, we report a high-quality chromosome-scale genome assembly of C. burmannii (chvar. Borneol) using Pacbio single-molecule sequencing and Hi-C technology. The assembled genome size was 1.14 GB with a scaffold N50 of 94.30 Mb, while 98.77% of the assembled sequences were anchored on 12 pseudochromosomes including 41549 protein-coding genes. Genomic evolution analysis revealed C. burmannii and C. micranthum shared two Lauraceae unique ancestral whole-genome duplication (WGD) events. Likewise, comparative genomic analysis showed strong collinearity between these two species. Besides, the analysis for Long Terminal Repeat Retrotransposons (LTR-RTs) indicated the outbreak of LTR-RTs insertion made a great contribution to the size difference of genomes between C. burmannii and C. micranthum. Furthermore, the candidate genes in pathway associated with natural borneol synthesis were identified on the genome and their differential expressions were analyzed in various biological tissues. We considered that several of genes in Mevalonate (MVA) Methylerythritol Phosphate (MEP) pathways or in downstream pathway have the potential to be the key factors in the biosynthesis of D-borneol. We also constructed the genome database (CAMD; http://www.cinnamomumdatabase.com/) of Cinnamomum species for a better data utilization in the future. All these results will enrich the genomic data of Lauraceae plants and facilitate genetic improvement of this commercially important plant.
High quality mapping of chromatin at or near the nuclear lamina for small numbers of cells The chromatin associated with the nuclear lamina (NL) is referred to as Lamina-Associated Domains (LADs). While mapping of this feature has been done using various technologies, technical limitations exist for each of the methods. Here, we present an adaptation of the Tyramide-Signal Amplification sequencing (TSA-seq) protocol, which we call chromatin pull down-based TSA-seq (cTSA-seq), that can be used to map chromatin regions at or near the NL from as little as 50,000 cells without using carriers. The cTSA-seq mapped regions are composed of LADs and smaller chromatin regions that fall within the chromatin B-compartment known to be enriched for heterochromatin and be present at the nuclear periphery. As a proof of principle, we used cTSA-seq to map chromatin at or near the assembling NL as cells exit mitosis and progress through early and later G1. Consistent with previous reports, lamin-B1 based cTSA-seq revealed that regions toward the distal ends of chromosomes are near or at the reassembling NL during early G1. The cTSA-seq mapping and analyses revealed similarity between the early G1 chromatin and oncogene-induced senescent cell populations. The cTSA-seq reported here represents a useful method for analyzing chromatin at or near the NL from small numbers of cells.
Small extracellular vesicles from young mice prevent frailty, improve healthspan and decrease epigenetic age in old mice. Aging is associated with an increased risk of frailty, disability, comorbidities, institutionalization, falls, fractures, hospitalization, and mortality. Searching for strategies to delay the degenerative changes associated with aging and frailty is interesting. We treated old animals intravenously with small extracellular vesicles (sEVs) derived from adipose mesenchymal stem cells (ADSCs) of young animals, and we found an improvement of several functional parameters usually altered with aging, such as motor coordination, grip strength, fatigue resistance, fur regeneration, and renal function. Frailty index analysis showed that 40% of old control mice were frail, whereas none of the old ADSCs-sEVs treated mice were. Molecular and structural benefits in muscle and kidney accompanied this functional improvement. ADSCs-sEVs induced pro-regenerative effects and a decrease in oxidative stress, inflammation, and senescence markers. Moreover, predicted epigenetic age was lower in tissues of old mice treated with ADSCs-sEVs and their metabolome changed to a youth-like pattern. Finally, we gained some insight into the miRNAs contained in sEVs that might be, at least in part, responsible for the effects observed. We propose that young sEVs treatment can be beneficial against frailty and therefore can promote healthy aging. Spindle position dictates division site during asymmetric cell division in moss Asymmetric cell division (ACD) underlies the development of multicellular organisms. It is considered that the division site in land plants is predetermined prior to mitosis and that the localization of the mitotic spindle does not govern the division plane. This contrasts with animal ACD, in which the division site is defined by active spindle-positioning mechanisms. Here, we isolated a hypomorphic mutant of the conserved microtubule-associated protein TPX2 in the moss Physcomitrium patens and observed abnormal spindle motility during cell division. This defect compromised the position of the division site and produced inverted daughter cell sizes in the first ACD of gametophore (leafy shoot) development. The phenotype was rescued by restoring endogenous TPX2 function and, unexpectedly, by depolymerizing actin filaments. Thus, we have identified an active spindle-positioning mechanism that, reminiscent of acentrosomal ACD in animals, involves microtubules and actin filaments and sets the division site in plants. 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. Decreased cell stiffness facilitates detachment and migration of breast cancer cells in 3D collagen matrices: An exploratory study BackgroundTumour-cell detachment is a critical early event in the metastatic cascade. However, the role of the cells mechanical properties in detachment and migration is not well understood. This exploratory study aimed to assess how intracellular stiffness changes these processes.

MethodsMDA-MB-231 cells were embedded as 10,000-cell spheroids in 2 and 4 mg/ml collagen matrices. Intracellular stiffness was assessed with mitochondria tracking microrheology of cells that migrated distances equivalent to four and six times the cell diameter (dC) from the spheroid and compared to cells at the spheroid surface (0dC), representing medium, high and no migration, respectively.

Results and discussionThe mitochondrial mean square displacement and intracellular stiffness decreased during detachment and migration for both collagen concentrations (i.e. rigidities). The mean square displacement of 4dC and 6dC cells was similar, whereas cell stiffness was lower for 4dC than for 6dC cells. With increasing matrix rigidity, the intracellular stiffness decreased for 0dC cells and did not change for 4dC and 6dC cells. It is proposed that decreased cell stiffness drives detachment and migration and increased matrix rigidity physically hinders migration, and cells need to become softer or remodel the environment to migrate. The independence of the stiffness of migrated cells from matrix rigidity suggests that cells remodel their environment through matrix proteins cleavage to migrate.

ConclusionsThe study revealed the collective effects of enhanced migratory conditions and increased matrix rigidity on the mechanical properties of the cells. The expression of matrix metalloproteinases and transforming growth factor {beta} and the role of cell volume on detachment and migration in matrices with varying pore sizes are proposed targets for further studies on metastatic cancer cells. An alternatively spliced TREM2 isoform lacking the ligand binding domain is expressed in human brain BackgroundGenetic variants in TREM2 are strongly associated with Alzheimers Disease (AD) risk but alternative splicing in TREM2 transcripts has not been comprehensively described.

ObjectiveRecognizing that alternative splice variants can result in reduced gene expression and/or altered function, we sought to fully characterize splice variation in TREM2.

MethodsHuman blood and anterior cingulate autopsy tissue from 61 donors were used for end-point and quantitative PCR and Western blotting to identify and quantify novel TREM2 isoforms.

ResultsIn addition to previously described transcripts lacking exon 3 or exon 4, or retaining part of intron 3, we identified novel isoforms lacking exon 2, along with isoforms lacking multiple exons. Isoforms lacking exon 2 were predominant at approximately 10% of TREM2 mRNA in the brain. Expression of TREM2 and frequency of exon 2 skipping did not differ between AD samples and non-AD controls (p = 0.1268 and p = 0.4909, respectively). Further, these novel splice isoforms were also observed across multiple tissues with similar frequency (range 5.3 - 13.0%). We found that the exon 2 skipped isoform D2-TREM2 is translated to protein and localizes similarly to full-length TREM2 protein, that both proteins are primarily retained in the Golgi complex, and that D2-TREM2 is expressed in AD and non-AD brain.

ConclusionSince the TREM2 ligand binding domain is encoded by exon 2, and skipping this exon retains reading frame while conserving localization, we hypothesize that D2-TREM2 acts as an inhibitor of TREM2 and targeting TREM2 splicing may be a novel therapeutic pathway for AD. 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. Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins - the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening. TRIM28-dependent SUMOylation protects the adult ovary from the testicular pathway Gonadal sexual fate in mammals is determined during embryonic development and must be actively maintained in adulthood. In the mouse ovary, oestrogen receptors and FOXL2 protect ovarian granulosa cells from transdifferentiation into Sertoli cells, their testicular counterpart. However, the mechanism underlying their protective effect is unknown. Here, we show that TRIM28 is required to prevent female-to-male sex reversal of the mouse ovary after birth. We found that upon loss of Trim28, ovarian granulosa cells transdifferentiate to Sertoli cells through an intermediate cell type, different from gonadal embryonic progenitors. TRIM28 is recruited on chromatin in the proximity of FOXL2 to maintain the ovarian pathway and to repress testicular-specific genes. The role of TRIM28 in ovarian maintenance depends on its E3-SUMO ligase activity that regulates the sex-specific SUMOylation profile of ovarian-specific genes. Our study identifies TRIM28 as a key factor in protecting the adult ovary from the testicular pathway. A computational modeling approach for predicting multicell patterns based on signaling-induced differential adhesion Differential adhesion within cell populations enables the emergence of unique patterns in heterogeneous multicellular systems. This process has been previously explored using synthetically engineered heterogenous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models may be leveraged to more systematically explore the key parameters that drive the emergence of different patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. The model was used to predict how varying initial cell seedings, cadherin induction probabilities, or homotypic adhesion strengths between cells impact spheroid patterning, and unsupervised machine learning techniques were used to map system parameters to unique spheroid patterns. The model was then deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern. A cross biomes bacterial diversity shed light on ocean-atmosphere microbial transmission. Microbes are ubiquitous in the oceans and the atmosphere, playing essential roles in biogeochemical processes. The bio-exchanges between the two environments can provide important insights into microbial distribution and diversity but are still not well understood. We simultaneously surveyed the genomic diversity of airborne and marine bacterial communities across 15 000 kilometers in the Atlantic and Pacific oceans. Higher variability of microbial community composition was observed in the atmosphere than in the ocean surface waters. In addition, a greater similarity was observed between oceans than their overlaying atmosphere, and between atmospheric samples than with the ocean beneath. We additionally detected a higher coverage rate and relative abundance of marine bacteria in the Pacific atmosphere as compared to the Atlantic, while the dominant fraction in the Atlantic atmosphere was annotated as soil-associated bacteria. This study advances our understanding of microbial dispersion in the ocean, the atmosphere, and the exchange between them, as well as their potential impact on microbial composition, ecology, and biogeochemistry. An isogenic collection of pluripotent stem cell lines with elevated α-synuclein expression -Synuclein (Syn) is a small, disordered protein that becomes aggregated in Lewy body diseases, such as Parkinsons disease (PD) and dementia with Lewy bodies (DLB). Human induced pluripotent stem cells (hiPSCs) potentially provide a tractable disease model to monitor early molecular changes associated with PD/DLB. We and others have previously derived hiPSC lines from patients with duplication and triplication of the SNCA gene, encoding for Syn. It is now recognised that to perform meaningful disease modelling with these hiPSC lines, it is critical to generate isogenic control cell lines that lack the disease causing mutations. In order to complement the existing and emerging hiPSC models for PD/DLB, we have generated an allelic series of Syn over-expressing hESC lines on the same isogenic background. An unresolved question is whether pluripotent stem cell lines, with elevated levels of Syn, can undergo efficient differentiation into dopaminergic and cortical neurons to model PD and DLB, respectively. We took advantage of our isogenic collection of hESC lines to determine if increased expression of Syn affects neural induction and neuronal differentiation. Clonal hESC lines with significantly different levels of Syn expression proliferated normally and maintained expression of pluripotent markers, such as OCT4. All cell lines efficiently produced PAX6+ neuroectoderm and there was no correlation between Syn expression and neural induction efficiency. Finally, global transcriptomic analysis of cortical differentiation of hESC lines with low or high levels of Syn expression demonstrated robust and similar induction of cortical neuronal expression profiles. Gene expression differences observed were unrelated to neural induction and neuronal differentiation. We conclude that elevated expression of Syn in human pluripotent stem cells does not adversely affect their neuronal differentiation potential and that collections of isogenic cell lines with differing levels of Syn expression are valid and suitable models to investigate synucleinopathies. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. The conserved zinc-finger protein GTSF1helps PIWI proteins achieve their full catalytic potential Argonaute proteins use nucleic acid guides to find and bind specific DNA or RNA target sequences. Argonaute proteins can be found in all kingdoms of life, and play diverse biological functions including genome defense, gene regulation, and chromosome partitioning. Many Argonautes retain their ancestral endoribonuclease activity, cleaving the phosphodiester bond between target nucleotides t10 and t11. In animals, a specialized class of Argonautes, the PIWI proteins, use 21-35 nt PIWI-interacting RNAs (piRNAs) to direct transposon silencing, protect the germline genome, and regulate gene expression during gametogenesis1. The piRNA pathway is required for fertility in one or both sexes of nearly all animals. Both piRNA production and function require RNA cleavage catalyzed by PIWI proteins. Spermatogenesis in mice and other placental mammals requires three distinct, developmentally regulated PIWI proteins: MIWI (PIWIL1), MILI (PIWIL2), and MIWI2 (PIWIL4)2-4. The piRNA-guided endoribonuclease activities of MIWI and MILI are essential to produce functional sperm5, 6. piRNA-directed silencing in mice and insects also requires Gametocyte-Specific Factor 1 (GTSF1), a PIWI-associated protein of unknown function7-12. Here, we report that GTSF1 potentiates the weak, intrinsic, piRNA-directed RNA cleavage activities of PIWI proteins, transforming them into efficient endoribonucleases. GTSF1 represents the first example of an auxiliary protein that potentiates the catalytic activity of an Argonaute protein. The combined treatment of Molnupiravir and Favipiravir results in a marked potentiation of efficacy in a SARS-CoV2 hamster infection model through an increased frequency of mutations in the viral genome. Favipiravir and Molnupiravir, orally available antivirals, have been reported to exert antiviral activity against SARS-CoV2. In recent days preliminary efficacy data have been reported in COVID-19 patients. We here studied the combined antiviral effect of the drugs in the SARS-CoV2 hamster infection model. We first demonstrate that Molnupiravir can reduce infectious virus titers in lungs of infected animals in a dose-dependent manner by up to 3.5 log10 which is associated with a marked improvement of virus-induced lung pathology. When animals are treated with a combination of suboptimal doses of Molnupiravir and Favipiravir (that each alone result in respectively a 1.3 log10 and 1.1 log10 reduction of infectious virus titers in the lungs), a marked combined potency is observed. Infectious virus titers in the lungs of animals treated with the combo are on average reduced by 4.5 log10 and infectious virus are no longer detected in the lungs of 60% of treated infected animals. Both drugs result in an increased mutation frequency of the remaining viral RNA recovered from the lungs. In the combo-treated hamsters an increased frequency of C-to-T and G-to-A mutations in the viral RNA is observed as compared to the single treatment groups which may explain the pronounced antiviral potency of the combination. Our findings may lay the basis for the design of clinical studies to test the efficacy of the combination of Molnupiravir and Favipiravir in the treatment of COVID-19. Distinct repertoire development history during affinity maturation of antigen-responding VHH antibody Time-resolved antibody repertoire analysis is essential to elucidate the antibody maturation process during immunization. However, most analyses are merely snapshot comparisons of unimmunized and fully immunized animals. In this study, an alpaca was immunized with multiple antigen injections and its repertoire development was traced by high-throughput next-generation sequencing every week for 3 months. The use of alpaca-derived, single-domain antibody circumvents the problems associated with conventional antibodies including heavy-chain and light-chain matching and labor-intensive production. The sequences were processed into clusters and the antibodies were generated to determine whether the clusters included antigen-binding antibodies. For most of the antigen-responding clusters, repertoire development exhibited distinct time-dependent patterns compared with clusters unresponsive to antigens. Moreover, the sequences at the tips and root of the phylogenetic tree had strong and weak antigen-binding affinity, respectively. The foregoing characteristics could be exploited in the prediction of antigen-responding clusters and the physical properties of antibodies.

TeaserClusters of antigen-responding, single-domain antibodies exhibited distinct time-dependent repertoire development history A vascularized tumoroid model for human glioblastoma angiogenesis Glioblastoma (GBM) angiogenesis is critical for tumor growth and recurrence, making it a compelling therapeutic target. Here, a disease-relevant, vascularized tumoroid in vitro model with stem-like features and stromal surrounds is reported. The model is used to recapitulate how individual components of the GBMs complex brain microenvironment such as hypoxia, vasculature-related stro-mal cells and growth factors support GBM angiogenesis. It is scalable, tractable, cost-effective and can be used with biologically-derived or biomimetic matrices. Patient-derived primary GBM cells are found to closely participate in blood vessel formation in contrast to a GBM cell line containing differentiated cells. Exogenous growth factors amplify this effect under normoxia but not at hypoxia suggesting that a significant amount of growth factors is already being produced under hypoxic conditions. Under hypoxia, primary GBM cells and umbilical vein endothelial cells are found to strongly co-localize in a mosaic pattern to form sprouting vascular networks, with GBM cells acquiring an endothelial-like behaviour, which has been reported to occur in vivo. These findings demonstrate that our 3D tumoroid in vitro model exhibits biomimetic attributes that may permit its use as a preclinical model in studying microenvironment cues of tumor angiogenesis. Tackling polypharmacology of kinase inhibitors and selecting chemical kinase probes by transcription factor activity profiling. Protein kinase inhibitors (PKI) evolved as promising drugs for multiple diseases. However, PKIs promiscuity causes polypharmacological effects challenging therapeutic development. We show that the polypharmacology of PKIs can be inferred from their impact on transcription factors (TF) lying at the apexes of signaling pathways. Using a high-content reporter assay, we evaluated TF activity profiles (TFAP) in cells treated with Akt, CDK, Aurora, Raf, MEK, and ERK inhibitors. For each kinase family, we found kinase-specific consensus TFAP signatures signifying the on-target PKI activity. Remarkably, proximal kinases had the highest similarity of their PKIs consensus signatures. Furthermore, we show that individual PKIs exhibited the consensus, "on-target signatures" within certain "specificity windows" and distinct "off-target signatures" at other concentrations. We also show that the TFAP signatures permit pinpointing the off-target PKI activity at irrelevant kinases and non-kinase effectors. Therefore, the TFAP approach provides clear-cut quantitative metrics for assessing PKI activity and concentration ranges where the on-target activity dominates cell response. By illuminating PKI biology invisible to target-based techniques, this effect-based approach streamlines the selection of kinase chemical probes and PKI drug prioritization. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Navigating the amino acid sequence space between functional proteins using a deep learning framework MotivationShedding light on the relationships between protein se-quences and functions is a challenging task with many implications in protein evolution, diseases understanding, and protein design. Protein sequence / function space is however hard to comprehend due to its com-plexity. Generative models help to decipher complex systems thanks to their abilities to learn and recreate data specificity. Applied to protein sequences, they can point out relationships between protein positions and functions capture the sequence patterns associated with functions or navigate through uncharted area of molecular evolution.

ResultsIn this study, an unsupervised generative approach based on adversarial auto-encoder (AAE) is proposed to generate and explore new sequences with respect to their functions thanks to the prior distribution allowing a continuous exploration of the latent space. AAEs are tested on three protein families known for their multiple functions. Clustering re-sults on the encoded sequences from the latent space computed by AAEs display high level of homogeneity regarding the protein sequence func-tions. The study also reports and analyzes for the first time two sampling strategies based on latent space interpolation and latent space arithmetic to generate intermediate protein sequences sharing sequential and functional properties of original sequences issued from different families and functions. Generated sequences by interpolation between latent space data points demonstrate the ability of the AAE to generalize and to pro-duce meaningful biological sequences from an evolutionary uncharted area of the biological sequence space. Finally, 3D structure models generated by comparative modelling between different combinations of structures of different sub-families and of generated sequences from latent space or sub-family sequences point out to the ability of the latent space arithmetic to successfully transfer functional properties between sub-families. All in all this study confirms the ability of deep learning frameworks to model biological complexity and bring new tools to explore amino acid sequence and functional spaces.

AvailabilityCode and data used for this study are freely available at https://github.com/T-B-F/aae4seq.

Contacttristan@bitardfeildel.fr

Supplementary informationSupplementary data are available at online. Genome-wide analysis of the non-coding RNA synthetic genetic network reveals extensive plasticity and unique phenotypes in yeast Recent studies identified non-coding RNAs (ncRNAs) with unknown function that are responsible for major fitness changes in yeast. To understand ncRNA interplay and aid their functional assignment, the synthetic genetic array (SGA) methodology was employed to create >15,000 double mutants and to score their epistasis in different environments. Unlike the protein network, ncRNAs mostly displayed positive epistasis in rich medium. Interestingly, the negative interactions significantly increased under stressors, showing environmental-dependent functions for ncRNAs. No correlation was found between the network of ncRNAs and that of their neighbouring genes, suggesting functional independence. The U3 paralogs, SNR17A and SNR17B, share the majority of genetic interactions in rich medium as expected. For example, SUT480 interacted with both paralogs and its function was linked to 18S rRNA processing. However, under stressors, a large number of unique epistatic interactions were observed, supporting the notion that SNR17A and SNR17B have diverged and sub-functionalised after genome duplication. Oscillations support short latency co-firing of neurons during human episodic memory formation. Theta and gamma oscillations in the medial temporal lobe are suggested to play a critical role for human memory formation via establishing synchrony in neural assemblies. Arguably, such synchrony facilitates efficient information transfer between neurons and enhances synaptic plasticity, both of which benefit episodic memory formation. However, to date little evidence exists from humans that would provide direct evidence for such a specific role of theta and gamma oscillations for episodic memory formation. Here we investigate how oscillations shape the temporal structure of neural firing during memory formation in the medial temporal lobe. We measured neural firing and local field potentials in human epilepsy patients via micro-wire electrode recordings to analyze how brain oscillations promote co-incidences of firing between neurons during successful and unsuccessful encoding of episodic memories. The results show that phase-coupling of neurons to faster theta and gamma oscillations correlates with co-firing at short latencies ([~]20-30 ms) and occurs during successful memory formation. Phase-coupling at slower oscillations in these same frequency bands, in contrast, correlates with longer co-firing latencies and occurs during memory failure. Thus, our findings provide evidence which implicate neural oscillations into the synchronization of neural firing in the medial temporal lobe during the encoding of episodic memories. Transcription factors in the fungus Aspergillus nidulans: Markers of genetic innovation, network rewiring and conflict between genomics and transcriptomics Gene Regulatory Networks (GRNs) are shaped by the democratic/hierarchical relationships among transcription factors (TFs) and associated proteins, together with the cis-regulatory sequences (CRSs) bound by these TFs at target promoters. GRNs control all cellular processes, including metabolism, stress-response, growth and development. Due to the ability to modify morphogenetic and developmental patterns, there is the consensus view that the reorganization of GRNs is a driving force of species evolution and differentiation. Duplication of TF-coding genes, their divergent sequence evolution and gain/loss/modification of CRSs are events rewiring GRNs. Fungi (mainly Saccharomycotina) have served as a reference kingdom for the study of GRN evolution. Here, I studied the genes predicted to encode TFs in the fungus Aspergillus nidulans (Pezizomycotina). The analysis of the expansion of different families of TFs suggests that duplication of TFs impacts the species level and that the expansion in Zn2Cys6 TFs is mainly due to dispersed duplication events. Comparison of genomic annotation and transcriptomic data suggests that a significant percentage of genes should be re-annotated while many others remain silent. Finally, a new regulator of growth and development is identified and characterized. Overall, this study establishes a novel theoretical framework in synthetic biology, since overexpression of silent TF forms would provide additional tools to assess how GRNs are rewired.

Impact StatementGRNs controlling fungal sexual and asexual development, stress response, carbon/nitrogen metabolism or secondary metabolism, have been widely characterized in the fungus A. nidulans and, in several cases, were characterized for the first time in this fungus. The functional and genetic relationships among the TFs of these networks are known. By comparing the genomes/proteomes of hundreds of fungal genomes available in public repositories, researchers can track the emergence of regulators of interest, and identify groups of paralog TFs. The present work uses basic bioinformatics tools and standard RNA-seq analysis tools to compare genomic annotation and transcriptomics data, find incorrectly annotated TFs, identify paralog clusters of TFs, track the trajectories of duplication events and expansion of bHLH-type TFs, and identify TFs with a role in the control of growth and development. The study of TFs and GRNs of A. nidulans offers the opportunity to dissect and analyze sequentially the contribution of each cis and trans regulatory element to network rewiring and the modification of developmental patterns.

Data SummaryThe data underlying this article are available in supplementary tables and supplementary figures. RNA-seq, genomic and proteomic data used is this article (and previously deposited by others; see references and accession numbers in Materials and Methods) were retrieved from sources (www.ebi.ac.uk/arrayexpress/; www.ncbi.nlm.nih.gov/bioproject/; https://www.ncbi.nlm.nih.gov/geo/) and databases (https://fungidb.org/fungidb/app) in the public domain. Aspergillus nidulans strains generated in this article will be shared on reasonable request to the corresponding author.

One-sentence summaryA study of the predicted set of transcriptional regulators in the fungus Aspergillus nidulans, their conservation patterns and the correlation between gene annotations and transcriptomic data Greater spear nosed bats in Panama do not use social proximity to improve foraging efficiency Group foraging animals may change their spatial relationships in response to the resource landscape. The strength of social network ties across foraging groups could then reflect the energy balance of individual group members through their effect on foraging efficiency. Greater spear-nosed bats (Phyllostomus hastatus) live in social groups where unrelated females develop group-specific social calls, defend each others offspring, and are hypothesized to coordinate social foraging to feed on balsa flowers. To test if social groups forage together, and if this increases foraging efficiency, we tracked 18 individuals from three social groups of P. hastatus on Isla Colon, Panama using high frequency GPS for 1 - 3 days each. We found that individuals commuted alone but then joined group members in small foraging patches with high densities of flowering balsas, 25 km from their roost on the mainland. Proximity-based social network measures did not predict foraging efficiency or energy expenditure, but individuals with higher social associations spent more time outside of the cave. Individuals did not coordinate movement, but close proximity among individuals increased resting bout duration and synchrony, particularly when in the same patch. These results show that stronger social proximity associations were more closely tied to resting behavior, and indicate that factors other than increased feeding efficiency may structure network relationships of group members while foraging.

Lay SummarySocial foraging can increase foraging efficiency on hard-to-find food sources. Miniaturized high-frequency GPS allowed us to test this on groups of greater spear-nosed bats. We estimated food intake rates, energy expenditure, and the strength of bat social relationships. In contrast to our expectations, we did not find a relationship between bat social relationships and energy expenditure or foraging efficiency. However, bats took longer resting breaks when group members were close. Why these bats commute 25 km individually across the ocean to then forage near each other thus requires further investigation. 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.