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.

Accounting for linkage disequilibrium in genome scans for selection without individual genotypes: The local score approach.

Detecting genomic footprints of selection is an important step in the understanding of evolution. Accounting for linkage disequilibrium in genome scans increases detection power, but haplotype-based methods require individual genotypes and are not applicable on pool-sequenced samples. We propose to take advantage of the local score approach to account for linkage disequilibrium in genome scans for selection, cumulating (possibly small) signals from single markers over a genomic segment, to clearly pinpoint a selection signal. Using computer simulations, we demonstrate that this approach detects selection with higher power than several state-of-the-art single-marker, windowing or haplotype-based approaches. We illustrate this on two benchmark data sets including individual genotypes, for which we obtain similar results with the local score and one haplotype-based approach. Finally, we apply the local score approach to Pool-Seq data obtained from a divergent selection experiment on behaviour in quail and obtain precise and biologically coherent selection signals: while competing methods fail to highlight any clear selection signature, our method detects several regions involving genes known to act on social responsiveness or autistic traits. Although we focus here on the detection of positive selection from multiple population data, the local score approach is general and can be applied to other genome scans for selection or other genomewide analyses such as GWAS.

Expressed alleles of imprinted IGF2, DLK1 and MEG3 colocalize in 3D-preserved nuclei of porcine fetal cells.

BACKGROUND: To explore the relationship between spatial genome organization and gene expression in the interphase nucleus, we used a genomic imprinting model, which offers parental-specific gene expression. Using 3D FISH in porcine fetal liver cells, we compared the nuclear organization of the two parental alleles (expressed or not) of insulin-like growth factor 2 (IGF2), a paternally imprinted gene located on chromosome 2. We investigated whether its nuclear positioning favors specific locus associations. We also tested whether IGF2 is implicated in long-range chromatin trans-associations as previously shown in the mouse model species for its reciprocal imprinted gene H19. RESULTS: We focused on the 3D position of IGF2 alleles, with respect to their individual chromosome 2 territories. The paternally expressed allele was tagged with nascent RNA. There were no significant differences in the position of the two alleles (p = 0.06). To determine long-range chromatin trans-interactions, we chose 12 genes, some of which are known to be imprinted in mammalian model species and belong to a network of imprinted genes (i.e. SLC38A4, DLK1, MEG3, and ZAC1). We screened them and ABCG2, OSBP2, OSBPL1, RPL32, NF1, ZAR1, SEP15, GPC3 for associations with IGF2 in liver cells. All imprinted genes tested showed an association with IGF2. The DLK1/MEG3 locus showed the highest rate of colocalization. This gene association was confirmed by 3D FISH (in 20 % of the nuclei analyzed), revealing also the close proximity of chromosomes 2 and 7 (in 60 % of nuclei). Furthermore, our observations showed that the expressed paternal IGF2 allele is involved in this association. This IGF2-(DLK1/MEG3) association also occurred in a high percentage of fetal muscle cells (36 % of nuclei). Finally, we showed that nascent IGF2, DLK1 and MEG3 RNAs can associate in pairs or in a three-way combination. CONCLUSION: Our results show that trans-associations occur between three imprinted genes IGF2, DLK1 and MEG3 both in fetal liver and muscle cells. All three expressed alleles associated in muscle cells. Our findings suggest that the 3D nuclear organization is linked to the transcriptional state of these genes.

3D organization of telomeres in porcine neutrophils and analysis of LPS-activation effect.

BACKGROUND: While the essential role of 3D nuclear architecture on nuclear functions has been demonstrated for various cell types, information available for neutrophils, essential components of the immune system, remains limited. In this study, we analysed the spatial arrangements of telomeres which play a central role in cell fate. Our studies were carried out in swine, which is an excellent model organism for both biomedical research and agronomic applications. We isolated bacterial artificial chromosome (BAC)-containing subtelomeric p and q sequences specific to each porcine chromosome. This allowed us to study the behaviour of p and q telomeres of homologous chromosomes for seven pairs chosen for their difference in length and morphology. This was performed using 3D-FISH on structurally preserved neutrophils, and confocal microscopy. Resting and lipopolysaccharide (LPS)-activated states were investigated to ascertain whether a response to a pathogen aggression modifies this organization. RESULTS: The positions of the p and q telomeres relative to the nuclear outer border were determined in the two states. All p telomeres changed their position significantly during the activation process, although the effect was less pronounced for the q telomeres. The patterns of telomeric associations between homologs and their frequencies were analysed for 7 pairs of chromosomes. This analysis revealed that the distribution of pp, qq and pq associations differs significantly among the 7 chromosomes. This distribution does not fit with the theoretical distribution for each chromosome, suggesting that preferential associations occur between subtelomeres. CONCLUSIONS: The percentage of nuclei harbouring at least one telomeric association between homologs varies significantly among the chromosomes, the smallest metacentric chromosome SSC12, which is also the richest in gene-density, harbouring the highest value. The distribution of types of telomeric associations is highly dependent on the chromosomes and is not affected by the activation process. The frequencies of telomeric associations are also highly dependent on the type of association and the type of chromosome. Overall, the LPS-activation process induces only minor changes in these patterns of associations. When telomeric associations occur, the associations of p and q arms from the same chromosome are the most frequent, suggesting that "chromosome bending" occurs in neutrophils as previously observed in gametes.

Transcriptomic and nuclear architecture of immune cells after LPS activation.

Changes in the nuclear positioning of specific genes, depending on their expression status, have been observed in a large diversity of physiological processes. However, gene position is poorly documented for immune cells which have been subjected to activation following bacterial infection. Using a pig model, we focused our study on monocyte-derived macrophages and neutrophils, as they are the first lines of defence against pathogens. We examined whether changes in gene expression due to LPS activation imply that genes have repositioned in the nuclear space. We first performed a transcriptomic analysis to identify the differentially expressed genes and then analysed the networks involved during lypopolysaccharide/interferon gamma activation in monocyte-derived macrophages. This allowed us to select four up-regulated (IL1ß, IL8, CXCL10 and TNFα) and four down-regulated (VIM, LGALS3, TUBA3 and IGF2) genes. Their expression statuses were verified by quantitative real-time RT-PCR before studying their behaviour in the nuclear space during macrophage activation by means of 3D fluorescence in situ hybridization. No global correlation was found between gene activity and radial positioning. Only TNFα belonging to the highly transcribed MHC region on chromosome 7 became more peripherally localized in relation to the less decondensed state of its chromosome territory (CT) in activated macrophages. The analysis of gene positioning towards their CT revealed that IL8 increases significantly its tendency to be outside its CT during the activation process. In addition, the gene to CT edge distances increase for the three up-regulated genes (IL8, CXCL10 and TNFα) among the four analysed. Contrarily, the four down-regulated genes did not change their position. The analysis of gene behaviour towards their CT was extended to include neutrophils for three (TNFα, IL8 and IL1ß) up- and two (IGF2 and TUBA3) down-regulated genes, and similar results were obtained. The analysis was completed by studying the four up-regulated genes in fibroblasts, not involved in immune response. Our data suggest that relocation in the nuclear space of genes that are differentially expressed in activated immune cells is gene and cell type specific but also closely linked to the entire up-regulation status of their chromosomal regions.

Major Reorganization of Chromosome Conformation During Muscle Development in Pig.

The spatial organization of the genome in the nucleus plays a crucial role in eukaryotic cell functions, yet little is known about chromatin structure variations during late fetal development in mammals. We performed in situ high-throughput chromosome conformation capture (Hi-C) sequencing of DNA from muscle samples of pig fetuses at two late stages of gestation. Comparative analysis of the resulting Hi-C interaction matrices between both groups showed widespread differences of different types. First, we discovered a complex landscape of stable and group-specific Topologically Associating Domains (TADs). Investigating the nuclear partition of the chromatin into transcriptionally active and inactive compartments, we observed a genome-wide fragmentation of these compartments between 90 and 110 days of gestation. Also, we identified and characterized the distribution of differential cis- and trans-pairwise interactions. In particular, trans-interactions at chromosome extremities revealed a mechanism of telomere clustering further confirmed by 3D Fluorescence in situ Hybridization (FISH). Altogether, we report major variations of the three-dimensional genome conformation during muscle development in pig, involving several levels of chromatin remodeling and structural regulation.

The 3D nuclear conformation of the major histocompatibility complex changes upon cell activation both in porcine and human macrophages.

BACKGROUND: The crucial role of the major histocompatibility complex (MHC) for the immune response to infectious diseases is well-known, but no information is available on the 3D nuclear organization of this gene-dense region in immune cells, whereas nuclear architecture is known to play an essential role on genome function regulation. We analyzed the spatial arrangement of the three MHC regions (class I, III and II) in macrophages using 3D-FISH. Since this complex presents major differences in humans and pigs with, notably, the presence of the centromere between class III and class II regions in pigs, the analysis was implemented in both species to determine the impact of this organization on the 3D conformation of the MHC. The expression level of the three genes selected to represent each MHC region was assessed by quantitative real-time PCR. Resting and lipopolysaccharide (LPS)-activated states were investigated to ascertain whether a response to a pathogen modifies their expression level and their 3D organization. RESULTS: While the three MHC regions occupy an intermediate radial position in porcine macrophages, the class I region was clearly more peripheral in humans. The BAC center-to-center distances allowed us to propose a 3D nuclear organization of the MHC in each species. LPS/IFNγ activation induces a significant decompaction of the chromatin between class I and class III regions in pigs and between class I and class II regions in humans. We detected a strong overexpression of TNFα (class III region) in both species. Moreover, a single nucleus analysis revealed that the two alleles can have either the same or a different compaction pattern. In addition, macrophage activation leads to an increase in alleles that present a decompacted pattern in humans and pigs. CONCLUSIONS: The data presented demonstrate that: (i) the MHC harbors a different 3D organization in humans and pigs; (ii) LPS/IFNγ activation induces chromatin decompaction, but it is not the same area affected in the two species. These findings were supported by the application of an original computation method based on the geometrical distribution of the three target genes. Finally, the position of the centromere inside the swine MHC could influence chromatin reorganization during the activation process.

Nuclear architecture of resting and LPS-stimulated porcine neutrophils by 3D FISH.

Neutrophils are essential components of the innate immune system due to their ability to kill and phagocytose invading microbes. They possess a lobulated nucleus and are capable of extensive and complex changes in response to bacterial stimulation. The aim of our study was to investigate whether the 3D nuclear organization of porcine neutrophils was modified upon stimulation. The organization of centromeres, telomeres, and chromosome territories (chromosomes 2, 3, 7, 8, 12, 13, and 17) was studied on structurally preserved nuclei using 3D fluorescence in situ hybridization, confocal microscopy, and image analysis. By differential labeling of centromeres of acrocentric and metacentric/submetacentric chromosomes, we showed that centromeres associated to form chromocenters but did so preferentially between chromosomes with the same morphology. Upon activation, some of these chromocenters dispersed. Telomeres were also found to form clusters, but their number remained unchanged in lipopolysaccharide-stimulated neutrophils. The analysis of the position of chromocenters and telomere clusters showed a more internal location of the latter compared to the former. The analysis of chromosome territories revealed that homologs were distributed randomly among lobes whatever the cell's status. The volume of these territories was not proportional to chromosome length, and some chromosomes (chr 3, 12, 13, and 17) were more prone to decondensation when neutrophils were stimulated. Thus, our study demonstrated that activation of neutrophils resulted in several modifications of their nuclear architecture: a decrease in the number of non-acrocentric chromocenters and the decondensation of several chromosomes.

Detecting local high-scoring segments: a first-stage approach for genome-wide association studies.

Genetic epidemiology aims at identifying biological mechanisms responsible for human diseases. Genome-wide association studies, made possible by recent improvements in genotyping technologies, are now promisingly investigated. In these studies, common first-stage strategies focus on marginal effects but lead to multiple-testing and are unable to capture the possibly complex interplay between genetic factors. We have adapted the use of the local score statistic, already successfully applied to analyse long molecular sequences. Via sum statistics, this method captures local and possible distant dependences between markers. Dedicated to genome-wide association studies, it is fast to compute, able to handle large datasets, circumvents the the multiple-testing problem and outlines a set of genomic regions (segments) for further analyses. Applied to simulated and real data, our approach outperforms classical Bonferroni and FDR corrections for multiple-testing. It is implemented in a software termed LHiSA for Local High-scoring Segments for Association and available at: http://stat.genopole.cnrs.fr/software/lhisa.

SIC: A tool to detect short inverted segments in a biological sequence.

The web software SIC provides a tool to search for short inverted segments (length 3-5000 bp) in a DNA sequence. The sequence is assumed to follow a Markov model. A statistic which is sensitive to inversion is presented. Searching inverted segments is done by a scanning approach after the user specifies the size of the scanning window and the order of the Markov chain. A list of the highest score segments is given with an assessment of the randomness of the result. SIC can be accessed via the URL: http://stat.genopole.cnrs.fr/SIC/.

Meiotic recombination analyses of individual chromosomes in male domestic pigs (Sus scrofa domestica).

For the first time in the domestic pig, meiotic recombination along the 18 porcine autosomes was directly studied by immunolocalization of MLH1 protein. In total, 7,848 synaptonemal complexes from 436 spermatocytes were analyzed, and 13,969 recombination sites were mapped. Individual chromosomes for 113 of the 436 cells (representing 2,034 synaptonemal complexes) were identified by immunostaining and fluorescence in situ hybridization (FISH). The average total length of autosomal synaptonemal complexes per cell was 190.3 µm, with 32.0 recombination sites (crossovers), on average, per cell. The number of crossovers and the lengths of the autosomal synaptonemal complexes showed significant intra- (i.e. between cells) and inter-individual variations. The distributions of recombination sites within each chromosomal category were similar: crossovers in metacentric and submetacentric chromosomes were concentrated in the telomeric regions of the p- and q-arms, whereas two hotspots were located near the centromere and in the telomeric region of acrocentrics. Lack of MLH1 foci was mainly observed in the smaller chromosomes, particularly chromosome 18 (SSC18) and the sex chromosomes. All autosomes displayed positive interference, with a large variability between the chromosomes.

The sheep genome illuminates biology of the rumen and lipid metabolism.

Sheep (Ovis aries) are a major source of meat, milk, and fiber in the form of wool and represent a distinct class of animals that have a specialized digestive organ, the rumen, that carries out the initial digestion of plant material. We have developed and analyzed a high-quality reference sheep genome and transcriptomes from 40 different tissues. We identified highly expressed genes encoding keratin cross-linking proteins associated with rumen evolution. We also identified genes involved in lipid metabolism that had been amplified and/or had altered tissue expression patterns. This may be in response to changes in the barrier lipids of the skin, an interaction between lipid metabolism and wool synthesis, and an increased role of volatile fatty acids in ruminants compared with nonruminant animals.

Pool-hmm: a Python program for estimating the allele frequency spectrum and detecting selective sweeps from next generation sequencing of pooled samples.

Due to its cost effectiveness, next generation sequencing of pools of individuals (Pool-Seq) is becoming a popular strategy for genome-wide estimation of allele frequencies in population samples. As the allele frequency spectrum provides information about past episodes of selection, Pool-seq is also a promising design for genomic scans for selection. However, no software tool has yet been developed for selection scans based on Pool-Seq data. We introduce Pool-hmm, a Python program for the estimation of allele frequencies and the detection of selective sweeps in a Pool-Seq sample. Pool-hmm includes several options that allow a flexible analysis of Pool-Seq data, and can be run in parallel on several processors. Source code and documentation for Pool-hmm is freely available at https://qgsp.jouy.inra.fr/.

seq++: analyzing biological sequences with a range of Markov-related models.

SUMMARY: The seq++ package offers a reference set of programs and an extensible library to biologists and developers working on sequence statistics. Its generality arises from the ability to handle sequences described with any alphabet (nucleotides, amino acids, codons and others). seq++ enables sequence modelling with various types of Markov models, including variable length Markov models and the newly developed parsimonious Markov models, all of them potentially phased. Simulation modules are supplied for Monte Carlo methods. Hence, this toolbox allows the study of any biological process which can be described by a series of states taken from a finite set.