Cell specification and functional interactions in the pig blastocyst inferred from single-cell transcriptomics and uterine fluids proteomics.

The embryonic development of the pig comprises a long in utero pre- and peri-implantation development, which dramatically differs from mice and humans. During this peri-implantation period, a complex series of paracrine signals establishes an intimate dialogue between the embryo and the uterus. To better understand the biology of the pig blastocyst during this period, we generated a large dataset of single-cell RNAseq from early and hatched blastocysts, spheroid and ovoid conceptus and proteomic datasets from corresponding uterine fluids. Our results confirm the molecular specificity and functionality of the three main cell populations. We also discovered two previously unknown subpopulations of the trophectoderm, one characterised by the expression of LRP2, which could represent progenitor cells, and the other, expressing pro-apoptotic markers, which could correspond to the Rauber's layer. Our work provides new insights into the biology of these populations, their reciprocal functional interactions, and the molecular dialogue with the maternal uterine environment.

Epithelial-mesenchymal transitions in development and disease.

The epithelial to mesenchymal transition (EMT) plays crucial roles in the formation of the body plan and in the differentiation of multiple tissues and organs. EMT also contributes to tissue repair, but it can adversely cause organ fibrosis and promote carcinoma progression through a variety of mechanisms. EMT endows cells with migratory and invasive properties, induces stem cell properties, prevents apoptosis and senescence, and contributes to immunosuppression. Thus, the mesenchymal state is associated with the capacity of cells to migrate to distant organs and maintain stemness, allowing their subsequent differentiation into multiple cell types during development and the initiation of metastasis.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates.

BACKGROUND: Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken). RESULTS: Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates. CONCLUSIONS: Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

BMAL1 Knockdown Leans Epithelial-Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells.

The circadian clock coordinates biological and physiological functions to day/night cycles. The perturbation of the circadian clock increases cancer risk and affects cancer progression. Here, we studied how BMAL1 knockdown (BMAL1-KD) by shRNA affects the epithelial-mesenchymal transition (EMT), a critical early event in the invasion and metastasis of colorectal carcinoma (CRC). In corresponding to a gene set enrichment analysis, which showed a significant enrichment of EMT and invasive signatures in BMAL1_high CRC patients as compared to BMAL1_low CRC patients, our results revealed that BMAL1 is implicated in keeping the epithelial-mesenchymal equilibrium of CRC cells and influences their capacity of adhesion, migration, invasion, and chemoresistance. Firstly, BMAL1-KD increased the expression of epithelial markers (E-cadherin, CK-20, and EpCAM) but decreased the expression of Twist and mesenchymal markers (N-cadherin and vimentin) in CRC cell lines. Finally, the molecular alterations after BMAL1-KD promoted mesenchymal-to-epithelial transition-like changes mostly appeared in two primary CRC cell lines (i.e., HCT116 and SW480) compared to the metastatic cell line SW620. As a consequence, migration/invasion and drug resistance capacities decreased in HCT116 and SW480 BMAL1-KD cells. Together, BMAL1-KD alerts the delicate equilibrium between epithelial and mesenchymal properties of CRC cell lines, which revealed the crucial role of BMAL1 in EMT-related CRC metastasis and chemoresistance.

Snail2 and Zeb2 repress P-cadherin to define embryonic territories in the chick embryo.

Snail and Zeb transcription factors induce epithelial-to-mesenchymal transition (EMT) in embryonic and adult tissues by direct repression of E-cadherin transcription. The repression of E-cadherin transcription by the EMT inducers Snail1 and Zeb2 plays a fundamental role in defining embryonic territories in the mouse, as E-cadherin needs to be downregulated in the primitive streak and in the epiblast, concomitant with the formation of mesendodermal precursors and the neural plate, respectively. Here, we show that in the chick embryo, E-cadherin is weakly expressed in the epiblast at pre-primitive streak stages where it is substituted for by P-cadherin We also show that Snail2 and Zeb2 repress P-cadherin transcription in the primitive streak and the neural plate, respectively. This indicates that E- and P-cadherin expression patterns evolved differently between chick and mouse. As such, the Snail1/E-cadherin axis described in the early mouse embryo corresponds to Snail2/P-cadherin in the chick, but both Snail factors and Zeb2 fulfil a similar role in chick and mouse in directly repressing ectodermal cadherin genes to contribute to the delamination of mesendodermal precursors at gastrulation and the proper specification of the neural ectoderm during neural induction.

Preimplantation development in ungulates: a 'ménage à quatre' scenario.

In ungulates, early embryonic development differs dramatically from that of mice and humans and is characterized by an extended period of pre- and peri-implantation development in utero. After hatching from the zona pellucida, the ungulate blastocyst will stay free in the uterus for many days before implanting within the uterine wall. During this protracted peri-implantation period, an intimate dialog between the embryo and the uterus is established through a complex series of paracrine signals. The blastocyst elongates, leading to extreme growth of extra-embryonic tissues, and at the same time, the inner cell mass moves up into the trophoblast and evolves into the embryonic disc, which is directly exposed to molecules present in the uterine fluids. In the peri-implantation period, uterine glands secrete a wide range of molecules, including enzymes, growth factors, adhesion proteins, cytokines, hormones, and nutrients like amino and fatty acids, which are collectively referred to as histotroph. The identification, role, and effects of these secretions on the biology of the conceptus are still being described; however, the studies that have been conducted to date have demonstrated that histotroph is essential for embryonic development and serves a critical function during the pre- and peri implantation periods. Here, we present an overview of current knowledge on the molecular dialogue among embryonic, extraembryonic, and maternal tissues prior to implantation. Taken together, the body of work described here demonstrates the extent to which this dialog enables the coordination of the development of the conceptus with respect to the establishment of embryonic and extra-embryonic tissues as well as in preparation for implantation.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

Meiotic pairing and gene expression disturbance in germ cells from an infertile boar with a balanced reciprocal autosome-autosome translocation.

Individuals carrying balanced constitutional reciprocal translocations generally have a normal phenotype, but often present reproductive disorders. The aim of our research was to analyze the meiotic process in an oligoasthenoteratospermic boar carrying an asymmetric reciprocal translocation involving chromosomes 1 and 14. Different multivalent structures (quadrivalent and trivalent plus univalent) were identified during chromosome pairing analysis. Some of these multivalents were characterized by the presence of unpaired autosomal segments with histone γH2AX accumulation sometimes associated with the XY body. Gene expression in spermatocytes was studied by RNA-DNA-FISH and microarray-based testis transcriptome analysis. Our results revealed a decrease in gene expression for chromosomes 1 and 14 and an up-regulated expression of X-chromosome genes for the translocated boar compared with normal individuals. We hypothesized that the observed meiotic arrest and reproductive failure in this boar might be due to silencing of crucial autosomal genes (MSUC) and disturbance of meiotic sex chromosome inactivation (MSCI). Further analysis revealed abnormal meiotic recombination (frequency and distribution) and the production of a high rate of unbalanced spermatozoa.

Long noncoding RNA repertoire in chicken liver and adipose tissue.

BACKGROUND: Improving functional annotation of the chicken genome is a key challenge in bridging the gap between genotype and phenotype. Among all transcribed regions, long noncoding RNAs (lncRNAs) are a major component of the transcriptome and its regulation, and whole-transcriptome sequencing (RNA-Seq) has greatly improved their identification and characterization. We performed an extensive profiling of the lncRNA transcriptome in the chicken liver and adipose tissue by RNA-Seq. We focused on these two tissues because of their importance in various economical traits for which energy storage and mobilization play key roles and also because of their high cell homogeneity. To predict lncRNAs, we used a recently developed tool called FEELnc, which also classifies them with respect to their distance and strand orientation to the closest protein-coding genes. Moreover, to confidently identify the genes/transcripts expressed in each tissue (a complex task for weakly expressed molecules such as lncRNAs), we probed a particularly large number of biological replicates (16 per tissue) compared to common multi-tissue studies with a larger set of tissues but less sampling. RESULTS: We predicted 2193 lncRNA genes, among which 1670 were robustly expressed across replicates in the liver and/or adipose tissue and which were classified into 1493 intergenic and 177 intragenic lncRNAs located between and within protein-coding genes, respectively. We observed similar structural features between chickens and mammals, with strong synteny conservation but without sequence conservation. As previously reported, we confirm that lncRNAs have a lower and more tissue-specific expression than mRNAs. Finally, we showed that adjacent lncRNA-mRNA genes in divergent orientation have a higher co-expression level when separated by less than 1 kb compared to more distant divergent pairs. Among these, we highlighted for the first time a novel lncRNA candidate involved in lipid metabolism, lnc_DHCR24, which is highly correlated with the DHCR24 gene that encodes a key enzyme of cholesterol biosynthesis. CONCLUSIONS: We provide a comprehensive lncRNA repertoire in the chicken liver and adipose tissue, which shows interesting patterns of co-expression between mRNAs and lncRNAs. It contributes to improving the structural and functional annotation of the chicken genome and provides a basis for further studies on energy storage and mobilization traits in the chicken.

GO-FAANG meeting: a Gathering On Functional Annotation of Animal Genomes.

The Functional Annotation of Animal Genomes (FAANG) Consortium recently held a Gathering On FAANG (GO-FAANG) Workshop in Washington, DC on October 7-8, 2015. This consortium is a grass-roots organization formed to advance the annotation of newly assembled genomes of domesticated and non-model organisms (www.faang.org). The workshop gathered together from around the world a group of 100+ genome scientists, administrators, representatives of funding agencies and commodity groups to discuss the latest advancements of the consortium, new perspectives, next steps and implementation plans. The workshop was streamed live and recorded, and all talks, along with speaker slide presentations, are available at www.faang.org. In this report, we describe the major activities and outcomes of this meeting. We also provide updates on ongoing efforts to implement discussions and decisions taken at GO-FAANG to guide future FAANG activities. In summary, reference datasets are being established under pilot projects; plans for tissue sets, morphological classification and methods of sample collection for different tissues were organized; and core assays and data and meta-data analysis standards were established.

Sperm DNA methylation analysis in swine reveals conserved and species-specific methylation patterns and highlights an altered methylation at the GNAS locus in infertile boars.

Male infertility is an increasing health issue in today's society for both human and livestock populations. In livestock, male infertility slows the improvement of animal selection programs and agricultural productivity. There is increasing evidence that epigenetic marks play an important role in the production of good-quality sperm. We therefore screened for specific or common epigenetic signatures of livestock infertility. To do so, we compared DNA methylation level in sperm DNA from fertile and infertile boars. We evaluated first the global level of sperm DNA methylation and found no difference between the two groups of boars. We then selected 42 loci of interest, most of them known to be imprinted in human or mice, and assessed the imprinting status of five of them not previously described in swine tissues: WT1, CNTN3, IMPACT, QPCT, and GRB10. DNA methylation level was then quantified in fertile and infertile boars at these 42 loci. Results from fertile boars indicated that the methylation level of the selected loci is highly conserved between pig, human, and mice, with a few exceptions, including the POU5F1 (OCT4) promoter and RTL1. Comparison between fertile and infertile boars revealed that one imprinted region, the GNAS locus, shows an increase in sperm DNA methylation in three out of eight infertile boars with low semen quality. This increase in DNA methylation is associated with an altered expression of the genes belonging to the GNAS locus, suggesting a new role for GNAS in the proper formation of functional gametes.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues.

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

Astacin-like metallo-endopeptidase is dynamically expressed in embryonic stem cells and embryonic epithelium during morphogenesis.

BACKGROUND: Astacin-like metallo-proteases are zinc endopeptidases conserved among vertebrates and invertebrates. First described as hatching gland enzymes, many members of the family possess other functions during embryonic development. In the chick, however, functions of Astacin-like proteins remain elusive. RESULTS: We report here that Astacin-like (ASTL) is strongly expressed in mouse and chicken embryonic stem (ES) cells and exhibits a very dynamic expression pattern during embryogenesis and organogenesis, mostly in remodeled epithelia. Consistent with its expression in ES cells, chick ASTL is detected in vivo in the pluripotent cells of the epiblast and then disappears from the newly induced neural plate. ASTL expression remains at the junction of non-neural and neural ectoderm, just before neural tube closure. At later stages, chick ASTL is detected in the ventral epidermis before ventral closure, in the intermediate mesoderm, in the gonads and in the forming nephric duct and tubules of the mesonephros and metanephros. CONCLUSIONS: ASTL is dynamically expressed in the embryonic epithelium and in embryonic stem cells, suggesting an important function for the control of epithelial cell behavior during early development.

Water contamination by delorazepam induces epigenetic defects in the embryos of the clawed frog Xenopus laevis.

Delorazepam, a derivative of diazepam, is a psychotropic drug belonging to the benzodiazepine class. Used as a nervous-system inhibitor, it treats anxiety, insomnia, and epilepsy, but is also associated with misuse and abuse. Nowadays benzodiazepines are considered emerging pollutants: conventional wastewater treatment plants indeed are unable to eliminate these compounds. Consequently, they persist in the environment and bioaccumulate in non-target aquatic organisms with consequences still not fully clear. To collect more information, we investigated the possible epigenetic activity of delorazepam, at three concentrations (1, 5 and 10 µg/L) using Xenopus laevis embryos as a model. Analyses demonstrated a significant increase in genomic DNA methylation and differential methylation of the promoters of some early developmental genes (otx2, sox3, sox9, pax6, rax1, foxf1, and myod1). Moreover, studies on gene expression highlighted an unbalancing in apoptosis/proliferation pathways and an aberrant expression of DNA-repair genes. Results are alarming considering the growing trend of benzodiazepine concentrations in superficial waters, especially after the peak occurred as a consequence of the pandemic COVID-19, and the fact that benzodiazepine GABA-A receptors are highly conserved and present in all aquatic organisms.

Modeling Global Genomic Instability in Chronic Myeloid Leukemia (CML) Using Patient-Derived Induced Pluripotent Stem Cells (iPSCs).

METHODS: We used a patient-specific induced pluripotent stem cell (iPSC) line treated with the mutagenic agent N-ethyl-N-nitrosourea (ENU). Genomic instability was validated using γ-H2AX and micronuclei assays and CGH array for genomic events. RESULTS: An increased number of progenitors (x5-Fold), which proliferated in liquid cultures with a blast cell morphology, was observed in the mutagenized condition as compared to the unmutagenized one. CGH array performed for both conditions in two different time points reveals several cancer genes in the ENU-treated condition, some known to be altered in leukemia (BLM, IKZF1, NCOA2, ALK, EP300, ERG, MKL1, PHF6 and TET1). Transcriptome GEO-dataset GSE4170 allowed us to associate 125 of 249 of the aberrations that we detected in CML-iPSC with the CML progression genes already described during progression from chronic and AP to BC. Among these candidates, eleven of them have been described in CML and related to tyrosine kinase inhibitor resistance and genomic instability. CONCLUSIONS: These results demonstrated that we have generated, for the first time to our knowledge, an in vitro genetic instability model, reproducing genomic events described in patients with BC.

The endogenous retrovirus ENS-1 provides active binding sites for transcription factors in embryonic stem cells that specify extra embryonic tissue.

BACKGROUND: Long terminal repeats (LTR) from endogenous retroviruses (ERV) are source of binding sites for transcription factors which affect the host regulatory networks in different cell types, including pluripotent cells. The embryonic epiblast is made of pluripotent cells that are subjected to opposite transcriptional regulatory networks to give rise to distinct embryonic and extraembryonic lineages. To assess the transcriptional contribution of ERV to early developmental processes, we have characterized in vitro and in vivo the regulation of ENS-1, a host adopted and developmentally regulated ERV that is expressed in chick embryonic stem cells. RESULTS: We show that Ens-1 LTR activity is controlled by two transcriptional pathways that drive pluripotent cells to alternative developmental fates. Indeed, both Nanog that maintains pluripotency and Gata4 that induces differentiation toward extraembryonic endoderm independently activate the LTR. Ets coactivators are required to support Gata factors' activity thus preventing inappropriate activation before epigenetic silencing occurs during differentiation. Consistent with their expression patterns during chick embryonic development, Gata4, Nanog and Ets1 are recruited on the LTR in embryonic stem cells; in the epiblast the complementary expression of Nanog and Gata/Ets correlates with the Ens-1 gene expression pattern; and Ens-1 transcripts are also detected in the hypoblast, an extraembryonic tissue expressing Gata4 and Ets2, but not Nanog. Accordingly, over expression of Gata4 in embryos induces an ectopic expression of Ens-1. CONCLUSION: Our results show that Ens-1 LTR have co-opted conditions required for the emergence of extraembryonic tissues from pluripotent epiblasts cells. By providing pluripotent cells with intact binding sites for Gata, Nanog, or both, Ens-1 LTR may promote distinct transcriptional networks in embryonic stem cells subpopulations and prime the separation between embryonic and extraembryonic fates.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues.

BACKGROUND: In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle. RESULTS: More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region. CONCLUSIONS: The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

Role of circadian rhythm disorders on EMT and tumour-immune interactions in endocrine-related cancers.

The circadian rhythm is a major environmental regulator of plants and animal physiology. The alternation of days and nights is translated at the cell and tissue level thanks to a molecular machinery, called the circadian clock. This clock controls in particular numerous endocrine functions, and its imbalances can have serious consequences on homeostasis. This is particularly true for the development of endocrine-related cancers, like breast, ovarian and prostate cancer. Circadian rhythm disorder (CRD) not only affects key hormone levels (including oestrogen, melatonin, insulin, glucagon, cortisol) but also favours a pro-inflammatory and immunosuppressive phenotype in the tumour microenvironment. This particular aspect is conducive to epithelial-mesenchymal transition (EMT) of solid epithelial tumours and cancer cell dissemination. It also favours resistance to chemo- and immunotherapy. Here, we discuss the current knowledge on this crosstalk between CRD, EMT and the immune microenvironment in endocrine-related cancers and its consequences for the development of efficient therapies.