An in-depth investigation of the microbiota and its virulence factors associated with severe udder cleft dermatitis lesions

Udder cleft dermatitis (UCD) is a skin condition affecting the anterior parts of the udder in dairy cattle. In the present study, we aimed to shed light on the microbiota in severe UCD lesions versus healthy udder skin by putting forward a taxonomic and functional profile based on a virulence factor analysis. Through shotgun metagenomic sequencing, we found a high proportion of bacteria in addition to a low abundance of archaea. A distinct clustering of healthy udder skin versus UCD lesion samples was shown by applying principal component analysis and (sparse) partial least squares analysis on the metagenomic data. Proteobacteria, Bacillota, and Actinomycetota were among the most abundant phyla in healthy udder skin samples. In UCD samples, Bacteroidota was the most abundant phylum. At genus level, Bifidobacterium spp. had the highest relative abundance in healthy skin samples, whereas Porphyromonas spp. and Corynebacterium spp. had the highest relative abundance in UCD samples. In the differential abundance analysis, Porphyromonas spp. and Bacteroides spp. were significantly differentially abundant in UCD samples, whereas Bifidobacterium spp., Staphylococcus sp. AntiMn-1, and Staphylococcus equorum were more commonly found in healthy samples. Moreover, the abundance of several treponeme phylotypes was significantly higher in lesion samples. The streptococcal cysteine protease speB was among the most abundant virulence factors present in severe UCD lesions, while a plethora of virulence factors such as the antitoxin relB were downregulated, possibly contributing to creating the ideal wound climate for the dysbiotic community. Network analysis showed healthy lesion samples had a large network ofpositive, correlations between the abundances of beneficial species such as Aerococcus urinaeequi and Bifidobacterium angulatum, indicating that the healthy skin microbiome forms an active protective bacterial network, which is disrupted in case of UCD. In UCD samples, a smaller microbial network mainly consisting of positive correlations between the abundances of Bacteroides fragilis and anaerobic Bacteroidota was exposed. Moreover, a high correlation between the taxonomic data and virulence factors was revealed, concurrently with 2 separate networks of microbes and virulence factors. One network, matching with the taxonomic findings in the healthy udder skin samples, showcased a community of harmless or beneficial bacteria, such as Bifidobacterium spp. and Butyrivibrio proteoclasticus, associated with hcnB, hcnC, relB, glyoxalase, and cupin 2. The other network, corresponding to UCD samples, consisted of pathogenic or facultative pathogenic and mainly anaerobic bacteria such as Treponema spp., Mycoplasmopsis spp., and bovine gammaherpesvirus 4, that correlated with virulence factors SpvB, fhaB, and haemagglutination activity domain-associated factor. Our results point toward a dysbiotic community with a notable decrease in diversity and evenness, with a loss of normal skin inhabitants and innocuous or useful species making way for predominantly anaerobic, facultative pathogens. The shift in the abundance of virulence factors such as fhaB and SpvB could play a role in the manifestation of a local micro-environment favorable to the microbiome associated with udder skin lesions. Lastly, the presence of specific networks between microbial species, and between microbes and virulence factors was shown.

Mycothione reductase as a potential target in the fight against Mycobacterium abscessus infections

While infections caused by Mycobacterium abscessus complex (MABC) are rising worldwide, the current treatment of these infections is far from ideal due to its numerous shortcomings thereby increasing the urge for novel drug targets. In this study, mycothione reductase (Mtr) was evaluated for its potential as a drug target for MABC infections since it is a key enzyme needed in the recycling of mycothiol, the main low-molecular-weight thiol protecting the bacteria against reactive oxygen species and other reactive intermediates. First, a Mab∆mtr mutant strain was generated, lacking mtr expression. Next, the in vitro sensitivity of Mab∆mtr to oxidative stress and antimycobacterial drugs was determined. Finally, we evaluated the intramacrophage survival and the virulence of Mab∆mtr in Galleria mellonella larvae. Mab∆mtr demonstrated a 39.5-fold reduction in IC90 when exposed to bedaquiline in vitro. Furthermore, the Mab∆mtr mutant showed a decreased ability to proliferate inside macrophages and larvae, suggesting that Mtr plays an important role during MABC infection. Altogether, these findings support the assumption of Mtr being a potential target for antimycobacterial drugs.IMPORTANCEMycobacterium abscessus complex (MABC) is a group of bacteria causing a serious public health problem worldwide due to its ability to cause progressive disease, its highly resistant profile against various antibiotics, and its lengthy treatment. Therefore, new drugs are needed to alleviate antibiotic resistance and reduce the length of the current treatment. A potential new target for new antibiotics is mycothione reductase (Mtr), an important enzyme belonging to a pathway that protects the bacteria against harmful conditions. Our research created a bacterium deficient of mtr by using advanced genetic techniques and demonstrated that mtr-deficient bacteria have a decreased ability to multiply during infection. Furthermore, we show evidence that currently used antibiotics combined with mtr deficiency can lead to a better treatment of MABC infection. Altogether, our results validate Mtr as a potential new target and suggest that Mtr plays a role during MABC infection.

Severe udder cleft dermatitis lesion transcriptomics points to an impaired skin barrier, defective wound repair and a dysregulated inflammatory response as key elements in the pathogenesis

This study is the first to investigate the transcriptomic changes occurring in severe udder cleft dermatitis lesions (UCD) in Holstein-Friesian cows. An examination of the gene expression levels in natural UCD lesions and healthy udder skin through RNA Seq-Technology provided a deeper insight into the inflammatory pathways associated with this disease. A clear distinction between the gene expression patterns of UCD lesions and healthy skin was shown in the principal component analysis. Genes coding for inflammatory molecules were upregulated such as the chemokines C-X-C motif ligand 2 (CXCL2), 5 (CXCL5) and 8 (CXCL8), and C-C motif ligand 11 (CCL11). Moreover, the genes coding for the multifunctional molecules ADAM12 and SLPI were amongst the highest upregulated ones, whereas the most downregulated genes included the ones coding for keratins and keratin-associated molecules. Predominantly inflammatory pathways such as the chemokine signaling, cytokine receptor interaction and IL-17 signaling pathway were significantly upregulated in the pathway analysis. These results point towards a fulminant, dysregulated inflammatory response concomitant with a disruption of the skin barrier integrity and a hampered wound repair mechanism in severe UCD lesions.

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.

O-099 TEAD4 regulates trophectoderm differentiation upstream of CDX2 in human preimplantation embryos

Study question

What is the main pathway regulating trophectoderm (TE) differentiation during pre-implantation development in mouse versus human embryos?

Summary answer

TEAD4 is acting upstream of CDX2 and is involved in TE differentiation, as TEAD4-null human embryos exhibit compromised TE lineage differentiation.

What is known already

TEAD4 is the earliest transcription factor during early embryo development, required for the expression of TE-associated genes leading to successful TE differentiation and subsequent blastocoel formation in mouse. Functional knock-out studies in mouse, inactivating Tead4 by site-specific recombination have shown that Tead4-null embryos do not express TE specific genes, including Caudal-Type Homeobox Protein 2 (Cdx2) and GATA Binding Protein 3 (Gata3), but expression of inner cell mass (ICM)-specific genes, remains unaffected. Furthermore, ablation of Tead4 compromises embryonic development and subsequent blastocoel formation in mouse. The role of TEAD4, during human pre-implantation development has not been functionally characterized yet.

Study design, size, duration

CRISPR-Cas9 was introduced in mouse zygotes and editing efficiency was evaluated by next-generation sequencing (NGS) on 4.5dpc embryos (n = 55). Developmental kinetics were monitored in CRISPR-Cas9 targeted (n = 83), sham-injected (n = 26) and non-injected media-control (n = 51) mouse embryos. Immunofluorescence analysis was performed in Tead4 targeted (n = 57) and non-injected media-control embryos (n = 94). The same methodology was applied in human donated in vitro matured (IVM) metaphase-II (MII) oocytes, which were CRISPR-Cas9 targeted (n = 74) during ICSI or used as media-Control (n = 33).

Participants/materials, setting, methods

A gRNA-Cas9 mixture targeting exon 2 of Tead4/TEAD4 was microinjected in respectively mouse 2PN (pronuclear) stage zygotes, or human IVM MII oocytes along with the sperm. Generated embryos were cultured in vitro for 4 days in mouse or 6.5 days in human. Embryonic development and morphology were assessed daily, followed by a detailed scoring at the late blastocyst stage. Successful targeting following CRISPR-Cas9 introduction was assessed by immunostaining and NGS analysis of the targeted locus.

Main results and the role of chance

In mouse, we confirmed previous findings, as the developmental capacity of Tead4 targeted embryos was significantly reduced starting from the morula stage and blastocyst formation rates were 8.97% in the targeted group, compared to 87.23% in the control and 87.50% in the sham group, respectively. Immunofluorescence analysis of late morula and blastocyst stage embryos confirmed the absence of Tead4, Cdx2 and Gata3, resulting from the successful interruption of the Tead4 locus (n = 57). Exon 2 of TEAD4 was successfully targeted in human. In total, 21 embryos from various developmental stages were successfully NGS analyzed and 90,48% (19 out of 21) of the embryos carried genetic modifications as a result of CRISPR-Cas9 genome editing and seven blastocysts were identified carrying exclusively frameshift mutations. In contrast to mouse, the developmental capacity of human targeted embryos (25%) did not differ significantly from the control group (23%). However, the blastocyst morphology and quality were compromised in the targeted group showing mostly grade C TE scores, containing very few cells. Immunofluorescence analysis of targeted blastocysts (n = 6) confirmed successful editing by complete absence of TEAD4 and its downstream TE marker CDX2.

Limitations, reasons for caution

CRISPR-Cas9 germline genome editing results in multiple editing outcomes with variable phenotypic penetrance, the mosaic nature of which complicates the phenotypic analysis and developmental behaviour of the injected embryos.

Wider implications of the findings

Elucidation of the evolutionary conserved molecular mechanisms that regulate self-renewal of the trophoblast lineage can give us fundamental insights on early implantation failure.

Trial registration number

Not Applicable

Comparative analysis of mouse and human preimplantation development following POU5F1 CRISPR/Cas9 targeting reveals interspecies differences

STUDY QUESTION

What is the role of POU class 5 homeobox 1 (POU5F1) in human preimplantation development and how does it compare with the mouse model?

SUMMARY ANSWER

POU5F1 is required for successful development of mouse and human embryos to the blastocyst stage as knockout embryos exhibited a significantly lower blastocyst formation rate, accompanied by lack of inner cell mass (ICM) formation.

WHAT IS KNOWN ALREADY

Clustered regularly interspaced short palindromic repeats-CRISPR associated genes (CRISPR-Cas9) has previously been used to examine the role of POU5F1 during human preimplantation development. The reported POU5F1-targeted blastocysts always retained POU5F1 expression in at least one cell, because of incomplete CRISPR-Cas9 editing. The question remains of whether the inability to obtain fully edited POU5F1-targeted blastocysts in human results from incomplete editing or the actual inability of these embryos to reach the blastocyst stage.

STUDY DESIGN, SIZE, DURATION

The efficiency of CRISPR-Cas9 to induce targeted gene mutations was first optimized in the mouse model. Two CRISPR-Cas9 delivery methods were compared in the B6D2F1 strain: S-phase injection (zygote stage) (n = 135) versus metaphase II-phase (M-phase) injection (oocyte stage) (n = 23). Four control groups were included: non-injected media-control zygotes (n = 43)/oocytes (n = 48); sham-injected zygotes (n = 45)/oocytes (n = 47); Cas9-protein injected zygotes (n = 23); and Cas9 protein and scrambled guide RNA (gRNA)-injected zygotes (n = 27). Immunofluorescence analysis was performed in Pou5f1-targeted zygotes (n = 37), media control zygotes (n = 19), and sham-injected zygotes (n = 15). To assess the capacity of Pou5f1-null embryos to develop further in vitro, additional groups of Pou5f1-targeted zygotes (n = 29) and media control zygotes (n = 30) were cultured to postimplantation stages (8.5 dpf). Aiming to identify differences in developmental capacity of Pou5f1-null embryos attributed to strain variation, zygotes from a second mouse strain-B6CBA (n = 52) were targeted. Overall, the optimized methodology was applied in human oocytes following IVM (metaphase II stage) (n = 101). The control group consisted of intracytoplasmically sperm injected (ICSI) IVM oocytes (n = 33). Immunofluorescence analysis was performed in human CRISPR-injected (n = 10) and media control (n = 9) human embryos.

PARTICIPANTS/MATERIALS, SETTING, METHODS

A gRNA-Cas9 protein mixture targeting exon 2 of Pou5f1/POU5F1 was microinjected in mouse oocytes/zygotes or human IVM oocytes. Reconstructed embryos were cultured for 4 days (mouse) or 6.5 days (human) in sequential culture media. An additional group of mouse-targeted zygotes was cultured to postimplantation stages. Embryonic development was assessed daily, with detailed scoring at late blastocyst stage. Genomic editing was assessed by immunofluorescence analysis and next-generation sequencing.

MAIN RESULTS AND THE ROLE OF CHANCE

Genomic analysis in mouse revealed very high editing efficiencies with 95% of the S-Phase and 100% of the M-Phase embryos containing genetic modifications, of which 89.47% in the S-Phase and 84.21% in the M-Phase group were fully edited. The developmental capacity was significantly compromised as only 46.88% embryos in the S-Phase and 19.05% in the M-Phase group reached the blastocyst stage, compared to 86.36% in control M-Phase and 90.24% in control S-Phase groups, respectively. Immunofluorescence analysis confirmed the loss of Pou5f1 expression and downregulation of the primitive marker SRY-Box transcription factor (Sox17). Our experiments confirmed the requirement of Pou5f1 expression for blastocyst development in the second B6CBA strain. Altogether, our data obtained in mouse reveal that Pou5f1 expression is essential for development to the blastocyst stage. M-Phase injection in human IVM oocytes (n = 101) similarly resulted in 88.37% of the POU5F1-targeted embryos being successfully edited. The developmental capacity of generated embryos was compromised from the eight-cell stage onwards. Only 4.55% of the microinjected embryos reached the late blastocyst stage and the embryos exhibited complete absence of ICM and an irregular trophectoderm cell layer. Loss of POU5F1 expression resulted in absence of SOX17 expression, as in mouse. Interestingly, genetic mosaicism was eliminated in a subset of targeted human embryos (9 out of 38), three of which developed into blastocysts.

LIMITATIONS, REASONS FOR CAUTION

One of the major hurdles of CRISPR-Cas9 germline genome editing is the occurrence of mosaicism, which may complicate phenotypic analysis and interpretation of developmental behavior of the injected embryos. Furthermore, in this study, spare IVM human oocytes were used, which may not recapitulate the developmental behavior of in vivo matured oocytes.

WIDER IMPLICATIONS OF THE FINDINGS

Comparison of developmental competency following CRISPR-Cas-mediated gene targeting in mouse and human may be influenced by the selected mouse strain. Gene targeting by CRISPR-Cas9 is subject to variable targeting efficiencies. Therefore, striving to reduce mosaicism can provide novel molecular insights into mouse and human embryogenesis.

STUDY FUNDING/COMPETING INTEREST(S)

The research was funded by the Ghent University Hospital and Ghent University and supported by the FWO-Vlaanderen (Flemish fund for scientific research, Grant no. G051516N), and Hercules funding (FWO.HMZ.2016.00.02.01). The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Effect of overconditioning on the hepatic global gene expression pattern of dairy cows at the end of pregnancy

Overconditioning is a risk factor for upregulated pre- and postpartum fat mobilization. Therefore, we hypothesized that overconditioning at the end of pregnancy leads to the accumulation of lipids in the liver and modifications of the hepatic gene expression pattern. The aim of this study was to evaluate the effect of normal- versus overconditioning on the hepatic transcriptomic profile of dairy cows at the end of pregnancy. Ten dry multiparous Holstein cows were killed 2 wk before expected calving. Body condition score (BCS) and backfat thickness (BFT) were evaluated, and blood samples for nonesterified fatty acids (NEFA) were taken before cows were killed. After cows were killed, liver biopsy samples were collected for further assessment of total lipids and RNA sequencing. Five cows were classified as normal-conditioned (median BCS = 3, range 2.75-3.5) and 5 as overconditioned (median BCS = 4, range 4-5). Regression models confirmed that normal-conditioned cows had lower BFT (1.29 ± 0.29 cm; least squares means ± standard error) and serum NEFA (0.16 ± 0.04 mmol/L) in comparison to overconditioned cows (3.14 ± 0.43 cm and 0.38 ± 0.07 mmol/L for BFT and NEFA, respectively). Total liver lipid percentage tended to be lower in normal- versus overconditioned cows (4.63 ± 0.40% and 6.06 ± 0.44%, respectively). In comparison to the mean liver lipid percentage of the normal- and overconditioned cows, 1 overconditioned cow had a relatively low (5.21%) and 1 normal-conditioned cow had a relatively high (6.07%) liver lipid percentage. Differentially expressed genes analysis (edgeR quasi-likelihood method) showed that normal-conditioned cows presented 11 upregulated and 12 downregulated genes in comparison to overconditioned cows. Linear discriminant analysis effects size revealed 133 differentially expressed genes between normal- versus overconditioned cows. Notably, the liver of normal-conditioned cows had upregulated genes associated with liver functionality (ALB, SELENOP, IGF1, and IGF2). On the other hand, overconditioned cows had upregulated genes associated with the acute-phase response (C3, HPX, and, LBP). High basal lipolysis in overconditioned cows at the end of pregnancy increased liver lipid content, and this may alter the hepatic gene expression pattern to a pro-inflammatory state.

Comparative analysis of different nuclear transfer techniques to prevent the transmission of mitochondrial DNA variants

Prevention of mitochondrial DNA (mtDNA) diseases may currently be possible using germline nuclear transfer (NT). However, scientific evidence to compare efficiency of different NT techniques to overcome mtDNA diseases is lacking. Here, we performed four types of NT, including first or second polar body transfer (PB1/2T), maternal spindle transfer (ST) and pronuclear transfer (PNT), using NZB/OlaHsd and B6D2F1 mouse models. Embryo development was assessed following NT, and mtDNA carry-over levels were measured by next generation sequencing (NGS). Moreover, we explored two novel protocols (PB2T-a and PB2T-b) to optimize PB2T using mouse and human oocytes. Chromosomal profiles of NT-generated blastocysts were evaluated using NGS. In mouse, our findings reveal that only PB2T-b successfully leads to blastocysts. There were comparable blastocyst rates among PB1T, PB2T-b, ST and PNT embryos. Furthermore, PB1T and PB2T-b had lower mtDNA carry-over levels than ST and PNT. After extrapolation of novel PB2T-b to human in vitro matured (IVM) oocytes and in vivo matured oocytes with smooth endoplasmic reticulum aggregate (SERa) oocytes, the reconstituted embryos successfully developed to blastocysts at a comparable rate to ICSI controls. PB2T-b embryos generated from IVM oocytes showed a similar euploidy rate to ICSI controls. Nevertheless, our mouse model with non-mutated mtDNAs is different from a mixture of pathogenic and non-pathogenic mtDNAs in a human scenario. Novel PB2T-b requires further optimization to improve blastocyst rates in human. Although more work is required to elucidate efficiency and safety of NT, our study suggests that PBT may have the potential to prevent mtDNA disease transmission.

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

Chromosomal mosaicism in human blastocysts: the ultimate challenge of preimplantation genetic testing?

STUDY QUESTION

To what extent does a trophectoderm (TE) biopsy reliably reflect the chromosomal constitution of the inner cell mass (ICM) in human blastocysts?

SUMMARY ANSWER

Concordance between TE and ICM was established in 62.1% of the embryos analysed.

WHAT IS KNOWN ALREADY

Next generation sequencing (NGS) platforms have recently been optimised for preimplantation genetic testing for aneuploidies (PGT-A). However, higher sensitivity has led to an increase in reports of chromosomal mosaicism within a single TE biopsy. This has raised substantial controversy surrounding the prevalence of mosaicism in human blastocysts and the clinical implications of heterogeneity between the TE and ICM.

STUDY DESIGN, SIZE, DURATION

To define the distribution and rate of mosaicism in human blastocysts, we assessed chromosomal profiles of the ICM and multiple TE portions obtained from the same embryo. We evaluated donated embryos with an unknown chromosomal profile (n = 34), as well as PGT-A blastocysts, previously diagnosed as abnormal or mosaic (n = 24). Our intra-embryo comparison included a total of 232 samples, obtained from 58 embryos.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Four embryo samples, including the ICM and three distinct TE portions, were acquired from good quality blastocysts by micromanipulation. Whole genome amplification (WGA), followed by NGS was performed on all embryo segments. Profiles were compared between samples from the same embryo, while the results from pretested blastocysts were further correlated to the original report. The embryos investigated in our untested group were obtained from good prognosis patients (n = 25), with maternal age ranging from 23 to 39 years. For the pretested embryo group, maternal age ranged from 23 to 40 years (n = 18).

MAIN RESULTS AND THE ROLE OF CHANCE

We uncover chromosomal mosaicism, involving both numerical and structural aberrations, in up to 37.9% of the blastocysts analysed. Within the untested group, the overall concordance between the ICM and all TE portions was 55.9%. A normal ICM was detected in 20.6% of blastocysts for which at least one TE portion showed a chromosomal aberration. Conversely, 17.6% of embryos presented with mosaic or uniform abnormalities within the ICM, while showing normal or mosaic TE profiles. For the pretested blastocysts, the overall concordance between the ICM and all TE samples was 70.8%. However, 50% of embryos previously diagnosed with mosaicism did not confirm the original diagnosis. Notably, 31.3% of embryos with a mosaic aberration reported in the original TE biopsy, revealed a euploid profile in the ICM and all three TE samples. Taken together, concordance between the ICM and all TE portions was established in 62.1% of blastocysts, across both embryo groups. Finally, we could not observe a significant effect of age on embryo mosaicism (P = 0.101 untested group; P = 0.7309 pretested group). Similarly, ICM and TE quality were not found to affect the occurrence of chromosomal mosaicism (P = 0.718 and P = 0.462 untested group; P = 1.000 and P = 0.2885 pretested group).

LARGE SCALE DATA

All data that support the findings of this study are available online in Vivar (http://cmgg.be/vivar) upon request.

LIMITATIONS, REASONS FOR CAUTION

Evaluating biological variation in some instances remains challenging. The technological limitations of sampling mitotic errors that lead to mosaicism, as well as WGA artefacts, warrant careful interpretation.

WIDER IMPLICATIONS OF THE FINDINGS

Our results highlight the complex nature of genetic (in)stability during early ontogenesis and indicate that blastocysts harbour a higher rate of chromosomal mosaicism than may have been anticipated. Moreover, our findings reveal an overall high diagnostic sensitivity and relatively low specificity in the context of PGT-A. This suggests that a considerable proportion of embryos are potentially being classified as clinically unsuitable. Ultimately, more precise quantification will benefit the clinical management of embryo mosaicism.

STUDY FUNDING/COMPETING INTEREST(S)

M.P. is supported by the Special Research Fund, Bijzonder Onderzoeksfonds (BOF01D08114). J.T. and L.D. are supported by the agency for innovation through science (131673, 141441). B.H. and this research are supported by the Special Research Fund, Bijzonder Onderzoeksfonds (BOF15/GOA/011). The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

Not applicable.

Toward the most ideal case-control design with related and unrelated dogs in whole-exome sequencing studies

With the recent development of whole-exome sequencing enrichment designs for the dog, a novel tool for disease-association studies became available. The aim of disease-association studies is to identify one or a very limited number of putative causal variants or genes from the large pool of genetic variation. To maximize the efficiency of these studies and to provide some directions of what to expect, we evaluated the effect on variant reduction for various combinations of cases and controls for both dominant and recessive types of inheritance assuming variable degrees of penetrance and detectance. In this study, variant data of 14 dogs (13 Labrador Retrievers and one Dogue de Bordeaux), obtained by whole-exome sequencing, were analyzed. In the filtering process, we found that unrelated dogs from the same breed share up to 70% of their variants, which is likely a consequence of the breeding history of the dog. For the designs tested with unrelated dogs, combining two cases and two controls gave the best result. These results were improved further by adding closely related dogs. Reduced penetrance and/or detectance has a drastic effect on the efficiency and is likely to have a profound effect on the sample size needed to elucidate the causal variant. Overall, we demonstrated that sequencing a small number of dogs results in a marked reduction of variants that are likely sufficient to pinpoint causal variants or genes.

Challenges and opportunities in the investigation of unexplained intellectual disability using family-based whole-exome sequencing

Intellectual disability (ID), characterized by an intellectual performance of at least 2 SD (standard deviations) below average is a frequent, lifelong disorder with a prevalence of 2-3%. Today, only for at most half of patients a diagnosis is made. Knowing the cause of the ID is important for patients and their relatives, as it allows for appropriate medical care, prognosis on further development of the disorder, familial counselling or access to support groups. Whole-exome sequencing (WES) now offers the possibility to identify the genetic cause for patients for which all previously available genetic tests, including karyotyping, specific gene analysis, or microarray analysis did not reveal causative abnormalities. However, data analysis of WES experiments is challenging. Here we present an analysis workflow implementable in any laboratory, requiring no bioinformatics knowledge. We demonstrated its feasibility on a cohort of 10 patients, in which we found a conclusive diagnosis in 3 and a likely diagnosis in 2 more patients. Of the three conclusive diagnoses, one was a clinically suspected mutation missed by Sanger sequencing, and one was an atypical presentation of a known monogenic disorder, highlighting two essential strengths of WES-based diagnostics.

Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation

Background

Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth.

Results

Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism.

Conclusions

Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation.

Impact of allelic dropout on evidential value of forensic DNA profiles using RMNE

Motivation: Two methods are commonly used to report on evidence carried by forensic DNA profiles: the ‘Random Man Not Excluded’ (RMNE) approach and the likelihood ratio (LR) approach. It is often claimed a major advantage of the LR method that dropout can be assessed probabilistically.

Results: In this article, a new RMNE measure is proposed that like-wise accounts for allelic dropout in an observed forensic DNA profile. We discuss the necessary calculations, underline their simplicity and provide a tool for performing the calculations.

Availability: An Excel file with preprogrammed calculations of RMNE probabilities for DNA profiles up to 16 loci and with a maximum of two dropouts is available at: http://www.labfbt.UGent.be/RMNE.php

Contact:dieter.deforce@ugent.be