Targeted haplotyping in pharmacogenomics using Oxford Nanopore Technologies' adaptive sampling

Pharmacogenomics (PGx) studies the impact of interindividual genomic variation on drug response, allowing the opportunity to tailor the dosing regimen for each patient. Current targeted PGx testing platforms are mainly based on microarray, polymerase chain reaction, or short-read sequencing. Despite demonstrating great value for the identification of single nucleotide variants (SNVs) and insertion/deletions (INDELs), these assays do not permit identification of large structural variants, nor do they allow unambiguous haplotype phasing for star-allele assignment. Here, we used Oxford Nanopore Technologies' adaptive sampling to enrich a panel of 1,036 genes with well-documented PGx relevance extracted from the Pharmacogenomics Knowledge Base (PharmGKB). By evaluating concordance with existing truth sets, we demonstrate accurate variant and star-allele calling for five Genome in a Bottle reference samples. We show that up to three samples can be multiplexed on one PromethION flow cell without a significant drop in variant calling performance, resulting in 99.35% and 99.84% recall and precision for the targeted variants, respectively. This work advances the use of nanopore sequencing in clinical PGx settings.

Post-feeding transcriptomics reveals essential genes expressed in the midgut of the desert locust

The digestive tract constitutes an important interface between an animal's internal and external environment. In insects, available gut transcriptome studies are mostly exploratory or look at changes upon infection or upon exposure to xenobiotics, mainly performed in species belonging to holometabolan orders, such as Diptera, Lepidoptera or Coleoptera. By contrast, studies focusing on gene expression changes after food uptake and during digestion are underrepresented. We have therefore compared the gene expression profiles in the midgut of the desert locust, Schistocerca gregaria, between three different time points after feeding, i.e., 24 h (no active digestion), 10 min (the initial stage of feeding), and 2 h (active food digestion). The observed gene expression profiles were consistent with the polyphagous herbivorous lifestyle of this hemimetabolan (orthopteran) species. Our study reveals the upregulation of 576 genes 2 h post-feeding. These are mostly predicted to be associated with digestive physiology, such as genes encoding putative digestive enzymes or nutrient transporters, as well as genes putatively involved in immunity or in xenobiotic metabolism. The 10 min time point represented an intermediate condition, suggesting that the S. gregaria midgut can react rapidly at the transcriptional level to the presence of food. Additionally, our study demonstrated the critical importance of two transcripts that exhibited a significant upregulation 2 h post-feeding: the vacuolar-type H(+)-ATPase and the sterol transporter Niemann-Pick 1b protein, which upon RNAi-induced knockdown resulted in a marked increase in mortality. Their vital role and accessibility via the midgut lumen may make the encoded proteins promising insecticidal target candidates, considering that the desert locust is infamous for its huge migrating swarms that can devastate the agricultural production in large areas of Northern Africa, the Middle East, and South Asia. In conclusion, the transcriptome datasets presented here will provide a useful and promising resource for studying the midgut physiology of S. gregaria, a socio-economically important pest species.

Haplotyping pharmacogenes using TLA combined with Illumina or Nanopore sequencing

The currently used pharmacogenetic genotyping assays offer limited haplotype information, which can potentially cause specific functional effects to be missed. This study tested if Targeted Locus Amplification (TLA), when using non-patient-specific primers combined with Illumina or Nanopore sequencing, can offer an advantage in terms of accurate phasing. The TLA method selectively amplifies and sequences entire genes based on crosslinking DNA in close physical proximity. This way, DNA fragments that were initially further apart in the genome are ligated into one molecule, making it possible to sequence distant variants within one short read. In this study, four pharmacogenes, CYP2D6, CYP2C19, CYP1A2 and BRCA1, were sequenced after enrichment using different primer pairs. Only 24% or 38% of the nucleotides mapped on target when using Illumina or Nanopore sequencing, respectively. With an average depth of more than 1000X for the regions of interest, none of the genes were entirely covered with either sequencing method. For three of the four genes, less than half of the variants were phased correctly compared to the reference. The Nanopore dataset with the optimized primer pair for CYP2D6 resulted in the correct haplotype, showing that this method can be used for reliable genotyping and phasing of pharmacogenes but does require patient-specific primer design and optimization to be effective.

Cas9 targeted nanopore sequencing with enhanced variant calling improves CYP2D6-CYP2D7 hybrid allele genotyping

CYP2D6 is a very important pharmacogene as it is responsible for the metabolization or bioactivation of 20 to 30% of the clinically used drugs. However, despite its relatively small length of only 4.4 kb, it is one of the most challenging pharmacogenes to genotype due to the high similarity with its neighboring pseudogenes and the frequent occurrence of CYP2D6-CYP2D7 hybrids. Unfortunately, most current genotyping methods are therefore not able to correctly determine the complete CYP2D6-CYP2D7 sequence. Therefore, we developed a genotyping assay to generate complete allele-specific consensus sequences of complex regions by optimizing the PCR-free nanopore Cas9-targeted sequencing (nCATS) method combined with adaptive sequencing, and developing a new comprehensive long read genotyping (CoLoRGen) pipeline. The CoLoRGen pipeline first generates consensus sequences of both alleles and subsequently determines both large structural and small variants to ultimately assign the correct star-alleles. In reference samples, our genotyping assay confirms the presence of CYP2D6-CYP2D7 large structural variants, single nucleotide variants (SNVs), and small insertions and deletions (INDELs) that go undetected by most current assays. Moreover, our results provide direct evidence that the CYP2D6 genotype of the NA12878 DNA should be updated to include the CYP2D6-CYP2D7 *68 hybrid and several additional single nucleotide variants compared to existing references. Ultimately, the nCATS-CoLoRGen genotyping assay additionally allows for more accurate gene function predictions by enabling the possibility to detect and phase de novo mutations in addition to known large structural and small variants.

During the last decades, the usefulness of personalized medicine has become increasingly apparent. Directly linked to that is the need for accurate genotyping assays to determine the pharmacogenetic profile of patients. Continuing research has led to the development of genotyping assays that perform quite robustly. However, complex genes remain an issue when it comes to determining the complete sequence correctly. An example of such a complex but very important pharmacogene is CYP2D6. Therefore, we developed a genotyping assay in an attempt to generate complete allele-specific consensus sequences of CYP2D6, by optimizing a targeted amplification-free long-read sequencing method and developing a new analysis pipeline. In reference samples, we showed that our genotyping assay performed accurately and confirmed the presence of variants that go undetected by most current assays. However, the implementation of this assay in practice is still hampered as the selected enrichment strategies inherently lead to a low percentage of on-target reads, resulting in low on-target sequencing depths. Further optimization and validation of the assay is thus needed, but definitely worth considering for follow-up research as we already demonstrated the added value for generating more complete genotypes, which on its turn will result in more accurate gene function predictions.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

Comparative study of preimplantation development following distinct assisted oocyte activation protocols in a PLC-zeta knockout mouse model

Mammalian fertilization encompasses a series of Ca2+ oscillations initiated by the sperm factor phospholipase C zeta (PLCζ). Some studies have shown that altering the Ca2+ oscillatory regime at fertilization affects preimplantation blastocyst development. However, assisted oocyte activation (AOA) protocols can induce oocyte activation in a manner that diverges profoundly from the physiological Ca2+ profiling. In our study, we used the newly developed PLCζ-null sperm to investigate the independent effect of AOA on mouse preimplantation embryogenesis. Based on previous findings, we hypothesized that AOA protocols with Ca2+ oscillatory responses might improve blastocyst formation rates and differing Ca2+ profiles might alter blastocyst transcriptomes. A total of 326 MII B6D2F1-oocytes were used to describe Ca2+ profiles and to compare embryonic development and individual blastocyst transcriptomes between four control conditions: C1 (in-vivo fertilization), C2 (ICSI control sperm), C3 (parthenogenesis) and C4 (ICSI-PLCζ-KO sperm) and four AOA groups: AOA1 (human recombinant PLCζ), AOA2 (Sr2+), AOA3 (ionomycin) and AOA4 (TPEN). All groups revealed remarkable variations in their Ca2+ profiles; however, oocyte activation rates were comparable between the controls (91.1% ± 13.8%) and AOA (86.9% ± 11.1%) groups. AOA methods which enable Ca2+ oscillatory responses (AOA1: 41% and AOA2: 75%) or single Ca2+ transients (AOA3: 50%) showed no significantly different blastocyst rates compared to ICSI control group (C2: 70%). In contrast, we observed a significant decrease in compaction (53% vs. 83%) and blastocyst rates (41% vs. 70%) in the absence of an initial Ca2+ trigger (AOA4) compared with the C2 group. Transcription profiles did not identify significant differences in gene expression levels between the ICSI control group (C2) and the four AOA groups.

MYCN-induced nucleolar stress drives an early senescence-like transcriptional program in hTERT-immortalized RPE cells

MYCN is an oncogenic driver in neural crest-derived neuroblastoma and medulloblastoma. To better understand the early effects of MYCN activation in a neural-crest lineage context, we profiled the transcriptome of immortalized human retina pigment epithelial cells with inducible MYCN activation. Gene signatures associated with elevated MYC/MYCN activity were induced after 24 h of MYCN activation, which attenuated but sustained at later time points. Unexpectedly, MYCN activation was accompanied by reduced cell growth. Gene set enrichment analysis revealed a senescence-like signature with strong induction of p53 and p21 but in the absence of canonical hallmarks of senescence such as β-galactosidase positivity, suggesting incomplete cell fate commitment. When scrutinizing the putative drivers of this growth attenuation, differential gene expression analysis identified several regulators of nucleolar stress. This process was also reflected by phenotypic correlates such as cytoplasmic granule accrual and nucleolar coalescence. Hence, we propose that the induction of MYCN congests the translational machinery, causing nucleolar stress and driving cells into a transient pre-senescent state. Our findings shed new light on the early events induced by MYCN activation and may help unravelling which factors are required for cells to tolerate unscheduled MYCN overexpression during early malignant transformation.

Endogenous suppression of WNT signalling in human embryonic stem cells leads to low differentiation propensity towards definitive endoderm

Low differentiation propensity towards a targeted lineage can significantly hamper the utility of individual human pluripotent stem cell (hPSC) lines in biomedical applications. Here, we use monolayer and micropatterned cell cultures, as well as transcriptomic profiling, to investigate how variability in signalling pathway activity between human embryonic stem cell lines affects their differentiation efficiency towards definitive endoderm (DE). We show that endogenous suppression of WNT signalling in hPSCs at the onset of differentiation prevents the switch from self-renewal to DE specification. Gene expression profiling reveals that this inefficient switch is reflected in NANOG expression dynamics. Importantly, we demonstrate that higher WNT stimulation or inhibition of the PI3K/AKT signalling can overcome the DE commitment blockage. Our findings highlight that redirection of the activity of Activin/NODAL pathway by WNT signalling towards mediating DE fate specification is a vulnerable spot, as disruption of this process can result in poor hPSC specification towards DE.

Metallothioneins alter macrophage phenotype and represent novel therapeutic targets for acetaminophen-induced liver injury

Acetaminophen (APAP) intoxication is the foremost cause of drug-induced liver failure in developed countries. The only pharmacologic treatment option, N-acetylcysteine (NAC), is not effective for patients who are admitted too late and/or who have excessive liver damage, emphasizing the need for alternative treatment options. APAP intoxication results in hepatocyte death and release of danger signals, which further contribute to liver injury, in part by hepatic monocyte/macrophage infiltration and activation. Metallothionein (MT) 1 and 2 have important danger signaling functions and might represent novel therapeutic targets in APAP overdose. Therefore, we evaluated hepatic MT expression and the effect of anti-MT antibodies on the transcriptional profile of the hepatic macrophage population and liver injury following APAP overdose in mice. Hepatic MT expression was significantly induced in APAP-intoxicated mice and abundantly present in human livers. APAP intoxication in mice resulted in increased serum transaminase levels, extended necrotic regions on liver histology and induced expression of proinflammatory markers, which was significantly less pronounced in mice treated with anti-MT antibodies. Anti-MT antibody therapy attenuated proinflammatory macrophage polarization, as demonstrated by RNA sequencing analyses of isolated liver macrophages and in LPS-stimulated bone marrow-derived macrophages. Importantly, NAC and anti-MT antibodies were equally effective whereas administration of anti-MT antibody in combination with NAC exceeded the efficiency of both monotherapies in APAP-induced liver injury (AILI). We conclude that the neutralization of secreted MTs using a monoclonal antibody is a novel therapeutic strategy as mono- or add-on therapy for AILI. In addition, we provide evidence suggesting that MTs in the extracellular environment are involved in macrophage polarization.

Comparative genomics of Flavobacterium columnare unveils novel insights in virulence and antimicrobial resistance mechanisms

This study reports the comparative analyses of four Flavobacterium columnare isolates that have different virulence and antimicrobial resistance patterns. The main research goal was to reveal new insights into possible virulence genes by comparing the genomes of bacterial isolates that could induce tissue damage and mortality versus the genome of a non-virulent isolate. The results indicated that only the genomes of the virulent isolates possessed unique genes encoding amongst others a methyl-accepting chemotaxis protein possibly involved in the initial colonization of tissue, and several VgrG proteins engaged in interbacterial competition. Furthermore, comparisons of genes unique for the genomes of the highly virulent (HV) carp and trout isolates versus the, respectively, low and non-virulent carp and trout isolates were performed. An important part of the identified unique virulence genes of the HV-trout isolate was located in one particular gene region identified as a genomic island. This region contained araC and nodT genes, both linked to pathogenic and multidrug-resistance, and a luxR-gene, functional in bacterial cell-to-cell communication. Furthermore, the genome of the HV-trout isolate possessed unique sugar-transferases possibly important in bacterial adhesion. The second research goal was to obtain insights into the genetic basis of acquired antimicrobial resistance. Several point-mutations were discovered in gyrase-genes of an isolate showing phenotypic resistance towards first and second-generation quinolones, which were absent in isolates susceptible to quinolones. Tetracycline-resistance gene tetA was found in an isolate displaying acquired phenotypic resistance towards oxytetracycline. Although not localized on a prophage, several flanking genes were indicative of the gene’s mobile character.

Activin A-derived human embryonic stem cells show increased competence to differentiate into primordial germ cell-like cells

Protocols for specifying human primordial germ cell‐like cells (hPGCLCs) from human embryonic stem cells (hESCs) remain hindered by differences between hESC lines, their derivation methods, and maintenance culture conditions. This poses significant challenges for establishing reproducible in vitro models of human gametogenesis. Here, we investigated the influence of activin A (ActA) during derivation and maintenance on the propensity of hESCs to differentiate into PGCLCs. We show that continuous ActA supplementation during hESC derivation (from blastocyst until the formation of the post‐inner cell mass intermediate [PICMI]) and supplementation (from the first passage of the PICMI onwards) is beneficial to differentiate hESCs to PGCLCs subsequently. Moreover, comparing isogenic primed and naïve states prior to differentiation, we showed that conversion of hESCs to the 4i‐state improves differentiation to (TNAP [tissue nonspecific alkaline phosphatase]+/PDPN [podoplanin]+) PGCLCs. Those PGCLCs expressed several germ cell markers, including TFAP2C (transcription factor AP‐2 gamma), SOX17 (SRY‐box transcription factor 17), and NANOS3 (nanos C2HC‐type zinc finger 3), and markers associated with germ cell migration, CXCR4 (C‐X‐C motif chemokine receptor 4), LAMA4 (laminin subunit alpha 4), ITGA6 (integrin subunit alpha 6), and CDH4 (cadherin 4), suggesting that the large numbers of PGCLCs obtained may be suitable to differentiate further into more mature germ cells. Finally, hESCs derived in the presence of ActA showed higher competence to differentiate to hPGCLC, in particular if transiently converted to the 4i‐state. Our work provides insights into the differences in differentiation propensity of hESCs and delivers an optimized protocol to support efficient human germ cell derivation.

HES1 and HES4 have non-redundant roles downstream of Notch during early human T-cell development

In both mouse and human, Notch1 activation is the main initial driver to induce T-cell development in hematopoietic progenitor cells. The initiation of this developmental process coincides with Notch1-dependent repression of differentiation towards other hematopoietic lineages. Although well described in mice, the role of the individual Notch1 target genes during these hematopoietic developmental choices is still unclear in human, particularly for HES4 since no orthologous gene is present in the mouse. Here, we investigated the functional capacity of the Notch1 target genes HES1 and HES4 to modulate human Notch1-dependent hematopoietic lineage decisions and their requirement during early T-cell development. We show that both genes are upregulated in a Notch-dependent manner during early T-cell development and that HES1 acts as a repressor of differentiation by maintaining a quiescent stem cell signature in CD34+ hematopoietic progenitor cells. While HES4 can also inhibit natural killer and myeloid cell development like HES1, it acts differently on the T- versus B-cell lineage choice. Surprisingly, HES4 is not capable of repressing B-cell development, the most sensitive hematopoietic lineage with respect to Notch-mediated repression. In contrast to HES1, HES4 promotes initiation of early T-cell development, but ectopic expression of HES4, or HES1 and HES4 combined, is not sufficient to induce T-lineage differentiation. Importantly, knockdown of HES1 or HES4 significantly reduces human T-cell development. Overall, we show that the Notch1 target genes HES1 and HES4 have non-redundant roles during early human T-cell development which may relate to differences in mediating Notch-dependent human hematopoietic lineage decisions.

Pan-cancer pharmacogenetics: targeted sequencing panels or exome sequencing?

Aim: This study provides clinicians and researchers with an informed choice between current commercially available targeted sequencing panels and exome sequencing panels in the context of pan-cancer pharmacogenetics. Materials & methods: Nine contemporary commercially available targeted pan-cancer panels and the xGen Exome Research Panel v2 were investigated to determine to what extent they cover the pharmacogenetic variant-drug interactions in five available cancer knowledgebases, and the driver mutations and fusion genes in the Cancer Genome Atlas. Results: xGen Exome Research Panel v2 and TrueSight Oncology 500 target 71.0 and 68.9% of the pharmacogenetic interactions in the available knowledgebases; and 93.7 and 86.0% of the driver mutations in the Cancer Genome Atlas, respectively. All other studied panels target lower percentages. Conclusion: Exome sequencing outperforms pan-cancer targeted sequencing panels in terms of covered cancer pharmacogenetic variant-drug interactions and pharmacogenetic cancer variants.

Kinship analysis on single cells after whole genome amplification

Short Tandem Repeat (STR-) and Single Nucleotide Polymorphism (SNP-) genotyping have been extensively studied within forensic kinship analysis. Nevertheless, no results have been reported on kinship analysis after whole genome amplification (WGA) of single cells. This WGA step is a necessary procedure in several applications, such as cell-based non-invasive prenatal testing (cbNIPT) and pre-implantation genetic diagnosis (PGD). In cbNIPT, all putative fetal cells must be discriminated from maternal cells after enrichment from whole blood. This study investigates the efficacy and evidential value of STR- and SNP-genotyping methods for the discrimination of 24 single cells after WGA, within three families. Formaldehyde-fixed and unfixed cells are assessed in offspring-parent duos and offspring-mother-father trios. Results demonstrate that both genotyping methods can be used in all tested conditions and scenarios with 100% sensitivity and 100% specificity, with a similar evidential value for fixed and unfixed cells. Moreover, sequence-based SNP-genotyping results in a higher evidential value than length-based STR-genotyping after WGA, which is not observed using high-quality offspring bulk DNA samples. Finally, it is also demonstrated that the availability of the DNA genotypes of both parents strongly increases the evidential value of the results.

Myeloid-specific IRE1alpha deletion reduces tumour development in a diabetic, non-alcoholic steatohepatitis-induced hepatocellular carcinoma mouse model

BACKGROUND AND AIMS

Obesity, diabetes and associated non-alcoholic steatohepatitis (NASH) are rising risk factors for hepatocellular carcinoma (HCC). Macrophages are important immune cells involved in inflammation and tumour development. Macrophage inositol-requiring enzyme 1 alpha (IRE1α), an ER-stress protein, has been shown to be involved in macrophage cytokine production, and myeloid-specific IRE1α knock-out (myeloid IRE1α-KO) mice showed reduced weight gain during high-fat diet feeding. However, the effect of myeloid IRE1α on NASH and subsequent HCC development has not been examined. Here, we characterized the transcriptional profile of the hepatic macrophage population in a diabetes-NASH-HCC mouse model, and investigated the effect of myeloid-specific IRE1α deletion on the phenotype of hepatic macrophage subsets and experimental NASH-HCC development.

METHODS

Mice with non-functional myeloid IRE1α were created by crossing Ire1a floxed mice with Lysm-Cre mice. Two-day old myeloid IRE1α-KO and wild type (WT) mice were subcutaneously injected with streptozotocin (STZ), and male mice were fed a high-fat, -sucrose, -cholesterol diet (Western diet, WD) from the age of 4 weeks until 21 weeks. Control myeloid IRE1α-KO and WT mice received a PBS injection and were fed a matched control diet. These mice were evaluated for obesity, diabetes, NASH and HCC. The hepatic macrophage population was evaluated by flow cytometry and RNA sequencing on FACS-isolated macrophage subsets.

RESULTS

STZ-injection and WD feeding resulted in an impaired glucose tolerance, advanced NASH with fibrosis, and HCC development. Myeloid IRE1α-KO STZ mice showed lower fasting glucose levels at the start of WD feeding, and an improved glucose tolerance and attenuated HCC development after 17 weeks of WD feeding despite a similar degree of liver steatosis and inflammation compared to WT mice. Transcriptomic analysis of WT liver Kupffer cells, macrophages and monocytes revealed phenotypical changes in those cell subsets during NASH-HCC development. Isolated liver Kupffer cells and macrophages from mice with a myeloid IRE1α deletion showed downregulated pathways involved in immune system activation and metabolic pathways (only in Kupffer cells), whereas pathways involved in cell division and metabolism were upregulated in monocytes. These transcriptional differences were attenuated during NASH-HCC development.

CONCLUSION

Our results show that myeloid-specific IRE1α deletion results in an altered transcriptional profile of hepatic macrophages and dampens diabetes-induced NASH-HCC development, possibly by attenuated diabetes induction.

Maternal Recognition of Pregnancy in the Horse: Are MicroRNAs the Secret Messengers?

The signal for maternal recognition of pregnancy (MRP) has still not been identified in the horse. High-throughput molecular biology at the embryo-maternal interface has substantially contributed to the knowledge on pathways affected during MRP, but an integrated study in which proteomics, transcriptomics and miRNA expression can be linked directly is currently lacking. The aim of this study was to provide such analysis. Endometrial biopsies, uterine fluid, embryonic tissues, and yolk sac fluid were collected 13 days after ovulation during pregnant and control cycles from the same mares. Micro-RNA-Sequencing was performed on all collected samples, mRNA-Sequencing on the same tissue samples and mass spectrometry was conducted previously on the same fluid samples. Differential expression of miRNA, mRNA and proteins showed high conformity with literature and confirmed involvement in pregnancy establishment, embryo quality, steroid synthesis and prostaglandin regulation, but the link between differential miRNAs and their targets was limited and did not indicate the identity of an unequivocal signal for MRP in the horse. Differential expression at the embryo-maternal interface was prominent, highlighting a potential role of miRNAs in embryo-maternal communication during early pregnancy in the horse. These data provide a strong basis for future targeted studies.

Dual RNA Sequencing of Vitis vinifera during Lasiodiplodia theobromae Infection Unveils Host-Pathogen Interactions

Lasiodiplodia theobromae is one of the most aggressive agents of the grapevine trunk disease Botryosphaeria dieback. Through a dual RNA-sequencing approach, this study aimed to give a broader perspective on the infection strategy deployed by L. theobromae, while understanding grapevine response. Approximately 0.05% and 90% of the reads were mapped to the genomes of L. theobromae and Vitis vinifera, respectively. Over 2500 genes were significantly differentially expressed in infected plants after 10 dpi, many of which are involved in the inducible defense mechanisms of grapevines. Gene expression analysis showed changes in the fungal metabolism of phenolic compounds, carbohydrate metabolism, transmembrane transport, and toxin synthesis. These functions are related to the pathogenicity mechanisms involved in plant cell wall degradation and fungal defense against antimicrobial substances produced by the host. Genes encoding for the degradation of plant phenylpropanoid precursors were up-regulated, suggesting that the fungus could evade the host defense response using the phenylpropanoid pathway. The up-regulation of many distinct components of the phenylpropanoid pathway in plants supports this hypothesis. Moreover, genes related to phytoalexin biosynthesis, hormone metabolism, cell wall modification enzymes, and pathogenesis-related proteins seem to be involved in the host responses observed. This study provides additional insights into the molecular mechanisms of L. theobromae and V. vinifera interactions.

YtrASa, a GntR-Family Transcription Factor, Represses Two Genetic Loci Encoding Membrane Proteins in Sulfolobus acidocaldarius

In bacteria, the GntR family is a widespread family of transcription factors responsible for the regulation of a myriad of biological processes. In contrast, despite their occurrence in archaea only a little information is available on the function of GntR-like transcription factors in this domain of life. The thermoacidophilic crenarchaeon Sulfolobus acidocaldarius harbors a GntR-like regulator belonging to the YtrA subfamily, encoded as the first gene in an operon with a second gene encoding a putative membrane protein. Here, we present a detailed characterization of this regulator, named YtrA_Sa, with a focus on regulon determination and mechanistic analysis with regards to DNA binding. Genome-wide chromatin immunoprecipitation and transcriptome experiments, the latter employing a ytrA_Sa overexpression strain, demonstrate that the regulator acts as a repressor on a very restricted regulon, consisting of only two targets including the operon encoding its own gene and a distinct genetic locus encoding another putative membrane protein. For both targets, a conserved 14-bp semi-palindromic binding motif was delineated that covers the transcriptional start site and that is surrounded by additional half-site motifs. The crystallographic structure of YtrA_Sa was determined, revealing a compact dimeric structure in which the DNA-binding motifs are oriented ideally to enable a specific high-affinity interaction with the core binding motif. This study provides new insights into the functioning of a YtrA-like regulator in the archaeal domain of life.

Transcriptional landscape changes during human embryonic stem cell derivation

STUDY QUESTION

What are the transcriptional changes occurring during the human embryonic stem cell (hESC) derivation process, from the inner cell mass (ICM) to post-ICM intermediate stage (PICMI) to hESC stage, that have downstream effects on pluripotency states and differentiation?

SUMMARY ANSWER

We reveal that although the PICMI is transcriptionally similar to the hESC profile and distinct from ICM, it exhibits upregulation of primordial germ cell (PGC) markers, dependence on leukemia inhibitory factor (LIF) signaling, upregulation of naïve pluripotency-specific signaling networks and appears to be an intermediate switching point from naïve to primed pluripotency.

WHAT IS KNOWN ALREADY

It is currently known that the PICMI exhibits markers of early and late-epiblast stage. It is suggested that hESCs acquire primed pluripotency features due to the upregulation of post-implantation genes in the PICMI which renders them predisposed towards differentiation cues. Despite this current knowledge, the transcriptional landscape changes during hESC derivation from ICM to hESC and the effect of PICMI on pluripotent state is still not well defined.

STUDY DESIGN, SIZE, DURATION

To gain insight into the signaling mechanisms that may govern the ICM to PICMI to hESC transition, comparative RNA sequencing (RNA-seq) analysis was performed on preimplantation ICMs, PICMIs and hESCs in biological and technical triplicates (n = 3).

PARTICIPANTS/MATERIALS, SETTING, AND METHODS

Primed hESCs (XX) were maintained in feeder-free culture conditions on Matrigel for two passages and approximately 50 cells were collected in biological and technical triplicates (n = 3). For ICM sample collection, Day 3, frozen-thawed human embryos were cultured up to day five blastocyst stage and only good quality blastocysts were subjected to laser-assisted micromanipulation for ICM collection (n = 3). Next, day six expanded blastocysts were cultured on mouse embryonic fibroblasts and manual dissection was performed on the PICMI outgrowths between post-plating Day 6 and Day 10 (n = 3). Sequencing of these samples was performed on NextSeq500 and statistical analysis was performed using edgeR (false discovery rate (FDR) < 0.05).

MAIN RESULTS AND THE ROLE OF CHANCE

Comparative RNA-seq data analysis revealed that 634 and 560 protein-coding genes were significantly up and downregulated in hESCs compared to ICM (FDR < 0.05), respectively. Upon ICM to PICMI transition, 471 genes were expressed significantly higher in the PICMI compared to ICM, while 296 genes were elevated in the ICM alone (FDR < 0.05). Principle component analysis showed that the ICM was completely distinct from the PICMI and hESCs while the latter two clustered in close proximity to each other. Increased expression of E-CADHERIN1 (CDH1) in ICM and intermediate levels in the PICMI was observed, while CDH2 was higher in hESCs, suggesting a role of extracellular matrix components in facilitating pluripotency transition during hESC derivation. The PICMI also showed regulation of naïve-specific LIF and bone morphogenetic protein signaling, differential regulation of primed pluripotency-specific fibroblast growth factor and NODAL signaling pathway components, upregulation of phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway (PI3K/AKT/mTORC), as well as predisposition towards the germ cell lineage, further confirmed by gene ontology analysis. Hence, the data suggest that the PICMI may serve as an intermediate pluripotency stage which, when subjected to an appropriate culture niche, could aid in enhancing naïve hESC derivation and germ cell differentiation efficiency.

LARGE-SCALE DATA

Gene Expression Omnibus (GEO) Accession number GSE119378.

LIMITATIONS, REASONS FOR CAUTION

Owing to the limitation in sample availability, the sex of ICM and PICMI have not been taken into consideration. Obtaining cells from the ICM and maintaining them in culture is not feasible as it will hamper the formation of PICMI and hESC derivation. Single-cell quantitative real-time PCR on low ICM and PICMI cell numbers, although challenging due to limited availability of human embryos, will be advantageous to further corroborate the RNA-seq data on transcriptional changes during hESC derivation process.

WIDER IMPLICATIONS OF THE FINDINGS

We elucidate the dynamics of transcriptional network changes from the naïve ICM to the intermediate PICMI stage and finally the primed hESC lines. We provide an in-depth understanding of the PICMI and its role in conferring the type of pluripotent state which may have important downstream effects on differentiation, specifically towards the PGC lineage. This knowledge contributes to our limited understanding of the true nature of the human pluripotent state in vitro.

STUDY FUNDING/COMPETING INTEREST(S)

This research is supported by the Concerted Research Actions funding from Bijzonder Onderzoeksfonds University Ghent (BOF GOA 01G01112).The authors declare no conflict of interest.

Gain of 20q11.21 in Human Pluripotent Stem Cells Impairs TGF-β-Dependent Neuroectodermal Commitment

Gain of 20q11.21 is one of the most common recurrent genomic aberrations in human pluripotent stem cells. Although it is known that overexpression of the antiapoptotic gene Bcl-xL confers a survival advantage to the abnormal cells, their differentiation capacity has not been fully investigated. RNA sequencing of mutant and control hESC lines, and a line transgenically overexpressing Bcl-xL, shows that overexpression of Bcl-xL is sufficient to cause most transcriptional changes induced by the gain of 20q11.21. Moreover, the differentially expressed genes in mutant and Bcl-xL overexpressing lines are enriched for genes involved in TGF-β- and SMAD-mediated signaling, and neuron differentiation. Finally, we show that this altered signaling has a dramatic negative effect on neuroectodermal differentiation, while the cells maintain their ability to differentiate to mesendoderm derivatives. These findings stress the importance of thorough genetic testing of the lines before their use in research or the clinic.

•

Bcl-xL overexpression drives the transcriptomic profile of 20q11.21 mutant lines

•

20q11.21 mutant lines downregulate CHCHD2, a known TGF-β pathway modulator

•

Mutant lines differentially express genes involved in TGF-β and SMAD signaling

•

Mutant lines show impaired ectoderm commitment due to TGF-β signaling deregulation

WNT Inhibition and Increased FGF Signaling Promotes Derivation of Less Heterogeneous Primed Human Embryonic Stem Cells, Compatible with Differentiation

Human embryonic stem cells (hESCs) hold great value for future clinical applications. However, standard culture conditions maintain hESCs in a primed state, which bears heterogeneity in pluripotency and a tendency for spontaneous differentiation. To counter these drawbacks, primed hESCs have been converted to a naive state, but this has restricted the efficiency of existing directed differentiation protocols. In mouse, WNT inhibition by inhibitor of WNT production-2, together with a higher dose of fibroblast growth factor 2 (12 ng/mL) in DMEM/F12 basal medium (DhiFI), markedly improved derivation and maintenance of primed mouse epiblast stem cells. In this study, we show that DhiFI conditions similarly improved primed hESC traits, such as conferring a primed transcriptional signature with high levels of pluripotency markers and reduced levels of differentiation markers. When triggered to differentiate to neuronal and cardiac lineages, DhiFI hESCs and isogenic primed hESCs progressed similarly. Moreover, DhiFI conditions supported the derivation of hESC lines from a post-inner cell mass intermediate (PICMI). DhiFI-derived hESCs showed less spontaneous differentiation and expressed significantly lower levels of lineage-specific markers, compared to primed-derived lines from the same PICMI. Overall, DhiFI hESCs retained advantages of both primed and naive pluripotency and may ultimately represent a more favorable starting point for differentiation toward clinically desired cell types.

Contemporary pharmacogenetic assays in view of the PharmGKB database

AIM

Six modern PGx assays were compared with the Pharmacogenomics Knowledge Base (PharmGKB) to determine the proportion of the currently known PGx genotypes that are assessed by these assays.

MATERIALS & METHODS

Investigated assays were 'Ion AmpliSeq Pharmacogenomics', 'iPLEX PGx Pro', 'DMET Plus,' 'PharmcoScan,' 'Living DNA' and '23andMe.'

RESULTS

PharmGKB contains 3474 clinical annotations of which 75, 70 and 45% can be determined by PharmacoScan, Living DNA and 23andMe, respectively. The other assays are designed to test a specific subset of PGx variants.

CONCLUSION

Assaying all known PGx variants would only comprise a minor fraction of the current assays' capacity. Unfortunately, this is not achieved. Moreover, not necessarily the variants with the highest effects or the highest evidence are selected.

The transcriptome of Mycobacterium tuberculosis in a lipid-rich dormancy model through RNAseq analysis

Tuberculosis (TB) is currently the number one killer among infectious diseases worldwide. Lipids are abundant molecules during the infectious cycle of Mycobacterium tuberculosis (Mtb) and studies better mimicking its actual metabolic state during pathogenesis are needed. Though most studies have focused on the mycobacterial lipid metabolism under standard culture conditions, little is known about the transcriptome of Mtb in a lipid environment. Here we determined the transcriptome of Mtb H37Rv in a lipid-rich environment (cholesterol and fatty acid) under aerobic and hypoxic conditions, using RNAseq. Lipids significantly induced the expression of 368 genes. A main core lipid response was observed involving efflux systems, iron caption and sulfur reduction. In co-expression with ncRNAs and other genes discussed below, may act coordinately to prepare the machinery conferring drug tolerance and increasing a persistent population. Our findings could be useful to tag relevant pathways for the development of new drugs, vaccines and new strategies to control TB.

STR profiling and Copy Number Variation analysis on single, preserved cells using current Whole Genome Amplification methods

The growing interest in liquid biopsies for cancer research and cell-based non-invasive prenatal testing (NIPT) invigorates the need for improved single cell analysis. In these applications, target cells are extremely rare and fragile in peripheral circulation, which makes the genetic analysis very challenging. To overcome these challenges, cell stabilization and unbiased whole genome amplification are required. This study investigates the performance of four WGA methods on single or a limited number of cells after 24 hour of Streck Cell-Free DNA BCT preservation. The suitability of the DNA, amplified with Ampli1, DOPlify, PicoPLEX and REPLI-g, was assessed for both short tandem repeat (STR) profiling and copy number variant (CNV) analysis after shallow whole genome massively parallel sequencing (MPS). Results demonstrate that Ampli1, DOPlify and PicoPLEX perform well for both applications, with some differences between the methods. Samples amplified with REPLI-g did not result in suitable STR or CNV profiles, indicating that this WGA method is not able to generate high quality DNA after Streck Cell-Free DNA BCT stabilization of the cells.

Performance of four modern whole genome amplification methods for copy number variant detection in single cells

Whole genome amplification (WGA) has become an invaluable tool to perform copy number variation (CNV) detection in single, or a limited number of cells. Unfortunately, current WGA methods introduce representation bias that limits the detection of small CNVs. New WGA methods have been introduced that might have the potential to reduce this bias. We compared the performance of PicoPLEX DNA-Seq (Picoseq), DOPlify, REPLI-g and Ampli-1 WGA for aneuploidy screening and copy number analysis using shallow whole genome massively parallel sequencing (MPS), starting from single or a limited number of cells. Although the four WGA methods perform differently, they are all suited for this application.