title: Temporal characteristics of functional connectome reconfigurations are heritable and related to cognition.

abstract: The brains functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly fractional occupancy (FO) and transition probability (TP), describing the trajectory of state transitions. Genetic effects explained a substantial proportion of phenotypic variance of these features (h2=.39, 95% CI= [.24.54] for FO; h2=.43, 95% CI=[.29.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the duration spent in each connectivity configuration and the frequency of shifting between configurations, and such temporal features may act as endophenotypes for cognitive abilities.


title: Dynamic trajectories of connectome state transitions are heritable

category: neuroscience


title: P53 is a direct regulator of the immune co-stimulatory molecule CD80

abstract: Increasing evidence indicates oncogenes and tumor suppressors not only influence cell fitness but can also control the immunophenotype of cells. Here, we examined how 34 commonly mutated genes in colorectal cancer (CRC) may influence the expression of 8 key immunomodulatory proteins. To do this, we employed a functional genomics approach utilizing Pro-Code/CRISPR libraries for high-dimensional analysis. We introduced a library of 102 Pro-Code/gRNA combinations, targeting each of the 34 genes, in CT26 cells, a CRC cell model, and measured the expression of each of the immunomodulatory proteins by CyTOF mass cytometry. Notably, cells carrying a Pro-Code/CRISPR targeting the Trp53 lost expression of the immune co-stimulatory molecule CD80. Validation confirmed that Trp53 knockout resulted in the loss of CD80 and that activation of P53, through DNA damage or stabilization, resulted in CD80 upregulation. P53 ChIP-seq identified the CD80 promoter as a direct target of P53. CD80 regulation by P53 was identified in other cells, including normal epithelial cells and macrophages. Functionally, CD80 reduction caused by P53 loss led to a reduced capacity for CRC to prime antigen-specific T cells. These studies establish CD80, a canonical co-stimulatory molecule, as a direct target of the tumor suppressor and DNA damage response gene, P53.

category: immunology


title: Microbial community structure in rice, crops, and pastures rotation systems with different intensification levels in the temperate region of Uruguay

abstract: Rice is an important crop in Uruguay associated mostly with livestock production in a rice and pasture rotation system since the 1920s. However, in recent years there has been interest in intensifying the production in some of these systems to satisfy market demands and increase income. Intensification occurs by augmenting the rice frequency in the rotation, including new crops like sorghum and soybean, or shortening the pasture phase. A long-term experiment was established in 2012 in the main rice producing area of Uruguay with the objective to study the impact of intensification in rice rotations. After the first cycle of rotation soils from seven rotation phases were sampled and microbial communities were studied by means of high-throughput sequencing of Illumina NovaSeq 6000. Archaeal/bacterial and fungal community composition were studied (16S rRNA and 18S gene regions) detecting 3662 and 807 bacterial and fungal Operational Taxonomic Units (OTUs), respectively. Actinobacteria, Firmicutes and Proteobacteria were the most common bacterial phyla. Among them, only Proteobacteria differed significantly between rotations. Although most fungal OTUs were unidentified, Ascomycota, Basidiomycota and Mucoromycota were the most abundant fungal classes within identified taxa. Bacterial communities differed between rotations forming three groups according to the percentage of rice in the system. Fungal communities clustered in four groups, although not well differentiated, and mostly associated with the antecessor crop. Only P and C:N varied between rotations among soil physicochemical variables after six years, and individual bacterial OTUs appeared weakly influenced by P, pH, Mg and fungal OTUs by P. The results suggest that after six years, bacteria/archaeal communities were influenced by the time with rice in the rotation, and fungal communities were more influenced by the antecessor crop. More studies are needed to associate fungal communities with certain rotational or environmental variables. Some taxa were associated with a particular rotation, and some bacterial taxa were identified as biomarkers. Fungal indicator taxa were not identified at the species level for any rotation.

category: microbiology


title: Sex matters: effects of sex and mating in the presence and absence of a protective microbe

abstract: In most animals, female investment in offspring production is greater than for males. Lifetime reproductive success (LRS) is predicted to be optimised in females through extended lifespans to maximise reproductive events. Extended lifespan can be achieved through increased investment in immunity. Males, however, maximise lifetime reproductive success by obtaining as many matings as possible. Microbe-mediated protection (MMP) is known to affect both immunity and reproduction, but whether the two sexes respond differently to the provision of MMP remains to be explored. Here, we investigated the sex-specific differences in host life history traits between female and male Caenorhabditis elegans following pathogenic infection with Staphylococcus aureus with or without MMP by Enterococcus faecalis. Overall, female survival decreased with increased mating. With MMP, females increased investment into offspring production, while males displayed higher behavioural activity. These results highlight the different strategies employed by the two sexes under pathogen infection with and without MMP.


title: The gut bacterial natural product colibactin triggers induction of latent viruses in diverse bacteria

abstract: Colibactin is a chemically unstable small molecule genotoxin produced by multiple different bacteria, including members of the human gut microbiome.1,2 While the biological activity of colibactin has been extensively investigated in mammalian systems,3 little is known about its effects on other microorganisms. Here, we discover that colibactin targets bacteria carrying prophages, inducing lytic development via the bacterial SOS response. DNA, added exogenously, protects bacteria from colibactin, as does expressing a colibactin resistance protein (ClbS) in non-colibactin-producing cells. The prophage-inducing effects we observe apply broadly across taxonomically diverse phage-bacteria systems. Finally, we identify bacteria that possess colibactin resistance genes but lack colibactin biosynthetic genes. Many of these bacteria are infected with predicted prophages, and we show that the expression of their ClbS homologs provides immunity from colibactin-triggered induction. Our study reveals a mechanism by which colibactin production could impact microbiomes and highlights an underappreciated role for microbial natural products in influencing population-level events such as phage outbreaks.


title: Cryo-ET reveals two major tubulin-based cytoskeleton structures in Toxoplasma gondii

abstract: In the obligate intracellular parasite, Toxoplasma gondii, the subpellicular microtubules (SPMTs) help maintain shape, while the apical conoid (also tubulin-based) is implicated in invasion. Here, we use cryo-electron tomography to determine the molecular structures of the SPMTs and the conoid-fibrils (CFs) in vitrified and detergent-lysed parasites. Subvolume densities from detergent-extracted parasites yielded averaged density maps at subnanometer resolutions, and these were related back to their architecture in situ. An intraluminal spiral (IS) lines the interior of the 13-protofilament SPMTs, revealing a preferred orientation of these microtubules relative to the parasites long axis. Each CF is composed of 9 tubulin protofilaments, that produce a comma-shaped cross-section, plus additional associated components. Conoid protrusion, a crucial step in invasion, is associated with an altered pitch of each CF. The use of basic building blocks of protofilaments and different accessory proteins in one organism, illustrates the versatility of these critical structures.

category: molecular biology


title: ZNF146/OZF and ZNF507 target LINE-1 sequences

abstract: Repetitive sequences including transposable elements (TEs) and transposon-derived fragments account for nearly half of the human genome. While transposition-competent TEs must be repressed to maintain genomic stability, mutated and fragmented TEs comprising the bulk of repetitive sequences can also contribute to regulation of host gene expression and broader genome organization. Here we analyzed published ChIP-seq data sets to identify proteins broadly enriched on TEs in the human genome. We show two of the proteins identified, C2H2 zinc finger-containing proteins ZNF146 (also known as OZF) and ZNF507, are targeted to distinct sites within LINE-1 ORF2 at thousands of locations in the genome. ZNF146 binding sites are found at old and young LINE-1 elements. In contrast, ZNF507 preferentially binds at young LINE-1 sequences correlated to sequence changes in LINE-1 elements at ZNF507s binding site. To gain further insight into ZNF146 and ZNF507 function, we disrupt their expression in HEK293 cells using CRISPR/Cas9 and perform RNA sequencing, finding modest gene expression changes in cells where ZNF507 has been disrupted. We further identify a physical interaction between ZNF507 and PRMT5, suggesting ZNF507 may target arginine methylation activity to LINE-1 sequences.

category: cell biology


title: Preservation of HIV-1 Gag helical bundle symmetry by bevirimat is central tomaturation inhibition

abstract: The assembly and maturation of human immunodeficiency virus type-1 (HIV-1) requires proteolytic cleavage of the Gag polyprotein. The rate-limiting step resides at the junction between the capsid protein CA and spacer peptide 1, which assembles as a six-helix bundle (6HB). bevirimat (BVM), the first-in-class maturation inhibitor drug, targets the 6HB and impedes proteolytic cleavage, yet the molecular mechanisms of its activity, and relatedly, the escape mechanisms of mutant viruses, remain unclear. Here, we employed extensive molecular dynamics (MD) simulations and free energy calculations to quantitatively investigate molecular structure-activity relationships, comparing wild-type and mutant viruses in the presence and absence of BVM and inositol hexakisphosphate (IP6), an assembly cofactor. Our analysis shows that the efficacy of BVM is directly correlated with preservation of six-fold symmetry in the 6HB, which exists as an ensemble of structural states. We identified two primary escape mechanisms, and both lead to loss of symmetry, thereby facilitating helix uncoiling to aid access of protease. Our findings also highlight specific interactions that can be targeted for improved inhibitor activity and support the use of MD simulations for future inhibitor design.

category: biophysics


title: Unsupervised explainable AI for the collective analysis of a massive number of genome sequences: various examples from the small genome of pandemic SARS-CoV-2 to the human genome

abstract: In genetics and related fields, huge amounts of data, such as genome sequences, are accumulating, and the use of artificial intelligence (AI) suitable for big data analysis has become increasingly important. Unsupervised AI that can reveal novel knowledge from big data without prior knowledge or particular models is highly desirable for analyses of genome sequences, particularly for obtaining unexpected insights. We have developed a batch-learning self-organizing map (BLSOM) for oligonucleotide compositions that can reveal various novel genome characteristics. Here, we explain the data mining by the BLSOM: unsupervised and explainable AI. As a specific target, we first selected SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) because a large number of the viral genome sequences have been accumulated via worldwide efforts. We analyzed more than 0.6 million sequences collected primarily in the first year of the pandemic. BLSOMs for short oligonucleotides (e.g., 4~6-mers) allowed separation into known clades, but longer oligonucleotides further increased the separation ability and revealed subgrouping within known clades. In the case of 15-mers, there is mostly one copy in the genome; thus, 15-mers appeared after the epidemic start could be connected to mutations. Because BLSOM is an explainable AI, BLSOM for 15-mers revealed the mutations that contributed to separation into known clades and their subgroups. After introducing the detailed methodological strategies, we explained BLSOMs for various topics. The tetranucleotide BLSOM for over 5 million 5-kb fragment sequences derived from almost all microorganisms currently available and its use in metagenome studies. We also explained BLSOMs for various eukaryotes, such as fishes, frogs and Drosophila species, and found a high separation ability among closely related species. When analyzing the human genome, we found evident enrichments in transcription factor-binding sequences (TFBSs) in centromeric and pericentromeric heterochromatin regions. The tDNAs (tRNA genes) were separated by the corresponding amino acid.

category: bioinformatics


title: Interplay between protein stability, binding to ACE2 and escape from neutralizing antibodies determines the natural selection of SARS-CoV-2 receptor binding domain variants

abstract: Emergence of new SARS-CoV-2 variants has raised concerns at the effectiveness of vaccines and antibody therapeutics developed against the unmutated wild-type virus. We examined the effect of 12 most commonly occurring mutations in the receptor binding domain on its expression, stability, activity, and antibody escape potential-some of the factors that may influence the natural selection of mutants. Recombinant proteins were expressed in human cells. Stability was measured using thermal denaturation melts. Activity and antibody escape potential were measured using isothermal titration calorimetry in terms of binding to ACE2 and to a neutralizing human antibody CC12.1, respectively. Our results show that variants differ in their expression levels with the two least stable variants showing lesser expression. Out of the 8 well-expressed mutants, only 2 (N501Y and K417T/E484K/N501Y) showed stronger affinity to ACE2, 4 (Y453F, S477N, T478I and S494P) have similar affinity, whereas the other 2 (K417N and E484K) have weaker affinity when compared to the wild-type. In terms of CC12.1 binding, when compared to the wild-type, 4 variants (K417N, Y453F, N501Y and K417T/E484K/N501Y) have weaker affinity, 2 (S477N and S494P) have similar affinity, and 2 (T478I and E484K) have stronger affinity. Taken together, these results indicate that multiple factors contribute towards the natural selection of variants, and all these factors need be considered to understand the evolution of the virus. In addition, since not all variants can escape a given neutralizing antibody, antibodies to treat new variants can be chosen based on the specific mutations in that variant.


title: Interplay between receptor binding, immune escape and protein stability determines the natural selection of SARS-CoV-2 variants


title: Evaluation of mRNA-1273 against SARS-CoV-2 B.1.351 Infection in Nonhuman Primates

abstract:
BackgroundVaccine efficacy against the B.1.351 variant following mRNA-1273 vaccination in humans has not been determined. Nonhuman primates (NHP) are a useful model for demonstrating whether mRNA-1273 mediates protection against B.1.351.

MethodsNonhuman primates received 30 or 100 {micro}g of mRNA-1273 as a prime-boost vaccine at 0 and 4 weeks, a single immunization of 30 {micro}g at week 0, or no vaccine. Antibody and T cell responses were assessed in blood, bronchioalveolar lavages (BAL), and nasal washes. Viral replication in BAL and nasal swabs were determined by qRT-PCR for sgRNA, and histopathology and viral antigen quantification were performed on lung tissue post-challenge.

ResultsEight weeks post-boost, 100 {micro}g x2 of mRNA-1273 induced reciprocal ID50 neutralizing geometric mean titers against live SARS-CoV-2 D614G and B.1.351 of 3300 and 240, respectively, and 430 and 84 for the 30 {micro}g x2 group. There were no detectable neutralizing antibodies against B.1351 after the single immunization of 30 {micro}g. On day 2 following B.1.351 challenge, sgRNA in BAL was undetectable in 6 of 8 NHP that received 100 {micro}g x2 of mRNA-1273, and there was a [~]2-log reduction in sgRNA in NHP that received two doses of 30 {micro}g compared to controls. In nasal swabs, there was a 1-log10 reduction observed in the 100 {micro}g x2 group. There was limited inflammation or viral antigen in lungs of vaccinated NHP post-challenge.

ConclusionsImmunization with two doses of mRNA-1273 achieves effective immunity that rapidly controls lower and upper airway viral replication against the B.1.351 variant in NHP.


title: Evidence for Deleterious Original Antigenic Sin in SARS-CoV-2 Immune Response

abstract:
A previous report demonstrated the strong association between the presence of antibodies binding to an epitope region from SARS-CoV-2 nucleocapsid, termed Ep9, and COVID-19 disease severity. Patients with anti-Ep9 antibodies (Abs) had hallmarks of antigenic imprinting (AIM), including early IgG upregulation and cytokine-associated injury. Thus, the immunological memory of a previous infection was hypothesized to drive formation of suboptimal anti-Ep9 Abs in severe COVID-19 infections. This study identifies a putative primary antigen capable of stimulating production of cross-reactive, anti-Ep9 Abs. Binding assays with patient blood samples directly show cross-reactivity between Abs binding to Ep9 and only one bioinformatics-derived, homologous potential antigen, a sequence derived from the neuraminidase protein of H3N2 Influenza A virus. This cross-reactive binding is highly influenza strain specific and sensitive to even single amino acid changes in epitope sequence. The neuraminidase protein is not present in the influenza vaccine, and the anti-Ep9 Abs likely resulted from the widespread influenza infection in 2014. Therefore, AIM from a previous infection could underlie some cases of COVID-19 disease severity.

ImportanceInfections with SARS-COV-2 result in diverse disease outcomes, ranging from asymptomatic to fatal. The mechanisms underlying different disease outcomes remain largely unexplained. Previously, our laboratory identified a strong association between the presence of an antibody and increased disease severity in a subset of COVID-19 patients. Here, we report that this severity-associated antibody cross-reacts with viral proteins from an influenza A viral strain from 2014. Therefore, we speculate that antibodies generated against previous infections, like the 2014 influenza A, play a significant role in directing some peoples immune responses against SARS-COV-2. Such understanding of the sources and drivers of COVID-19 disease severity can help early identification and pre-emptive treatment.


title: Evidence for Deleterious Antigenic Imprinting in SARS-CoV-2 Immune Response


title: Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core

abstract: The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind, fuse, and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes the preparation of the fusion machinery by dissociation of the S1 domains. We report the development of bottom-up coarse-grained (CG) models validated with cryo-electron tomography (cryo-ET) data, and the use of CG molecular dynamics simulations to investigate the dynamical mechanisms involved in viral binding and exposure of the S2 trimeric core. We show that spike trimers cooperatively bind to multiple ACE2 dimers at virion-cell interfaces. The multivalent interaction cyclically and processively induces S1 dissociation, thereby exposing the S2 core containing the fusion machinery. Our simulations thus reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.


title: Standard Sample Storage Conditions Impact on Inferred Microbiome Composition and Antimicrobial Resistance Patterns

abstract:
Storage of biological specimens is crucial in the life and medical sciences. The storage conditions for samples can be different for a number of reasons, and it is unclear which effect this can have on the inferred microbiome composition in metagenomics analyses. Here, we assess the effect of common storage temperatures (deep freezer: -80{degrees}C, freezer: -20{degrees}C, fridge: 5{degrees}C, room temperature: 22{degrees}C) and storage times, (immediate sample processing: 0h, next day: 16h, over weekend: 64h, and longer term: 4, 8, 12 months), as well as repeated sample freezing and thawing (2-4 freeze-thaw cycles). We examine two different pig feces and sewage samples, unspiked and spiked with a mock community, and in triplicates, respectively, amounting to a total of 438 samples (777 Gbp; 5.1 billion reads). Storage conditions had a significant and systematic effect on the taxonomic and functional composition of microbiomes. Distinct microbial taxa and antimicrobial resistance classes were in some situations similarly effected across samples, while others were not, suggesting an impact of individual inherent sample characteristics. With an increasing number of freeze-thaw cycles, an increasing abundance of Firmicutes, Actinobacteria, and eukaryotic microorganisms was observed. We include recommendations for sample storage, and strongly suggest including more detailed information in the metadata together with the DNA sequencing data in public repositories to better facilitate meta-analyses and reproducibility of findings.

IMPORTANCEPrevious research has reported effects of DNA isolation, library preparation, and sequencing technology on metagenomics-based microbiome composition; however, the effect of biospecimen storage conditions has not been thoroughly assessed. We examined the effect of common sample storage conditions on metagenomics-based microbiome composition and find significant and, in part, systematic effects. Repeated freeze-thaw cycles could be used to improve the detection of microorganisms with more rigid cell walls, including parasites. We provide a dataset that could also be used for benchmarking algorithms to identify and correct for batch effects. Overall, the findings suggest that all samples of a microbiome study should be stored in the same way. Furthermore, there is a need to mandate more detailed information about sample storage and processing published together with the DNA sequencing data at INSDC (ENA/EBI, NCBI, DDBJ) or other repositories.


title: Selection for cell yield does not reduce overflow metabolism in E. coli

abstract:
Overflow metabolism is ubiquitous in nature, and it is often considered inefficient because it leads to a relatively low biomass yield per consumed carbon. This metabolic strategy has been described as advantageous because it supports high growth rates during nutrient competition.

Here we experimentally evolved bacteria without nutrient competition by repeatedly growing and mixing millions of parallel batch cultures of E. coli. Each culture originated from a water-in-oil emulsion droplet seeded with a single cell. Unexpectedly we found that overflow metabolism (acetate production) did not change. Instead the numerical cell yield during the consumption of the accumulated acetate increased as a consequence of a reduction in cell size. Our experiments and a mathematical model show that fast growth and overflow metabolism followed by the consumption of the overflow metabolite, leads to a higher numerical cell yield and therefore a higher fitness compared to full respiration of the substrate. This provides an evolutionary scenario where overflow metabolism can be favourable even in the absence of nutrient competition.


title: Intimate genetic relationships and fungicide resistance in multiple strains of human pathogenic fungus Aspergillus fumigatus isolated from a plant bulb

abstract: Fungal infections are increasingly dangerous because of environmentally-dispersed resistance to antifungal drugs. Azoles are commonly used antifungal drugs, but they are also used as fungicides in agriculture, which may enable enrichment of azole-resistant strains of the human pathogen Aspergillus fumigatus in the environment. Understanding of environmental dissemination and enrichment of genetic variation associated with azole resistance in A. fumigatus is required to suppress resistant strains. Here, we focused on eight strains of azole-resistant A. fumigatus isolated from a single tulip bulb for sale in Japan. This set includes strains with TR34/L98H/T289A/I364V/G448S and TR46/Y121F/T289A/S363P/I364V/G448S mutations in the cyp51A gene, which showed higher tolerance to several azoles than strains harboring TR46/Y121F/T289A mutation. The strains were typed by microsatellite typing, single nucleotide polymorphism profiles, and mitochondrial and nuclear genome analyses. The strains grouped differently using each typing method, suggesting historical genetic recombination among the strains. Our data also revealed that some strains isolated from the tulip bulb showed tolerance to other classes of fungicide, such as QoI and carbendazim, followed by related amino acid alterations in the target proteins. Considering spatial-temporal factors, plant bulbs are an excellent environmental niche for fungal strains to encounter partners, and to obtain and spread resistance-associated mutations.


title: A defined medium based on R2A for cultivation and exometabolite profiling of soil bacteria

abstract:
Exometabolomics is an approach to assess how microorganisms alter their environments through the depletion and secretion of chemical compounds. Comparisons of inoculated with uninoculated media can be used to provide direct biochemical observations on depleted and secreted metabolites which can be used to predict resource competition, cross-feeding and secondary metabolite production in microbial isolates and communities. This approach is most powerful when used with defined media that enable tracking of all depleted metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through LC-MS/MS. Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse bacteria but is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 53 phylogenetically diverse soil bacterial isolates and all of its metabolites were trackable through LC-MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for cultivating and characterizing diverse microbial isolates and communities.

Originality-Significance StatementWe build a defined medium based on the metabolite composition of R2A medium and soil, elemental stoichiometry requirements, and knowledge of metabolite usage by different bacteria. The newly formulated defined medium was evaluated on its ability to support the growth of soil isolates and its application for metabolite utilization profiling. We found that of 53 phylogenetically diverse soil bacterial isolates grew on the defined medium and all of its metabolites were trackable through LC-MS/MS analysis. This demonstrates the viability and utility of the constructed defined medium for cultivating and characterizing diverse microbial isolates and communities.


title: Isolation of SARS-CoV-2 B.1.1.28.2 P2 variant and pathogenicity comparison with D614G variant in hamster model

abstract:
BackgroundConsidering the potential threat from emerging SARS-CoV-2 variants and the rising COVID-19 cases, SARS-CoV-2 genomic surveillance is ongoing in India. We report herewith the isolation of the P.2 variant (B.1.1.28.2) from international travelers and further its pathogenicity evaluation and comparison with D614G variant (B.1) in hamster model.

MethodsVirus isolation was performed in Vero CCL81 cells and genomic characterization by next generation sequencing. The pathogenicity of the isolate was assessed in Syrian hamster model and compared with B.1 variant.

ResultsB.1.1.28.2 variant was isolated from nasal/throat swabs of international travelers returned to India from United Kingdom and Brazil. The B.1.1.28.2 variant induced body weight loss, viral replication in the respiratory tract, lung lesions and caused severe lung pathology in infected Syrian hamster model in comparison, with B.1 variant infected hamsters. The sera from B.1.1.28.2 infected hamsters efficiently neutralized the D614G variant virus whereas 6-fold reduction in the neutralization was seen in case of D614G variant infected hamsters sera with the B.1.1.28.2 variant.

ConclusionsB.1.1.28.2 lineage variant could be successfully isolated and characterization could be performed. Pathogenicity of the isolate was demonstrated in Syrian hamster model and in comparison, with B.1 variant was found more pathogenic. The findings of increased disease severity and neutralization reduction is of great concern and point towards the need for screening the vaccines for efficacy.


title: Epigenomic signature of the progeroid Cockayne syndrome exposes distinct and common features with physiological ageing

abstract: Cockayne syndrome (CS) and UV-sensitivity syndrome (UVSS) are rare genetic disorders caused by mutation of the DNA repair and chromatin remodelling proteins CSA or CSB, but only CS patients display a progeroid and neurodegenerative phenotype. As epigenetic modifications constitute a well-established hallmark of ageing, we characterized genome-wide DNA methylation (DNAm) of fibroblasts from CS versus UVSS patients and healthy donors. The analysis of differentially methylated positions and regions revealed a CS-specific epigenetic signature, enriched in developmental transcription factors, transmembrane transporters, and cell adhesion factors. The CS-specific signature compared to DNAm changes in other progeroid diseases and regular ageing, identifyied commonalities and differences in epigenetic remodelling. CS shares DNAm changes with normal ageing more than other progeroid diseases do, and according to the methylation clock CS samples show up to 13-fold accelerated ageing. Thus, CS is characterized by a specific epigenomic signature that partially overlaps with and exacerbates DNAm changes occurring in physiological aging. Our results unveil new genes and pathways that are potentially relevant for the progeroid/degenerative CS phenotype.


title: Dengue transmission risk in a changing climate: Bangladesh could experience a longer dengue fever season in the future

abstract: Our changing climate is already affecting the transmission of vector borne diseases such as dengue fever. This issue presents a significant public health concern for some nations, such as Bangladesh, which already experience regular seasonal outbreaks of dengue fever under present day conditions. To provide guidance for proactive public health planning to potentially mitigate future infections, we explore the impact of climate change on dengue infections by calculating the change in vectorial capacity of Aedes aegypti mosquito at a seasonal level for all regions in Bangladesh under two scenarios for future atmospheric greenhouse gas concentrations. For each of the four climate models used, and for both scenarios, our analysis reveals that the annual vectorial capacity remains at a level that would enable potential dengue epidemic transmission in all regions in Bangladesh. We found a slight decline in vectorial capacity in half of the regions examined during the last two decades of the 21st Century for the lower-concentration scenario, with a pronounced decline in vectorial capacity in all geographic regions beginning in 2060 for the higher-concentration scenario. The likely reason is that in many regions, warming is leading to sub-optimal mosquito breeding temperatures. However, seasonal differences in vectorial capacity dissipates as the climate warms, to the point that there is almost no observable seasonality for the higher-concentration scenario during the last two decades of this century. This finding suggests there is the potential for the dengue season to extend all year, with outbreaks occurring at any time. Consequently, disease surveillance and control activities would need to be geographically and temporally adapted to mitigate dengue epidemic risk in response to climate change.

category: ecology


title: De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

abstract: While it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human {beta}-globin (HBB) gene and to the identical, homologous{delta} -globin (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than Europeans in the HBB region studied. Finally, the rate of the 20A[->]T mutation, called the "HbS mutation" when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations, and no instances were observed in the European HBB ROI or in the European or African HBD ROI. Further studies will be needed to examine de novo mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

category: evolutionary biology


title: Nucleocapsid mutation R203K/G204R increases the infectivity, fitness and virulence of SARS-CoV-2

abstract: In addition to the mutations on the spike protein (S), co-occurring mutations on nucleocapsid (N) protein are also emerging in SARS-CoV-2 world widely. Mutations R203K/G204R on N, carried by high transmissibility SARS-CoV-2 lineages including B.1.1.7 and P.1, has a rapid spread in the pandemic during the past year. In this study, we performed comprehensive population genomic analyses and virology experiment concerning on the evolution, causation and virology consequence of R203K/G204R mutations. The global incidence frequency (IF) of 203K/204R has rose up from nearly zero to 76% to date with a shrinking from August to November in 2020 but bounced later. Our results show that the emergence of B.1.1.7 is associated with the second growth of R203K/G204R mutants. We identified positive selection evidences that support the adaptiveness of 203K/204R variants. The R203K/G204R mutant virus was created and compared with the native virus. The virus competition experiments show that 203K/204R variants possess a replication advantage over the preceding R203/G204 variants, possibly in relation to the ribonucleocapsid (RNP) assemble during the virus replication. Moreover, the 203K/204R virus increased the infectivity in a human lung cell line and induced an enhanced damage to blood vessel of infected hamsters lungs. In consistence, we observed a positive association between the increased severity of COVID-19 and the IF of 203K/204R from in silicon analysis of global clinical and epidemic data. In combination with the informatics and virology experiment, our work suggested the contribution of 203K/204R to the increased transmission and virulence of the SARS-CoV-2. In addition to mutations on the S protein, the mutations on the N protein are also important to virus spread during the pandemic.


title: Severe Acute Respiratory Syndrome Coronavirus-2 genome sequence variations relate to morbidity and mortality in Coronavirus Disease-19

abstract: Outcome of infection with Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) may depend on the host, virus or the host-virus interaction-related factors. Complete SARS-CoV-2 genome was sequenced using Illumina and Nanopore platforms from naso-/oro-pharyngeal ribonucleic acid (RNA) specimens from COVID-19 patients of varying severity and outcomes, including patients with mild upper respiratory symptoms (n=35), severe disease ad-mitted to intensive care with respiratory and gastrointestinal symptoms (n=21), fatal COVID-19 outcome (n=17) and asymptomatic (n=42). Of a number of genome variants observed, p.16L>L (Nsp1), p.39C>C (Nsp3), p.57Q>H (ORF3a), p.71Y>Y (Membrane glycoprotein), p.194S>L (Nucleocapsid protein) were observed in similar frequencies in different patient subgroups. However, seventeen other variants were observed only in symptomatic patients with severe and fatal COVID-19. Out of the latter, one was in the 5UTR (g.241C>T), eight were synonymous (p.14V>V and p.92L>L in Nsp1 protein, p.226D>D, p.253V>V, and p.305N>N in Nsp3, p.34G>G and p.79C>C in Nsp10 protein, p.789Y>Y in Spike protein), and eight were non-synonymous (p.106P>S, p.157V>F and p.159A>V in Nsp2, p.1197S>R and p.1198T>K in Nsp3, p.97A>V in RdRp, p.614D>G in Spike protein, p.13P>L in nucleocapsid). These were completely absent in the asymptomatic group. SARS-CoV-2 genome variations have a significant impact on COVID-19 presentation, severity and outcome.

category: genomics


title: Investigating the possible origin and transmission routes of SARS-CoV-2 genomes and variants of concern in Bangladesh

abstract: The COVID-19 pandemic induced by the SARS-CoV-2 virus and its variants has ravaged most countries around the world including Bangladesh. We have analyzed publicly available genomic data to understand the current COVID-19 outbreak scenario as well as the evolutionary origin and transmission routes of SARS-CoV-2 isolates in Bangladesh. All the early isolates as well as recent B.1.1.7 and B.1.351 variants had already spread across the major divisional cities of Bangladesh. A sex biasness towards male COVID-19 patient samples sequencing has observed over female in all age-group, that could be the trend in infection rate. Phylogenetic analysis indicated a total of 13 estimated countries, including Italy, India, United Kingdom, Saudi Arabia, United Arab Emirates, Germany, Australia, New Zealand, South Africa, Democratic Republic of the Congo, United States, Russia, and Denmark, could be the possible origin introduced SARS-CoV-2 isolates in Bangladesh due to regional and intercontinental travel. Recent, B.1.1.7 variant could be imported from a total of 7 estimated countries including UK, India, Nigeria, Spain, Ireland, Australia, and Indonesia, while South Africa and the United States are the most likely sources of B.1351 variant in Bangladesh. Based on these findings, public health strategies could be designed and implemented to reduce the local transmission of the virus.


title: Structure-function relationship for a divergent Atg8 protein required for a non-autophagic function in malaria parasites

abstract:
Atg8 family proteins are highly-conserved eukaryotic proteins with diverse autophagy and non-autophagic functions in eukaryotes. While the structural features required for conserved autophagy functions of Atg8 are well-established, little is known about the molecular changes that facilitated acquisition of divergent, non-autophagic functions of Atg8. The malaria parasite Plasmodium falciparum offers a unique opportunity to study non-autophagic functions of Atg8 family proteins because it encodes a single Atg8 homolog whose only essential function is in the inheritance of an unusual secondary plastid called the apicoplast. Here we used functional complementation to investigate the structure-function relationship for this divergent Atg8 protein. We showed that the LC3-interacting region (LIR) docking site (LDS), the major interaction interface of Atg8 protein family, is not sufficient for PfAtg8 apicoplast function. Other regions previously implicated in canonical Atg8 interactions, the ubiquitin-interacting motif (UIM) docking site (UDS) and the N-terminal helix are not required for PfAtg8 function. Finally, the unique Apicomplexan-specific loop previously implicated in interaction with membrane conjugation machinery in vitro, is not required in vivo neither for membrane conjugation nor for the effector function of PfAtg8. These results suggest that the effector function of PfAtg8 is mediated by structural features distinct from those previously identified for macroautophagy and selective autophagy functions.

ImportanceThe most extensively studied role of Atg8 proteins is in autophagy. However, it is clear that they have other non-autophagic functions critical to cell function and disease pathogenesis yet understudied compared to their canonical role in autophagy. Mammalian cells contain multiple Atg8 paralogs that have diverse, specialized functions. Gaining molecular insight into their non-autophagic functions is difficult because of redundancy between the homologs and their role in both autophagy and non-autophagic pathways. Malaria parasites such as Plasmodium falciparum are a unique system to study a novel, non-autophagic function of Atg8 separate from its role in autophagy: They have only one Atg8 protein whose only essential function is in the inheritance of the apicoplast, a unique secondary plastid organelle. Insights into the molecular basis of PfAtg8s function in apicoplast biogenesis will have important implications for the evolution of diverse non-autophagic functions of the Atg8 protein family.

category: biochemistry


title: Controllable drug release of pH-sensitive liposomes encapsulating artificial cytosol system

abstract: The fabrication of cell model containing artificial cytosol is challenging. Herein we constructed an artificial cytosol contained cell model by electroformation method. Agarose was selected as the main component of the artificial cytosol. Sucrose was added into agarose to regulate the sol viscosity and phase transition temperature (Tm). The viscosity of the sol with the mass ratio (agarose-sucrose) 1:9 was closest to the natural cytosolic. DSPC/20 mol% Chol was used to form large unilamellar vesicle (LUV) as cell model compartment. The rhodamine release experiment confirmed that the release of rhodamine from LUVs containing artificial cytosol took more time than that from LUVs containing pure water. The unique release profile makes agarose-sucrose@LUVs suitable as a drug carrier. Doxorubicin (DOX) is loaded in the agarose-sucrose@LUVs, and their half maximum inhibition concentration on HeLa cells is 0.015 mol L-1, which means 31.7 times increase in inhibition efficiency over free DOX.


title: Genome-wide identification and prediction of SARS-CoV-2 mutations show an abundance of variants: Integrated study of bioinformatics and deep neural learning.

abstract: Genomic data analysis is a fundamental system for monitoring pathogen evolution and the outbreak of infectious diseases. Based on bioinformatics and deep learning, this study was designed to identify the genomic variability of SARS-CoV-2 worldwide and predict the impending mutation rate. Analysis of 259044 SARS-CoV-2 isolates identify 3334545 mutations (14.01 mutations per isolate), suggesting a high mutation rate. Strains from India showed the highest no. of mutations (48) followed by Scotland, USA, Netherlands, Norway, and France having up to 36 mutations. Besides the most prominently occurring mutations (D416G, F106F, P314L, and UTR:C241T), we identify L93L, A222V, A199A, V30L, and A220V mutations which are in the top 10 most frequent mutations. Multi-nucleotide mutations GGG>AAC, CC>TT, TG>CA, and AT>TA have come up in our analysis which are in the top 20 mutational cohort. Future mutation rate analysis predicts a 17%, 7%, and 3% increment of C>T, A>G, and A>T, respectively in the future. Conversely, 7%, 7%, and 6% decrement is estimated for T>C, G>A, and G>T mutations, respectively. T>G\A, C>G\A, and A>T\C are not anticipated in the future. Since SARS-CoV-2 is evolving continuously, our findings will facilitate the tracking of mutations and help to map the progression of the COVID-19 intensity worldwide.


title: Protein structure prediction using AI and quantum computers

abstract: This work seeks to combine the combined advantage of leveraging these emerging areas of Artificial Intelligence and quantum computing in applying it to solve the specific biological problem of protein structure prediction using Quantum Machine Learning algorithms. The CASP dataset from ProteinNet was downloaded which is a standardized data set for machine learning of protein structure. Its large and standardized dataset of PDB entries contains the coordinates of the backbone atoms, corresponding to the sequential chain of N, C_alpha, and C atoms. This dataset was used to train a quantum-classical hybrid Keras deep neural network model to predict the structure of the proteins. To visually qualify the quality of the predicted versus the actual protein structure, protein contact maps were generated with the experimental and predicted protein structure data and qualified. Therefore this model is recommended for the use of protein structure prediction using AI leveraging the power of quantum computers.


title: Entry Pathway for the Inverse Agonist Ligand in the G Protein-Coupled Receptor Rhodopsin

abstract: While the number of high-resolution structures of ligand-bound G protein-coupled receptors (GPCRs) has been steadily climbing, ligand binding and unbinding pathways remain largely undefined. The visual photoreceptor rhodopsin (Rho) represents a curious case among GPCRs because its primary ligand 11-cis-retinal (11CR) is an inverse agonist, which partitions into the bilayer and is likely to enter its orthosteric binding pocket through an intermembranous pathway. Light activates Rho by converting 11CR to all-trans-retinal (ATR), which serves as an agonist ligand. The light-triggered switch from the inactive to the active conformation creates two openings in the transmembrane region, suggesting pathways for exit of ATR and subsequent entry of 11CR to regenerate Rho. However, stabilization of an active ligand-free opsin conformation has been found to inhibit 11CR binding. Here we address this paradox of opsin regeneration with 11CR. We used genetic code expansion to engineer Rho mutants that serve as fluorescence sensors for measuring 11CR binding kinetics and energetics. We found mutations that alter a channel between transmembrane helices 5 and 6 (TM5/6) dramatically affect 11CR binding kinetics, but not ATR release kinetics. Our data provide direct experimental evidence for 11CR entry between TM5/6 in Rho that involves dynamic allosteric control of the ligand entry channel. Our findings can be extended to other visual pigments and a wide range of GPCRs with hydrophobic ligands that are hypothesized to enter their binding pockets through transmembrane pores.


title: Transcriptional overlap links DNA hypomethylation with DNA hypermethylation at adjacent promoters in cancer

abstract: DNA methylation is an epigenetic mark associated with gene repression. It is now well established that tumor development involves alterations in DNA methylation patterns, which include both gains (hypermethylation) and losses (hypomethylation) of methylation marks in different genomic regions. The mechanisms underlying these two opposite, yet co-existing, alterations in tumors remain unclear. While studying the human MAGEA6/GABRA3 gene locus, we observed that DNA hypomethylation in tumor cells can lead to the activation of a long transcript (CT-GABRA3) that overlaps downstream promoters (GABRQ and GABRA3) and triggers their hypermethylation. Overlapped promoters displayed increases in H3K36me3, a histone mark known to be deposited during progression of the transcription machinery and to stimulate de novo DNA methylation. Consistent with such a processive mechanism, increases in H3K36me3 and DNA methylation were observed over the entire region covered by the CT-GABRA3 overlapping transcript. Importantly, experimental induction of CT-GABRA3 by depletion of DNMT1 DNA methyltransferase, resulted in a similar pattern of increased DNA methylation in the MAGEA6/GABRA3 locus. Bioinformatics analyses in lung cancer datasets identified other genomic loci displaying this process of coupled DNA hypo- and hypermethylation. In several of these loci, DNA hypermethylation affected tumor suppressor genes, e.g. RERG and PTPRO. Together, our work reveals that focal DNA hypomethylation in tumors can indirectly contribute to hypermethylation of nearby promoters through activation of overlapping transcription, and establishes therefore an unsuspected connection between these two opposite epigenetic alterations.

category: cancer biology


title: A UAV-based high-throughput phenotyping approach to assess time-series nitrogen responses and identify traits associated genetic components in maize

abstract: Advancements in the use of genome-wide markers have provided new opportunities for dissecting the genetic components that control phenotypic trait variation. However, cost-effectively characterizing agronomically important phenotypic traits on a large scale remains a bottleneck. Unmanned aerial vehicle (UAV)-based high-throughput phenotyping has recently become a prominent method, as it allows large numbers of plants to be analyzed in a time-series manner. In this experiment, 233 inbred lines from the maize diversity panel were grown in a replicated incomplete block under both nitrogen-limited conditions and following conventional agronomic practices. UAV images were collected during different plant developmental stages throughout the growing season. A pipeline for extracting plot-level images, filtering images to remove non-foliage elements, and calculating canopy coverage and greenness ratings based on vegetation indices (VIs) was developed. After applying the pipeline, about half a million plot-level image clips were obtained for 12 different time points. High correlations were detected between VIs and ground truth physiological and yield-related traits collected from the same plots, i.e., Vegetative Index (VEG) vs. leaf nitrogen levels (Pearson correlation coefficient, R = 0.73), Woebbecke index vs. leaf area (R = -0.52), and Visible Atmospherically Resistant Index (VARI) vs. 20 kernel weight - a yield component trait (R = 0.40). The genome-wide association study was performed using canopy coverage and each of the VIs at each date, resulting in N = 29 unique genomic regions associated with image extracted traits from three or more of the 12 total time points. A candidate gene Zm00001d031997, a maize homolog of the Arabidopsis HCF244 (high chlorophyll fluorescence 244), located underneath the leading SNPs of the canopy coverage associated signals that were repeatedly detected under both nitrogen conditions. The plot-level time-series phenotypic data and the trait-associated genes provide great opportunities to advance plant science and to facilitate plant breeding.

category: plant biology


title: Reprogrammed tracrRNAs enable repurposing RNAs as crRNAs and detecting RNAs

abstract: In type II CRISPR systems, the guide RNA (gRNA) consists of a CRISPR RNA (crRNA) and a hybridized trans-acting CRISPR RNA (tracrRNA) which interacts directly with Cas9 and is essential to its guided DNA targeting function. Though tracrRNAs are diverse in sequences and structures across type II CRISPR systems, the programmability of crRNA-tracrRNA hybridization for particular Cas9 has not been studied adequately. Here, we revealed the high programmability of crRNA-tracrRNA hybridization for Streptococcus pyogenes Cas9. By reprogramming the crRNA-tracrRNA hybridized sequence, reprogrammed tracrRNAs can repurpose various RNAs as crRNAs to trigger CRISPR function. We showed that the engineered crRNA-tracrRNA pairs enable design of orthogonal cellular computing devices and hijacking of endogenous RNAs as crRNAs. We next designed novel RNA sensors that can monitor the transcriptional activity of specific genes on the host genome and detect SARS-CoV-2 RNA in vitro. The engineering potential of crRNA-tracrRNA interaction has therefore redefined the capabilities of CRISPR/Cas9 system.

category: synthetic biology


title: Adaptation to chronic ER stress enforces pancreatic β-cell plasticity

abstract: Pancreatic {beta}-cells undergo high levels of endoplasmic reticulum (ER) stress due to their role in insulin secretion. Hence, they require sustainable and efficient adaptive stress responses to cope with the stress. Whether duration and episodes of chronic ER stress directly compromises {beta}-cell identity is largely unknown. We show that under reversible, chronic ER stress, {beta}-cells undergo a distinct transcriptional and translational reprogramming. During reprogramming, expression of master regulators of {beta}-cell function and identity and proinsulin processing is impaired. Upon recovery from stress, {beta}-cells regain their identity, highlighting a high-degree of adaptive {beta}-cell plasticity. Remarkably, when stress episodes exceed a certain threshold, {beta}-cell identity is gradually lost. Single cell RNA-seq analysis of islets from type 1 diabetes (T1D) patients, identifies the severe deregulation of the chronic stress-adaptation program, and reveals novel biomarkers for progression of T1D. Our results suggest {beta}-cell adaptive exhaustion ({beta}EAR) is a significant component of the pathogenesis of T1D.

category: systems biology


title: Non-neutralizing antibodies targeting the immunogenic regions of HIV-1 envelope reduce mucosal infection and virus burden in humanized mice

abstract: Antibodies are principal immune components elicited by vaccines to induce protection from microbial pathogens. In the Thai RV144 HIV-1 vaccine trial, vaccine efficacy was 31% and the sole primary correlate of reduced risk was shown to be vigorous antibody response targeting the V1V2 region of HIV-1 envelope. Antibodies against V3 also were inversely correlated with infection risk in subsets of vaccinees. Antibodies recognizing these regions, however, do not exhibit potent neutralizing activity. Therefore, we examined the antiviral potential of poorly neutralizing monoclonal antibodies (mAbs) against immunodominant V1V2 and V3 sites by passive administration of human mAbs to humanized mice engrafted with CD34+ hematopoietic stem cells, followed by mucosal challenge with an HIV-1 infectious molecular clone (IMC) expressing the envelope of a tier 2 resistant HIV-1 strain. Treatment with anti-V1V2 mAb 2158 or anti-V3 mAb 2219 did not prevent infection, but both reduced the virus burden, and V3 mAb 2219 displayed a superior potency compared to V1V2 mAb 2158. While these mAbs had no or weak neutralizing activity and elicited undetectable levels of antibody-dependent cellular cytotoxicity (ADCC), V3 mAb 2219 displayed a greater capacity to bind virus- and cell-associated HIV-1 envelope and to mediate antibody-dependent cellular phagocytosis (ADCP) and C1q complement binding as compared to V1V2 mAb 2158. Mutations in the Fc region of 2219 abolished these effector activities and abrogated virus control in humanized mice. These results demonstrate the importance of Fc functions other than ADCC for antibodies without potent neutralizing activity.


title: HIV-1 Nef and CycK:CDK13 antagonize SERINC5 for optimal viral infectivity

abstract: HIV-1 Nef antagonizes SERINC5 by redirecting this potent restriction factor to the endosomes and lysosomes for degradation. However, the precise mechanism remains unclear. Using affinity purification/mass spectrometry, we identified cyclin K and cyclin-dependent kinase 13 (CycK:CDK13) as a new Nef-associated kinase complex. CycK:CDK13 phosphorylates the serine at position 360 (S360) in SERINC5, which is required for Nef downregulation of SERINC5 from the cell surface and its counter activity of the SERINC5 antiviral activity. To understand the role of S360 phosphorylation, we created chimeric proteins between CD8 and SERINC5. Nef not only downregulates, but importantly, also binds to this chimera in a S360-dependent manner. Thus, S360 phosphorylation increases interactions between Nef and SERINC5 and initiates the destruction of SERINC5 by the endocytic machinery.


title: Origins of Scaling Laws in Microbial Dynamics

abstract: Analysis of high-resolution time series data from the human and mouse gut microbiomes revealed that the gut microbial dynamics can be characterized by several robust and simple scaling laws. It is still unknown if those scaling laws are universal across different body sites, host species, or even free-living microbial communities. Moreover, the underlying mechanisms responsible for those scaling laws remain poorly understood. Here, we demonstrate that those scaling laws are not unique to gut microbiome, but universal across different habitats, from human skin and oral microbiome to marine plankton bacteria and eukarya communities. Since completely shuffled time series yield very similar scaling laws, we conjecture that the universal scaling laws in various microbiomes are largely driven by temporal stochasticity of the host or environmental factors. We leverage a simple population dynamics model with both deterministic inter-species interactions and stochastic noise to confirm our conjecture. In particular, we find that those scaling laws are jointly determined by inter-species interactions and linear multiplicative noises. The presented results deepen our understanding of the nature of scaling laws in microbial dynamics.


title: The mechanism underlying the organization of Borna disease virus inclusion bodies is unique among mononegaviruses

abstract: Inclusion bodies (IBs) are characteristic biomolecular condensates organized by mononegaviruses. Here, we characterize the IBs of Borna disease virus 1 (BoDV-1), a unique mononegavirus that forms IBs in the nucleus, in terms of liquid-liquid phase separation (LLPS). The BoDV-1 phosphoprotein (P) alone induces LLPS and the nucleoprotein (N) is incorporated into the P droplet in vitro. In contrast, co-expression of N and P is required for the formation of IB-like structure in cells. Furthermore, while BoDV-1 P binds to RNA, an excess amount of RNA dissolves the liquid droplets formed by N and P. Notably, the N-terminal intrinsically disordered region of BoDV-1 P is essential to drive LLPS and bind to RNA, suggesting that both abilities could compete with one another. These features are unique among mononegaviruses, and thus this study will contribute to a deeper understanding of LLPS-driven organization and RNA-mediated regulation of biomolecular condensates.


title: CZON-cutter: a CRISPR-Cas9 system with multiplexed organelle imaging in a simple unicellular alga

abstract:
The simple cellular structure of the unicellular alga Cyanidioschyzon merolae consists of one nucleus, one mitochondrion, one chloroplast, and one peroxisome per cell and offers unique advantages to investigate mechanisms of organellar proliferation and the cell cycle. Here, we describe an engineered clustered, regularly interspaced, short palindromic repeats (CRISPR)-associated protein 9 (Cas9) system, CZON-cutter, for simultaneous genome editing and organellar visualization. We engineered a C. merolae strain expressing a nuclear-localized Cas9-Venus nuclease to target editing at a locus defined by a single-guide RNA (sgRNA). We then successfully edited the algal genome and visualized the mitochondrion and peroxisome in transformants by fluorescent protein reporters with different excitation wavelengths. Fluorescent protein labeling of organelles in living transformants allows validation of phenotypes associated with organellar proliferation and the cell cycle, even when the edited gene is essential. Combined with the exceptional biological features of C. merolae, CZON-cutter will be instrumental for investigating cellular and organellar division in a high-throughput manner.

SummaryAn engineered CRISPR-Cas9 system, named CZON-cutter, for simultaneous genome editing and fluorescent protein labeling of organelles in Cyanidioschyzon merolae can be used to validate intracellular function of a particular gene, even if it is essential.


title: Rab3a regulates melanin exocytosis induced by keratinocyte-conditioned medium

abstract: Skin pigmentation relies on melanin and is crucial for photoprotection against ultraviolet radiation-induced toxicity. Melanin is synthesized and stored in melanosomes, within melanocytes and then transferred to keratinocytes. While the molecular players involved in melanogenesis have been extensively studied, those underlying melanin transfer remain poorly characterized. Previously, our group proposed that coupled exo/phagocytosis is the predominant mechanism of melanin transfer in human skin and showed an essential role for Rab11b and the exocyst tethering complex in this process. Using a fluorescence-based assay, we show here that keratinocyte-conditioned medium (KCM) specifically induces melanin exocytosis from melanocytes. Moreover, we found that Rab3a, but not Rab11b, regulates melanin exocytosis upon KCM stimulation. In fact, melanosomes accumulate in melanocyte dendrites upon KCM stimulation, co-localizing with Rab3a mainly in the vicinity of the plasma membrane. Additionally, Rab3a silencing does not affect melanin transfer in melanocyte/keratinocyte co-cultures, in contrast with Rab11b depletion, indicating that Rab11b regulates non-KCM-stimulated melanin exocytosis. Thus, our results suggest the existence of at least two distinct routes of melanin exocytosis: one controlled by Rab11b and another Rab3a-dependent, stimulated by KCM. Furthermore, these results provide evidence that soluble factors secreted by keratinocytes can control skin pigmentation via induction of melanocyte signaling pathways that promote peripheral transport of melanosomes and a Rab3a-mediated exocytosis mechanism.


title: The RNA chaperone protein CspA stimulates translation during cold acclimation by promoting the progression of the ribosomes

abstract: CspA is an RNA binding protein expressed during cold-shock in Escherichia coli, capable of stimulating translation of several mRNAs - including its own - at low temperature. We used reconstituted translation systems to monitor the effects of CspA on the different steps of the translation process and probing experiments to analyze the interactions with its target mRNAs. We specifically focused on cspA mRNA which adopts a cold-induced secondary structure at temperatures below 20{degrees}C and a more closed conformation at 37{degrees}C. We show that at low temperature CspA specifically promotes the translation of the mRNA folded in the conformation less accessible to the ribosome (37{degrees}C form). CspA interacts with its mRNA without inducing large structural rearrangement, does not bind the ribosomal subunits and is not able to stimulate the formation of the translation initiation complexes. On the other hand, CspA promotes the progression of the ribosomes during translation of its mRNA at low temperature and this stimulation is mRNA structure-dependent. A similar structure-dependent mechanism may be responsible for the CspA- dependent translation stimulation observed with other probed mRNAs, for which the transition to the elongation phase is progressively facilitated during cold acclimation with the accumulation of CspA.


title: Island-specific evolution of a sex-primed autosome in the planarian Schmidtea mediterranea

abstract: The sexual biotype of the planarian Schmidtea mediterranea is a hermaphrodite indigenous to Tunisia and several Mediterranean islands. Here, we isolated individual chromosomes and used sequencing, Hi-C and linkage mapping to assemble a chromosome-scale genome reference. The linkage map revealed an extremely low rate of recombination on chromosome 1. We confirmed suppression of recombination on chromosome 1 by quantifying recombination events in individual sperm and oocytes. We showed that the extensive heterozygous regions, previously designated as J and V haplotypes, comprise essentially all of chromosome 1. Genome sequencing of individuals isolated in the wild indicated that this heterozygosity has evolved specifically in populations from Sardinia and Corsica. We found that chromosome 1 acquired many genes that determine the development of female and male reproductive systems, along with haplotype-specific expression of such sex specific genes. These molecular signatures and restricted recombination in a hermaphrodite led us to propose that chromosome 1 is a sex-primed autosome, providing direct molecular evidence for the traditional model of how sex chromosomes may have evolved from autosomes.

category: genetics


title: CREB-mediated transcriptional activation of NRMT1 drives muscle differentiation

abstract: The N-terminal methyltransferase NRMT1 is an important regulator of protein-DNA interactions and plays a role in many cellular processes, including mitosis, cell cycle progression, chromatin organization, DNA damage repair, and transcriptional regulation. Accordingly, loss of NRMT1 results in both developmental pathologies and oncogenic phenotypes. Though NRMT1 plays such important and diverse roles in the cell, little is known about its own regulation. To better understand the mechanisms governing NRMT1 expression, we first identified its predominant transcriptional start site and minimal promoter region with predicted transcription factor motifs. We then used a combination of luciferase and binding assays to confirm CREB1 as the major regulator of NRMT1 transcription. We tested which conditions known to activate CREB1 also activated NRMT1 transcription, and found CREB1-mediated NRMT1 expression was increased during recovery from serum starvation and muscle cell differentiation. To determine how NRMT1 expression affects myoblast differentiation, we used CRISPR/Cas9 technology to knock out NRMT1 expression in immortalized C2C12 mouse myoblasts. C2C12 cells depleted of NRMT1 lacked Pax7 expression and were unable to proceed down the muscle differentiation pathway. Instead, they took on characteristics of C2C12 cells that have transdifferentiated into osteoblasts, including increased alkaline phosphatase and type I collagen expression and decreased proliferation. These data implicate NRMT1 as an important downstream target of CREB1 during muscle cell differentiation.


title: Dimerization of SARS-CoV-2 nucleocapsid protein affects sensitivity of ELISA based diagnostics of COVID-19

abstract:
Diagnostics has played a significant role in effective management of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Nucleocapsid protein (N protein) is the primary antigen of the virus for development of sensitive diagnostic assays. Thus far, limited knowledge exists about the antigenic properties of the N protein. In this paper, we demonstrate the significant impact of dimerization of SARS-CoV-2 nucleocapsid protein on sensitivity of enzyme-linked immunosorbent assay (ELISA) based diagnostics of COVID-19. The expressed purified protein from E.coli consists of two forms, dimeric and monomeric forms, which have been further characterized by biophysical and immunological means. Indirect ELISA indicated elevated susceptibility of the dimeric form of the nucleocapsid protein for identification of protein-specific monoclonal antibody as compared to the monomeric form of the protein. These findings have also been confirmed with the modelled structure of monomeric and dimeric nucleocapsid protein via HHPred software and its solvent accessible surface area, which indicates higher stability and antigenicity of the dimeric type as compared to the monomeric form. It is evident that use of the dimeric form will increase the sensitivity of the current nucleocapsid dependent ELISA for rapid COVID-19 diagnostic. Further, the results indicate that monitoring and maintaining of the monomerdimer composition is critical for accurate and robust diagnostics.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/445305v2_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@16bf153org.highwire.dtl.DTLVardef@1b0be6corg.highwire.dtl.DTLVardef@a5b3aeorg.highwire.dtl.DTLVardef@b784a3_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: bioengineering


title: SARS-CoV-2 spike protein unlikely to bind to integrins via the Arg-Gly-Asp (RGD) motif of the Receptor Binding Domain: evidence from structural analysis and microscale accelerated molecular dynamics

abstract: The Receptor Binding Domain (RBD) of SARS-CoV-2 virus harbors a sequence of Arg-Gly-Asp tripeptide named RGD motif, which has also been identified in extracellular matrix proteins that bind integrins as well as other disintegrins and viruses. Accordingly, integrins have been proposed as host receptors for SARS-CoV-2. The hypothesis was supported by sequence and structural analysis. However, given that the microenvironment of the RGD motif imposes structural hindrance to the protein-protein association, the validity of this hypothesis is still uncertain. Here, we used normal mode analysis, accelerated molecular dynamics microscale simulation, and protein-protein docking to investigate the putative role of RGD motif of SARS-CoV-2 RBD for interacting with integrins. We found, by molecular dynamics, that neither RGD motif nore its microenvironment show any significant conformational shift in the RBD structure. Highly populated clusters were used to run a protein-protein docking against three RGD-binding integrin types, showing no capability of the RBD domain to interact with the RGD binding site. Moreover, the free energy landscape revealed that the RGD conformation within RBD could not acquire an optimal geometry to allow the interaction with integrins. Our results highlighted different structural features of the RGD motif that may prevent its involvement in the interaction with integrins. We, therefore, suggest, in the case where integrins are confirmed to be the direct host receptors for SARS-CoV-2, a possible involvement of other residues to stabilize the interaction.


title: Natural variation in a short region of the Acidovorax citrulli type III-secreted effector AopW1 is associated with differences in cytotoxicity

abstract: Bacterial fruit blotch (BFB) is a serious disease of melon and watermelon caused by the Gram-negative bacterium Acidovorax citrulli. The strains of the pathogen can be divided into two major genetic groups, I and II. While group I strains are strongly associated with melon, group II strains are more aggressive on watermelon. Like many pathogenic bacteria, A. citrulli secretes a variety of protein effectors to the host cell via the type III secretion system. In the present study, we characterized AopW1, an A. citrulli type III-secreted effector that shares similarity with the actin cytoskeletondisrupting effector HopW1 of Pseudomonas syringae and with effectors from other plant-pathogenic bacterial species. aopW1 is present in group I and II strains, encoding products of 485 amino acids. Although highly conserved in most of the sequence, AopW1 has a highly variable region (HVR) within amino acid positions 147 to 192, including 14 amino acid differences between groups. Here we show that group I AopW1 is more toxic to yeast and plant cells than group II AopW1, having a stronger actin filament disruption activity, and increased ability to reduce plant callose deposition. We demonstrate the role of some of these 14 amino acid positions in determining the phenotypic differences between the two versions of the effector. Moreover, cellular analyses revealed that in addition to the interaction with actin filaments, AopW1 is localized to the endoplasmic reticulum, chloroplasts, and early and recycling plant endosomes, with differences observed between the two AopW1 versions. Finally, we show that overexpression of the endosome-associated protein EHD1 that increases cellular recycling, attenuates the toxic effects exerted by AopW1 and increases defence responses. This study provides insights into the HopW1 family of bacterial effectors and their interactions with the plant cell and provides first evidence on the involvement of EHD1 in response to biotic stress.


title: The Nicotiana tabacum Mediator subunit MED25 regulates nicotine biosynthesis through interacting with the basic helix-loop-helix (bHLH) transcription factor NtMYC2s

abstract: Nicotine is one of the most important secondary metabolites in tobacco, and its biosynthesis can be induced by topping and jasmonic acid treatment. NtMYC2s play pivotal roles in the regulation of nicotine. The mediator server as a bridge betwen the transcription factors and RNA polymerase in order to facilitates transcription and functions in plants. However, the role of mediator in the regulation of nicotine biosynthesis remains unknown. In this study, we firstly identify the NtMED25 through homologous analysis. NtMED25 interacts with NtMYC2s through the MD region. Interestingly, the nicotine content is decreased in the the knock-down transgenic lines of NtMED25, and the expression levels of two nicotine biosynthesis genes, NtQPT2 and NtPMT2, are also reduced when compared with that in the wild-type plants. Furthermore, NtMED25 enhances the binding of NtMYC2a/ NtMYC2b to the promoter of NtPMT2 and NtQPT2, and then facilitates the nicotine biosynthesis. Therefore, our study revealed the function of mediator in the regulation of nicotine, and provide the insight role on the transcriptional regulation of plant secondary metabolites.


title: Three classes of epigenomic regulators converge to hyperactivate the essential maternal gene deadhead within a heterochromatin mini-domain.

abstract:
The formation of a diploid zygote is a highly complex cellular process that is entirely controlled by maternal gene products stored in the egg cytoplasm. This highly specialized transcriptional program is tightly controlled at the chromatin level in the female germline. As an extreme case in point, the massive and specific ovarian expression of the essential thioredoxin Deadhead (DHD) is critically regulated in Drosophila by the histone demethylase Lid and its partner, the histone deacetylase complex scaffold Sin3A, via yet unknown mechanisms. Here, we identified the Brahma chromatin remodeler sub-unit Snr1 and the insulator component Mod(mdg4) as essential for dhd expression and investigated how these epigenomic effectors act with Lid and Sin3A to hyperactivate dhd. Using Cut&Run chromatin profiling with a dedicated data analysis procedure, we found that dhd is intriguingly embedded in an H3K27me3/H3K9me3-enriched mini-domain flanked by DNA regulatory elements, including a dhd promoter-proximal element essential for its expression. Surprisingly, Lid, Sin3A, Snr1 and Mod(mdg4) impact H3K27me3 and this regulatory element in distinct manners. However, we show that these effectors activate dhd independently of H3K27me3/H3K9me3, and that these marks are not required to repress dhd. Together, our study demonstrates an atypical and critical role for chromatin regulators Lid, Sin3A, Snr1 and Mod(mdg4) to trigger tissue-specific hyperactivation within a unique heterochromatin mini-domain.

AUTHOR SUMMARYGene expression is tightly regulated by conserved protein complexes that act at the chromatin level to allow or restrict transcription. Such epigenetic control of gene activity defines the identity of different cell types during development, as well as their response to environmental cues. Yet, how multiple chromatin factors converge to achieve precise gene regulation remains difficult to address, partly due to the lack of biological situations where these intricate relationships can be studied. In this paper, we have addressed this issue by dissecting the regulation of deadhead, an essential gene specifically and massively expressed in the Drosophila germline. Unexpectedly, we found that its hyperactivation occurs despite deadhead being embedded in an apparently unfavorable chromatin mini-domain, notably featuring repressive histone modifications. We further demonstrate that four chromatin effectors, Lid, Sin3A, Snr1 and Mod(mdg4), have distinct, atypical and essential roles to ensure deadhead expression within this chromatin environment. Together, our findings put into perspective our understanding on these regulatory factors by illustrating how they can exert a biologically essential function via non-canonical mechanisms.


title: Pathway analysis in metabolomics: pitfalls and best practice for the use of over-representation analysis

abstract:
Over-representation analysis (ORA) is one of the commonest pathway analysis approaches used for the functional interpretation of metabolomics datasets. Despite the widespread use of ORA in metabolomics, the community lacks guidelines detailing its best-practice use. Many factors have a pronounced impact on the results, but to date their effects have received little systematic attention in the field. We developed in-silico simulations using five publicly available datasets and illustrated that changes in parameters, such as the background set, differential metabolite selection methods, and pathway database choice, could all lead to profoundly different ORA results. The use of a non-assay-specific background set, for example, resulted in large numbers of false-positive pathways. Pathway database choice, evaluated using three of the most popular metabolic pathway databases: KEGG, Reactome, and BioCyc, led to vastly different results in both the number and function of significantly enriched pathways. Metabolomics data specific factors, such as reliability of compound identification and assay chemical bias also impacted ORA results. Simulated metabolite misidentification rates as low as 4% resulted in both gain of false-positive pathways and loss of truly significant pathways across all datasets. Our results have several practical implications for ORA users, as well as those using alternative pathway analysis methods. We offer a set of recommendations for the use of ORA in metabolomics, alongside a set of minimal reporting guidelines, as a first step towards the standardisation of pathway analysis in metabolomics.

Author summaryMetabolomics is a rapidly growing field of study involving the profiling of small molecules within an organism. It allows researchers to understand the effects of biological status (such as health or disease) on cellular biochemistry, and has wide-ranging applications, from biomarker discovery and personalised medicine in healthcare to crop protection and food security in agriculture. Pathway analysis helps to understand which biological pathways, representing collections of molecules performing a particular function, are involved in response to a disease phenotype, or drug treatment, for example. Over-representation analysis (ORA) is perhaps the most common pathway analysis method used in the metabolomics community. However, ORA can give drastically different results depending on the input data and parameters used. In this work, we have established the effects of these factors on ORA results using computational simulations applied to five real-world datasets. Based on our results, we offer the research community a set of best-practice recommendations applicable not only to ORA but also to other pathway analysis methods to help ensure the reliability and reproducibility of results.


title: RNA secondary structure prediction with Convolutional Neural Networks

abstract:
Predicting the secondary, i.e. base-pairing structure of a folded RNA strand is an important problem in synthetic and computational biology. First-principle algorithmic approaches to this task are challenging because existing models of the folding process are inaccurate, and even if a perfect model existed, finding an optimal solution would be in general NP-complete. In this paper, we propose a simple, yet extremely effective data-driven approach. We represent RNA sequences in the form of three-dimensional tensors in which we encode possible relations between all pairs of bases in a given sequence. We then use a convolutional neural network to predict a two-dimensional map which represents the correct pairings between the bases. Our model achieves significant accuracy improvements over existing methods on two standard datasets. Our experiments show excellent performance of the model across a wide range of sequence lengths and RNA families. We also observe considerable improvements in predicting complex pseudoknotted RNA structures, as compared to previous approaches.

Author summaryStructure prediction for RNA sequences is a computationally difficult task that is of increasing importance in applications such as medical diagnostics and drug design; this is because the structure of a folded RNA strand to a large extent defines its function. An open RNA strand can fold to many different structures of varying thermal stability, and the goal of structure prediction is to determine a most stable one among these. There are two main difficulties to this task. Firstly, a given RNA sequence can fold into an enormous number of alternative structures, and a computational search for a most stable one in this huge space can be very demanding. The search can however be facilitated by using heuristics that take into account some underlying principles of the folding process. Here is where machine learning methods come into play: they are suitable for discovering patterns in data, and can thus predict features of the desired structure based on previously learned patterns. Secondly, there do not yet exist fully satisfactory coarse-grained models for the most popular metric for stability, the free energy of the folded structure. Although in principle a minimum free energy (MFE) structure should be a good candidate for a most stable one, MFE structures determined according to current energy models do not match experimental data on native RNA conformations very well. We show how to use an artificial neural network design to predict the structure for a given RNA sequence with high accuracy only by learning from samples whose native structures have been experimentally characterized, independent of any stability metric or energy model.


title: Insights into substrate recognition and specificity for IgG by Endoglycosidase S2

abstract:
Antibodies bind foreign antigens with high affinity and specificity leading to their neutralization and/or clearance by the immune system. The conserved N-glycan on IgG has significant impact on antibody effector function, with the endoglycosidases of Streptococcus pyogenes deglycosylating the IgG to evade the immune system, a process catalyzed by the endoglycosidase EndoS2. Studies have shown that two of the four domains of EndoS2, the carbohydrate binding module (CBM) and the glycoside hydrolase (GH) domain site are critical for catalytic activity. To yield structural insights into contributions of the CBM and the GH domains, as well as the overall flexibility of EndoS2, to the proteins catalytic activity, models of EndoS2-Fc complexes were generated through enhanced-sampling molecular-dynamics (MD) simulation and site-identification by ligand competitive saturation (SILCS) docking followed by reconstruction and long-time MD simulations. Modeling results predict that EndoS2 initially interacts with the IgG through its CBM followed by interactions with the GH yielding catalytically competent states. These may involve the CBM and GH of EndoS2 simultaneously interacting with either the same Fc CH2/CH3 domain or individually with the two Fc CH2/CH3 domains, with EndoS2 predicted to assume closed conformations in the former case and open conformations in the latter. Apo EndoS2 is predicted to sample both the open and closed states, suggesting that either complex can directly form following initial IgG-EndoS2 encounter. Interactions of the CBM and GH domains with the IgG are predicted to occur through both its glycan and protein regions. Simulations also predict that the Fc glycan can directly transfer from the CBM to the GH, facilitating formation of catalytically competent complexes. The predicted models are compared and consistent with Hydrogen/Deuterium Exchange data. In addition, the complex models are consistent with the high specificity of EndoS2 for the glycans on IgG supporting the validity of the predicted models.

Author SummaryThe pathogen Streptococcus pyogenes uses the endoglycosidases S and S2 to cleave the glycans on the Fc portion of IgG antibodies, leading to a decreased cytotoxicity of the antibodies, thereby evading the host immune response. To identify potential structures of the complex of EndoS2 with IgG that could lead to the catalytic hydrolysis of the IgG glycan, molecular modeling and molecular dynamics simulations were applied. The resulting structural models predict that EndoS2 initially interacts through its carbohydrate binding module (CBM) with the IgG with subsequent interactions with the catalytic glycoside hydrolase (GH) domain yielding stable complexes. In the modeled complexes the CBM and the GH interact either simultaneously with the same Fc CH2/CH3 domain or with the two individual Fc CH2/CH3 domains separately to yield potentially catalytically competent species. In addition, apo EndoS2 is shown to assume both open and closed conformations allowing it to directly form either type of complex from which deglycosylation of either mono- or diglycosylated IgG species may occur.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=99 SRC="FIGDIR/small/445403v1_ufig1.gif" ALT="Figure 1">
View larger version (35K):
org.highwire.dtl.DTLVardef@15ffa95org.highwire.dtl.DTLVardef@3c8032org.highwire.dtl.DTLVardef@452bdorg.highwire.dtl.DTLVardef@1d3f85a_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Modelling, Analysis, and Optimization of Three-Dimensional Restricted Visual Field Metric-Free Swarms

abstract: Models of collective behaviour have been proved helpful in revealing what mechanism may underlie characteristics of a flock of birds, a school of fish, and a swarm of herds. Recently, the metric-free model gradually occupies a dominant position in the research field of collective intelligence. Most of these models endow every single individual with the ability of a global visual field, which can offer each particle sufficient external information. In this paper, we mainly focus on whether the global visual field is necessary to form a consistent and cohesive group or not. Inspired by the biological characteristic of starlings, we develop a three-dimensional restricted visual field metric-free(RVFMF) model based on Pearce and Turners previous work. We further investigate several vital factors governing the convergent consistency of the RVFMF model with the assistance of extensive numerical simulations. According to the simulation results, we conclude that the best view angle of each particle in a swarm increases with the expansion of the population size. Besides, the best view angle gradually becomes stable around 155 degrees when the population size is larger than 1000. We also offer quantitative analysis data to prove that a flock of birds could obtain better consistency under optimal restricted visual field than under global visual field.


title: Genetically engineered coral: A mixed-methods analysis of initial public opinion

abstract: Rising sea surface water temperatures is contributing to coral degradation in the Great Barrier Reef. Synthetic biology technologies offer the potential to enhance coral resilience to higher water temperatures. To explore what the public think of genetically engineered coral, qualitative responses to an open-ended question in a survey of 1,148 of the Australian public were analysed. More respondents supported the technology (59%) than did not (11%). However, a considerable proportion indicated moderate or neutral support (29%). Participants commented about the (moral) right to interfere with nature and uncertainty regarding the consequences of implementing the technology. Participants also mentioned the need to take responsibility and act to save the reef, as well as the benefits likely to result from implementing the technology. Other themes included a desire for further testing and proof, more information, and tight regulation and controls when introducing the technology.


title: In vivo complex haploinsufficiency-based genetic analysis identifies a transcription factor circuit regulating Candida albicans oropharyngeal infection and epithelial cell endocytosis

abstract:
Oropharyngeal candidiasis (OPC) is a common infection that complicates a wide range of medical conditions which can cause either mild or severe disease depending on the patient. The pathobiology of OPC shares many features with candidal biofilms of abiotic surfaces. The transcriptional regulation of C. albicans formation of biofilms on abiotic surfaces has been extensively characterized and involves six key transcription factors (Efg1, Ndt80, Rob1, Bcr1, Brg1, and Tec1). To determine whether this same in vitro biofilm transcriptional regulatory network played a role in OPC, we have carried out the first systematic genetic interaction analysis in a mouse model of C. albicans infection. Whereas all six transcription factors are required for in vitro biofilm formation, only three homozygous deletion mutants (tec1{Delta}{Delta}, bcr1{Delta}{Delta}, and rob1{Delta}{Delta}) and one heterozygous mutant (tec1{Delta}/TEC1) have reduced infectivity in a mouse model of OPC, indicating the network is more robust in vivo than in vitro. Although single mutants (heterozygous or homozygous) of BRG1 and EFG1 have no effect on fungal burden, the double heterozygous and homozygous mutants have dramatically reduced infectivity, indicating a critical genetic interaction between these two transcription factors. Using epistasis analysis, we have formulated a genetic circuit [EFG1+BRG1][->]TEC1[->]BCR1 that is required for OPC infectivity and oral epithelial cell endocytosis. Surprisingly, we also found transcription factor mutants with in vitro defects in filamentation such as efg1{Delta}{Delta} and brg1{Delta}{Delta} filament during oral infection and that decreased filamentation did not correlate with decreased infectivity. Taken together, these data indicate that key in vitro biofilm transcription factors are involved in OPC but that the network characteristics and functional connections are remodeled significantly during interactions with tissues.

Author SummaryThe pathology of oral candidiasis has features of biofilm formation, a well-studied process in vitro. Based on that analogy, we hypothesized that network of transcription factors that regulates in vitro biofilm formation might have similarities and differences in during oral infection. To test this, we employed the first systematic genetic interaction analysis of C. albicans in a mouse model of oropharyngeal infection. This revealed that the six regulators involved in in vitro biofilm formation played roles in vivo but that the functional connections between factors were quite distinct. Surprisingly, we also found that, while many of the factors are required for filamentation in vitro, none of the transcription factor deletion mutants was deficient for this key virulence trait in vivo. These observations clearly demonstrate that C. albicans regulates key aspects of its biology differently in vitro and in vivo.


title: Knockdown of a novel ATPase domain of capsid protein inhibits genome packaging in potato leaf roll virus

abstract: Potato leaf roll virus (PLRV) uses powerful molecular machines to package its genome into a viral capsid employing ATP as fuel. Although, recent bioinformatics and structural studies have revealed detailed mechanism of DNA packaging, little is known about the mechanochemistry of genome packaging in small plant viruses such as PLRV. We have identified a novel P-loop-containing ATPase domain with two Walker A-like motifs, two arginine fingers, and two sensor motifs distributed throughout the polypeptide chain of PLRV capsid protein (CP). The composition and arrangement of the ATP binding and hydrolysis domain of PLRV CP is unique and rarely reported. The discovery of the system sheds new light on the mechanism of viral genome packaging, regulation of viral assembly process, and evolution of plant viruses. Here, we used the RNAi approach to suppress CP gene expression, which in turn prevented PLRV genome packaging and assembly in Solanum tuberosum cv. Khufri Ashoka. Potato plants agroinfiltrated with siRNA constructs against the ATPase domain of CP exhibited no rolling symptoms upon PLRV infection, indicating that the silencing of CP gene expression is an efficient method for generating PLRV-resistant potato plants. Moreover, our findings provide a robust approach to generate PLRV-resistant potato plants, which can be further extended to other species. Finally, we propose a new mechanism of genome packaging and assembly in plant viruses.


title: Quantifying the Inhibitory Impact of Soluble Phenolics on Carbon Mineralization from Sphagnum-rich Peatlands

abstract:
The mechanisms controlling the extraordinarily slow carbon (C) mineralization rates characteristic of Sphagnum-rich peatlands ("bogs") remain somewhat elusive, despite decades of research on this topic. Soluble phenolic compounds have been invoked as potentially significant contributors to bog peat recalcitrance due to their affinity to slow microbial metabolism and cell growth. Despite this potentially significant role, the effects of soluble phenolic compounds on bog peat C mineralization remain unclear.

We analyzed this effect by manipulating the concentration of free soluble phenolics in anaerobic bog peat incubations using water-soluble polyvinylpyrrolidone (PVP), a compound that binds with and inactivates phenolics, preventing phenolic-enzyme interactions. CO2 and CH4 production rates (end-products of C mineralization) correlated positively with PVP concentration following Michaelis-Menten (M.M.) saturation functions. Using M.M. parameters, we determined that soluble phenolics inhibit, at minimum, 57 {+/-} 16% of total C (CO2+CH4) mineralization in the anaerobic incubation conditions studied. These findings are consistent with other studies that have indicated that soluble phenolics play a significant role in regulating bog peat stability in the face of decomposition.