Accurate prediction of protein structures and interactions using a three-track neural network

DeepMind presented notably accurate predictions at the recent 14th Critical Assessment of Structure Prediction (CASP14) conference. We explored network architectures that incorporate related ideas and obtained the best performance with a three-track network in which information at the one-dimensional (1D) sequence level, the 2D distance map level, and the 3D coordinate level is successively transformed and integrated. The three-track network produces structure predictions with accuracies approaching those of DeepMind in CASP14, enables the rapid solution of challenging x-ray crystallography and cryo-electron microscopy structure modeling problems, and provides insights into the functions of proteins of currently unknown structure. The network also enables rapid generation of accurate protein-protein complex models from sequence information alone, short-circuiting traditional approaches that require modeling of individual subunits followed by docking. We make the method available to the scientific community to speed biological research.

Phosphorylation of innate immune adaptor proteins MAVS, STING, and TRIF induces IRF3 activation

During virus infection, the adaptor proteins MAVS and STING transduce signals from the cytosolic nucleic acid sensors RIG-I and cGAS, respectively, to induce type I interferons (IFNs) and other antiviral molecules. Here we show that MAVS and STING harbor two conserved serine and threonine clusters that are phosphorylated by the kinases IKK and/or TBK1 in response to stimulation. Phosphorylated MAVS and STING then bind to a positively charged surface of interferon regulatory factor 3 (IRF3) and thereby recruit IRF3 for its phosphorylation and activation by TBK1. We further show that TRIF, an adaptor protein in Toll-like receptor signaling, activates IRF3 through a similar phosphorylation-dependent mechanism. These results reveal that phosphorylation of innate adaptor proteins is an essential and conserved mechanism that selectively recruits IRF3 to activate the type I IFN pathway.

Computed structures of core eukaryotic protein complexes

Protein-protein interactions play critical roles in biology, but the structures of many eukaryotic protein complexes are unknown, and there are likely many interactions not yet identified. We take advantage of advances in proteome-wide amino acid coevolution analysis and deep-learning–based structure modeling to systematically identify and build accurate models of core eukaryotic protein complexes within the Saccharomyces cerevisiae proteome. We use a combination of RoseTTAFold and AlphaFold to screen through paired multiple sequence alignments for 8.3 million pairs of yeast proteins, identify 1505 likely to interact, and build structure models for 106 previously unidentified assemblies and 806 that have not been structurally characterized. These complexes, which have as many as five subunits, play roles in almost all key processes in eukaryotic cells and provide broad insights into biological function.

Piezo1/2 mediate mechanotransduction essential for bone formation through concerted activation of NFAT-YAP1-ß-catenin

Mechanical forces are fundamental regulators of cell behaviors. However, molecular regulation of mechanotransduction remain poorly understood. Here, we identified the mechanosensitive channels Piezo1 and Piezo2 as key force sensors required for bone development and osteoblast differentiation. Loss of Piezo1, or more severely Piezo1/2, in mesenchymal or osteoblast progenitor cells, led to multiple spontaneous bone fractures in newborn mice due to inhibition of osteoblast differentiation and increased bone resorption. In addition, loss of Piezo1/2 rendered resistant to further bone loss caused by unloading in both bone development and homeostasis. Mechanistically, Piezo1/2 relayed fluid shear stress and extracellular matrix stiffness signals to activate Ca²⁺ influx to stimulate Calcineurin, which promotes concerted activation of NFATc1, YAP1 and ß-catenin transcription factors by inducing their dephosphorylation as well as NFAT/YAP1/ß-catenin complex formation. Yap1 and ß-catenin activities were reduced in the Piezo1 and Piezo1/2 mutant bones and such defects were partially rescued by enhanced ß-catenin activities.

Protein interaction networks revealed by proteome coevolution

Residue-residue coevolution has been observed across a number of protein-protein interfaces, but the extent of residue coevolution between protein families on the whole-proteome scale has not been systematically studied. We investigate coevolution between 5.4 million pairs of proteins in Escherichia coli and between 3.9 millions pairs in Mycobacterium tuberculosis We find strong coevolution for binary complexes involved in metabolism and weaker coevolution for larger complexes playing roles in genetic information processing. We take advantage of this coevolution, in combination with structure modeling, to predict protein-protein interactions (PPIs) with an accuracy that benchmark studies suggest is considerably higher than that of proteome-wide two-hybrid and mass spectrometry screens. We identify hundreds of previously uncharacterized PPIs in E. coli and M. tuberculosis that both add components to known protein complexes and networks and establish the existence of new ones.

CASP9 assessment of free modeling target predictions

We present an overview of the ninth round of Critical Assessment of Protein Structure Prediction (CASP9) "Template free modeling" category (FM). Prediction models were evaluated using a combination of established structural and sequence comparison measures and a novel automated method designed to mimic manual inspection by capturing both global and local structural features. These scores were compared to those assigned manually over a diverse subset of target domains. Scores were combined to compare overall performance of participating groups and to estimate rank significance. Moreover, we discuss a few examples of free modeling targets to highlight the progress and bottlenecks of current prediction methods. Notably, a server prediction model for a single target (T0581) improved significantly over the closest structure template (44% GDT increase). This accomplishment represents the "winner" of the CASP9 FM category. A number of human expert groups submitted slight variations of this model, highlighting a trend for human experts to act as "meta predictors" by correctly selecting among models produced by the top-performing automated servers. The details of evaluation are available at http://prodata.swmed.edu/CASP9/ .

Genomes of skipper butterflies reveal extensive convergence of wing patterns

For centuries, biologists have used phenotypes to infer evolution. For decades, a handful of gene markers have given us a glimpse of the genotype to combine with phenotypic traits. Today, we can sequence entire genomes from hundreds of species and gain yet closer scrutiny. To illustrate the power of genomics, we have chosen skipper butterflies (Hesperiidae). The genomes of 250 representative species of skippers reveal rampant inconsistencies between their current classification and a genome-based phylogeny. We use a dated genomic tree to define tribes (six new) and subtribes (six new), to overhaul genera (nine new) and subgenera (three new), and to display convergence in wing patterns that fooled researchers for decades. We find that many skippers with similar appearance are distantly related, and several skippers with distinct morphology are close relatives. These conclusions are strongly supported by different genomic regions and are consistent with some morphological traits. Our work is a forerunner to genomic biology shaping biodiversity research.

p38α MAPK regulates proliferation and differentiation of osteoclast progenitors and bone remodeling in an aging-dependent manner

Bone mass is determined by the balance between bone formation, carried out by mesenchymal stem cell-derived osteoblasts, and bone resorption, carried out by monocyte-derived osteoclasts. Here we investigated the potential roles of p38 MAPKs, which are activated by growth factors and cytokines including RANKL and BMPs, in osteoclastogenesis and bone resorption by ablating p38α MAPK in LysM+monocytes. p38α deficiency promoted monocyte proliferation but regulated monocyte osteoclastic differentiation in a cell-density dependent manner, with proliferating p38α^−/− cultures showing increased differentiation. While young mutant mice showed minor increase in bone mass, 6-month-old mutant mice developed osteoporosis, associated with an increase in osteoclastogenesis and bone resorption and an increase in the pool of monocytes. Moreover, monocyte-specific p38α ablation resulted in a decrease in bone formation and the number of bone marrow mesenchymal stem/stromal cells, likely due to decreased expression of PDGF-AA and BMP2. The expression of PDGF-AA and BMP2 was positively regulated by the p38 MAPK-Creb axis in osteoclasts, with the promoters of PDGF-AA and BMP2 having Creb binding sites. These findings uncovered the molecular mechanisms by which p38α MAPK regulates osteoclastogenesis and coordinates osteoclastogenesis and osteoblastogenesis.

CASP9 target classification

The Critical assessment of protein structure prediction round 9 (CASP9) aimed to evaluate predictions for 129 experimentally determined protein structures. To assess tertiary structure predictions, these target structures were divided into domain-based evaluation units that were then classified into two assessment categories: template based modeling (TBM) and template free modeling (FM). CASP9 targets were split into domains of structurally compact evolutionary modules. For the targets with more than one defined domain, the decision to split structures into domains for evaluation was based on server performance. Target domains were categorized based on their evolutionary relatedness to existing templates as well as their difficulty levels indicated by server performance. Those target domains with sequence-related templates and high server prediction performance were classified as TMB, whereas those targets without identifiable templates and low server performance were classified as FM. However, using these generalizations for classification resulted in a blurred boundary between CASP9 assessment categories. Thus, the FM category included those domains without sequence detectable templates (25 target domains) as well as some domains with difficult to detect templates whose predictions were as poor as those without templates (five target domains). Several interesting examples are discussed, including targets with sequence related templates that exhibit unusual structural differences, targets with homologous or analogous structure templates that are not detectable by sequence, and targets with new folds.

RN1, a novel galectin-3 inhibitor, inhibits pancreatic cancer cell growth in vitro and in vivo via blocking galectin-3 associated signaling pathways

Galectin-3 (Gal-3) has been implicated in pancreatic ductal adenocarcinoma (PDAC), and its candidacy as a therapeutic target has been evaluated. Gal-3 is widely upregulated in tumors, and its expression is associated with the development and malignancy of PDAC. In the present study, we demonstrate that a polysaccharide, RN1, purified from the flower of Panax notoginseng binds to Gal-3 and suppresses its expression. In addition, RN1 markedly inhibits PDAC cells growth in vitro, in vivo and in patient-derived xenografts. Mechanistically, RN1 binds to epidermal growth factor receptor (EGFR) and Gal-3, thereby disrupting the interaction between Gal-3 and EGFR and downregulating extracellular-related kinase (ERK) phosphorylation and the transcription factor of Gal-3, Runx1 expression. Inhibiting the expression of Runx1 by RN1, suppresses Gal-3 expression and inactivates Gal-3-associated signaling pathways, including the EGFR/ERK/Runx1, BMP/smad/Id-3 and integrin/FAK/JNK signaling pathways. In addition, RN1 can also bind to bone morphogenetic protein receptors (BMPR1A and BMPR2) and block the interaction between Gal-3 and the BMPRs. Thus, our results suggest that a novel Gal-3 inhibitor RN1 may be a potential candidate for human PDAC treatment via multiple targets and multiple signaling pathways.

Protein contact prediction using metagenome sequence data and residual neural networks

MOTIVATION

Almost all protein residue contact prediction methods rely on the availability of deep multiple sequence alignments (MSAs). However, many proteins from the poorly populated families do not have sufficient number of homologs in the conventional UniProt database. Here we aim to solve this issue by exploring the rich sequence data from the metagenome sequencing projects.

RESULTS

Based on the improved MSA constructed from the metagenome sequence data, we developed MapPred, a new deep learning-based contact prediction method. MapPred consists of two component methods, DeepMSA and DeepMeta, both trained with the residual neural networks. DeepMSA was inspired by the recent method DeepCov, which was trained on 441 matrices of covariance features. By considering the symmetry of contact map, we reduced the number of matrices to 231, which makes the training more efficient in DeepMSA. Experiments show that DeepMSA outperforms DeepCov by 10-13% in precision. DeepMeta works by combining predicted contacts and other sequence profile features. Experiments on three benchmark datasets suggest that the contribution from the metagenome sequence data is significant with P-values less than 4.04E-17. MapPred is shown to be complementary and comparable the state-of-the-art methods. The success of MapPred is attributed to three factors: the deeper MSA from the metagenome sequence data, improved feature design in DeepMSA and optimized training by the residual neural networks.

AVAILABILITY AND IMPLEMENTATION

http://yanglab.nankai.edu.cn/mappred/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Complete genomes of Hairstreak butterflies, their speciation, and nucleo-mitochondrial incongruence

Comparison of complete genomes of closely related species enables research on speciation and how phenotype is determined by genotype. Lepidoptera, an insect order of 150,000 species with diverse phenotypes, is well-suited for such comparative genomics studies if new genomes, which cover additional Lepidoptera families are acquired. We report a 729 Mbp genome assembly of the Calycopis cecrops, the first genome from the family Lycaenidae and the largest available Lepidoptera genome. As detritivore, Calycopis shows expansion in detoxification and digestion enzymes. We further obtained complete genomes of 8 Calycopis specimens: 3 C. cecrops and 5 C. isobeon, including a dry specimen stored in the museum for 30 years. The two species differ subtly in phenotype and cannot be differentiated by mitochondrial DNA. However, nuclear genomes revealed a deep split between them. Genes that can clearly separate the two species (speciation hotspots) mostly pertain to circadian clock, mating behavior, transcription regulation, development and cytoskeleton. The speciation hotspots and their function significantly overlap with those we previously found in Pterourus, suggesting common speciation mechanisms in these butterflies.

Quantitative analysis of routine chemical constituents in tobacco by near-infrared spectroscopy and support vector machine

It is important to monitor quality of tobacco during the production of cigarette. Therefore, in order to scientifically control the tobacco raw material and guarantee the cigarette quality, fast and accurate determination routine chemical of constituents of tobacco, including the total sugar, reducing sugar, Nicotine, the total nitrogen and so on, is needed. In this study, 50 samples of tobacco from different cultivation areas were surveyed by near-infrared (NIR) spectroscopy, and the spectral differences provided enough quantitative analysis information for the tobacco. Partial least squares regression (PLSR), artificial neural network (ANN), and support vector machine (SVM), were applied. The quantitative analysis models of 50 tobacco samples were studied comparatively in this experiment using PLSR, ANN, radial basis function (RBF) SVM regression, and the parameters of the models were also discussed. The spectrum variables of 50 samples had been compressed through the wavelet transformation technology before the models were established. The best experimental results were obtained using the (RBF) SVM regression with gamma=1.5, 1.3, 0.9, and 0.1, separately corresponds to total sugar, reducing sugar, Nicotine, and total nitrogen, respectively. Finally, compared with the back propagation (BP-ANN) and PLSR approach, SVM algorithm showed its excellent generalization for quantitative analysis results, while the number of samples for establishing the model is smaller. The overall results show that NIR spectroscopy combined with SVM can be efficiently utilized for rapid and accurate analysis of routine chemical compositions in tobacco. Simultaneously, the research can serve as the technical support and the foundation of quantitative analysis of other NIR applications.

Complete genome of Pieris rapae, a resilient alien, a cabbage pest, and a source of anti-cancer proteins

The Small Cabbage White ( Pieris rapae) is originally a Eurasian butterfly. Being accidentally introduced into North America, Australia, and New Zealand a century or more ago, it spread throughout the continents and rapidly established as one of the most abundant butterfly species. Although it is a serious pest of cabbage and other mustard family plants with its caterpillars reducing crops to stems, it is also a source of pierisin, a protein unique to the Whites that shows cytotoxicity to cancer cells. To better understand the unusual biology of this omnipresent agriculturally and medically important butterfly, we sequenced and annotated the complete genome from USA specimens. At 246 Mbp, it is among the smallest Lepidoptera genomes reported to date. While 1.5% positions in the genome are heterozygous, they are distributed highly non-randomly along the scaffolds, and nearly 20% of longer than 1000 base-pair segments are SNP-free (median length: 38000 bp). Computational simulations of population evolutionary history suggest that American populations started from a very small number of introduced individuals, possibly a single fertilized female, which is in agreement with historical literature. Comparison to other Lepidoptera genomes reveals several unique families of proteins that may contribute to the unusual resilience of Pieris. The nitrile-specifier proteins divert the plant defense chemicals to non-toxic products. The apoptosis-inducing pierisins could offer a defense mechanism against parasitic wasps. While only two pierisins from Pieris rapae were characterized before, the genome sequence revealed eight, offering additional candidates as anti-cancer drugs. The reference genome we obtained lays the foundation for future studies of the Cabbage White and other Pieridae species.

A self-amplifying loop of YAP and SHH drives formation and expansion of heterotopic ossification

Heterotopic ossification (HO) occurs as a common complication after injury or in genetic disorders. The mechanisms underlying HO remain incompletely understood, and there are no approved prophylactic or secondary treatments available. Here, we identify a self-amplifying, self-propagating loop of Yes-associated protein (YAP)-Sonic hedgehog (SHH) as a core molecular mechanism underlying diverse forms of HO. In mouse models of progressive osseous heteroplasia (POH), a disease caused by null mutations in GNAS, we found that Gnas^-/- mesenchymal cells secreted SHH, which induced osteoblast differentiation of the surrounding wild-type cells. We further showed that loss of Gnas led to activation of YAP transcription activity, which directly drove Shh expression. Secreted SHH further induced YAP activation, Shh expression, and osteoblast differentiation in surrounding wild-type cells. This self-propagating positive feedback loop was both necessary and sufficient for HO expansion and could act independently of Gnas in fibrodysplasia ossificans progressiva (FOP), another genetic HO, and nonhereditary HO mouse models. Genetic or pharmacological inhibition of YAP or SHH abolished HO in POH and FOP and acquired HO mouse models without affecting normal bone homeostasis, providing a previously unrecognized therapeutic rationale to prevent, reduce, and shrink HO.

Speciation in Cloudless Sulphurs Gleaned from Complete Genomes

For 200 years, zoologists have relied on phenotypes to learn about the evolution of animals. A glance at the genotype, even through several gene markers, revolutionized our understanding of animal phylogeny. Recent advances in sequencing techniques allow researchers to study speciation mechanisms and the link between genotype and phenotype using complete genomes. We sequenced and assembled a complete genome of the Cloudless Sulphur (Phoebis sennae) from a single wild-caught specimen. This genome was used as reference to compare genomes of six specimens, three from the eastern populations (Oklahoma and north Texas), referred to as a subspeciesPhoebis sennae eubule, and three from the southwestern populations (south Texas) known as a subspeciesPhoebis sennae marcellina While the two subspecies differ only subtly in phenotype and mitochondrial DNA, comparison of their complete genomes revealed consistent and significant differences, which are more prominent than those between tiger swallowtailsPterourus canadensisandPterourus glaucus The two sulphur taxa differed in histone methylation regulators, chromatin-associated proteins, circadian clock, and early development proteins. Despite being well separated on the whole-genome level, the two taxa show introgression, with gene flow mainly fromP. s. marcellinatoP. s. eubule Functional analysis of introgressed genes reveals enrichment in transmembrane transporters. Many transporters are responsible for nutrient uptake, and their introgression may be of selective advantage for caterpillars to feed on more diverse food resources. Phylogenetically, complete genomes place family Pieridae away from Papilionidae, which is consistent with previous analyses based on several gene markers.

The ABC transporters in Candidatus Liberibacter asiaticus

Candidatus Liberibacter asiaticus (Ca. L. asiaticus) is a Gram-negative bacterium and the pathogen of Citrus Greening disease (Huanglongbing, HLB). As a parasitic bacterium, Ca. L. asiaticus harbors ABC transporters that play important roles in exchanging chemical compounds between Ca. L. asiaticus and its host. Here, we analyzed all the ABC transporter-related proteins in Ca. L. asiaticus. We identified 14 ABC transporter systems and predicted their structures and substrate specificities. In-depth sequence and structure analysis including multiple sequence alignment, phylogenetic tree reconstruction, and structure comparison further support their function predictions. Our study shows that this bacterium could use these ABC transporters to import metabolites (amino acids and phosphates) and enzyme cofactors (choline, thiamine, iron, manganese, and zinc), resist to organic solvent, heavy metal, and lipid-like drugs, maintain the composition of the outer membrane (OM), and secrete virulence factors. Although the features of most ABC systems could be deduced from the abundant experimental data on their orthologs, we reported several novel observations within ABC system proteins. Moreover, we identified seven nontransport ABC systems that are likely involved in virulence gene expression regulation, transposon excision regulation, and DNA repair. Our analysis reveals several candidates for further studies to understand and control the disease, including the type I virulence factor secretion system and its substrate that are likely related to Ca. L. asiaticus pathogenicity and the ABC transporter systems responsible for bacterial OM biosynthesis that are good drug targets.

Skipper genome sheds light on unique phenotypic traits and phylogeny

BACKGROUND

Butterflies and moths are emerging as model organisms in genetics and evolutionary studies. The family Hesperiidae (skippers) was traditionally viewed as a sister to other butterflies based on its moth-like morphology and darting flight habits with fast wing beats. However, DNA studies suggest that the family Papilionidae (swallowtails) may be the sister to other butterflies including skippers. The moth-like features and the controversial position of skippers in Lepidoptera phylogeny make them valuable targets for comparative genomics.

RESULTS

We obtained the 310 Mb draft genome of the Clouded Skipper (Lerema accius) from a wild-caught specimen using a cost-effective strategy that overcomes the high (1.6 %) heterozygosity problem. Comparative analysis of Lerema accius and the highly heterozygous genome of Papilio glaucus revealed differences in patterns of SNP distribution, but similarities in functions of genes that are enriched in non-synonymous SNPs. Comparison of Lepidoptera genomes revealed possible molecular bases for unique traits of skippers: a duplication of electron transport chain components could result in efficient energy supply for their rapid flight; a diversified family of predicted cellulases might allow them to feed on cellulose-enriched grasses; an expansion of pheromone-binding proteins and enzymes for pheromone synthesis implies a more efficient mate-recognition system, which compensates for the lack of clear visual cues due to the similarities in wing colors and patterns of many species of skippers. Phylogenetic analysis of several Lepidoptera genomes suggested that the position of Hesperiidae remains uncertain as the tree topology varied depending on the evolutionary model.

CONCLUSION

Completion of the first genome from the family Hesperiidae allowed comparative analyses with other Lepidoptera that revealed potential genetic bases for the unique phenotypic traits of skippers. This work lays the foundation for future experimental studies of skippers and provides a rich dataset for comparative genomics and phylogenetic studies of Lepidoptera.

Bone Size and Quality Regulation: Concerted Actions of mTOR in Mesenchymal Stromal Cells and Osteoclasts

The bone size and quality, acquired during adolescent growth under the influence of anabolic hormones, growth factors, and nutrients, determine the height and bone stability and forecast osteoporosis risks in late life. Yet bone size and quality control mechanisms remain enigmatic. To study the roles of mammalian target of rapamycin (mTOR) signaling, sensor of growth factors and nutrients, in bone size and quality regulation, we ablated Tsc1, a suppressor of mTOR, in mesenchymal stromal cells (MSCs), monocytes, or their progenies osteoblasts and osteoclasts. mTOR activation in MSCs, but much less in osteoblasts, increased bone width and mass due to MSC hyperproliferation, but decreased bone length and mineral contents due to defective MSC differentiation. mTOR activation promotes bone mineral accretion by inhibiting osteoclast differentiation and activity directly or via coupling with MSCs. Tuberous sclerosis complex patient studies confirmed these findings. Thus, mTOR regulates bone size via MSCs and bone quality by suppressing catabolic activities of osteoclasts.

High-accuracy refinement using Rosetta in CASP13

Because proteins generally fold to their lowest free energy states, energy-guided refinement in principle should be able to systematically improve the quality of protein structure models generated using homologous structure or co-evolution derived information. However, because of the high dimensionality of the search space, there are far more ways to degrade the quality of a near native model than to improve it, and hence, refinement methods are very sensitive to energy function errors. In the 13th Critial Assessment of techniques for protein Structure Prediction (CASP13), we sought to carry out a thorough search for low energy states in the neighborhood of a starting model using restraints to avoid straying too far. The approach was reasonably successful in improving both regions largely incorrect in the starting models as well as core regions that started out closer to the correct structure. Models with GDT-HA over 70 were obtained for five targets and for one of those, an accuracy of 0.5 å backbone root-mean-square deviation (RMSD) was achieved. An important current challenge is to improve performance in refining oligomers and larger proteins, for which the search problem remains extremely difficult.

When COI barcodes deceive: complete genomes reveal introgression in hairstreaks

Two species of hairstreak butterflies from the genus Calycopis are known in the United States: C. cecrops and C. isobeon Analysis of mitochondrial COI barcodes of Calycopis revealed cecrops-like specimens from the eastern US with atypical barcodes that were 2.6% different from either USA species, but similar to Central American Calycopis species. To address the possibility that the specimens with atypical barcodes represent an undescribed cryptic species, we sequenced complete genomes of 27 Calycopis specimens of four species: C. cecrops, C. isobeon, C. quintana and C. bactra Some of these specimens were collected up to 60 years ago and preserved dry in museum collections, but nonetheless produced genomes as complete as fresh samples. Phylogenetic trees reconstructed using the whole mitochondrial and nuclear genomes were incongruent. While USA Calycopis with atypical barcodes grouped with Central American species C. quintana by mitochondria, nuclear genome trees placed them within typical USA C. cecrops in agreement with morphology, suggesting mitochondrial introgression. Nuclear genomes also show introgression, especially between C. cecrops and C. isobeon About 2.3% of each C. cecrops genome has probably (p-value < 0.01, FDR < 0.1) introgressed from C. isobeon and about 3.4% of each C. isobeon genome may have come from C. cecrops. The introgressed regions are enriched in genes encoding transmembrane proteins, mitochondria-targeting proteins and components of the larval cuticle. This study provides the first example of mitochondrial introgression in Lepidoptera supported by complete genome sequencing. Our results caution about relying solely on COI barcodes and mitochondrial DNA for species identification or discovery.

Preparation of a Superhydrophobic and Peroxidase-like Activity Array Chip for H2O2 Sensing by Surface-Enhanced Raman Scattering

In this paper, we propose a novel and simple method for preparing a dual-biomimetic functional array possessing both superhydrophobic and peroxidase-like activity that can be used for hydrogen peroxide (H2O2) sensing. The proposed method is an integration innovation that combines the above two properties and surface-enhanced Raman scattering (SERS). We integrated a series of well-ordered arrays of Au points (d = 1 mm) onto a superhydrophobic copper (Cu)/silver (Ag) surface by replicating an arrayed molybdenum template. Instead of using photoresists and the traditional lithography method, we utilized a chemical etching method (a substitution reaction between Cu and HAuCl4) with a Cu/Ag superhydrophobic surface as the barrier layer, which has the benefit of water repellency. The as-prepared Au points were observed to possess peroxidase-like activity, allowing for catalytic oxidation of the chromogenic molecule o-phenylenediamine dihydrochloride (OPD). Oxidation was evidenced by a color change in the presence of H2O2, which allows the array chip to act as an H2O2 sensor. In this study, the water repellency of the superhydrophobic surface was used to fabricate the array chip and increase the local reactant concentration during the catalytic reaction. As a result, the catalytic reaction occurred when only 2 μL of an aqueous sample (OPD/H2O2) was placed onto the Au point, and the enzymatic product, 2,3-diaminophenazine, showed a SERS signal distinguishable from that of OPD after mixing with 2 μL of colloidal Au. Using the dual-biomimetic functional array chip, quantitative analysis of H2O2 was performed by observing the change in the SERS spectra, which showed a concentration-dependent behavior for H2O2. This method allows for the detection of H2O2 at concentrations as low as 3 pmol per 2 μL of sample, which is a considerable advantage in H2O2 analysis. The as-prepared substrate was convenient for H2O2 detection because only a small amount of sample was required in each analysis. Highly sensitive detection was realized using SERS. Therefore, this chip was shown to exhibit significant potential for applications in bioanalysis.

A tissue injury sensing and repair pathway distinct from host pathogen defense

Pathogen infection and tissue injury are universal insults that disrupt homeostasis. Innate immunity senses microbial infections and induces cytokines/chemokines to activate resistance mechanisms. Here, we show that, in contrast to most pathogen-induced cytokines, interleukin-24 (IL-24) is predominately induced by barrier epithelial progenitors after tissue injury and is independent of microbiome or adaptive immunity. Moreover, Il24 ablation in mice impedes not only epidermal proliferation and re-epithelialization but also capillary and fibroblast regeneration within the dermal wound bed. Conversely, ectopic IL-24 induction in the homeostatic epidermis triggers global epithelial-mesenchymal tissue repair responses. Mechanistically, Il24 expression depends upon both epithelial IL24-receptor/STAT3 signaling and hypoxia-stabilized HIF1α, which converge following injury to trigger autocrine and paracrine signaling involving IL-24-mediated receptor signaling and metabolic regulation. Thus, parallel to innate immune sensing of pathogens to resolve infections, epithelial stem cells sense injury signals to orchestrate IL-24-mediated tissue repair.

Gαs signaling in skeletal development, homeostasis and diseases

Skeletal development is exquisitely controlled both spatially and temporally by cell signaling networks. Gα_s is the stimulatory α-subunit in a heterotrimeric G protein complex transducing the signaling of G-protein-coupled receptors (GPCRs), responsible for controlling both skeletal development and homeostasis. Gα_s, encoded by the GNAS gene in humans, plays critical roles in skeletal development and homeostasis by regulating commitment, differentiation and maturation of skeletal cells. Gα_s-mediated signaling interacts with the Wnt and Hedgehog signaling pathways, both crucial regulators of skeletal development, remodeling and injury repair. Genetic mutations that disrupt Gα_s functions cause human disorders with severe skeletal defects, such as fibrous dysplasia of bone and heterotopic bone formation. This chapter focuses on the crucial roles of Gα_s signaling during skeletal development and homeostasis, and the pathological mechanisms underlying skeletal diseases caused by GNAS mutations.