Bone morphogenetic protein (BMP) signaling determines neuroblastoma cell fate and sensitivity to retinoic acid.

Retinoic acid (RA) is a standard-of-care neuroblastoma drug thought to be effective by inducing differentiation. Curiously, RA has little effect on primary human tumors during upfront treatment but can eliminate neuroblastoma cells from the bone marrow during post-chemo consolidation therapy-a discrepancy that has never been explained. To investigate this, we treated a large cohort of neuroblastoma cell lines with RA and observed that the most RA-sensitive cells predominantly undergo apoptosis or senescence, rather than differentiation. We conducted genome-wide CRISPR knockout screens under RA treatment, which identified BMP signaling as controlling the apoptosis/senescence vs differentiation cell fate decision and determining RA's overall potency. We then discovered that BMP signaling activity is markedly higher in neuroblastoma patient samples at bone marrow metastatic sites, providing a plausible explanation for RA's ability to clear neuroblastoma cells specifically from the bone marrow, seemingly mimicking interactions between BMP and RA during normal development.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Dietary Chinese herbal mixture supplementation improves production performance by regulating reproductive hormones, antioxidant capacity, immunity, and intestinal health of broiler breeders.

Chinese herbs have been used as feed additives and are commonly utilized in domestic intensive livestock farming. However, their impact on the production performance and intestinal health of broiler breeders has yet to be thoroughly explored. This study aimed to evaluate the effects of a Chinese herbal mixture (CHM) on the production performance of broiler breeders in terms of reproductive hormones, antioxidant capacity, immunity, and intestinal health of broiler breeders. A total of 336 thirty-wk-old hens were randomly allotted to 4 groups with 6 replicates of fourteen hens each, which fed a basal diet supplemented with 0 (CON), 500 (CHM500), 1,000 (CHM1000), and 1,500 (CHM1500) mg/kg CHM for 56 days, respectively. Our results showed that dietary supplementation with CHM1000 increased the laying rate and number of SYF and decreased the feed conversion ratio (P < 0.05). All CHM groups increased oviduct and ovarian indexes, serum E2 and T-AOC levels, and decreased serum TG and MDA levels compared with CON (P < 0.05). In comparison to the CON group, the CHM1000 and CHM1500 groups increased serum ALB, IgM, and IL-10 levels, whereas the CHM1000 group also increased serum TP and SOD levels, and the CHM1500 group increased serum P and decreased serum TNF-α (P < 0.05). The addition of CHM increased FSHR expressions in the ovary, Claudin-1 expressions in the jejunum, and SOD1 expressions in the liver and ovary, but decreased the mRNA expressions of INH in the ovary as well as IL-2 and IL-6 expressions in the jejunum (P < 0.05). Moreover, CHM500 and CHM1000 groups increased CAT, GPx, and HO-1 expression in the ovary, and SOD1 and GPx expression in the jejunum, while decreasing IL-17A expression in the jejunum (P < 0.05). In addition, CHM1000 and CHM1500 groups increased villus height, VCR, and the mRNA expressions of Nrf2, HO-1, Occludin, and MUC2 in the jejunum, and IL-10 expression in the ovary, while decreasing IL-2 and IL-17A expression in the ovary, in addition to increasing GPx, Nrf2, HO-1, NQO1, and IL-10 expression in the liver (P < 0.05). Supplementation with CHM1000 increased ESR-α, ESR-ß, GnRH, Nrf2, and NQO1 expression in the ovary, but decreased IFN-γ expression in the ovary as well as crypt depth in the jejunum (P < 0.05). Supplementing CHM1500 increased NQO1 and ZO-1 expression in the jejunum and decreased IL-2 in the liver (P < 0.05). The high-throughput sequencing results showed that dietary CHM1000 supplementation altered the composition of the intestinal microbiota, as evidenced by the regulation of the genera Lactobacillus, Faecalibacterium, and Phascolarctobacterium. PICRUSt analysis revealed that metabolic pathways of bacterial chemotaxis, butanoate metabolism, and synthesis and degradation of ketone bodies were enriched in the CHM1000 group. Spearman's correlation analysis indicated that the differentiated genera were significantly associated with the production performance, serum hormone, and gut barrier-related genes. Taken together, supplementation of CHM, especially at 1,000 mg/kg, could improve production performance by regulating reproductive hormones, antioxidant capacity, immunity, and intestinal health of broiler breeders, and maybe provide insights into its application as a potential feed additive to promote the performance of broiler breeders.

Depth estimation from a single-shot fringe pattern based on DD-Inceptionv2-UNet.

The quick and accurate retrieval of an object's depth from a single-shot fringe pattern in fringe projection profilometry has been a topic of ongoing research. In recent years, with the development of deep learning, a deep learning technique to FPP for single-shot 3D measurement is being used. To improve the accuracy of depth estimation from a single-shot fringe pattern, we propose the depthwise separable Dilation Inceptionv2-UNet (DD-Inceptionv2-UNet) by adjusting the depth and width of the network model simultaneously. And we evaluate the model on both simulated and experimental datasets. The experimental results show that the error between the depth map predicted by the proposed method and the label is smaller, and the depth curve map is closer to the ground truth. And on the simulated dataset, the MAE of the proposed method decreased by 35.22%, compared to UNet. On the experimental dataset, the MAE of the proposed method decreased by 34.62%, compared to UNet. The proposed method is relatively outstanding in both quantitative and qualitative evaluations, effectively improving the accuracy of 3D measurement results from a single-shot fringe pattern.

A CRISPR-drug perturbational map for identifying compounds to combine with commonly used chemotherapeutics.

Combination chemotherapy is crucial for successfully treating cancer. However, the enormous number of possible drug combinations means discovering safe and effective combinations remains a significant challenge. To improve this process, we conduct large-scale targeted CRISPR knockout screens in drug-treated cells, creating a genetic map of druggable genes that sensitize cells to commonly used chemotherapeutics. We prioritize neuroblastoma, the most common extracranial pediatric solid tumor, where ~50% of high-risk patients do not survive. Our screen examines all druggable gene knockouts in 18 cell lines (10 neuroblastoma, 8 others) treated with 8 widely used drugs, resulting in 94,320 unique combination-cell line perturbations, which is comparable to the largest existing drug combination screens. Using dense drug-drug rescreening, we find that the top CRISPR-nominated drug combinations are more synergistic than standard-of-care combinations, suggesting existing combinations could be improved. As proof of principle, we discover that inhibition of PRKDC, a component of the non-homologous end-joining pathway, sensitizes high-risk neuroblastoma cells to the standard-of-care drug doxorubicin in vitro and in vivo using patient-derived xenograft (PDX) models. Our findings provide a valuable resource and demonstrate the feasibility of using targeted CRISPR knockout to discover combinations with common chemotherapeutics, a methodology with application across all cancers.

Absence of terminal deoxynucleotidyl transferase expression in T-ALL/LBL accumulates chromosomal abnormalities to induce drug resistance.

T-acute lymphoblastic leukemia/lymphoma (T-ALL/LBL) is a malignant neoplasm of immature lymphoblasts. Terminal deoxynucleotidyl transferase (TDT) is a template-independent DNA polymerase that plays an essential role in generating diversity for immunoglobulin genes. T-ALL/LBL patients with TDT- have a worse prognosis. However, how TDT- promotes the disease progression of T-ALL/LBL remains unknown. Here we analyzed the prognosis of T-ALL/LBL patients in Shanghai Children's Medical Center (SCMC) and confirmed that TDT- patients had a higher rate of recurrence and remission failure and worse outcomes. Cellular experiments demonstrated that TDT was involved in DNA damage repair. TDT knockout delayed DNA repair, arrested the cell cycle and decreased apoptosis to induce the accumulation of chromosomal abnormalities and tolerance to abnormal karyotypes. Our study demonstrated that the poor outcomes in TDT- T-ALL/LBL might be due to the drug resistance (VP16 and MTX) induced by chromosomal abnormalities. Our findings revealed novel functions and mechanisms of TDT in T-ALL/LBL and supported that hematopoietic stem cell transplantation (HSCT) might be a better choice for these patients.

Dynamic genome evolution in a model fern.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Molecular properties of epimutation hotspots.

Mistakes in the maintenance of CG methylation are a source of heritable epimutations in plants. Multigenerational surveys indicate that the rate of these stochastic events varies substantially across the genome, with some regions harbouring localized 'epimutation hotspots'. Using Arabidopsis as a model, we show that epimutation hotspots are indexed by a specific set of chromatin states that map to subregions of gene body methylation genes. Although these regions comprise only ~12% of all CGs in the genome, they account for ~63% of all epimutation events per unit time. Molecular profiling revealed that these regions contain unique sequence features, harbour steady-state intermediate methylation levels and act as putative targets of antagonistic DNA methylation pathways. We further demonstrate that experimentally induced shifts in steady-state methylation in these hotspot regions are sufficient to significantly alter local epimutation intensities. Our work lays the foundation for dissecting the molecular mechanisms and evolutionary consequences of epimutation hotspots in plants.

Heterochromatin is a quantitative trait associated with spontaneous epiallele formation.

Epialleles are meiotically heritable variations in expression states that are independent from changes in DNA sequence. Although they are common in plant genomes, their molecular origins are unknown. Here we show, using mutant and experimental populations, that epialleles in Arabidopsis thaliana that result from ectopic hypermethylation are due to feedback regulation of pathways that primarily function to maintain DNA methylation at heterochromatin. Perturbations to maintenance of heterochromatin methylation leads to feedback regulation of DNA methylation in genes. Using single base resolution methylomes from epigenetic recombinant inbred lines (epiRIL), we show that epiallelic variation is abundant in euchromatin, yet, associates with QTL primarily in heterochromatin regions. Mapping three-dimensional chromatin contacts shows that genes that are hotspots for ectopic hypermethylation have increases in contact frequencies with regions possessing H3K9me2. Altogether, these data show that feedback regulation of pathways that have evolved to maintain heterochromatin silencing leads to the origins of spontaneous hypermethylated epialleles.

Natural variation in DNA methylation homeostasis and the emergence of epialleles.

In plants and mammals, DNA methylation plays a critical role in transcriptional silencing by delineating heterochromatin from transcriptionally active euchromatin. A homeostatic balance between heterochromatin and euchromatin is essential to genomic stability. This is evident in many diseases and mutants for heterochromatin maintenance, which are characterized by global losses of DNA methylation coupled with localized ectopic gains of DNA methylation that alter transcription. Furthermore, we have shown that genome-wide methylation patterns in Arabidopsis thaliana are highly stable over generations, with the exception of rare epialleles. However, the extent to which natural variation in the robustness of targeting DNA methylation to heterochromatin exists, and the phenotypic consequences of such variation, remain to be fully explored. Here we describe the finding that heterochromatin and genic DNA methylation are highly variable among 725 A. thaliana accessions. We found that genic DNA methylation is inversely correlated with that in heterochromatin, suggesting that certain methylation pathway(s) may be redirected to genes upon the loss of heterochromatin. This redistribution likely involves a feedback loop involving the DNA methyltransferase, CHROMOMETHYLASE 3 (CMT3), H3K9me2, and histone turnover, as highly expressed, long genes with a high density of CMT3-preferred CWG sites are more likely to be methylated. Importantly, although the presence of CG methylation in genes alone may not affect transcription, genes containing CG methylation are more likely to become methylated at non-CG sites and silenced. These findings are consistent with the hypothesis that natural variation in DNA methylation homeostasis may underlie the evolution of epialleles that alter phenotypes.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation.

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Evolutionary and Experimental Loss of Gene Body Methylation and Its Consequence to Gene Expression.

In flowering plants, gene body methylation (gbM) is associated with a subset of constitutively expressed genes. It has been proposed that gbM modulates gene expression. Here, we show that there are no consistent and direct differences to expression following the loss of gbM. By comparing expression of gbM genes in Arabidopsis thaliana accessions to orthologous genes in two Eutrema salsugineum genotypes, we identified both positive and negative expression differences associated with gbM loss. However, expression is largely unaffected by gbM loss in E. salsugineum Expression differences between species were within the variation of expression observed within A. thaliana accessions that displayed variation in gbM. Furthermore, experimentally induced loss of gbM did not consistently lead to differences in expression compared to wild type. To date, there is no convincing data to support a direct causal link between the presence/absence of gbM and the modulation of expression in flowering plants.

Pericentromeric hypomethylation elicits an interferon response in an animal model of ICF syndrome.

Pericentromeric satellite repeats are enriched in 5-methylcytosine (5mC). Loss of 5mC at these sequences is common in cancer and is a hallmark of Immunodeficiency, Centromere and Facial abnormalities (ICF) syndrome. While the general importance of 5mC is well-established, the specific functions of 5mC at pericentromeres are less clear. To address this deficiency, we generated a viable animal model of pericentromeric hypomethylation through mutation of the ICF-gene ZBTB24. Deletion of zebrafish zbtb24 caused a progressive loss of 5mC at pericentromeres and ICF-like phenotypes. Hypomethylation of these repeats triggered derepression of pericentromeric transcripts and activation of an interferon-based innate immune response. Injection of pericentromeric RNA is sufficient to elicit this response in wild-type embryos, and mutation of the MDA5-MAVS dsRNA-sensing machinery blocks the response in mutants. These findings identify activation of the innate immune system as an early consequence of pericentromeric hypomethylation, implicating derepression of pericentromeric transcripts as a trigger of autoimmunity. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Genome-wide association study of 8 carcass traits in Jinghai Yellow chickens using specific-locus amplified fragment sequencing technology.

Carcass traits are important to the commercial chicken industry, and understanding the genetics of these traits will be useful in the development of commercially viable varieties of chickens. We conducted a genome-wide association study based on 8 carcass trait phenotypes in a population of 400 43-week-old Jinghai Yellow chickens. Specific-locus amplified fragment sequencing technology was used to identify 90,961 single nucleotide polymorphisms (SNP) distributed among 29 chromosomes and the mitochondrial genome. SNP that were significantly associated with phenotypic traits were identified by a simple general linear model. Fifteen SNP attained genome-wide significance (P < 1.87E−6) and were associated with 5 of the 8 carcass traits; only one SNP was significantly associated with 2 traits (foot weight and wing weight). Twelve genes were associated with these 15 SNP. A region of chromosome 4 between 75.5 and 76.1 Mb was associated with carcass weight, foot weight, and wing weight. An 84-kb region on chromosome 3 (51.2 Mb) was associated with eviscerated weight and semi-eviscerated weight.

Genome-wide association study of antibody level response to NDV and IBV in Jinghai yellow chicken based on SLAF-seq technology.

Newcastle disease (ND) and avian infectious bronchitis (IB) are contagious diseases of chickens. To identify genes associated with antibody levels against ND and IB, a genome-wide association study was performed using specific-locus amplified fragment sequencing (SLAF-seq) technology in Jinghai yellow chickens. This determined six single-nucleotide polymorphisms (SNPs) that were associated with antibody levels against Newcastle disease virus (NDV): rsZ2494661, rsZ2494710, rs1211307701, rs1211307711, rs1218289310 and rs420701988. Of these, rsZ2494661 and rsZ2494710 reached the 5 % Bonferroni genome-wide significance level (5.5E-07) and they were both 134.7 kb downstream of the SETBP1 gene. The remaining four SNPs had 'suggestive' genome-wide significance levels (1.1E-05) and they were within or near the Plexin B1, LRRN1 and PDGFC genes. IB had two SNPs associated with antibody levels: rs149988433 and rs16170823; both reached chromosome-wide significance levels and they were near the USP7 and TRIM27 genes, respectively. Bioinformatics, GO annotation and pathway analysis indicated that five of these genes (Plexin B1, TRIM27, PDGFC, SETBP1 and USP7) may be important for the generation of protective antibodies against NDV and infectious bronchitis virus (IBV). This paves the way for further research on host immune responses against NDV.

A new computational strategy for predicting essential genes.

BACKGROUND: Determination of the minimum gene set for cellular life is one of the central goals in biology. Genome-wide essential gene identification has progressed rapidly in certain bacterial species; however, it remains difficult to achieve in most eukaryotic species. Several computational models have recently been developed to integrate gene features and used as alternatives to transfer gene essentiality annotations between organisms. RESULTS: We first collected features that were widely used by previous predictive models and assessed the relationships between gene features and gene essentiality using a stepwise regression model. We found two issues that could significantly reduce model accuracy: (i) the effect of multicollinearity among gene features and (ii) the diverse and even contrasting correlations between gene features and gene essentiality existing within and among different species. To address these issues, we developed a novel model called feature-based weighted Naïve Bayes model (FWM), which is based on Naïve Bayes classifiers, logistic regression, and genetic algorithm. The proposed model assesses features and filters out the effects of multicollinearity and diversity. The performance of FWM was compared with other popular models, such as support vector machine, Naïve Bayes model, and logistic regression model, by applying FWM to reciprocally predict essential genes among and within 21 species. Our results showed that FWM significantly improves the accuracy and robustness of essential gene prediction. CONCLUSIONS: FWM can remarkably improve the accuracy of essential gene prediction and may be used as an alternative method for other classification work. This method can contribute substantially to the knowledge of the minimum gene sets required for living organisms and the discovery of new drug targets.

Genome-wide analysis of selective constraints on high stability regions of mRNA reveals multiple compensatory mutations in Escherichia coli.

Message RNA (mRNA) carries a large number of local secondary structures, with structural stability to participate in the regulations of gene expression. A worthy question is how the local structural stability is maintained under the constraint that multiple selective pressures are imposed on mRNA local regions. Here, we performed the first genome-wide study of natural selection operating on high structural stability regions (HSRs) of mRNAs in Escherichia coli. We found that HSR tends to adjust the folded conformation to reduce the harm of mutations, showing a high level of mutational robustness. Moreover, guanine preference in HSR was observed, supporting the hypothesis that the selective constraint for high structural stability may partly account for the high percentage of G content in Escherichia coli genome. Notably, we found a substantially reduced synonymous substitution rate in HSRs compared with that in their adjacent regions. Surprisingly and interestingly, the non-key sites in HSRs, which have slight effect on structural stability, have synonymous substitution rate equivalent to background regions. To explain this result, we identified compensatory mutations in HSRs based on structural stability, and found that a considerable number of synonymous mutations occur to restore the structural stability decreased heavily by the mutations on key sites. Overall, these results suggest a significant role of local structural stability as a selective force operating on mRNA, which furthers our understanding of the constraints imposed on protein-coding RNAs.