Understanding the genetic mechanisms underlying particular adaptations/phenotypes of organisms is one of the core issues of evolutionary biology. The use of genomic data has greatly advanced our understandings on this issue, as well as other aspects of evolutionary biology, including molecular adaptation, speciation, and even conservation of endangered species. Despite the well-recognized advantages, usages of genomic data are still limited to non-mammal vertebrate groups, partly due to the difficulties in assembling large or highly heterozygous genomes. Although this is particularly the case for amphibians, nonetheless, several comparative and population genomic analyses have shed lights into the speciation and adaptation processes of amphibians in a complex landscape, giving a promising hope for a wider application of genomics in the previously believed challenging groups of organisms. At the same time, these pioneer studies also allow us to realize numerous challenges in studying the molecular adaptations and/or phenotypic evolutionary mechanisms of amphibians. In this review, we first summarize the recent progresses in the study of adaptive evolution of amphibians based on genomic data, and then we give perspectives regarding how to effectively identify key pathways underlying the evolution of complex traits in the genomic era, as well as directions for future research.

As they belong to the most species-rich class of tetrapod vertebrates, birds have long been believed to possess an inferior taste system. However, the bitter taste is fundamental in birds to recognize dietary toxins (which are typically bitter) in potential food sources. To characterize the evolution of avian bitter taste receptor genes (Tas2rs) and to test whether dietary toxins have shaped the repertoire size of avian Tas2rs, we examined 48 genomes representing all but 3 avian orders. The total number of Tas2r genes was found to range from 1 in the domestic pigeon to 12 in the bar-tailed trogon, with an average of 4, which suggested that a much smaller Tas2r gene repertoire exists in birds than in other vertebrates. Furthermore, we uncovered a positive correlation between the number of putatively functional Tas2rs and the abundance of potential toxins in avian diets. Because plant products contain more toxins than animal tissues and insects release poisonous defensive secretions, we hypothesized that herbivorous and insectivorous birds may demand more functional Tas2rs than carnivorous birds feeding on noninsect animals. Our analyses appear to support this hypothesis and highlight the critical role of taste perception in birds.

With the development of high throughput technologies, there are more and more protein-protein interaction (PPI) networks available, which provide a need for efficient computational tools for network alignment. Network alignment is widely used to predict functions of certain proteins, identify conserved network modules, and study the evolutionary relationship across species or biological entities. However, network alignment is an NP-complete problem, and previous algorithms are usually slow or less accurate in aligning big networks like human vs. yeast. In this study, we proposed a fast yet accurate algorithm called Network Alignment by Integrating Biological Process (NAIGO). Specifically, we first divided the networks into subnets taking the advantage of known prior knowledge, such as gene ontology. For each subnet pair, we then developed a novel method to align them by considering both protein orthologous information and their local structural information. After that, we expanded the obtained local network alignments in a greedy manner. Taking the aligned pairs as seeds, we formulated the global network alignment problem as an assignment problem based on similarity matrix, which was solved by the Hungarian method. We applied NAIGO to align human and Saccharomyces cerevisiae S288c PPI network and compared the results with other popular methods like IsoRank, GRAAL, SANA, and NABEECO. As a result, our method outperformed the competitors by aligning more orthologous proteins or matched interactions. In addition, we found a few potential functional orthologous proteins such as RRM2B in human and DNA2 in S. cerevisiae S288c, which are related to DNA repair. We also identified a conserved subnet with six orthologous proteins EXO1, MSH3, MSH2, MLH1, MLH3, and MSH6, and six aligned interactions. All these proteins are associated with mismatch repair. Finally, we predicted a few proteins of S. cerevisiae S288c potentially involving in certain biological processes like autophagosome assembly.

Genome-wide association studies (GWAS) have hitherto identified several germline variants associated with cancer susceptibility, but the molecular functions of these risk modulators remain largely uncharacterized. Recent studies have begun to uncover the regulatory potential of noncoding GWAS SNPs using epigenetic information in corresponding cancer cell types and matched normal tissues. However, this approach does not explore the potential effect of risk germline variants on other important cell types that constitute the microenvironment of tumor or its precursor. This paper presents evidence that the breast-cancer-associated variant rs3903072 may regulate the expression of CTSW in tumor-infiltrating lymphocytes. CTSW is a candidate tumor-suppressor gene, with expression highly specific to immune cells and also positively correlated with breast cancer patient survival. Integrative analyses suggest a putative causative variant in a GWAS-linked enhancer in lymphocytes that loops to the 3' end of CTSW through three-dimensional chromatin interaction. Our work thus poses the possibility that a cancer-associated genetic variant could regulate a gene not only in the cell of cancer origin but also in immune cells in the microenvironment, thereby modulating the immune surveillance by T lymphocytes and natural killer cells and affecting the clearing of early cancer initiating cells.

MicroRNAs (miRNAs) play a critical role in the development of cancers. However, the role of miRNAs in glioma is still poorly understood. In this study, we demonstrate that microRNA-10a (miR-10a) promotes cell migration and invasion by negatively regulating the expression of Eph tyrosine kinase receptor A8 (EphA8). Ectopic expression of EphA8 counteracts the promotion of migration and invasion induced by miR-10a. We further demonstrate that miR-10a and EphA8 regulate epithelial-mesenchymal transition (EMT) to affect cell migration and invasion. Collectively, we unveil a branch of the miR-10a/EphA8 pathway that regulates the progression of glioma.

Increasing evidence has lent support to the notion that miRNAs regulate hepatocellular carcinoma (HCC) cell proliferation by directly targeting cell cycle-related genes. Among these genes, we identified PRPF4B, a CDK-like kinase, as a new target of miR-371-5p. Over-expression of miR-371-5p and knockdown of PRPF4B promotes cell growth by accelerating the G1/S transition in HCC cell lines. Moreover, miR-371-5p promotes tumor growth of QGY-7703 cells in vivo. Conversely, inhibition of miR-371-5p yields an opposing effect. Ectopic expression of PFPF4B abolishes the malignant phenotypes caused by miR-371-5p. Furthermore, contrary to PRPF4B, miR-371 was up-regulated in HCC tissues. Collectively, we highlight the significance of miR-371-5p and PRPF4B in cell cycle progression and hepatocarcinogenesis.

1α,25(OH)2D3 (vitamin D3) is crucial for mineral homeostasis in mammals, but the precise effects of 1α,25(OH)2D3 in adipose tissue remain to be clarified in vivo. The initial 25-hydroxylation is catalyzed by liver microsomal cytochrome P450 2R1 (CYP2R1), which is conserved in vertebrates. To probe the physiological function(s) of 1α,25(OH)2D3 in teleosts, we generated two independent cyp2r1-deficient zebrafish lines. These mutants exhibit retarded growth and increased obesity, especially in the visceral adipose tissue (VAT). These defects could be rescued with 25(OH)D3 treatments. ChIP-PCR analyses demonstrated that pgc1a is the target of the vitamin D receptor in the liver and VAT of zebrafish. Significantly decreased protein levels of Pgc1a, impaired mitochondrial biogenesis, and free fatty acid oxidation are also observed in the cyp2r1 mutant VAT. Our results demonstrate that regulation of 1α,25(OH)2D3 during lipid metabolism occurs through the regulation of Pgc1a for mitochondrial biogenesis and oxidative metabolism within zebrafish VAT.

Adult neurogenesis, the production and integration of new neurons into circuits in the brains of adult animals, is a common feature of a variety of organisms, ranging from insects and crustaceans to birds and mammals. In the mammalian brain the 1st-generation neuronal precursors, the astrocytic stem cells, reside in neurogenic niches and are reported to undergo self-renewing divisions, thereby providing a source of new neurons throughout an animal's life. In contrast, our work shows that the 1st-generation neuronal precursors in the crayfish (Procambarus clarkii) brain, which also have glial properties and lie in a neurogenic niche resembling that of vertebrates, undergo geometrically symmetrical divisions and both daughters appear to migrate away from the niche. However, in spite of this continuous efflux of cells, the number of neuronal precursors in the crayfish niche continues to expand as the animals grow and age. Based on these observations we have hypothesized that (1) the neuronal stem cells in the crayfish brain are not self-renewing, and (2) a source external to the neurogenic niche must provide cells that replenish the stem cell pool.

Piglet splay leg syndrome (PSL) is one of the most frequent genetic defects, and can cause considerable economic loss in pig production. The present understanding of etiology and pathogenesis of PSL is poor. The current study focused on identifying loci associated with PSL through a genome-wide association study (GWAS) performed with the Illumina Porcine60 SNP Beadchip v2.0. The study was a case/control design with four pig populations (Duroc, Landrace, Yorkshire and one crossbred of Landrace × Yorkshire).

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

Chimeric antigen receptor T (CAR-T) cell therapy is not satisfying in solid tumors. PD-1-mediated suppression greatly hinders CAR-T cells in the microenvironment. It has been shown that PD-1 blockade improves the effectiveness of CAR-T cells. Herein, we designed CAR-T cells than could secret α-PD-1 scFv by themselves. To obtain optimal secretions of scFv, we screened several signal peptides. And the segment from human increased the extracellular production of PD-1-neutralizing proteins. The secreted neutralizing scFv efficiently blocked PD-1 and enhanced T cell activation when PD-L1 was present. Further analysis showed that CAR-T cells themselves could secret α-PD-1 scFv with bioactivity. In contrast to the prototype, the scFv-producing CAR-T cells demonstrated decreased PD-1 but increases expansion and toxicity against solid tumor cells. In the subcutaneous and orthotopic xenograft models, the self-delivered α-PD-1 scFv increased CAR-T cell functionalities and tumor-suppressions. Our work suggested that engineering T cells to co-express antigen-responsive receptors and checkpoint inhibitors is effective to optimize CAR-T cell therapy for solid tumors.

Pre-eclampsia is a leading cause of maternal and foetal morbidity and mortality worldwide. Insufficient uteroplacental oxygenation is believed to be responsible for the disease. However, what molecular events involve in hypoxic responses and how they affect placental development remain unclear. Recently, miRNAs have emerged as a new class of molecules in response to hypoxia. We show here that the expression of microRNA-210 (mir-210) is up-regulated in patients with pre-eclampsia, as well as in trophoblast cells cultured under hypoxic conditions. Ectopic expression of mir-210 inhibited the migration and invasion capability of trophoblast cells. Ephrin-A3 and Homeobox-A9, which related with cell migration and vascular remodelling, were then experimentally validated as the functional targets of mir-210 both in vivo and in vitro. Using luciferase reporter, chromatin immunoprecipitation (ChIP) and small interfering RNA (siRNA) experiments, we finally identified a new transcriptional mechanism that the overexpression of mir-210 under hypoxia was regulated by NF-κB transcriptional factor p50, apart from the well-known HIF 1α. Taken together, our study implicates an important role for mir-210 in the molecular mechanism of pre-eclampsia.

MicroRNAs (miRNAs) are non-coding small RNAs that participate in the regulation of a variety of biological processes. Muscle fiber types were very important to meat quality traits, however, the molecular mechanism by which miRNAs regulate the muscle fiber type composition is not fully understood. The aim of this study was to investigate whether miRNA-23a can affect muscle fiber type composition. Luciferase reporter assays proved that miRNA-23a directly targets the 3' untranslated region (UTRs) of MEF2c. Overexpression of miRNA-23a significantly suppressed the expression of MEF2c both in mRNA and protein levels, thus caused down-regulation of the expression of some key downstream genes of MEF2c (PGC1-α, NRF1 and mtTFA). More interestingly, overexpression of miRNA-23a significantly restrained the myogenic differentiation and decreased the ratio of slow myosin heavy chain in myoblasts (p<0.05). Our findings hinted a novel role of miRNA-23a in the epigenetic regulation of meat quality via decreasing the ratio of slow myosin heavy chain isoforms.

Tle6 (Transducin-like enhancer of split 6) is a member of the Tle co-repressor superfamily, which is expressed in various tissues of invertebrates and vertebrates and participates in the developmental process. However, the current research has only found that the TLE6 mutation is related to infertility, and the key regulatory mechanism of TLE6 remains to be explored. In this study, we combined Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas9 and the Tet-on system to construct mouse spermatogonia cell lines that induced TLE6 protein knockout (KO), and studied the effect of Tle6 on mouse spermatogonia proliferation and the cell cycle. The results showed that, after drug induction, the Tle6 gene in mouse spermatogonia was successfully knocked out at the genome and protein levels, and the Tle6 gene knockout efficiency was confirmed to be 87.5% with gene-cloning technology. At the same time, we also found that the mouse spermatogonia proliferated slowly after the Tle6 knockout. Using flow cytometry, we found that the cells did not undergo significant apoptosis, and the number of cells in the S phase decreased. After real-time quantity PCR (qRT-PCR) analysis, we found that the expression of cell-proliferation-related genes, CCAAT enhancer-binding protein α(C/ebp α), granulocyte-colony stimulating factor(G-csf), cyclin-dependent kinases 4(Cdk 4), Cyclin E, proliferating cell nuclear antigen(Pcna), and S-phase kinase-associated protein 2 (Skp2) was significantly reduced, which further affected cell growth. In summary, Tle6 can regulate spermatogonia cell proliferation and the cell cycle and provide a scientific basis for studying the role of TLE6 on spermatogenesis.

The emergence of single-cell RNA-seq (scRNA-seq) technology has made it possible to measure gene expression variations at cellular level. This breakthrough enables the investigation of a wider range of problems including analysis of splicing heterogeneity among individual cells. However, compared to bulk RNA-seq, scRNA-seq data are much noisier due to high technical variability and low sequencing depth. Here we propose SCATS (Single-Cell Analysis of Transcript Splicing) for differential splicing analysis in scRNA-seq, which achieves high sensitivity at low coverage by accounting for technical noise. SCATS models scRNA-seq data either with or without Unique Molecular Identifiers (UMIs). For non-UMI data, SCATS explicitly models technical noise by accounting for capture efficiency and amplification bias through the use of external spike-ins; for UMI data, SCATS models capture efficiency and further accounts for transcriptional burstiness. A key aspect of SCATS lies in its ability to group "exons" that originate from the same isoform(s). Grouping exons is essential in splicing analysis of scRNA-seq data as it naturally aggregates spliced reads across different exons, making it possible to detect splicing events even when sequencing depth is low. To evaluate the performance of SCATS, we analyzed both simulated and real scRNA-seq datasets and compared with existing methods including Census and DEXSeq. We show that SCATS has well controlled type I error rate, and is more powerful than existing methods, especially when splicing difference is small. In contrast, Census suffers from severe type I error inflation, whereas DEXSeq is more conservative. When applied to mouse brain scRNA-seq datasets, SCATS identified more differential splicing events with subtle difference across cell types compared to Census and DEXSeq. With the increasing adoption of scRNA-seq, we believe SCATS will be well-suited for various splicing studies. The implementation of SCATS can be downloaded from https://github.com/huyustats/SCATS.

Neuropeptide F (NPF) signaling systems are widespread and highly evolutionarily conserved from vertebrates to invertebrates. In fact, NPF has been identified in many insect species and plays regulatory roles in diverse physiological processes, such as feeding, learning, reproduction and stress responses. NPF operates by interacting with the NPF receptor (NPFR). Here, we characterized and determined the presumed role of NPF signaling in the wingless parthenogenetic pea aphid, Acyrthosiphon pisum. Quantitative real-time reverse transcription-PCR (qRT-PCR) revealed that the expression levels of both NPF and NPFR transcripts varied across developmental stages, which implies that the NPF signaling system might participate in the developmental regulation of aphid physiological processes or behaviors. The NPF transcript was mainly detected in the head but not in the gut, whereas the NPFR transcript was mainly detected in both the gut and head. In addition, the NPF transcript levels were markedly up-regulated in starved aphids compared with satiated aphids, and the transcript levels recovered after re-feeding. In contrast, the NPFR transcript levels remained stable in starved and re-fed aphids. Furthermore, RNAi knockdown by the injection of NPF dsRNA into wingless adult aphids significantly reduced their food intake. Further analysis of the modification of aphid feeding behavior on broad bean plants using electrical penetration graphs (EPGs) revealed that both the probing time and the total duration of phloem activity decreased significantly in the NPF treatment group. These results indicated a lower appetite for food after NPF knockdown, which could explain the reduction in aphid food intake. NPF silencing was also shown to reduce reproduction but not survival in aphids. Overall, the results of these experiments suggest that NPF plays an important role in regulation of feeding in A. pisum.

Multidonor antibodies are of interest for vaccine design because they can in principle be elicited in the general population by a common set of immunogens. For influenza, multidonor antibodies have been observed against the hemagglutinin (HA) stem, but not the immunodominant HA head. Here, we identify and characterize a multidonor antibody class (LPAF-a class) targeting the HA head. This class exhibits potent viral entry inhibition against H1N1 A/California/04/2009 (CA09) virus. LPAF-a class antibodies derive from the HV2-70 gene and contain a "Tyr-Gly-Asp"-motif, which occludes the HA-sialic acid binding site as revealed by a co-crystal structure with HA. Both germline-reverted and mature LPAF antibodies potently neutralize CA09 virus and have nanomolar affinities for CA09 HA. Moreover, increased frequencies for LPFA-a class antibodies are observed in humans after a single vaccination. Overall, this work highlights the identification of a multidonor class of head-directed influenza-neutralizing antibodies and delineates the mechanism of their recurrent elicitation in humans.

MicroRNAs (miRNAs) have a critical role in tumorigenesis and metastasis, which are major obstacles of cancer therapy. However, the role of miRNAs in colorectal cancer (CRC) metastasis remains poorly understood. Here, we found that miRNA-10a (miR-10a) was upregulated in primary CRC tissues and cell line (SW480) derived from primary CRC compared with metastatic cancer tissues in lymph node and cell line (SW620). The differential expression of miR-10a was inversely correlated with distant metastasis and invasion depth. miR-10a promoted migration and invasion in vitro but inhibited metastasis in vivo by regulating the epithelial-to-mesenchymal transition and anoikis. Furthermore, matrix metalloproteinase 14 (MMP14) and actin gamma 1 (ACTG1) were validated as target genes of miR-10a in CRC cells. Ectopic expression of MMP14 and ACTG1 counteracted the decreased cell adhesion and anoikis resistance activities induced by miR-10a. These findings not only describe the mechanism by which miR-10a suppresses CRC metastasis but also suggest the potential prognostic and therapeutic value of miR-10a in CRC patients.

MicroRNAs (miRs) are involved in several physiological processes, including chondrogenic differentiation, however, their expression and roles in the chondrogenic differentiation of human adipose‑derived stem cells (hADSCs) remain to be fully elucidated to date. Our previous study showed that miR‑1307‑3p was significantly downregulated during chondrogenic differentiation by microarray and northern blot analysis. The present study aimed to investigate the effects of miR‑1307‑3p on chondrogenic differentiation and the underlying mechanisms. First, the decreased expression of miR‑1307‑3p was confirmed by reverse transcription‑quantitative polymerase chain reaction analysis. Subsequently, gain‑ and loss‑of‑function of miR‑1307‑3p experiments showed that the overexpression of miR‑1307‑3p suppressed the deposition of cartilage matrix proteoglycans and decreased the expression of cartilage‑related markers, including sex determining region Y‑box 9, collagen type II α1 chain and aggrecan, whereas the knockdown of miR‑1307‑3p had the opposite effect. In addition, bone morphogenetic protein receptor type 2 (BMPR2) was identified as a target of miR‑1307‑3p. Further mechanistic investigations showed that miR‑1307‑3p attenuated the chondrogenic differentiation of hADSCs at least partly by inhibiting BMPR2‑mothers against decapentaplegic signaling pathways. In conclusion, the findings revealed that miR‑1307‑3p inhibited the chondrogenic differentiation of hADSCs by targeting BMPR2 and its downstream signaling pathway, which may provide novel therapeutic clues for the treatment of cartilage injury.

Socially affected traits are affected by direct genetic effects (DGE) and social genetic effects (SGE). DGE and SGE of an individual directly quantify the genetic influence of its own phenotypes and the phenotypes of other individuals, respectively. In the current study, a total of 3,276 Large White pigs from different pens were used, and each pen contained 10 piglets. DGE and SGE were estimated for six socially affected traits, and then a GWAS was conducted to identify SNPs associated with DGE and SGE. Based on the whole-genome re-sequencing, 40 Large White pigs were genotyped and 10,501,384 high quality SNPs were retained for single-locus and multi-locus GWAS. For single-locus GWAS, a total of 54 SNPs associated with DGE and 33 SNPs with SGE exceeded the threshold (P < 5.00E-07) were detected for six growth traits. Of these, 22 SNPs with pleiotropic effects were shared by DGE and SGE. For multi-locus GWAS, a total of 72 and 110 putative QTNs were detected for DGE and SGE, respectively. Of these, 5 SNPs with pleiotropic effects were shared by DGE and SGE. It is noteworthy that 2 SNPs (SSC8: 16438396 for DGE and SSC17: 9697454 for SGE) were detected in single-locus and multi-locus GWAS. Furthermore, 15 positional candidate genes shared by SGE and DGE were identified because of their roles in behaviour, health and disease. Identification of genetic variants and candidate genes for DGE and SGE for socially affected traits will provide a new insight to understand the genetic architecture of socially affected traits in pigs.