Falco: high-speed FastQC emulation for quality control of sequencing data

Study of pathogenic genes in a pedigree with familial dilated cardiomyopathy

BACKGROUND

Dilated cardiomyopathy (DCM) is a genetically heterogeneous cardiac disorder characterized by left ventricular dilation and contractile dysfunction. The substantial genetic heterogeneity evident in patients with DCM contributes to variable disease severity and complicates overall prognosis, which can be very poor.

AIM

To identify pathogenic genes in DCM through pedigree analysis.

METHODS

Our research team identified a patient with DCM in the clinic. Through investigation, we found that the family of this patient has a typical DCM pedigree. High-throughput sequencing technology, next-generation sequencing, was used to sequence the whole exomes of seven samples in the pedigree.

RESULTS

A novel and potentially pathogenic gene mutation-ANK2p.F3067L-was discovered. The mutation was completely consistent with the clinical information for this DCM pedigree. Sanger sequencing was used to further verify the locus of the mutation in pedigree samples. These results were consistent with those of high-throughput sequencing.

CONCLUSIONS

ANK2p.F3067L is considered a novel and potentially pathogenic gene mutation in DCM.

Whole Transcriptome Analysis of Differentially Expressed Genes in Cultured Nile Tilapia (O. niloticus) Subjected to Chronic Stress Reveals Signaling Pathways Associated with Depressed Growth

Chronic stress is a serious threat to aquaculture as it lowers fish growth performance and compromises fish welfare. The exact mechanism by which growth is retarded is, however, not clearly understood. This study sought to elucidate the gene expression profiles associated with chronic stress in cultured Nile tilapia (Oreochromis niloticus) reared for 70 days at different ammonia concentrations and stocking densities. Fish in the treatment groups showed negative growth, while the controls showed positive allometric growth. The specific condition factor (Kn) ranged from 1.17 for the controls to 0.93 for the ammonia and 0.91 for the stocking density treatments. RNA was extracted from muscle tissue using TRIzol followed by library construction and Illumina sequencing. Comparative transcriptome analysis revealed 209 differentially expressed genes (DEGs) (156 up- and 53 down-regulated) in the ammonia and 252 DEGs (175 up- and 77 down-regulated) in the stocking density treatment. In both treatments, 24 and 17 common DEGs were up- and down-regulated, respectively. DEGs were significantly enriched in six pathways associated with muscle activity, energy mobilization and immunity. The heightened muscular activity consumes energy which would otherwise have been utilized for growth. These results bring to fore the molecular mechanisms underlying chronic stress’ suppression of growth in cultured Nile tilapia.

Genome Characterization of Bacteriophage KPP-1, a Novel Member in the Subfamily Vequintavirinae, and Use of Its Endolysin for the Lysis of Multidrug-Resistant Klebsiella variicola In Vitro

Multidrug-resistant members of the Klebsiella pneumoniae complex have become a threat to human lives and animals, including aquatic animals, owing to the limited choice of antimicrobial treatments. Bacteriophages are effective natural tools available to fight against multidrug-resistant bacteria. The bacteriophage KPP-1 was found to be strictly lytic against K. variicola, a multidrug-resistant isolate, producing clear plaques. The genome sequence analysis of KPP-1 revealed that it comprised 143,369 base pairs with 47% overall GC content. A total of 272 genes (forward 161, complementary 111) encode for 17 tRNAs and 255 open reading frames (ORFs). Among them, 32 ORFs could be functionally annotated using the National Center for Biotechnology Information (NCBI) Protein Basic Local Alignment Search Tool (BLASTp) algorithm while 223 were found to code for hypothetical proteins. Comparative genomic analysis revealed that the closest neighbor of KPP-1 can be found in the genus Mydovirus of the subfamily Vequintavirinae. KPP-1 not only markedly suppressed the growth of the host but also worked synergistically with ampicillin. Useful genes for pathogen control such as endolysin (locus tag: KPP_11591) were found to have activity against multidrug-resistant isolate of K. variicola. Further studies are necessary to develop a strategy to control the emerging pathogen K. variicola using bacteriophages such as KPP-1.

Praliciguat Promotes Ischemic Leg Reperfusion in Leptin Receptor-Deficient Mice

BACKGROUND

Lower-limb peripheral artery disease is one of the major complications of diabetes. Peripheral artery disease is associated with poor limb and cardiovascular prognoses, along with a dramatic decrease in life expectancy. Despite major medical advances in the treatment of diabetes, a substantial therapeutic gap remains in the peripheral artery disease population. Praliciguat is an orally available sGC (soluble guanylate cyclase) stimulator that has been reported both preclinically and in early stage clinical trials to have favorable effects in metabolic and hemodynamic outcomes, suggesting that it may have a potential beneficial effect in peripheral artery disease.

METHODS

We evaluated the effect of praliciguat on hind limb ischemia recovery in a mouse model of type 2 diabetes. Hind limb ischemia was induced in leptin receptor-deficient (Lepr^db/db) mice by ligation and excision of the left femoral artery. Praliciguat (10 mg/kg/day) was administered in the diet starting 3 days before surgery.

RESULTS

Twenty-eight days after surgery, ischemic foot perfusion and function parameters were better in praliciguat-treated mice than in vehicle controls. Improved ischemic foot perfusion was not associated with either improved traditional cardiovascular risk factors (ie, weight, glycemia) or increased angiogenesis. However, treatment with praliciguat significantly increased arteriole diameter, decreased ICAM1 (intercellular adhesion molecule 1) expression, and prevented the accumulation of oxidative proangiogenic and proinflammatory muscle fibers. While investigating the mechanism underlying the beneficial effects of praliciguat therapy, we found that praliciguat significantly downregulated Myh2 and Cxcl12 mRNA expression in cultured myoblasts and that conditioned medium form praliciguat-treated myoblast decreased ICAM1 mRNA expression in endothelial cells. These results suggest that praliciguat therapy may decrease ICAM1 expression in endothelial cells by downregulating Cxcl12 in myocytes.

CONCLUSIONS

Our results demonstrated that praliciguat promotes blood flow recovery in the ischemic muscle of mice with type 2 diabetes, at least in part by increasing arteriole diameter and by downregulating ICAM1 expression.

The Histone Chaperone Network Is Highly Conserved in Physarum polycephalum

The nucleosome is composed of histones and DNA. Prior to their deposition on chromatin, histones are shielded by specialized and diverse proteins known as histone chaperones. They escort histones during their entire cellular life and ensure their proper incorporation in chromatin. Physarum polycephalum is a Mycetozoan, a clade located at the crown of the eukaryotic tree. We previously found that histones, which are highly conserved between plants and animals, are also highly conserved in Physarum. However, histone chaperones differ significantly between animal and plant kingdoms, and this thus probed us to further study the conservation of histone chaperones in Physarum and their evolution relative to animal and plants. Most of the known histone chaperones and their functional domains are conserved as well as key residues required for histone and chaperone interactions. Physarum is divergent from yeast, plants and animals, but PpHIRA, PpCABIN1 and PpSPT6 are similar in structure to plant orthologues. PpFACT is closely related to the yeast complex, and the Physarum genome encodes the animal-specific APFL chaperone. Furthermore, we performed RNA sequencing to monitor chaperone expression during the cell cycle and uncovered two distinct patterns during S-phase. In summary, our study demonstrates the conserved role of histone chaperones in handling histones in an early-branching eukaryote.

Nitrogen-dependent binding of the transcription factor PBF1 contributes to the balance of protein and carbohydrate storage in maize endosperm

Fluctuations in nitrogen (N) availability influence protein and starch levels in maize (Zea mays) seeds, yet the underlying mechanism is not well understood. Here, we report that N limitation impacted the expression of many key genes in N and carbon (C) metabolism in the developing endosperm of maize. Notably, the promoter regions of those genes were enriched for P-box sequences, the binding motif of the transcription factor prolamin-box binding factor 1 (PBF1). Loss of PBF1 altered accumulation of starch and proteins in endosperm. Under different N conditions, PBF1 protein levels remained stable but PBF1 bound different sets of target genes, especially genes related to the biosynthesis and accumulation of N and C storage products. Upon N-starvation, the absence of PBF1 from the promoters of some zein genes coincided with their reduced expression, suggesting that PBF1 promotes zein accumulation in the endosperm. In addition, PBF1 repressed the expression of sugary1 (Su1) and starch branching enzyme 2b (Sbe2b) under normal N supply, suggesting that, under N-deficiency, PBF1 redirects the flow of C skeletons for zein toward the formation of C compounds. Overall, our study demonstrates that PBF1 modulates C and N metabolism during endosperm development in an N-dependent manner.

A Novel Hypomorphic Apex1 Mouse Model Implicates Apurinic/Apyrimidinic Endonuclease 1 in Oxidative DNA Damage Repair in Gastric Epithelial Cells

Aims: Though best known for its role in oxidative DNA damage repair, apurinic/apyrimidinic endonuclease 1 (APE1) is a multifunctional protein that regulates multiple host responses during oxidative stress, including the reductive activation of transcription factors. As knockout of the APE1-encoding gene, Apex1, is embryonically lethal, we sought to create a viable model with generalized inhibition of APE1 expression. Results: A hypomorphic (HM) mouse with decreased APE1 expression throughout the body was generated using a construct containing a neomycin resistance (NeoR) cassette knocked into the Apex1 site. Offspring were assessed for APE1 expression, breeding efficiency, and morphology with a focused examination of DNA damage in the stomach. Heterozygotic breeding pairs yielded 50% fewer HM mice than predicted by Mendelian genetics. APE1 expression was reduced up to 90% in the lungs, heart, stomach, and spleen. The HM offspring were typically smaller, and most had a malformed tail. Oxidative DNA damage was increased spontaneously in the stomachs of HM mice. Further, all changes were reversed when the NeoR cassette was removed. Primary gastric epithelial cells from HM mice differentiated more quickly and had more evidence of oxidative DNA damage after stimulation with Helicobacter pylori or a chemical carcinogen than control lines from wildtype mice. Innovation: A HM mouse with decreased APE1 expression throughout the body was generated and extensively characterized. Conclusion: The results suggest that HM mice enable studies of APE1's multiple functions throughout the body. The detailed characterization of the stomach showed that gastric epithelial cells from HM were more susceptible to DNA damage. Antioxid. Redox Signal. 38, 183-197.

Genetic Variants Associated with Bronchial Asthma Specific to the Population of the Russian Federation

Bronchial asthma (BA) is a disease that still lacks an exhaustive treatment protocol. In this regard, the global medical community pays special attention to the genetic prerequisites for the occurrence of this disease. Therefore, the search for the genetic polymorphisms underlying bronchial asthma has expanded considerably. As the present study progressed, a significant amount of scientific medical literature was analyzed and 167 genes reported to be associated with the development of bronchial asthma were identified. A group of participants (n = 7,303) who had voluntarily provided their biomaterial (venous blood) to be used in the research conducted by the Federal Medical Biological Agency of Russia was formed to subsequently perform a bioinformatic verification of known associations and search for new ones. This group of participants was divided into four cohorts, including two sex-distinct cohorts of individuals with a history of asthma and two sex-distinct cohorts of apparently healthy individuals. A search for polymorphisms was made in each cohort among the selected genes, and genetic variants were identified whose difference in occurrence in the different cohorts was statistically significant (significance level less than 0.0001). The study revealed 11 polymorphisms that affect the development of asthma: four genetic variants (rs869106717, rs1461555098, rs189649077, and rs1199362453), which are more common in men with bronchial asthma compared to apparently healthy men; five genetic variants (rs1923038536, rs181066119, rs143247175, rs140597386, and rs762042586), which are more common in women with bronchial asthma compared to apparently healthy women; and two genetic variants (rs1219244986 and rs2291651) that are rare in women with a history of asthma.

Chronic exposure to ampicillin alters lung microbial composition in laboratory rat

PURPOSE

High-throughput sequencing technologies have revealed that the lungs contain a variety of low biomass microbiota associated with various lung diseases. Rat model is an important tool to understand the possible causal relationship between pulmonary microbiota and diseases. Antibiotic exposure can alter the microbiota, however, a direct influence of long-term ampicillin exposure on commensal bacteria of healthy lungs has not been investigated, which could be useful in the study of the relation between microbiome and long-term lung diseases, especially in animal model-making of lung diseases.

METHODS

The rats were aerosolized ampicillin of different concentrations for five months, and then the effect on the lung microbiota was investigated using 16S rRNA gene sequencing.

RESULTS

The ampicillin treatment by a certain concentration (LA5, 0.2 ml of 5 mg/ml ampicillin) administration leads to profound changes in the rat lung microbiota but not in the low critical ampicillin concentration (LA01 and LA1, 0.1 and 1 mg/ml ampicillin), when compared to the untreated group (LC). The genus Acidobacteria_Gp16 dominated the ampicillin treated lung microbiota while the genera Brucella, Acinetobacter, Acidobacteria_Gp14, Sphingomonas, and Tumebacillus dominated the untreated lung microbiota. The predicted KEGG pathway analysis profile revealed some difference in the ampicillin treated group.

CONCLUSIONS

The study demonstrated the effects of different concentrations of ampicillin treatment on lung microbiota of rats in a relatively long term. It could serve as a basis for the clinical use of antibiotic and the use of ampicillin to control certain bacteria in the animal model-making of respiratory diseases such as chronic obstructive pulmonary disease.

Analysis of the response regulatory network of pepper genes under hydrogen peroxide stress

Hydrogen peroxide (H₂O₂) is a regulatory component related to plant signal transduction. To better understand the genome-wide gene expression response to H₂O₂ stress in pepper plants, a regulatory network of H₂O₂ stress-gene expression in pepper leaves and roots was constructed in the present study. We collected the normal tissues of leaves and roots of pepper plants after 40 days of H₂O₂ treatment and obtained the RNA-seq data of leaves and roots exposed to H₂O₂ for 0.5-24 h. By comparing the gene responses of pepper leaves and roots exposed to H₂O₂ stress for different time periods, we found that the response in roots reached the peak at 3 h, whereas the response in leaves reached the peak at 24 h after treatment, and the response degree in the roots was higher than that in the leaves. We used all datasets for K-means analysis and network analysis identified the clusters related to stress response and related genes. In addition, CaEBS1, CaRAP2, and CabHLH029 were identified through a co-expression analysis and were found to be strongly related to several reactive oxygen species-scavenging enzyme genes; their homologous genes in Arabidopsis showed important functions in response to hypoxia or iron uptake. This study provides a theoretical basis for determining the dynamic response process of pepper plants to H₂O₂ stress in leaves and roots, as well as for determining the critical time and the molecular mechanism of H₂O₂ stress response in leaves and roots. The candidate transcription factors identified in this study can be used as a reference for further experimental verification.

Postmortem submersion interval estimation of cadavers recovered from freshwater based on gut microbial community succession

Microbial community succession during decomposition has been proven to be a useful tool for postmortem interval (PMI) estimation. Numerous studies have shown that the intestinal microbial community presented chronological changes after death and was stable in terrestrial corpses with different causes of death. However, the postmortem pattern of intestinal microbial community succession in cadavers retrieved from water remains unclear. For immersed corpses, the postmortem submersion interval (PMSI) is a useful indicator of PMI. To provide reliable estimates of PMSI in forensic investigations, we investigated the gut microbial community succession of corpses submersed in freshwater and explored its potential application in forensic investigation. In this study, the intestinal microbial community of mouse submersed in freshwater that died of drowning or CO2 asphyxia (i.e., postmortem submersion) were characterized by 16S rDNA amplification and high-throughput sequencing, followed by bioinformatic analyses. The results demonstrated that the chronological changes in intestinal bacterial communities were not different between the drowning and postmortem submersion groups. α-diversity decreased significantly within 14 days of decomposition in both groups, and the β-diversity bacterial community structure ordinated chronologically, inferring the functional pathway and phenotype. To estimate PMSI, a regression model was established by random forest (RF) algorithm based on the succession of postmortem microbiota. Furthermore, 15 genera, including Proteus, Enterococcus, and others, were selected as candidate biomarkers to set up a concise predicted model, which provided a prediction of PMSI [MAE (± SE) = 0.818 (± 0.165) d]. Overall, our present study provides evidence that intestinal microbial community succession would be a valuable marker to estimate the PMSI of corpses submerged in an aquatic habitat.

RvD1n-3 DPA Downregulates the Transcription of Pro-Inflammatory Genes in Oral Epithelial Cells and Reverses Nuclear Translocation of Transcription Factor p65 after TNF-α Stimulation

Specialized pro-resolving mediators (SPMs) are multifunctional lipid mediators that participate in the resolution of inflammation. We have recently described that oral epithelial cells (OECs) express receptors of the SPM resolvin RvD1_{n-3 DPA} and that cultured OECs respond to RvD1_{n-3 DPA} addition by intracellular calcium release, nuclear receptor translocation and transcription of genes coding for antimicrobial peptides. The aim of the present study was to assess the functional outcome of RvD1_{n-3 DPA}-signaling in OECs under inflammatory conditions. To this end, we performed transcriptomic analyses of TNF-α-stimulated cells that were subsequently treated with RvD1_{n-3 DPA} and found significant downregulation of pro-inflammatory nuclear factor kappa B (NF-κB) target genes. Further bioinformatics analyses showed that RvD1_{n-3 DPA} inhibited the expression of several genes involved in the NF-κB activation pathway. Confocal microscopy revealed that addition of RvD1_{n-3 DPA} to OECs reversed TNF-α-induced nuclear translocation of NF-κB p65. Co-treatment of the cells with the exportin 1 inhibitor leptomycin B indicated that RvD1_{n-3 DPA} increases nuclear export of p65. Taken together, our observations suggest that SPMs also have the potential to be used as a therapeutic aid when inflammation is established.

Microbial communities in the liver and brain are informative for postmortem submersion interval estimation in the late phase of decomposition: A study in mouse cadavers recovered from freshwater

INTRODUCTION

Bodies recovered from water, especially in the late phase of decomposition, pose difficulties to the investigating authorities. Various methods have been proposed for postmortem submersion interval (PMSI) estimation and drowning identification, but some limitations remain. Many recent studies have proved the value of microbiota succession in viscera for postmortem interval estimation. Nevertheless, the visceral microbiota succession and its application for PMSI estimation and drowning identification require further investigation.

METHODS

In the current study, mouse drowning and CO2 asphyxia models were developed, and cadavers were immersed in freshwater for 0 to 14 days. Microbial communities in the liver and brain were characterized via 16S rDNA high-throughput sequencing.

RESULTS

Only livers and brains collected from 5 to 14 days postmortem were qualified for sequencing. There was significant variation between microbiota from liver and brain. Differences in microbiota between the cadavers of mice that had drowned and those only subjected to postmortem submersion decreased over the PMSI. Significant successions in microbial communities were observed among the different subgroups within the late phase of the PMSI in livers and brains. Eighteen taxa in the liver which were mainly related to Clostridium_sensu_stricto and Aeromonas, and 26 taxa in the brain which were mainly belonged to Clostridium_sensu_stricto, Acetobacteroides, and Limnochorda, were selected as potential biomarkers for PMSI estimation based on a random forest algorithm. The PMSI estimation models established yielded accurate prediction results with mean absolute errors ± the standard error of 1.282 ± 0.189 d for the liver and 0.989 ± 0.237 d for the brain.

CONCLUSIONS

The present study provides novel information on visceral postmortem microbiota succession in corpses submerged in freshwater which sheds new light on PMSI estimation based on the liver and brain in forensic practice.

Type I interferon receptor (IFNAR2) deficiency reveals Zika virus cytopathicity in human macrophages and microglia

Macrophages are key target cells of Zika virus (ZIKV) infection, implicated as a viral reservoir seeding sanctuary sites such as the central nervous system and testes. This rests on the apparent ability of macrophages to sustain ZIKV replication without experiencing cytopathic effects. ZIKV infection of macrophages triggers an innate immune response involving type I interferons (IFN-I), key antiviral cytokines that play a complex role in ZIKV pathogenesis in animal models. To investigate the functional role of the IFN-I response we generated human induced pluripotent stem cell (iPSC)-derived macrophages from a patient with complete deficiency of IFNAR2, the high affinity IFN-I receptor subunit. Accompanying the profound defect of IFN-I signalling in IFNAR2 deficient iPS-macrophages we observed significantly enhanced ZIKV replication and cell death, revealing the inherent cytopathicity of ZIKV towards macrophages. These observations were recapitulated by genetic and pharmacological ablation of IFN-I signalling in control iPS-macrophages and extended to a model of iPS-microglia. Thus, the capacity of macrophages to support noncytolytic ZIKV replication depends on an equilibrium set by IFN-I, suggesting that innate antiviral responses might counterintuitively promote ZIKV persistence via the maintenance of tissue viral reservoirs relevant to pathogenesis.

Transcriptomic and Proteomic Profiles for Elucidating Cisplatin Resistance in Head-and-Neck Squamous Cell Carcinoma

To identify the novel genes involved in chemoresistance in head and neck squamous cell carcinoma (HNSCC), we explored the expression profiles of the following cisplatin (CDDP) resistant (R) versus parental (sensitive) cell lines by RNA-sequencing (RNA-seq): JHU029, HTB-43 and CCL-138. Using the parental condition as a control, 30 upregulated and 85 downregulated genes were identified for JHU029-R cells; 263 upregulated and 392 downregulated genes for HTB-43-R cells, and 154 upregulated and 68 downregulated genes for CCL-138-R cells. Moreover, we crossed-checked the RNA-seq results with the proteomic profiles of HTB-43-R (versus HTB-43) and CCL-138-R (versus CCL-138) cell lines. For the HTB-43-R cells, 21 upregulated and 72 downregulated targets overlapped between the proteomic and transcriptomic data; whereas in CCL-138-R cells, four upregulated and three downregulated targets matched. Following an extensive literature search, six genes from the RNA-seq (CLDN1, MAGEB2, CD24, CEACAM6, IL1B and ISG15) and six genes from the RNA-seq and proteomics crossover (AKR1C3, TNFAIP2, RAB7A, LGALS3BP, PSCA and SSRP1) were selected to be studied by qRT-PCR in 11 HNSCC patients: six resistant and five sensitive to conventional therapy. Interestingly, the high MAGEB2 expression was associated with resistant tumours and is revealed as a novel target to sensitise resistant cells to therapy in HNSCC patients.

Increased AR expression in castration-resistant prostate cancer rapidly induces AR signaling reprogramming with the collaboration of EZH2

Elevated androgen receptor (AR) expression is a hallmark of castration-resistant prostate cancer (CRPC) and contributes to the restoration of AR signaling under the conditions of androgen deprivation. However, whether overexpressed AR alone with the stimulation of castrate levels of androgens can be sufficient to induce the reprogramming of AR signaling for the adaptation of prostate cancer (PCa) cells remains unclear. In this study, we used a PCa model with inducible overexpression of AR to examine the acute effects of AR overexpression on its cistrome and transcriptome. Our results show that overexpression of AR alone in conjunction with lower androgen levels can rapidly redistribute AR chromatin binding and activates a distinct transcription program that is enriched for DNA damage repair pathways. Moreover, using a recently developed bioinformatic tool, we predicted the involvement of EZH2 in this AR reprogramming and subsequently identified a subset of AR/EZH2 co-targeting genes, which are overexpressed in CRPC and associated with worse patient outcomes. Mechanistically, we found that AR-EZH2 interaction is impaired by the pre-castration level of androgens but can be recovered by the post-castration level of androgens. Overall, our study provides new molecular insights into AR signaling reprogramming with the engagement of specific epigenetic factors.

Complete Genome Sequence of Pseudomonas Phage Motto

We describe the complete genome sequence of bacteriophage Motto, which infects clinical strains of Pseudomonas aeruginosa. Motto is a T1-like siphovirus related to members of the family Drexlerviridae and has a capsid width of ~57 nm and a tail length of ~255 nm. The 49.9-kb genome contains 84 protein-coding genes.

Activation function 1 of progesterone receptor is required for progesterone antagonism of oestrogen action in the uterus

BACKGROUND

Progesterone receptor (PGR) is a master regulator of uterine function through antagonistic and synergistic interplays with oestrogen receptors. PGR action is primarily mediated by activation functions AF1 and AF2, but their physiological significance is unknown.

RESULTS

We report the first study of AF1 function in mice. The AF1 mutant mice are infertile with impaired implantation and decidualization. This is associated with a delay in the cessation of epithelial proliferation and in the initiation of stromal proliferation at preimplantation. Despite tissue selective effect on PGR target genes, AF1 mutations caused global loss of the antioestrogenic activity of progesterone in both pregnant and ovariectomized models. Importantly, the study provides evidence that PGR can exert an antioestrogenic effect by genomic inhibition of Esr1 and Greb1 expression. ChIP-Seq data mining reveals intermingled PGR and ESR1 binding on Esr1 and Greb1 gene enhancers. Chromatin conformation analysis shows reduced interactions in these genes' loci in the mutant, coinciding with their upregulations.

CONCLUSION

AF1 mediates genomic inhibition of ESR1 action globally whilst it also has tissue-selective effect on PGR target genes.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

Transcriptome-wide m6A methylome analysis uncovered the changes of m6A modification in oral pre-malignant cells compared with normal oral epithelial cells

As the most common post-transcriptional RNA modification, m⁶A methylation extensively regulates the structure and function of RNA. The dynamic and reversible modification of m⁶A is coordinated by m⁶A writers and erasers. m⁶A reader proteins recognize m⁶A modification on RNA, mediating different downstream biological functions. mRNA m⁶A modification and its corresponding regulators play an important role in cancers, but its characteristics in the precancerous stage are still unclear. In this study, we used oral precancerous DOK cells as a model to explore the characteristics of transcriptome-wide m⁶A modification and major m⁶A regulator expression in the precancerous stage compared with normal oral epithelial cell HOEC and oral cancer cell SCC-9 through MeRIP-seq and RT-PCR. Compared with HOEC cells, we found 1180 hyper-methylated and 1606 hypo-methylated m⁶A peaks and 354 differentially expressed mRNAs with differential m⁶A peaks in DOK cells. Although the change of m⁶A modification in DOK cells was less than that in SCC-9 cells, mRNAs with differential m⁶A in both cell lines were enriched into many identical GO terms and KEGG pathways. Among the 20 known m⁶A regulatory genes, FTO, ALKBH5, METTL3 and VIRMA were upregulated or downregulated in DOK cells, and the expression levels of 10 genes such as METTL14/16, FTO and IGF2BP2/3 were significantly changed in SCC-9 cells. Our data suggest that precancerous cells showed, to some extent, changes of m⁶A modification. Identifying some key m⁶A targets and corresponding regulators in precancerous stage may provide potential intervention targets for the prevention of cancer development through epigenetic modification in the future.

Decoding brain memory formation by single-cell RNA sequencing

To understand how distinct memories are formed and stored in the brain is an important and fundamental question in neuroscience and computational biology. A population of neurons, termed engram cells, represents the physiological manifestation of a specific memory trace and is characterized by dynamic changes in gene expression, which in turn alters the synaptic connectivity and excitability of these cells. Recent applications of single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) are promising approaches for delineating the dynamic expression profiles in these subsets of neurons, and thus understanding memory-specific genes, their combinatorial patterns and regulatory networks. The aim of this article is to review and discuss the experimental and computational procedures of sc/snRNA-seq, new studies of molecular mechanisms of memory aided by sc/snRNA-seq in human brain diseases and related mouse models, and computational challenges in understanding the regulatory mechanisms underlying long-term memory formation.

Crosstalk between mucosal microbiota, host gene expression, and sociomedical factors in the progression of colorectal cancer

Various omics-based biomarkers related to the occurrence, progression, and prognosis of colorectal cancer (CRC) have been identified. In this study, we attempted to identify gut microbiome-based biomarkers and detect their association with host gene expression in the initiation and progression of CRC by integrating analysis of the gut mucosal metagenome, RNA sequencing, and sociomedical factors. We performed metagenome and RNA sequencing on colonic mucosa samples from 13 patients with advanced CRC (ACRC), 10 patients with high-risk adenoma (HRA), and 7 normal control (NC) individuals. All participants completed a questionnaire on sociomedical factors. The interaction and correlation between changes in the microbiome and gene expression were assessed using bioinformatic analysis. When comparing HRA and NC samples, which can be considered to represent the process of tumor initiation, 28 genes and five microbiome species were analyzed with correlation plots. When comparing ACRC and HRA samples, which can be considered to represent the progression of CRC, seven bacterial species and 21 genes were analyzed. When comparing ACRC and NC samples, 16 genes and five bacterial species were analyzed, and four correlation plots were generated. A network visualizing the relationship between bacterial and host gene expression in the initiation and progression of CRC indicated that Clostridium spiroforme and Tyzzerella nexilis were hub bacteria in the development and progression of CRC. Our study revealed the interactions of and correlation between the colonic mucosal microbiome and host gene expression to identify potential roles of the microbiome in the initiation and progression of CRC. Our results provide gut microbiome-based biomarkers that may be potential diagnostic markers and therapeutic targets in patients with CRC.

RNA Sequencing Reveals the Expression Profiles of circRNAs and Indicates Hsa_circ_0070562 as a Pro-osteogenic Factor in Bone Marrow-Derived Mesenchymal Stem Cells of Patients With Ankylosing Spondylitis

Recent studies have reported that circular RNAs (circRNAs) play a crucial regulatory role in a variety of human diseases. However, the roles of circRNAs in pathological osteogenesis in ankylosing spondylitis (AS) remain unclear. We conducted circRNA and miRNA expression profiling of osteogenically differentiated bone marrow-derived mesenchymal stem cells (BMSCs) of patients with AS compared with those of healthy donors (HDs) by RNA sequencing (RNA-seq). Results showed that a total of 31806 circRNAs were detected in the BMSC samples, of which 418 circRNAs were significantly differentially expressed (DE) with a fold change ≥2 and p value <0.05. Among these, 204 circRNAs were upregulated, and 214 were downregulated. GO and KEGG analyses demonstrated that the DE circRNAs were mainly involved in the regulation of biological processes of the cell matrix adhesion and the TGF-beta signaling pathway, which are closely related to AS. circRNA-miRNA interaction networks related to the TGF-beta signaling pathway were established. The results of qRT-PCR showed that has_circ_0070562 was significantly up-regulated in AS-MSCs. In vitro experiments showed that silencing of has_circ_0070562 weakened osteogenesis of AS-BMSCs. In conclusion, we identified numerous circRNAs that were dysregulated in AS-BMSCs compared with HD-BMSCs. Bioinformatic analyses suggested that these dysregulated circRNAs might play important functional roles in AS-BMSCs osteogenesis. Circ_0070562 functioned as a pro-ostegenic factor and might serve as a potential biomarker and a therapeutic target for AS.

MIR133A regulates cell proliferation, migration, and apoptosis by targeting SOX9 in human colorectal cancer cells

The human microRNA 133A (MIR133A) was identified as a CRC-associated miRNA. It was down-regulated in human CRC tissues. We identified the putative MIR133A1 and A2 target genes by comparing the transcriptome analysis data of MIR133A1 and A2 knock-in cells with the candidate MIR133A target genes predicted by bioinformatics tools. We identified 29 and 33 putative MIR133A and A2 direct target genes, respectively. Among them, we focused on the master transcription regulator gene SRY-box transcription factor 9 (SOX9), which exhibits a pleiotropic role in cancer. We confirmed that SOX9 is a direct target gene of MIR133A by luciferase reporter assay, quantitative RT-PCR, and western blot analysis. Overexpression of MIR133A in CRC cell lines significantly decreased SOX9 and its downstream PIK3CA-AKT1-GSK3B-CTNNB1 and KRAS-BRAF-MAP2K1-MAPK1/3 pathways and increased apoptosis. Furthermore, functional studies reveal that cell proliferation, colony formation, and migration ability were significantly decreased by MIR133A-overexpressed CRC cell lines. Knockdown of SOX9 in CRC cell lines by SOX9 gene silencing showed similar results. We also used a xenograft model to show that MIR133A overexpression suppresses tumor growth and proliferation. Our results suggest that MIR133A regulates cell proliferation, migration, and apoptosis by targeting SOX9 in human colorectal cancer.

Blood-derived lncRNAs as biomarkers for cancer diagnosis: the Good, the Bad and the Beauty

Cancer ranks as one of the deadliest diseases worldwide. The high mortality rate associated with cancer is partially due to the lack of reliable early detection methods and/or inaccurate diagnostic tools such as certain protein biomarkers. Cell-free nucleic acids (cfNA) such as circulating long noncoding RNAs (lncRNAs) have been proposed as a new class of potential biomarkers for cancer diagnosis. The reported correlation between the presence of tumors and abnormal levels of lncRNAs in the blood of cancer patients has notably triggered a worldwide interest among clinicians and oncologists who have been actively investigating their potentials as reliable cancer biomarkers. In this report, we review the progress achieved ("the Good") and challenges encountered ("the Bad") in the development of circulating lncRNAs as potential biomarkers for early cancer diagnosis. We report and discuss the diagnostic performance of more than 50 different circulating lncRNAs and emphasize their numerous potential clinical applications ("the Beauty") including therapeutic targets and agents, on top of diagnostic and prognostic capabilities. This review also summarizes the best methods of investigation and provides useful guidelines for clinicians and scientists who desire conducting their own clinical studies on circulating lncRNAs in cancer patients via RT-qPCR or Next Generation Sequencing (NGS).

Genetic Profiling of Glucocorticoid (NR3C1) and Mineralocorticoid (NR3C2) Receptor Polymorphisms before Starting Therapy with Androgen Receptor Inhibitors: A Study of a Patient Who Developed Toxic Myocarditis after Enzalutamide Treatment

Enzalutamide is a nonsteroidal inhibitor of the androgen receptor (AR) signaling pathway and is used to treat patients with metastatic castration-resistant prostate cancer. However, the risk of cardiovascular-related hospitalization in patients with no contraindications for the use of enzalutamide is about 1–2%. To date, the underlying molecular basis of this has not been established. The androgen receptor, glucocorticoid receptor (GR) and mineralocorticoid receptor (MR) are nuclear receptors that share structural similarities and have closely related DNA-binding sites and coregulators. In non-epithelial cells, a fine balance of the activities of these receptors is essential to ensure correct cellular function. In this study, we present a molecular characterization of these nuclear receptors in a prostate cancer patient who developed congestive heart failure after enzalutamide treatment. White cell RNAseq revealed a homozygous rs5522 MR polymorphism and both the rs143711342 and rs56149945 GR polymorphisms, carried in different alleles. No different specific splice isoforms were detected. Recent research suggests that AR inhibition by enzalutamide makes available a coregulator that specifically interacts with the rs5522-mutated MR, increasing its activity and producing adverse effects on cardiovascular health. We suggest an evaluation of the MR rs5522 polymorphism before starting therapy with AR inhibitors.

MRG Proteins Are Shared by Multiple Protein Complexes With Distinct Functions

MRG15/MORF4L1 is a highly conserved protein in eukaryotes that contains a chromodomain (CHD) recognizing methylation of lysine 36 on histone H3 (H3K36me3) in chromatin. Intriguingly, it has been reported in the literature to interact with several different factors involved in chromatin modifications, gene regulation, alternative mRNA splicing, and DNA repair by homologous recombination. To get a complete and reliable picture of associations in physiological conditions, we used genome editing and tandem affinity purification to analyze the stable native interactome of human MRG15, its paralog MRGX/MORF4L2 that lacks the CHD, and MRGBP (MRG-binding protein) in isogenic K562 cells. We found stable interchangeable association of MRG15 and MRGX with the NuA4/TIP60 histone acetyltransferase/chromatin remodeler, Sin3B histone deacetylase/demethylase, ASH1L histone methyltransferase, and PALB2-BRCA2 DNA repair protein complexes. These associations were further confirmed and analyzed by CRISPR tagging of endogenous proteins and comparison of expressed isoforms. Importantly, based on structural information, point mutations could be introduced that specifically disrupt MRG15 association with some complexes but not others. Most interestingly, we also identified a new abundant native complex formed by MRG15/X-MRGBP-BRD8-EP400NL (EP400 N-terminal like) that is functionally similar to the yeast TINTIN (Trimer Independent of NuA4 for Transcription Interactions with Nucleosomes) complex. Our results show that EP400NL, being homologous to the N-terminal region of NuA4/TIP60 subunit EP400, creates TINTIN by competing for BRD8 association. Functional genomics indicate that human TINTIN plays a role in transcription of specific genes. This is most likely linked to the H4ac-binding bromodomain of BRD8 along the H3K36me3-binding CHD of MRG15 on the coding region of transcribed genes. Taken together, our data provide a complete detailed picture of human MRG proteins-associated protein complexes, which are essential to understand and correlate their diverse biological functions in chromatin-based nuclear processes.

A Novel Proline-Rich Cathelicidin from the Alpaca Vicugna pacos with Potency to Combat Antibiotic-Resistant Bacteria: Mechanism of Action and the Functional Role of the C-Terminal Region

Over recent years, a growing number of bacterial species have become resistant to clinically relevant antibiotics. Proline-rich antimicrobial peptides (PrAMPs) having a potent antimicrobial activity and a negligible toxicity toward mammalian cells attract attention as new templates for the development of antibiotic drugs. Here, we mined genomes of all living Camelidae species and found a novel family of Bac7-like proline-rich cathelicidins which inhibited bacterial protein synthesis. The N-terminal region of a novel peptide from the alpaca Vicugna pacos named VicBac is responsible for inhibition of bacterial protein synthesis with an IC₅₀ value of 0.5 µM in the E. coli cell-free system whereas the C-terminal region allows the peptide to penetrate bacterial membranes effectively. We also found that the full-length VicBac did not induce bacterial resistance after a two-week selection experiment, unlike the N-terminal truncated analog, which depended on the SbmA transport system. Both pro- and anti-inflammatory action of VicBac and its N-terminal truncated variant on various human cell types was found by multiplex immunoassay. The presence of the C-terminal tail in the natural VicBac does not provide for specific immune-modulatory effects in vitro but enhances the observed impact compared with the truncated analog. The pronounced antibacterial activity of VicBac, along with its moderate adverse effects on mammalian cells, make this molecule a promising scaffold for the development of peptide antibiotics.

Integrative Analysis of lncRNA-miRNA-mRNA Regulatory Network Reveals the Key lncRNAs Implicated Potentially in the Differentiation of Adipocyte in Goats

Goats are popular in China because of their superior meat quality, delicate flesh, and unique flavor. Long noncoding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation of gene expression. However, the effects of lncRNAs on adipocyte differentiation in goat has not been fully elucidated yet. In this investigation, we performed RNA-Seq analysis of intramuscular and subcutaneous adipocytes from Jianzhou Daer goat before and after differentiation, including both intramuscular preadipocytes (IMPA) vs. intramuscular adipocytes (IMA) and subcutaneous preadipocytes (SPA) vs. subcutaneous adipocytes (SA). A total of 289.49 G clean reads and 12,519 lncRNAs were obtained from 20 samples. In total, 3,733 differentially expressed RNAs (182 lncRNAs and 3,551 mRNAs) were identified by pairwise comparison. There were 135 differentially expressed lncRNAs (DELs) specific to intramuscular adipocytes, 39 DELs specific to subcutaneous adipocytes, and 8 DELs common to both adipocytes in these 182 DELs. Some well-known and novel pathways associated with preadipocyte differentiation were identified: fat acid metabolism, TGF-beta signaling pathway and PI3K-Akt signaling pathway. By integrating miRNA-seq data from another study, we also identified hub miRNAs in both types of fat cells. Our analysis revealed the unique and common lncRNA-miRNA-mRNA networks of two kinds of adipocytes. Several lncRNAs that regulate potentially goat preadipocyte differentiation were identified, such as XR_001918 647.1, XR_001917728.1, XR_001297263.2 and LNC_004191. Furthermore, our findings from the present study may contribute to a better understanding of the molecular mechanisms underlying in goat meat quality and provide a theoretical basis for further goat molecular breeding.

ASCL1 activates neuronal stem cell-like lineage programming through remodeling of the chromatin landscape in prostate cancer

Treatment with androgen receptor pathway inhibitors (ARPIs) in prostate cancer leads to the emergence of resistant tumors characterized by lineage plasticity and differentiation toward neuroendocrine lineage. Here, we find that ARPIs induce a rapid epigenetic alteration mediated by large-scale chromatin remodeling to support activation of stem/neuronal transcriptional programs. We identify the proneuronal transcription factor ASCL1 motif to be enriched in hyper-accessible regions. ASCL1 acts as a driver of the lineage plastic, neuronal transcriptional program to support treatment resistance and neuroendocrine phenotype. Targeting ASCL1 switches the neuroendocrine lineage back to the luminal epithelial state. This effect is modulated by disruption of the polycomb repressive complex-2 through UHRF1/AMPK axis and change the chromatin architecture in favor of luminal phenotype. Our study provides insights into the epigenetic alterations induced by ARPIs, governed by ASCL1, provides a proof of principle of targeting ASCL1 to reverse neuroendocrine phenotype, support luminal conversion and re-addiction to ARPIs.

Following androgen receptor pathway inhibition prostate cancers can differentiate towards the neuroendocrine lineage. Here, the authors identify epigenetic alterations regulated by ASCL1 and suggest targeting ASCL1 to reverse the neuroendocrine phenotype.

Bioluminescent Zebrafish Transplantation Model for Drug Discovery

In the last decade, zebrafish have accompanied the mouse as a robust animal model for cancer research. The possibility of screening small-molecule inhibitors in a large number of zebrafish embryos makes this model particularly valuable. However, the dynamic visualization of fluorescently labeled tumor cells needs to be complemented by a more sensitive, easy, and rapid mode for evaluating tumor growth in vivo to enable high-throughput screening of clinically relevant drugs. In this study we proposed and validated a pre-clinical screening model for drug discovery by utilizing bioluminescence as our readout for the determination of transplanted cancer cell growth and inhibition in zebrafish embryos. For this purpose, we used NanoLuc luciferase, which ensured rapid cancer cell growth quantification in vivo with high sensitivity and low background when compared to conventional fluorescence measurements. This allowed us large-scale evaluation of in vivo drug responses of 180 kinase inhibitors in zebrafish. Our bioluminescent screening platform could facilitate identification of new small-molecules for targeted cancer therapy as well as for drug repurposing.

Establishment of DNA Methylation Profile Associated with TCM Syndrome in Endometriosis

Objectives

To screen the potential epigenetic biomarkers associated with endometriosis (EMS) and traditional Chinese medicine (TCM) syndrome EMS types.

Methods

A cohort of 99 participants comprising 42 EMS patients with cold coagulation blood stasis (CCBS) syndrome, 35 EMS patients with Qi stagnation blood stasis (QSBS) syndrome, and 22 women of childbearing age without EMS were recruited. Reduced representation bisulfite sequencing (RRBS) was used to establish the differential DNA methylation profiles in human peripheral blood samples obtained from four non-EMS and four EMS patients with CCBS or QSBS syndrome, respectively. Differentially expressed genes (DEGs) were verified in 18 non-EMS, 38 CCBS-EMS, and 31 QSBS-EMS using pyrosequencing.

Results

Methylation sites of 123942, 127229, and 115961 were found in peripheral blood DNA of non-EMS, CCBS-EMS, and QSBS-EMS patients, respectively. GO and KEGG analyses showed that the pathological process of EMS may be closely related to the nervous system development, cell junctions, GABA-gated chloride ion channel activity, nicotine addiction, Hippo signaling pathway, mRNA surveillance pathway, and Wnt signaling pathway. The methylation level at CpG site within HDAC6 gene in QSBS-EMS patients was significantly different from that in control women.

Conclusions

The changes in DNA methylation in peripheral blood samples may be associated with EMS and TCM syndrome EMS types. The methylation level of HDAC6 gene may be used to distinguish QSBS-EMS patients from women without EMS.

Lineage Conversion in Pediatric B-Cell Precursor Acute Leukemia under Blinatumomab Therapy

We report incidence and deep molecular characteristics of lineage switch in 182 pediatric patients affected by B-cell precursor acute lymphoblastic leukemia (BCP-ALL), who were treated with blinatumomab. We documented six cases of lineage switch that occurred after or during blinatumomab exposure. Therefore, lineage conversion was found in 17.4% of all resistance cases (4/27) and 3.2% of relapses (2/63). Half of patients switched completely from BCP-ALL to CD19-negative acute myeloid leukemia, others retained CD19-positive B-blasts and acquired an additional CD19-negative blast population: myeloid or unclassifiable. Five patients had KMT2A gene rearrangements; one had TCF3::ZNF384 translocation. The presented cases showed consistency of gene rearrangements and fusion transcripts across initially diagnosed leukemia and lineage switch. In two of six patients, the clonal architecture assessed by IG/TR gene rearrangements was stable, while in others, loss of clones or gain of new clones was noted. KMT2A-r patients demonstrated very few additional mutations, while in the TCF3::ZNF384 case, lineage switch was accompanied by a large set of additional mutations. The immunophenotype of an existing leukemia sometimes changes via different mechanisms and with different additional molecular changes. Careful investigation of all BM compartments together with all molecular -minimal residual disease studies can lead to reliable identification of lineage switch.

A Feedback Regulatory Loop Containing McdR and WhiB2 Controls Cell Division and DNA Repair in Mycobacteria

Cell division must be coordinated with DNA repair, which is strictly regulated in response to different drugs and environmental stresses in bacteria. However, the mechanisms by which mycobacteria orchestrate these two processes remain largely uncharacterized. Here, we report a regulatory loop between two essential mycobacterial regulators, McdR (Rv1830) and WhiB2, in coordinating the processes of cell division and DNA repair. McdR inhibits cell division-associated whiB2 expression by binding to the AATnACAnnnnTGTnATT motif in the promoter region. Furthermore, McdR overexpression simultaneously activates imuAB and dnaE2 expression to promote error-prone DNA repair, which facilitates genetic adaptation to stress conditions. Through a feedback mechanism, WhiB2 activates mcdR expression by binding to the cGACACGc motif in the promoter region. Importantly, analyses of mutations in clinical Mycobacterium tuberculosis strains indicate that disruption of this McdR-WhiB2 feedback regulatory loop influences expression of both cell growth- and DNA repair-associated genes, which further supports the contribution of McdR-WhiB2 regulatory loop in regulating mycobacterial cell growth and drug resistance. This highly conserved feedback regulatory loop provides fresh insight into the link between mycobacterial cell growth control and stress responses. IMPORTANCE Drug-resistant M. tuberculosis poses a threat to the control and prevention of tuberculosis (TB) worldwide. Thus, there is a need to identify the mechanisms enabling M. tuberculosis to adapt and grow under drug-induced stress. Rv1830 has been shown to be associated with drug resistance in M. tuberculosis, but its mechanisms have not yet been elucidated. Here, we reveal a regulatory role of Rv1830, which coordinates cell division and DNA repair in mycobacteria, and rename it McdR (mycobacterial cell division regulator). An increase in McdR levels represses the expression of cell division-associated whiB2 but activates the DNA repair-associated, error-prone enzymes ImuA/B and DnaE2, which in turn facilitates adaptation to stress responses and drug resistance. Furthermore, WhiB2 activates the transcription of mcdR to form a conserved regulatory loop. These data provide new insights into the mechanisms controlling mycobacterial cell growth and stress responses.

Single-Cell Analysis Uncovers Osteoblast Factor Growth Differentiation Factor 10 as Mediator of Vascular Smooth Muscle Cell Phenotypic Modulation Associated with Plaque Rupture in Human Carotid Artery Disease

(1) Background: Vascular smooth muscle cells (VSMCs) undergo a complex phenotypic switch in response to atherosclerosis environmental triggers, contributing to atherosclerosis disease progression. However, the complex heterogeneity of VSMCs and how VSMC dedifferentiation affects human carotid artery disease (CAD) risk has not been clearly established. (2) Method: A single-cell RNA sequencing analysis of CD45- cells derived from the atherosclerotic aorta of Apolipoprotein E-deficient (Apoe-/-) mice on a normal cholesterol diet (NCD) or a high cholesterol diet (HCD), respecting the site-specific predisposition to atherosclerosis was performed. Growth Differentiation Factor 10 (GDF10) role in VSMCs phenotypic switch was investigated via flow cytometry, immunofluorescence in human atherosclerotic plaques. (3) Results: scRNAseq analysis revealed the transcriptomic profile of seven clusters, five of which showed disease-relevant gene signature of VSMC macrophagic calcific phenotype, VSMC mesenchymal chondrogenic phenotype, VSMC inflammatory and fibro-phenotype and VSMC inflammatory phenotype. Osteoblast factor GDF10 involved in ossification and osteoblast differentiation emerged as a hallmark of VSMCs undergoing phenotypic switch. Under hypercholesteremia, GDF10 triggered VSMC osteogenic switch in vitro. The abundance of GDF10 expressing osteogenic-like VSMCs cells was linked to the occurrence of carotid artery disease (CAD) events. (4) Conclusions: Taken together, these results provide evidence about GDF10-mediated VSMC osteogenic switch, with a likely detrimental role in atherosclerotic plaque stability.

Human pluripotent stem-cell-derived islets ameliorate diabetes in non-human primates

Human pluripotent stem-cell-derived islets (hPSC-islets) are a promising cell resource for diabetes treatment^1,2. However, this therapeutic strategy has not been systematically assessed in large animal models physiologically similar to humans, such as non-human primates³. In this study, we generated islets from human chemically induced pluripotent stem cells (hCiPSC-islets) and show that a one-dose intraportal infusion of hCiPSC-islets into diabetic non-human primates effectively restored endogenous insulin secretion and improved glycemic control. Fasting and average pre-prandial blood glucose levels significantly decreased in all recipients, accompanied by meal or glucose-responsive C-peptide release and overall increase in body weight. Notably, in the four long-term follow-up macaques, average hemoglobin A1c dropped by over 2% compared with peak values, whereas the average exogenous insulin requirement reduced by 49% 15 weeks after transplantation. Collectively, our findings show the feasibility of hPSC-islets for diabetic treatment in a preclinical context, marking a substantial step forward in clinical translation of hPSC-islets.

Isolation of AmpC- and extended spectrum β-lactamase-producing Enterobacterales from fresh vegetables in the United States

Vegetables may serve as a reservoir for antibiotic resistant bacteria and resistance genes. AmpC β-lactamases and extended spectrum beta-lactamases (ESBL) inactivate commonly used β-lactam antibiotics, including penicillins and cephalosporins. In this study, we determined the prevalence of AmpC and ESBL-producing Enterobacterales in retail vegetables in the United States. A total of 88 vegetable samples were collected for the screening of AmpC and ESBL-producing Enterobacterales using CHROMagar ESBL agar. These vegetables included washed ready-to-eat salad (23), microgreens/sprouts (13), lettuce (11), herbs (11), spinach (5), mushrooms (5), brussels sprouts (4), kale (3), and other vegetable samples (13). AmpC and ESBL activity in these isolates were determined using double disk combination tests. Two vegetable samples (2.27%), organic basil and brussels sprouts, were positive for AmpC-producing Enterobacterales and eight samples (9.09%), including bean sprouts, organic parsley, organic baby spinach, and several mixed salads, were positive for ESBL-producing Enterobacterales. Whole genome sequencing was used to identify the bacterial species and resistance genes in these isolates. Genes encoding AmpC β-lactamases were found in Enterobacter hormaechei strains S43-1 and 74-2, which were consistent with AmpC production phenotypes. Multidrug-resistant E. hormaechei strains S11-1, S17-1, and S45-4 possess an ESBL gene, blaSHV66 , whereas five Serratia fonticola isolates contain genes encoding a minor ESBL, FONA-5. In addition, we used shotgun metagenomic sequencing approach to examine the microbiome and resistome profiles of three spinach samples. We found that Pseudomonas was the most prevalent bacteria genus in the spinach samples. Within the Enterobacteriaceae family, Enterobacter was the most abundant genus in the spinach samples. Moreover, antibiotic resistance genes encoding 12 major classes of antibiotics, including β-lactam antibiotics, aminoglycoside, macrolide, fluoroquinolone, and others, were found in these spinach samples. Therefore, vegetables can serve as an important vehicle for transmitting antibiotic resistance. The study highlights the need for antibiotic resistance surveillance in vegetable products.

Molecular diagnosis of Kallmann syndrome with diabetes by whole exome sequencing and bioinformatic approaches

BACKGROUND

Kallmann syndrome (KS) is a hypogonadotropic hypogonadism accompanied by anosmia or hyposmia. It is associated with the low secretion of gonadotropins which can lead to other abnormal endocrine metabolism disorders such as diabetes. Through genetic and molecular biological methods, more than 10 KS pathogenic genes have been found.

AIM

To identify the existing mutation sites of KS with diabetes and reveal the relationship between genotype and phenotype.

METHODS

We studied KS pathogenesis through high-throughput exome sequencing on four diabetes’ patients with KS for screening the potential pathogenic sites and exploring the genotype-phenotype correlation. Clinical data and peripheral blood samples were collected from the patients. White blood cells were separated and genomic DNA was extracted. High-throughput sequencing of all exons in the candidate pathogenic genes of probands was performed, and the results obtained were analyzed.

RESULTS

Sequencing revealed mutations in the KLB p.T313M, ANOS1 p.C172F, and IGSF10 gene (p.Lys1819Arg and p.Arg1035Thr) at different sites, which may have been associated with disease onset.

CONCLUSION

The diagnosis of KS is challenging, especially in early puberty, and the clinical manifestations reflect physical delays in development and puberty. Timely diagnosis and treatment can induce puberty, thereby improving sexual, bone, metabolic and mental health.

Analysis of the long noncoding RNA profiles of RD and SH-SY5Y cells infected with coxsackievirus B5, using RNA sequencing

Hand, foot, and mouth disease caused by coxsackievirus B5 (CV-B5) is a considerable threat to infant health, especially with regard to neurological damage. Long noncoding RNAs (lncRNAs) are known to play pivotal roles in virus-host interactions. However, the roles of lncRNAs in CV-B5-host interactions have not yet been elucidated. In the current study, we used RNA sequencing to determine the expression profiles of lncRNAs in CV-B5-infected human rhabdomyosarcoma (RD) and SH-SY5Y cells. Our results showed that, of the differentially expressed lncRNAs, 508 were upregulated and 760 were downregulated in RD cells. Of these, 46.2% were long noncoding intergenic RNAs (lincRNAs), 28.6% were antisense lncRNAs, 24.1% were sense overlapping lncRNAs, and 1.0% were sense intronic lncRNAs. Moreover, 792 lncRNAs were upregulated and 811 lncRNAs were downregulated in SH-SY5Y cells, 48.6% of which were lincRNAs, 34.7% were antisense lncRNAs, 16.0% were sense overlapping lncRNAs, and 0.8% were sense intronic lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that lncRNAs that were differentially expressed in CV-B5-infected RD cells were associated with disease, and those differentially expressed in SH-SY5Y cells were involved in signaling pathways. RT-qPCR analysis of seven lncRNAs supported these results. Moreover, our study revealed that lncRNA-IL12A inhibits viral replication. We conclude that lncRNAs constitute potential novel molecular targets for the prevention and treatment of CV-B5 infection and also may serve to distinguish neurogenic diseases caused by CV-B5 infection.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

Molecular characterization of porcine epidemic diarrhoea virus (PEDV) in Poland reveals the presence of swine enteric coronavirus (SeCoV) sequence in S gene

Porcine epidemic diarrhoea (PED) is a highly contagious enteric viral disease of pigs with a high morbidity and mortality rate, which ultimately results in huge economic losses in the pig production sector. The etiological agent of this disease is the porcine epidemic diarrhoea virus (PEDV) which is an enveloped, positive single-stranded RNA virus. The aim of this study was to perform molecular characterization of PEDV to identify the strains circulating in Poland. In this study, 662 faecal samples from 2015 to 2021 were tested with reverse transcription quantitative real-time PCR (RT-qPCR) and the results showed that 3.8% of the tested samples revealed a positive result for PEDV. A phylogenetic analysis of the complete genome and complete S gene sequences showed that Polish PEDV strains belonged to the G1b (S-INDEL) subgroup and were closely related to the European PEDV strains isolated from 2014 to 2019. Furthermore, RDP4 analysis revealed that the Polish PEDV strains harboured a recombinant fragment of ~400 nt in the 5' end of S gene with PEDV and swine enteric coronavirus (SeCoV) being the major and minor parents, respectively. Antigenic analysis showed that the aa sequences of neutralizing epitopes were conserved among the Polish PEDV strains. Only one strain, #0100/5P, had a unique substitution in the COE epitope. However, Polish PEDV strains showed several substitutions, especially in the COE antigen, as compared to the classical strain CV777. To the best of our knowledge, this is the first report concerning the molecular characterization of porcine epidemic diarrhoea virus strains, as well as the first phylogenetic analysis for PEDV in Poland.

Regulation of the Muscarinic M3 Receptor by Myocardin-Related Transcription Factors

Myocardin-related transcription factors (MRTFs: myocardin/MYOCD, MRTF-A/MRTFA, and MRTF-B/MRTFB) are co-factors of serum response factor (SRF) that activate the smooth muscle cell (SMC) gene program and that play roles in cardiovascular development and mechanobiology. Gain and loss of function experiments have defined the SMC gene program under control of MRTFs, yet full understanding of their impact is lacking. In the present study, we tested the hypothesis that the muscarinic M₃ receptor (CHRM3) is regulated by MRTFs together with SRF. Forced expression of MYOCD (8d) in human coronary artery (SMC) followed by RNA-sequencing showed increased levels of M₂, M₃, and M₅ receptors (CHRM2: 2-fold, CHRM3: 16-fold, and CHRM5: 2-fold). The effect of MYOCD on M₃ was confirmed by RT-qPCR using both coronary artery and urinary bladder SMCs, and correlation analyses using human transcriptomic datasets suggested that M₃ may also be regulated by MRTF-B. Head-to-head comparisons of MYOCD, MRTF-A and MRTF-B, argued that while all MRTFs are effective, MRTF-B is the most powerful transactivator of CHRM3, causing a 600-fold increase at 120h. Accordingly, MRTF-B conferred responsiveness to the muscarinic agonist carbachol in Ca²⁺ imaging experiments. M₃ was suppressed on treatment with the MRTF-SRF inhibitor CCG-1423 using SMCs transduced with either MRTF-A or MRTF-B and using intact mouse esophagus in culture (by 92±2%). Moreover, silencing of SRF with a short hairpin reduced CHRM3 (by >60%) in parallel with α-actin (ACTA2). Tamoxifen inducible knockout of Srf in smooth muscle reduced Srf (by 54±4%) and Chrm3 (by 41±6%) in the urinary bladder at 10days, but Srf was much less reduced or unchanged in aorta, ileum, colon, trachea, and esophagus. Longer induction (21d) further accentuated the reduction of Chrm3 in the bladder and ileum, but no change was seen in the aorta. Single cell RNA-sequencing revealed that Mrtfb dominates in ECs, while Myocd dominates in SMCs, raising the possibility that Chrm3 may be driven by Mrtfb-Srf in the endothelium and by Myocd-Srf in SMCs. These findings define a novel transcriptional control mechanism for muscarinic M₃ receptors in human cells, and in mice, that could be targeted for therapy.

The epigenetic role of HTR1A antagonist in facilitaing GnRH expression for pubertal initiation control

Serotonin (5-hydroxytryptamine [5-HT]), a metabolite of tryptophan, acts on the components of the hypothalamus-hypophysis-gonad axis and induces puberty delay in mammals via 5-HT receptor 1A (HTR1A). However, the roles of HTR1A in the hypothalamus in pubertal regulation of gene expression are not fully understood. In the current study, the upregulated gonadotropin-releasing hormone (GnRH) expression in GT1-7 GnRH neuronal cells induced by the HTR1A antagonist WAY-100635 maleate was observed in vitro. Furthermore, RNA sequencing (RNA-seq) showed decreased expression of chromobox 4 (CBX4), a member of the polycomb-repressive complex 1 (PRC1), and the loss of RING2 and YY1 interaction with CBX4, suggesting the degradation of the PRC1 in GT1-7 cells treated with maleate. Chromatin immunoprecipitation sequencing (ChIP-seq) showed that the genome-wide occupancy of CBX4 and histone H2A lysine-119 ubiquitination (H2AK119ub) was compromised, especially on the promoter of GnRH. Finally, we determined that inactivation of phosphatidylinositol 3-kinase (PI3K)/Akt and mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) contributed to CBX4 downregulation. Taken together, we concluded that HTR1A antagonists could enhance GnRH transcription via PRC1 degradation and H2AK119ub loss driven by reduced CBX4 expression through PI3K/Akt and MAPK/ERK pathway suppression in GT1-7 cells and provided a potential epigenetic mechanism of action of HTR1A on GnRH gene expression for mammalian puberty onset.

Transcriptome Analysis Reveals the Important Role of WRKY28 in Fusarium oxysporum Resistance

Root rot of Populus davidiana × P. alba var. pyramidalis Louche (Pdpap) is caused by Fusarium oxysporum. We used RNA sequencing to study the molecular mechanisms and response pattern of Pdpap infected by F. oxysporum CFCC86068. We cloned the PdpapWRKY28 transcription factor gene and transformed the recombinant vector pBI121-PdpapWRKY28 into Pdpap. The resistance function of PdpapWRKY28 was verified using physiological and biochemical methods. By means of RNA sequencing, we detected 1,403 differentially expressed genes (DEGs) that are common in the different treatments by F. oxysporum. Furthermore, we found that overexpression of the PdpapWRKY28 gene may significantly improve the resistance of Pdpap plants to F. oxysporum. Our research reveals a key role for PdpapWRKY28 in the resistance response of Pdpap to F. oxysporum. Additionally, our results provide a theoretical basis for in-depth research on resistance breeding to combat root rot.

Genome Analysis of Phytophthora nicotianae JM01 Provides Insights into Its Pathogenicity Mechanisms

Phytophthora nicotianae is a widely distributed plant pathogen that can cause serious disease and cause significant economic losses to various crops, including tomatoes, tobacco, onions, and strawberries. To understand its pathogenic mechanisms and explore strategies for controlling diseases caused by this pathogen, we sequenced and analyzed the whole genome of Ph. nicotianae JM01. The Ph. nicotianae JM01 genome was assembled using a combination of approaches including shotgun sequencing, single-molecule sequencing, and the Hi-C technique. The assembled Ph. nicotianae JM01 genome is about 95.32 Mb, with contig and scaffold N50 54.23 kb and 113.15 kb, respectively. The average GC content of the whole-genome is about 49.02%, encoding 23,275 genes. In addition, we identified 19.15% of interspersed elements and 0.95% of tandem elements in the whole genome. A genome-wide phylogenetic tree indicated that Phytophthora diverged from Pythium approximately 156.32 Ma. Meanwhile, we found that 252 and 285 gene families showed expansion and contraction in Phytophthora when compared to gene families in Pythium. To determine the pathogenic mechanisms Ph. nicotianae JM01, we analyzed a suite of proteins involved in plant-pathogen interactions. The results revealed that gene duplication contributed to the expansion of Cell Wall Degrading Enzymes (CWDEs) such as glycoside hydrolases, and effectors such as Arg-Xaa-Leu-Arg (RXLR) effectors. In addition, transient expression was performed on Nicotiana benthamiana by infiltrating with Agrobacterium tumefaciens cells containing a cysteine-rich (SCR) protein. The results indicated that SCR can cause symptoms of hypersensitive response. Moreover, we also conducted comparative genome analysis among four Ph. nicotianae genomes. The completion of the Ph. nicotianae JM01 genome can not only help us understand its genomic characteristics, but also help us discover genes involved in infection and then help us understand its pathogenic mechanisms.

Protein Arginine Methyltransferase 1 Is Essential for the Meiosis of Male Germ Cells

Protein arginine methyltransferase 1 (PRMT1) is a major enzyme responsible for the formation of methylarginine in mammalian cells; however, its function in vivo is not well understood due to its early embryonic lethality in null mice exhibiting spontaneous DNA damage, cell cycle delays, and defects in check point activation. Here, we generated germ cell-specific Prmt1 knock-out (KO) mice to evaluate the function of PRMT1 in spermatogenesis. Our findings demonstrate that PRMT1 is vital for male fertility in mice. Spermatogenesis in Prmt1 KO mice was arrested at the zygotene-like stage of the first meiotic division due to an elevated number of DNA double-strand breaks (DSBs). There was a loss of methylation in meiotic recombination 11 (MRE11), the key endonuclease in MRE11/RAD50/NBS 1 (MRN) complex, resulting in the accumulation of SPO11 protein in DSBs. The ATM-mediated negative feedback control over SPO11 was lost and, consequently, the repair pathway of DSBs was highly affected in PRMT1 deficient male germ cells. Our findings provide a novel insight into the role of PRMT1-mediated asymmetric demethylation in mouse spermatogenesis.

Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed

Toll-like receptors (TLRs) 7 and 8 are important in single-stranded viral RNA recognition, so genetic variation of these genes may play a role in SRLVs infection and disease progression. Present study aimed to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs provirus concentration as index of disease progression. A total of 14 SNPs were detected, 6 SNPs in the TLR7 gene locus and 8 SNPs in the TLR8 gene. Nine of the 14 identified polymorphisms, 4 in the TLR7 gene and 5 in TLR8 gene, were significantly associated with the SRLVs proviral concentration. These SNPs were located in 3'UTR, 5'UTR and intron sequences as well as in the coding sequences, but they led to silent changes. Homozygous genotypes of three TLR7 SNPs (synonymous variant 1:50703293, 3'UTR variant 1:50701297 and 5'UTR variant 1:50718645) were observed in goats with lower provirus copy number as well as in seronegative animals. The results obtained in this study suggest that SNPs of TLR7/TLR8 genes may induce differential innate immune response towards SRLVs affecting proviral concentration and thereby disease pathogenesis and progression. These findings support a role for genetic variations of TLR7 and TLR8 in SRLVs infection and warrants further studies on the effect of TLR7/TLR8 polymorphisms on SRLVs infection in different populations.

Deciphering the diet of a wandering spider (Phoneutria boliviensis; Araneae: Ctenidae) by DNA metabarcoding of gut contents

Arachnids are the most abundant land predators. Despite the importance of their functional roles as predators and the necessity to understand their diet for conservation, the trophic ecology of many arachnid species has not been sufficiently studied. In the case of the wandering spider, Phoneutria boliviensis F. O. Pickard-Cambridge, 1897, only field and laboratory observational studies on their diet exist. By using a DNA metabarcoding approach, we compared the prey found in the gut content of males and females from three distant Colombian populations of P. boliviensis. By DNA metabarcoding of the cytochrome c oxidase subunit I (COI), we detected and identified 234 prey items (individual captured by the spider) belonging to 96 operational taxonomic units (OTUs), as prey for this wandering predator. Our results broaden the known diet of P. boliviensis with at least 75 prey taxa not previously registered in fieldwork or laboratory experimental trials. These results suggest that P. boliviensis feeds predominantly on invertebrates (Diptera, Lepidoptera, Coleoptera, and Orthoptera) and opportunistically on small squamates. Intersex and interpopulation differences were also observed. Assuming that prey preference does not vary between populations, these differences are likely associated with a higher local prey availability. Finally, we suggest that DNA metabarcoding can be used for evaluating subtle differences in the diet of distinct populations of P. boliviensis, particularly when predation records in the field cannot be established or quantified using direct observation.

An Investigation of the Male Genitourinary Abscess Originated from Urinary Tract in a Tertiary Hospital, Shanghai, China, from 2004 to 2019

Background

Male genitourinary abscess is one of the serious complications of urinary tract infections (UTIs). There were few researches on the clinical and pathogenic characteristics of male genitourinary abscess.

Patients and Methods

A retrospective observational study was conducted between January 2004 and April 2019. Male patients with genitourinary abscess originated from urinary tract, including sites of scrotum, testis, epididymis, spermatic cord, and prostate, were enrolled. Clinical and microbial records were collected and analyzed, and antimicrobial susceptibility testings were performed according to CLSI standard. Whole-genome sequencing was applied to detect the β-lactamase genes and virulence genes, as well as to determine the multilocus-sequence typing (MLST) of the collected Klebsiella pneumoniae (K. pneumoniae) isolates.

Results

A total of 22 male patients were included. The main clinical symptoms were fever (86.4%), scrotal swelling (68.2%), local skin warmth (59.1%) and ulceration (45.5%). Urinary irritation was often presented in prostate involved abscess. Ultrasound features had a 94.7% positive rate. Surgical treatment, including abscess drainage, was helpful to the prognosis. No matter where the specimens obtained from, including blood, urine or pus, multidrug-resistant K. pneumoniae was the dominant (11 cases, 50.0%) microorganism in positive cultures. Nine of eleven K. pneumoniae isolates had been preserved and recovered. As for MLST typing, all the nine available isolates of K. pneumoniae belonged to the ST11 type and characterized with bla_KPC-2 carbapenemase gene. Virulence genes rmpA2, ybtS, kfuC, wzi, aerobactin genes (iucABCD and iutA) and type 3 fimbriae genes (mrkAD) were identified in all the K. pneumoniae isolates.

Conclusion

It seemed that more patients under 35 years old were vulnerable to genitourinary abscess. There was an increasing trend that multidrug-resistant K. pneumoniae isolates with multiple virulence genes were involved in male genitourinary abscess. Prompt and proper antibiotic use, combined with adequate drainage of the abscess, was important to prognosis.