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.

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.

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.

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.

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.

In Utero Exposure to Metformin Reduces the Fertility of Male Offspring in Adulthood.

Metformin is a drug used for the treatment of type 2 diabetes and disorders associated with insulin resistance. Metformin is also used in the treatment of pregnancy disorders such as gestational diabetes. However, the consequences of foetal exposure to metformin on the fertility of exposed offspring remain poorly documented. In this study, we investigated the effect of in utero metformin exposure on the fertility of female and male offspring. We observed that metformin is detectable in the blood of the mother and in amniotic fluid and blood of the umbilical cord. Metformin was not measurable in any tissues of the embryo, including the gonads. The effect of metformin exposure on offspring was sex specific. The adult females that had been exposed to metformin in utero presented no clear reduction in fertility. However, the adult males that had been exposed to metformin during foetal life exhibited a 30% reduction in litter size compared with controls. The lower fertility was not due to a change in sperm production or the motility of sperm. Rather, the phenotype was due to lower sperm head quality - significantly increased spermatozoa head abnormality with greater DNA damage - and hypermethylation of the genomic DNA in the spermatozoa associated with lower expression of the ten-eleven translocation methylcytosine dioxygenase 1 (TET1) protein. In conclusion, while foetal metformin exposure did not dramatically alter gonad development, these results suggest that metabolic modification by metformin during the foetal period could change the expression of epigenetic regulators such as Tet1 and perturb the genomic DNA in germ cells, changes that might contribute to a reduced fertility.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species.

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

Tissue Resources for the Functional Annotation of Animal Genomes.

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

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.

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.

Inducing Sequential Cycles of Epithelial-Mesenchymal and Mesenchymal-Epithelial Transitions in Mammary Epithelial Cells.

Epithelial-Mesenchymal Transition (EMT) and its reciprocal Mesenchymal-Epithelial Transition (MET) occur naturally as a cycling process during embryonic and foetal development. The capacity of such iterative cycles to drive cell fate and cellular and molecular behaviour in physiology or pathology remains unclear. We describe here a protocol to induce successive cycles of EMT/MET in an untransformed human mammary epithelial cell line (MCF10A) as well as the necessary controls for cycle validation.

An integrative atlas of chicken long non-coding genes and their annotations across 25 tissues.

Long non-coding RNAs (LNC) regulate numerous biological processes. In contrast to human, the identification of LNC in farm species, like chicken, is still lacunar. We propose a catalogue of 52,075 chicken genes enriched in LNC ( http://www.fragencode.org/ ), built from the Ensembl reference extended using novel LNC modelled here from 364 RNA-seq and LNC from four public databases. The Ensembl reference grew from 4,643 to 30,084 LNC, of which 59% and 41% with expression ≥ 0.5 and ≥ 1 TPM respectively. Characterization of these LNC relatively to the closest protein coding genes (PCG) revealed that 79% of LNC are in intergenic regions, as in other species. Expression analysis across 25 tissues revealed an enrichment of co-expressed LNC:PCG pairs, suggesting co-regulation and/or co-function. As expected LNC were more tissue-specific than PCG (25% vs. 10%). Similarly to human, 16% of chicken LNC hosted one or more miRNA. We highlighted a new chicken LNC, hosting miR155, conserved in human, highly expressed in immune tissues like miR155, and correlated with immunity-related PCG in both species. Among LNC:PCG pairs tissue-specific in the same tissue, we revealed an enrichment of divergent pairs with the PCG coding transcription factors, as for example LHX5, HXD3 and TBX4, in both human and chicken.

Chronic circadian disruption modulates breast cancer stemness and immune microenvironment to drive metastasis in mice.

Breast cancer is the most common type of cancer worldwide and one of the major causes of cancer death in women. Epidemiological studies have established a link between night-shift work and increased cancer risk, suggesting that circadian disruption may play a role in carcinogenesis. Here, we aim to shed light on the effect of chronic jetlag (JL) on mammary tumour development. To do this, we use a mouse model of spontaneous mammary tumourigenesis and subject it to chronic circadian disruption. We observe that circadian disruption significantly increases cancer-cell dissemination and lung metastasis. It also enhances the stemness and tumour-initiating potential of tumour cells and creates an immunosuppressive shift in the tumour microenvironment. Finally, our results suggest that the use of a CXCR2 inhibitor could correct the effect of JL on cancer-cell dissemination and metastasis. Altogether, our data provide a conceptual framework to better understand and manage the effects of chronic circadian disruption on breast cancer progression.

Single-Cell Transcriptome in Chronic Myeloid Leukemia: Pseudotime Analysis Reveals Evidence of Embryonic and Transitional Stem Cell States.

Recent experimental data suggest that the heterogeneity of chronic myeloid leukemia (CML) stem cells may be the result of the development of unique molecular events generating functional consequences in terms of the resistance and persistence of leukemic stem cells. To explore this phenomenon, we designed a single-cell transcriptome assay evaluating simultaneously the expression of 87 genes. Highly purified CD34+ cells from three CML patients at diagnosis were immobilized in microfluidic chips, and the expression of 87 genes was evaluated in each cell. This analysis identified a group of 13 highly connected genes including NANOG, POU5F1, LIN28A, and SOX2, representing on average 8.59% of the cell population analyzed. Bioinformatics analysis with the corrected matrix and t-distributed stochastic neighbor embedding (tSNE) algorithm identified four distinct clusters, and the pseudotime analysis confirmed the presence of seven stem cell states in the four clusters identified. ALOX5 expression was associated with the group of cells expressing the pluripotency markers. In in vitro analyses, two genes that were predicted to undergo similar regulation using pseudotime analysis (ALOX5 and FGFR) were found to be similarly inhibited by ponatinib, an FGFR inhibitor. Finally, in an independent cohort of CML patients, we found that pluripotency gene expression is a common feature of CD34+ CML cells at diagnosis. Overall, these experiments allowed identification of individual CD34+ cells expressing high levels of pluripotency genes at diagnosis, in which a continuum of transitional states were identified using pseudotime analysis. These results suggest that leukemic stem cell persistence in CML needs to be targeted simultaneously rather than using a single pathway.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens.

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

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.