Exome sequencing of deer mice on two California Channel Islands identifies potential adaptation to strongly contrasting ecological conditions

Understanding the forces that drive genotypic and phenotypic change in wild populations is a central goal of evolutionary biology. We examined exome variation in populations of deer mice from two of the California Channel Islands: Peromyscus maniculatus elusus from Santa Barbara Island and P. m. santacruzae from Santa Cruz Island exhibit significant differences in olfactory predator recognition, activity timing, aggressive behavior, morphology, prevalence of Sin Nombre virus, and population densities. We characterized variation in protein‐coding regions using exome capture and sequencing of 25 mice from Santa Barbara Island and 22 mice from Santa Cruz Island. We identified and examined 386,256 SNPs using three complementary methods (BayeScan, pcadapt, and LFMM). We found strong differences in molecular variation between the two populations and 710 outlier SNPs in protein‐coding genes that were detected by all three methods. We identified 35 candidate genes from this outlier set that were related to differences in phenotypes between island populations. Enrichment analyses demonstrated that patterns of molecular variation were associated with biological processes related to response to chemical stimuli and regulation of immune processes. Candidate genes associated with olfaction (Gfy, Tlr2, Vmn13r2, numerous olfactory receptor genes), circadian activity (Cry1), anxiety (Brca1), immunity (Cd28, Eif2ak4, Il12a, Syne1), aggression (Cyp19a, Lama2), and body size (Bc16, Syne1) exhibited non‐synonymous mutations predicted to have moderate to large effects. Variation in olfaction‐related genes, including a stop codon in the Santa Barbara Island population, suggests loss of predator‐recognition traits at the molecular level, consistent with a lack of behavioral aversion to fox feces. These findings also suggest that divergent pathogen prevalence and population density may have influenced adaptive immunity and behavioral phenotypes, such as reduced aggression. Overall, our study indicates that ecological differences between islands are associated with signatures of selection in protein‐coding genes underlying phenotypes that promote success in those environments.     

Identification of immune infiltration landscape on prognosis and therapy of the ferroptosis-related genes signature in breast cancer

Ferroptosis is a unique iron-dependent cell death mechanism characterized by the generation of lipid reactive oxygen species (ROS) in cancer cells, which leads to mitochondrial metabolic dysregulation. However, how could the tumor immune microenvironment (TIME) modulates ferroptosis remains unclear. Thus, by integrating multiple algorithms, we revealed the novel functional and immune patterns of the ferroptosis-related genes (FRGs) in breast cancer. Five prognostic FRGs were finally selected for the prognostic signature and four of which were identified as the independent biomarkers for immunotherapies. The consensus cluster analysis illustrated the FRGs were characterized by the metabolism dysfunction and immune infiltration cells, meanwhile, these FRGs have the same stem cell characteristics and response efficacy to the immunotherapies. In conclusion, a comprehensive analysis of the FRGs in breast cancer was conducted to develop a prognostic gene signature. Functional and immunological evidence of vulnerabilities in the interaction between ferroptosis and the TIME was also revealed. Further data and research are required.     

Molecular profile of the NF‐κB signalling pathway in human colorectal cancer

The development and progression of colorectal cancer (CRC) have been associated with inflammation processes that involve the overactivation of the NF‐κB signalling pathway. The characterization of the NF‐κB expression profile in CRC is an important topic since the suppression of NF‐κB represents a potential therapeutic approach. In this study, we assessed the expression levels of 84 NF‐κB‐related genes in paired tumoral (T) and peritumoral (PT) tissues from 18 CRC patients and 18 normal colonic mucosae, and the expression levels of three miRNAs targeting the most dysregulated genes revealed by the case–control analysis. Comparing the gene expression profile of T and controls, 60 genes were dysregulated. The comparison of T and PT revealed 17 dysregulated genes in the tumoral tissues, with IL1B, CXCL8, IL1A, and CSF2 being the most upregulated. Notably, through a bioinformatics analysis, the differential gene expression of 11 out of the 17 genes was validated on a larger cohort of 308 CRC patients compared with 41 controls. Moreover, a decrease in the levels of RELA, NOD1, CASP8, BCL2L1, ELK1, and IKBKB was identified in poorly differentiated tumours compared to moderately differentiated tumours. The analysis of the three miRNAs targeting IL1B, CXCL8, IL1A, and CSF2 showed that miR‐182‐5p was upregulated in T compared with PT, whereas miR‐10b‐5p was downregulated in T compared with PT and control tissues. Our results may contribute to the design of new experimental therapeutic strategies based on endogenous molecules, such as miRNAs, to target the genetic key players of the NF‐ κB pathway.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

Selection of a reference gene for studies on lipid‐related aquatic adaptations of toothed whales (Grampus griseus)

Toothed whales are one group of marine mammals that has developed special adaptations, such as echolocation for predation, to successfully live in a dynamic aquatic environment. Their fat metabolism may differ from that of other mammals because toothed whales have acoustic fats. Gene expression in the metabolic pathways of animals can change with respect to their evolution and environment. A real‐time quantitative polymerase chain reaction (RT‐qPCR) is a reliable technique for studying the relative expressions of genes. However, since the accuracy of RT‐qPCR data is totally dependent on the reference gene, the selection of the reference gene is an essential step. In this study, 10 candidate reference genes (ZC3H10, FTL, LGALS1, RPL27, GAPDH, FTH1, DCN, TCTP, NDUS5, and UBIM) were initially tested for amplification efficiency using RT‐qPCR. After excluding DCN, the remaining nine genes, which are nearly 100% efficient, were selected for the gene stability analysis. Stable reference genes across eight different fat tissue, liver, and muscle samples from Grampus griseus were identified by four algorithms, which were provided in Genorm, NormFinder, BestKeeper, and Delta CT. Finally, a RefFinder comprehensive ranking was performed based on the stability values, and the nine genes were ranked as follows: LGALS1 > FTL > GAPDH > ZC3H10 > FTH1 > NDUS5 > TCTP > RPL27 > UBIM. The LGALS1 and FTL genes were identified as the most stable novel reference genes. The third‐ranked gene, GAPDH, is a well‐known housekeeping gene for mammals. Ultimately, we suggest the use of LGALS1 as a reliable novel reference gene for genomics studies on the lipid‐related aquatic adaptations of toothed whales.     

Cytochrome P450 26A1 regulates the clusters and killing activity of NK cells during the peri‐implantation period

Cytochrome P450 26A1 (CYP26A1) plays a vital role in early pregnancy in mice. Our previous studies have found that CYP26A1 affects embryo implantation by modulating natural killer (NK) cells, and that there is a novel population of CYP26A1⁺ NK cells in the uteri of pregnant mice. The aim of this study was to investigate the effects of CYP26A1 on the subsets and killing activity of NK cells. Through single‐cell RNA sequencing (scRNA‐seq), we identified four NK cell subsets in the uterus, namely, conventional NK (cNK), tissue‐resident NK (trNK) 1 and 2, and proliferating trNK (trNKp). The two most variable subpopulations after uterine knockdown of CYP26A1 were trNKp and trNK2 cells. CYP26A1 knockdown significantly downregulated the expression of the NK cell function‐related genes Cd44, Cd160, Vegfc, and Slamf6 in trNK2 cells, and Klra17 and Ogn in trNKp cells. Both RNA‐seq and cytotoxicity assays confirmed that CYP26A1⁺ NK cells had low cytotoxicity. These results indicate that CYP26A1 may affect the immune microenvironment at the maternal‐foetal interface by regulating the activity of NK cells.     

Impacts of seasonality on gene expression in the Chinese horseshoe bat

Seasonality can cause changes in many environmental factors which potentially affects gene expression. Here, we used a bat species (Rhinolophus sinicus) from eastern China as a model to explore the molecular mechanisms of seasonal effects, in particular during phenological shifts in the spring and autumn. Based on the analysis of 45 RNA‐seq samples, we found strong seasonal effects on gene expression, with a large number of genes identified as either specific or biased to each season. Weighted gene co‐expression network analysis also identified multiple modules significantly associated with each season. These seasonal genes were further enriched into different functional categories. Consistent with effects of phenological shifts on bats, we found that genes related to promoting food intake were highly expressed in both autumn and spring. In addition, immunity genes were also highly expressed in both seasons although this seasonal immune response had tissue specificity in different seasons. In female bats, genes related to the delay of ovulation (e.g., NPPC, natriuretic peptide precursor type C) were highly expressed in October, while genes associated with the promotion of reproduction (e.g., DIO2, iodothyronine deiodinase 2) were biasedly expressed in April. Lastly, we found multiple known core clock genes in both October‐biased and April‐biased expressed genes, which may be involved in regulating the start and end of hibernation, respectively. Overall, together with studies in other land and aquatic animals, our work supports that seasonal gene expression variations may be a general evolutionary response to environmental changes in wild animals.     

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein‐coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.     

Rhizosphere soil bacterial communities and nitrogen cycling affected by deciduous and evergreen tree species

Deciduous and evergreen trees differ in their responses to drought and nitrogen (N) demand. Whether or not these functional types affect the role of the bacterial community in the N cycle during drought remains uncertain. Two deciduous tree species (Alnus cremastogyne, an N₂‐fixing species, and Liquidambar formosana) and two evergreen trees (Cunninghamia lanceolata and Pinus massoniana) were used to assess factors in controlling rhizosphere soil bacterial community and N cycling functions. Photosynthetic rates and biomass production of plants, 16S rRNA sequencing and N‐cycling‐related genes of rhizosphere soil were measured. The relative abundance of the phyla Actinobacteria and Firmicutes was higher, and that of Proteobacteria, Acidobacteria, and Gemmatimondaetes was lower in rhizosphere soil of deciduous trees than that of evergreen. Beta‐diversity of bacterial community also significantly differed between the two types of trees. Deciduous trees showed significantly higher net photosynthetic rates and biomass production than evergreen species both at well water condition and short‐term drought. Root biomass was the most important factor in driving soil bacterial community and N‐cycling functions than total biomass and aboveground biomass. Furthermore, 44 bacteria genera with a decreasing response and 46 taxa showed an increased response along the root biomass gradient. Regarding N‐cycle‐related functional genes, copy numbers of ammonia‐oxidizing bacteria (AOB) and autotrophic ammonia‐oxidizing archaea (AOA), N₂ fixation gene (nifH), and denitrification genes (nirK, nirS) were significantly higher in the soil of deciduous trees than in that of the evergreen. Structural equation models explained 50.2%, 47.6%, 48.6%, 49.4%, and 37.3% of the variability in copy numbers of nifH, AOB, AOA, nirK, and nirS, respectively, and revealed that root biomass had significant positive effects on copy numbers of all N‐cycle functional genes. In conclusion, root biomass played key roles in affecting bacterial community structure and soil N cycling. Our findings have important implications for our understanding of plants control over bacterial community and N‐cycling function in artificial forest ecosystems.     

Exercise reduces circulating biomarkers of cellular senescence in humans

Cellular senescence has emerged as a significant and potentially tractable mechanism of aging and multiple aging‐related conditions. Biomarkers of senescent cell burden, including molecular signals in circulating immune cells and the abundance of circulating senescence‐related proteins, have been associated with chronological age and clinical parameters of biological age in humans. The extent to which senescence biomarkers are affected by interventions that enhance health and function has not yet been examined. Here, we report that a 12‐week structured exercise program drives significant improvements in several performance‐based and self‐reported measures of physical function in older adults. Impressively, the expression of key markers of the senescence program, including p16, p21, cGAS, and TNFα, were significantly lowered in CD3⁺ T cells in response to the intervention, as were the circulating concentrations of multiple senescence‐related proteins. Moreover, partial least squares discriminant analysis showed levels of senescence‐related proteins at baseline were predictive of changes in physical function in response to the exercise intervention. Our study provides first‐in‐human evidence that biomarkers of senescent cell burden are significantly lowered by a structured exercise program and predictive of the adaptive response to exercise.     

PPARα alleviates iron overload‐induced ferroptosis in mouse liver

Ferroptosis is an iron‐dependent form of non‐apoptotic cell death implicated in liver, brain, kidney, and heart pathology. How ferroptosis is regulated remains poorly understood. Here, we show that PPARα suppresses ferroptosis by promoting the expression of glutathione peroxidase 4 (Gpx4) and by inhibiting the expression of the plasma iron carrier TRF. PPARα directly induces Gpx4 expression by binding to a PPRE element within intron 3. PPARα knockout mice develop more severe iron accumulation and ferroptosis in the liver when fed a high‐iron diet than wild‐type mice. Ferrous iron (Fe²⁺) triggers ferroptosis via Fenton reactions and ROS accumulation. We further find that a rhodamine‐based "turn‐on" fluorescent probe(probe1) is suitable for the in vivo detection of Fe²⁺. Probe1 displays high selectivity towards Fe²⁺, and exhibits a stable response for Fe²⁺ with a concentration of 20 μM in tissue. Our data thus show that PPARα activation alleviates iron overload‐induced ferroptosis in mouse livers through Gpx4 and TRF, suggesting that PPARα may be a promising therapeutic target for drug discovery in ferroptosis‐related tissue injuries. Moreover, we identified a fluorescent probe that specifically labels ferrous ions and can be used to monitor Fe²⁺ in vivo.     

Genetic prediction of ICU hospitalization and mortality in COVID‐19 patients using artificial neural networks

There is an unmet need of models for early prediction of morbidity and mortality of Coronavirus disease‐19 (COVID‐19). We aimed to a) identify complement‐related genetic variants associated with the clinical outcomes of ICU hospitalization and death, b) develop an artificial neural network (ANN) predicting these outcomes and c) validate whether complement‐related variants are associated with an impaired complement phenotype. We prospectively recruited consecutive adult patients of Caucasian origin, hospitalized due to COVID‐19. Through targeted next‐generation sequencing, we identified variants in complement factor H/CFH, CFB, CFH‐related, CFD, CD55, C3, C5, CFI, CD46, thrombomodulin/THBD, and A Disintegrin and Metalloproteinase with Thrombospondin motifs (ADAMTS13). Among 381 variants in 133 patients, we identified 5 critical variants associated with severe COVID‐19: rs2547438 (C3), rs2250656 (C3), rs1042580 (THBD), rs800292 (CFH) and rs414628 (CFHR1). Using age, gender and presence or absence of each variant, we developed an ANN predicting morbidity and mortality in 89.47% of the examined population. Furthermore, THBD and C3a levels were significantly increased in severe COVID‐19 patients and those harbouring relevant variants. Thus, we reveal for the first time an ANN accurately predicting ICU hospitalization and death in COVID‐19 patients, based on genetic variants in complement genes, age and gender. Importantly, we confirm that genetic dysregulation is associated with impaired complement phenotype.     

SOX2 plays a crucial role in cell proliferation and lineage segregation during porcine pre‐implantation embryo development

ObjectivesGene regulation in early embryos has been widely studied for a long time because lineage segregation gives rise to the formation of a pluripotent cell population, known as the inner cell mass (ICM), during pre‐implantation embryo development. The extraordinarily longer pre‐implantation embryo development in pigs leads to the distinct features of the pluripotency network compared with mice and humans. For these reasons, a comparative study using pre‐implantation pig embryos would provide new insights into the mammalian pluripotency network and help to understand differences in the roles and networks of genes in pre‐implantation embryos between species.Materials and methodsTo analyse the functions of SOX2 in lineage segregation and cell proliferation, loss‐ and gain‐of‐function studies were conducted in pig embryos using an overexpression vector and the CRISPR/Cas9 system. Then, we analysed the morphological features and examined the effect on the expression of downstream genes through immunocytochemistry and quantitative real‐time PCR.ResultsOur results showed that among the core pluripotent factors, only SOX2 was specifically expressed in the ICM. In SOX2‐disrupted blastocysts, the expression of the ICM‐related genes, but not OCT4, was suppressed, and the total cell number was also decreased. Likewise, according to real‐time PCR analysis, pluripotency‐related genes, excluding OCT4, and proliferation‐related genes were decreased in SOX2‐targeted blastocysts. In SOX2‐overexpressing embryos, the total blastocyst cell number was greatly increased but the ICM/TE ratio decreased.ConclusionsTaken together, our results demonstrated that SOX2 is essential for ICM formation and cell proliferation in porcine early‐stage embryogenesis.     

3,6'‐disinapoyl sucrose attenuates Aβ1‐42‐ induced neurotoxicity in Caenorhabditis elegans by enhancing antioxidation and regulating autophagy

The aggregation of β‐amyloid (Aβ) has the neurotoxicity, which is thought to play critical role in the pathogenesis of Alzheimer''s disease (AD). Inhibiting Aβ deposition and neurotoxicity has been considered as an important strategy for AD treatment. 3,6''‐Disinapoyl sucrose (DISS), one of the oligosaccharide esters derived from traditional Chinese medicine Polygalae Radix, possesses antioxidative activity, neuroprotective effect and anti‐depressive activity. This study was to explore whether DISS could attenuate the pathological changes of Aβ_1‐42 transgenic Caenorhabditis elegans (C. elegans). The results showed that DISS (5 and 50 μM) treatment significantly prolonged the life span, increased the number of egg‐laying, reduced paralysis rate, decreased the levels of lipofuscin and ROS and attenuated Aβ deposition in Aβ_1‐42 transgenic C. elegans. Gene analysis showed that DISS could up‐regulate the mRNA expression of sod‐3, gst‐4, daf‐16, bec‐1 and lgg‐1, while down‐regulate the mRNA expression of daf‐2 and daf‐15 in Aβ_1‐42 transgenic C. elegans. These results suggested that DISS has the protective effect against Aβ_1‐42‐induced pathological damages and prolongs the life span of C. elegans, which may be related to the reduction of Aβ deposition and neurotoxicity by regulating expression of genes related to antioxidation and autophagy.