Genome‐scale metabolic modeling reveals SARS‐CoV‐2‐induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID‐19 pandemic caused by the SARS‐CoV‐2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS‐CoV‐2 infection using genome‐scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS‐CoV‐2 infection. We next applied the GEM‐based metabolic transformation algorithm to predict anti‐SARS‐CoV‐2 targets that counteract the virus‐induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco‐2 cells. Further generating and analyzing RNA‐sequencing data of remdesivir‐treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti‐SARS‐CoV‐2 drug. Our study provides clinical data‐supported candidate anti‐SARS‐CoV‐2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.     

Genome‐scale metabolic modeling reveals key features of a minimal gene set

Mesoplasma florum, a fast‐growing near‐minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein‐coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species‐specific constraints to produce the functional iJL208 genome‐scale model (GEM) of metabolism. Genome‐wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI‐syn3.0 and its parent JCVI‐syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping‐stone toward model‐driven whole‐genome engineering.     

NGSpeciesID: DNA barcode and amplicon consensus generation from long‐read sequencing data

Third‐generation sequencing technologies, such as Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), have gained popularity over the last years. These platforms can generate millions of long‐read sequences. This is not only advantageous for genome sequencing projects, but also advantageous for amplicon‐based high‐throughput sequencing experiments, such as DNA barcoding. However, the relatively high error rates associated with these technologies still pose challenges for generating high‐quality consensus sequences. Here, we present NGSpeciesID, a program which can generate highly accurate consensus sequences from long‐read amplicon sequencing technologies, including ONT and PacBio. The tool includes clustering of the reads to help filter out contaminants or reads with high error rates and employs polishing strategies specific to the appropriate sequencing platform. We show that NGSpeciesID produces consensus sequences with improved usability by minimizing preprocessing and software installation and scalability by enabling rapid processing of hundreds to thousands of samples, while maintaining similar consensus accuracy as current pipelines.     

Genome‐wide association study of body morphological traits in Sudanese goats

Long‐term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome‐wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185‐scaffold1365‐620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482‐scaffold89‐467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1‐Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes.     

Genome‐wide SNP identification in Fraxinus linking genetic characteristics to tolerance of Agrilus planipennis

Ash (Fraxinus spp.) is one of the most widely distributed tree genera in North America. Populations of ash in the United States and Canada have been decimated by the introduced pest Agrilus planipennis (Coleoptera: Buprestidae; emerald ash borer), having negative impacts on both forest ecosystems and economic interests. The majority of trees succumb to attack by A. planipennis, but some trees have been found to be tolerant to infestation despite years of exposure. Restriction site‐associated DNA (RAD) sequencing was used to sequence ash individuals, both tolerant and susceptible to A. planipennis attack, in order to identify single nucleotide polymorphism (SNP) patterns related to tolerance and health declines. de novo SNPs were called using SAMtools and, after filtering criteria were implemented, a set of 17,807 SNPs were generated. Principal component analysis (PCA) of SNPs aligned individual trees into clusters related to geography; however, five tolerant trees clustered together despite geographic location. A subset of 32 outlier SNPs identified within this group, as well as a subset of 17 SNPs identified based on vigor rating, are potential candidates for the selection of host tolerance. Understanding the mechanisms of host tolerance through genome‐wide association has the potential to restore populations with cultivars that are able to withstand A. planipennis infestation. This study was successful in using RAD‐sequencing in order to identify SNPs that could contribute to tolerance of A. planipennis. This was a first step toward uncovering the genetic basis for host tolerance to A. planipennis. Future studies are needed to identify the functionality of the loci where these SNPs occur and how they may be related to tolerance of A. planipennis attack.     

Genome‐wide detection of CNVs associated with beak deformity in chickens using high‐density 600K SNP arrays

Beak deformity (crossed beaks) is found in several indigenous chicken breeds including Beijing‐You studied here. Birds with deformed beaks have reduced feed intake and poor production performance. Recently, copy number variation (CNV) has been examined in many species and is recognized as a source of genetic variation, especially for disease phenotypes. In this study, to unravel the genetic mechanisms underlying beak deformity, we performed genome‐wide CNV detection using Affymetrix chicken high‐density 600K data on 48 deformed‐beak and 48 normal birds using penncnv . As a result, two and eight CNV regions (CNVRs) covering 0.32 and 2.45 Mb respectively on autosomes were identified in deformed‐beak and normal birds respectively. Further RT‐qPCR studies validated nine of the 10 CNVRs. The ratios of six CNVRs were significantly different between deformed‐beak and normal birds (P < 0.01). Within these six regions, three and 21 known genes were identified in deformed‐beak and normal birds respectively. Bioinformatics analysis showed that these genes were enriched in six GO terms and one KEGG pathway. Five candidate genes in the CNVRs were further validated using RT‐qPCR. The expression of LRIG2 (leucine rich repeats and immunoglobulin like domains 2) was lower in birds with deformed beaks (P < 0.01). Therefore, the LRIG2 gene could be considered a key factor in view of its known functions and its potential roles in beak deformity. Overall, our results will be helpful for future investigations of the genomic structural variations underlying beak deformity in chickens.     

Effect of MiR‐100‐5p on proliferation and apoptosis of goat endometrial stromal cell in vitro and embryo implantation in vivo

The growth of endometrial stromal cells (ESCs) at implantation sites may be a potential factor affecting the success rate of embryo implantation. Incremental proofs demonstrated that ncRNAs (e.g. miRNAs, lncRNAs and circRNAs) were involved in various biological procedures, including proliferation and apoptosis. In this study, the role of miR‐100‐5p on proliferation and apoptosis of goat ESCs in vitro and embryo implantation in vivo was determined. The mRNA expression of miR‐100‐5p was significantly inhibited in the receptive phase (RE) rather than in the pre‐receptive phase (PE). Overexpression of miR‐100‐5p suppressed ESCs proliferation and induced apoptosis. The molecular target of MiR‐100‐5p, HOXA1, was confirmed by 3′‐UTR assays. Meanwhile, the product of HOXA1 mRNA RT‐PCR increased in the RE more than that in the PE. The HOXA1‐siRNA exerted significant negative effects on growth arrest. Instead, incubation of ESCs with miR‐100‐5p inhibitor or overexpressed HOXA1 promoted the cell proliferation. In addition, Circ‐9110 which acted as a sponge for miR‐100‐5p reversed the relevant biological effects of miR‐100‐5p. The intrinsic apoptosis pathway was suppressed in ESCs, revealing a crosstalk between Circ‐9110/miR‐100‐5p/HOXA1 axis, PI3K/AKT/mTOR, and ERK1/2 pathways. To further evaluate the progress in study on embryo implantation regulating mechanism of miR‐100‐5p in vivo, the pinopodes of two phases were observed and analysed, suggesting that, as similar as in situ, miR‐100‐5p was involved in significantly regulating embryo implantation in vivo. Mechanistically, miR‐100‐5p performed its embryo implantation function through regulation of PI3K/AKT/mTOR and ERK1/2 pathways by targeting Circ‐9110/miR‐100‐5p/HOXA1 axis in vivo.     

Population genetic structure and evolution of Batesian mimicry in Papilio polytes from the Ryukyu Islands,Japan, analyzed by genotyping‐by‐sequencing

Batesian mimicry is a striking example of Darwinian evolution, in which a mimetic species resembles toxic or unpalatable model species, thereby receiving protection from predators. In some species exhibiting Batesian mimicry, nonmimetic individuals coexist as polymorphism in the same population despite the benefits of mimicry. In a previous study, we proposed that the abundance of mimics is limited by that of the models, leading to polymorphic Batesian mimicry in the swallowtail butterfly, Papilio polytes, on the Ryukyu Islands in Japan. We found that their mimic ratios (MRs), which varied among the Islands, were explained by the model abundance of each habitat, rather than isolation by distance or phylogenetic constraint based on the mitochondrial DNA (mtDNA) analysis. In the present study, this possibility was reexamined based on hundreds of nuclear single nucleotide polymorphisms (SNPs) of 93 P. polytes individuals from five Islands of the Ryukyus. We found that the population genetic and phylogenetic structures of P. polytes largely corresponded to the geographic arrangement of the habitat Islands, and the genetic distances among island populations show significant correlation with the geographic distances, which was not evident by the mtDNA‐based analysis. A partial Mantel test controlling for the present SNP‐based genetic distances revealed that the MRs of P. polytes were strongly correlated with the model abundance of each island, implying that negative frequency‐dependent selection interacting with model species shaped and maintained the mimetic polymorphism. Taken together, our results support the possibility that predation pressure, not isolation by distance or other neutral factors, is a major driving force of evolution of the Batesian mimicry in P. polytes from the Ryukyus.     

Analysis of microbial diversity and succession during Xiaoqu Baijiu fermentation using high‐throughput sequencing technology

In this study, high‐throughput sequencing (HTS) was used to compare and analyze the microbial diversity and succession during the brewing process of xiaoqu Baijiu. A total of 34 phyla and 378 genera of bacteria, as well as four phyla, 32 genera of fungi were detected. At the phylum level, Firmicutes, Proteobacteria, Ascomycota, and Mucoromycota were the dominant groups. During the brewing process of xiaoqu Baijiu, the dominant bacteria were Weissella and unidentified Rickettsiales within the first 2 days of brewing, followed by Lactobacillus at 3 days until to the end of brewing. The dominant fungi were Rhizopus, Saccharomyces, and Issatchenkia. The relative abundance of Rhizopus decreased with the extension of brewing time, while the relative abundance of Saccharomyces increased, and Saccharomyces became the dominant species at the second day of brewing. This study revealed the diversity and changes of the microbial community during the brewing process of xiaoqu Baijiu, providing theoretical support and laying a foundation for future study on the contribution of microbial metabolism during brewing of xiaoqu Baijiu, thereby promoting the development of xiaoqu Baijiu industry.     

Genome‐wide analysis reveals demographic and life‐history patterns associated with habitat modification in landlocked,deep‐spawning sockeye salmon (Oncorhynchus nerka)

Human‐mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life‐history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome‐wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high‐confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep‐water spawning behavior associated with this life‐history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.     

High diversity in Delta variant across countries revealed by genome‐wide analysis of SARS‐CoV‐2 beyond the Spike protein

The highly contagious Delta variant of SARS‐CoV‐2 has become a prevalent strain globally and poses a public health challenge around the world. While there has been extensive focus on understanding the amino acid mutations in the Delta variant’s Spike protein, the mutational landscape of the rest of the SARS‐CoV‐2 proteome (25 proteins) remains poorly understood. To this end, we performed a systematic analysis of mutations in all the SARS‐CoV‐2 proteins from nearly 2 million SARS‐CoV‐2 genomes from 176 countries/territories. Six highly prevalent missense mutations in the viral life cycle‐associated Membrane (I82T), Nucleocapsid (R203M, D377Y), NS3 (S26L), and NS7a (V82A, T120I) proteins are almost exclusive to the Delta variant compared to other variants of concern (mean prevalence across genomes: Delta = 99.74%, Alpha = 0.06%, Beta = 0.09%, and Gamma = 0.22%). Furthermore, we find that the Delta variant harbors a more diverse repertoire of mutations across countries compared to the previously dominant Alpha variant. Overall, our study underscores the high diversity of the Delta variant between countries and identifies a list of amino acid mutations in the Delta variant’s proteome for probing the mechanistic basis of pathogenic features such as high viral loads, high transmissibility, and reduced susceptibility against neutralization by vaccines.     

A functional assay for serum detection of antibodies against SARS‐CoV‐2 nucleoprotein

The humoral immune response to SARS‐CoV‐2 results in antibodies against spike (S) and nucleoprotein (N). However, whilst there are widely available neutralization assays for S antibodies, there is no assay for N‐antibody activity. Here, we present a simple in vitro method called EDNA (electroporated‐antibody‐dependent neutralization assay) that provides a quantitative measure of N‐antibody activity in unpurified serum from SARS‐CoV‐2 convalescents. We show that N antibodies neutralize SARS‐CoV‐2 intracellularly and cell‐autonomously but require the cytosolic Fc receptor TRIM21. Using EDNA, we show that low N‐antibody titres can be neutralizing, whilst some convalescents possess serum with high titres but weak activity. N‐antibody and N‐specific T‐cell activity correlates within individuals, suggesting N antibodies may protect against SARS‐CoV‐2 by promoting antigen presentation. This work highlights the potential benefits of N‐based vaccines and provides an in vitro assay to allow the antibodies they induce to be tested.     

Crystallographic molecular replacement using an in silico‐generated search model of SARS‐CoV‐2 ORF8

The majority of crystal structures are determined by the method of molecular replacement (MR). The range of application of MR is limited mainly by the need for an accurate search model. In most cases, pre‐existing experimentally determined structures are used as search models. In favorable cases, ab initio predicted structures have yielded search models adequate for MR. The ORF8 protein of SARS‐CoV‐2 represents a challenging case for MR using an ab initio prediction because ORF8 has an all β‐sheet fold and few orthologs. We previously determined experimentally the structure of ORF8 using the single anomalous dispersion (SAD) phasing method, having been unable to find an MR solution to the crystallographic phase problem. Following a report of an accurate prediction of the ORF8 structure, we assessed whether the predicted model would have succeeded as an MR search model. A phase problem solution was found, and the resulting structure was refined, yielding structural parameters equivalent to the original experimental solution.     

Asprin‐loaded strontium‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite promotes regeneration of critical bone defects

Our laboratory originally synthesized strontium(Sr)‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite (Sr‐α‐CSH/n‐HA) and demonstrated its ability to repair critical bone defects. This study attempted to incorporate aspirin into it to produce a better bone graft material for critical bone defects. After 5% Sr‐α‐CSH was prepared by coprecipitation and hydrothermal methods, it was mixed with aspirin solution of different concentrations (50 μg/ml, 200 μg/ml, 800 μg/ml and 3200 μg/ml) at a fixed liquid‐solid ratio (0.54 v/w) to obtain aspirin‐loaded Sr‐α‐CSH/n‐HA composite. In vitro experiments were performed on the composite extracts. The tibial defects (3 mm*5 mm) in SD rat model were filled with the composite for 4 weeks and 12 weeks to evaluate its osteogenic capacity in vivo. Our results showed its capability of proliferation, migration and osteogenesis of BMSCs in vitro got improved. In vivo treatment with 800 μg/ml aspirin–loaded Sr‐α‐CSH/n‐HA composite led to significantly more new bone formation in the defects compared with Sr‐α‐CSH/n‐HA composite and significantly promoted the expression of osteogenic‐related genes and inhibited osteoclast activity. In general, our research suggests that aspirin‐loaded Sr‐α‐CSH/n‐HA composite may have a greater capacity of repairing tibial defects in SD rats than simple Sr‐α‐CSH/n‐HA composite.     

Live‐cell super‐resolution imaging of actin using LifeAct‐14 with a PAINT‐based approach

We present direct‐LIVE‐PAINT, an easy‐to‐implement approach for the nanoscopic imaging of protein structures in live cells using labeled binding peptides. We demonstrate the feasibility of direct‐LIVE‐PAINT with an actin‐binding peptide fused to EGFP, the location of which can be accurately determined as it transiently binds to actin filaments. We show that direct‐LIVE‐PAINT can be used to image actin structures below the diffraction‐limit of light and have used it to observe the dynamic nature of actin in live cells. We envisage a similar approach could be applied to imaging other proteins within live mammalian cells.     

Stromal cell‐derived factor‐1/Exendin‐4 cotherapy facilitates the proliferation,migration and osteogenic differentiation of human periodontal ligament stem cells in vitro and promotes periodontal bone regeneration in vivo

ObjectivesStromal cell‐derived factor‐1 (SDF‐1) actively directs endogenous cell homing. Exendin‐4 (EX‐4) promotes stem cell osteogenic differentiation. Studies revealed that EX‐4 strengthened SDF‐1‐mediated stem cell migration. However, the effects of SDF‐1 and EX‐4 on periodontal ligament stem cells (PDLSCs) and bone regeneration have not been investigated. In this study, we aimed to evaluate the effects of SDF‐1/EX‐4 cotherapy on PDLSCs in vitro and periodontal bone regeneration in vivo.MethodsCell‐counting kit‐8 (CCK8), transwell assay, qRT‐PCR and western blot were used to determine the effects and mechanism of SDF‐1/EX‐4 cotherapy on PDLSCs in vitro. A rat periodontal bone defect model was developed to evaluate the effects of topical application of SDF‐1 and systemic injection of EX‐4 on endogenous cell recruitment, osteoclastogenesis and bone regeneration in vivo.ResultsSDF‐1/EX‐4 cotherapy had additive effects on PDLSC proliferation, migration, alkaline phosphatase (ALP) activity, mineral deposition and osteogenesis‐related gene expression compared to SDF‐1 or EX‐4 in vitro. Pretreatment with ERK inhibitor U0126 blocked SDF‐1/EX‐4 cotherapy induced ERK signal activation and PDLSC proliferation. SDF‐1/EX‐4 cotherapy significantly promoted new bone formation, recruited more CXCR4⁺ cells and CD90⁺/CD34^‐ stromal cells to the defects, enhanced early‐stage osteoclastogenesis and osteogenesis‐related markers expression in regenerated bone compared to control, SDF‐1 or EX‐4 in vivo.ConclusionsSDF‐1/EX‐4 cotherapy synergistically regulated PDLSC activities, promoted periodontal bone formation, thereby providing a new strategy for periodontal bone regeneration.     

Genome‐wide DNA methylation analysis using next‐generation sequencing to reveal candidate genes responsible for boar taint in pigs

Boar taint (BT) is an offensive flavor observed in non‐castrated male pigs that reduces the carcass price. Surgical castration effectively avoids the taint but is associated with animal welfare concerns. The functional annotation of farm animal genomes for understanding the biology of complex traits can be used in the selection of breeding animals to achieve favorable phenotypic outcomes. The characterization of pig epigenomes/methylation changes between animals with high and low BT and genome‐wide epigenetic markers that can predict BT are lacking. Reduced representation bisulfite sequencing of DNA methylation patterns based on next‐generation sequencing is an efficient technology to identify candidate epigenetic biomarkers associated with BT. Three different BT levels were analyzed using reduced representation bisulfite sequencing data to calculate the methylation levels of cytosine and guanine dinucleotide (CpG) sites. The co‐analysis of differentially methylated CpG sites identified by this study and differentially expressed genes identified by a previous study found 32 significant co‐located genes. The joint analysis of GO terms and pathways revealed that methylation and gene expression of seven candidate genes were associated with BT; in particular, FASN plays a key role in fatty acid biosynthesis, and PEMT might be involved in estrogen regulation and the development of BT. This study is the first to report the genome‐wide DNA methylation profiles of BT in pigs using next‐generation sequencing and summarize candidate genes associated with epigenetic markers of BT, which could contribute to the understanding of the functional biology of BT traits and selective breeding of pigs against BT based on epigenetic biomarkers.     

Identification of liver‐derived bone morphogenetic protein (BMP)‐9 as a potential new candidate for treatment of colorectal cancer

Colorectal cancer (CRC) is a high‐incidence malignancy worldwide which still needs better therapy options. Therefore, the aim of the present study was to investigate the responses of normal or malignant human intestinal epithelium to bone morphogenetic protein (BMP)‐9 and to find out whether the application of BMP‐9 to patients with CRC or the enhancement of its synthesis in the liver could be useful strategies for new therapy approaches. In silico analyses of CRC patient cohorts (TCGA database) revealed that high expression of the BMP‐target gene ID1, especially in combination with low expression of the BMP‐inhibitor noggin, is significantly associated with better patient survival. Organoid lines were generated from human biopsies of colon cancer (T‐Orgs) and corresponding non‐malignant areas (N‐Orgs) of three patients. The N‐Orgs represented tumours belonging to three different consensus molecular subtypes (CMS) of CRC. Overall, BMP‐9 stimulation of organoids promoted an enrichment of tumour‐suppressive gene expression signatures, whereas the stimulation with noggin had the opposite effects. Furthermore, treatment of organoids with BMP‐9 induced ID1 expression (independently of high noggin levels), while treatment with noggin reduced ID1.In summary, our data identify the ratio between ID1 and noggin as a new prognostic value for CRC patient outcome. We further show that by inducing ID1, BMP‐9 enhances this ratio, even in the presence of noggin. Thus, BMP‐9 is identified as a novel target for the development of improved anti‐cancer therapies of patients with CRC.