Profiling the landscape of transcription, chromatin accessibility and chromosome conformation of cattle, pig, chicken and goat genomes [FAANG pilot project]

Functional annotation of livestock genomes is a critical and obvious next step to derive maximum benefit for agriculture, animal science, animal welfare and human health. The aim of the Fr-AgENCODE project is to generate multi-species functional genome annotations by applying high-throughput molecular assays on three target tissues/cells relevant to the study of immune and metabolic traits. An extensive collection of stored samples from other tissues is available for further use (FAANG Biosamples ‘FR-AGENCODE’). From each of two males and two females per species (pig, cattle, goat, chicken), strand-oriented RNA-seq and chromatin accessibility ATAC-seq assays were performed on liver tissue and on two T-cell types (CD3+CD4+&CD3+CD8+) sorted from blood (mammals) or spleen (chicken). Chromosome Conformation Capture (in situ Hi-C) was also carried out on liver. Sequencing reads from the 3 assays were processed using standard processing pipelines. While most (50–70%) RNA-seq reads mapped to annotated exons, thousands of novel transcripts and genes were found, including extensions of annotated protein-coding genes and new lncRNAs (see abstract #69857). Consistency of ATAC-seq results was confirmed by the significant proportion of called peaks in promoter regions (36–66%) and by the specific accumulation pattern of peaks around gene starts (TSS) v. gene ends (TTS). Principal Component Analyses for RNA-seq (based on quantified gene expression) and ATAC-seq (based on quantified chromatin accessibility) highlighted clusters characterised by cell type and sex in all species. From Hi-C data, we generated 40kb-resolution interaction maps, profiled a genome-wide Directionality Index and identified from 4,100 (chicken) to 12,100 (pig) topologically-associating do- mains (TADs). Correlations were reported between RNA-seq and ATAC-seq results (see abstract #71581). In summary, we present here an overview of the first multi-species and -tissue annotations of chromatin accessibility and genome architecture related to gene expression for farm animals

RH mapping by sequencing: chromosome-scale assembly of the duck genome

Like many other species, the duck genome has been sequenced thanks to the technological breakthrough provided by the emergence of Next Generation Sequencing (NGS). The resulting de novo assemblies are however made of thousands of scattered scaffolds. To achieve chromosome-scale contiguity, long-range intermediate genome maps remain indispensable. Radiation Hybrid (RH) maps have been used to assist the generation of chromosome-scale genome assemblies by taking advantage of the high density SNP chips that provide a large number of markers that can be efficiently genotyped on the panel. In the absence of such a resource in duck, we sequenced 100 hybrid clones of a duck RH panel enabling direct genotyping of the assembly scaffolds on the panel. The rationale is to use scaffolds as markers and to genotype the scaffolds by sequencing the clones: the presence/absence of a scaffold in a particular sequenced hybrid is attested by the presence/absence of reads mapping specifically to this scaffold. The detection of scaffolds exhibiting a chromosomal breakage resulting from the irradiation process revealed itself to be a critical issue of this genotyping by sequencing process. This process resulted in the construction of RH vectors for 2,027 scaffolds, representing a total of about 1 Gb of sequences (95% of the current Duck genome assembly). The subsequent linkage analysis enabled the construction of RH maps and therefore to organize, i.e. order and orient, the scaffolds into pseudomolecules associated to the corresponding duck chromosomes. We describe here the whole mapping process, from sequence-based genotyping to the construction of comparative maps, as well as few examples of intra-chromosomal rearrangements that have been identified by the comparison with the chicken, turkey and zebra finch genomes and subsequently confirmed by FISH. We describe a method to order and orient sequence scaffolds into super-scaffolds spanning entire chromosomes. The method, which requires a pre-existing RH panel and sequence scaffolds from an NGS assembly, relies on a shallow sequencing of the RH clones. This approach was applied to the duck genome and produced chromosome-scale scaffolds for 29 out of the 41 duck chromosomes

Identification of a t(3;4)(p1.3;q1.5) translocation breakpoint in pigs using somatic cell hybrid mapping and high-resolution mate-pair sequencing

<div><p>Reciprocal translocations are the most frequently occurring constitutional structural rearrangements in mammalian genomes. In phenotypically normal pigs, an incidence of 1/200 is estimated for such rearrangements. Even if constitutional translocations do not necessarily induce defects and diseases, they are responsible for significant economic losses in domestic animals due to reproduction failures. Over the last 30 years, advances in molecular and cytogenetic technologies have led to major improvements in the resolution of the characterization of translocation events. Characterization of translocation breakpoints helps to decipher the mechanisms that lead to such rearrangements and the functions of the genes that are involved in the translocation. Here, we describe the fine characterization of a reciprocal translocation t(3;4) (p1.3;q1.5) detected in a pig line. The breakpoint was identified at the base-pair level using a positional cloning and chromosome walking strategy in somatic cell hybrids that were generated from an animal that carries this translocation. We show that this translocation occurs within the <i>ADAMTSL4</i> gene and results in a loss of expression in homozygous carriers. In addition, by taking this translocation as a model, we used a whole-genome next-generation mate-pair sequencing approach on pooled individuals to evaluate this strategy for high-throughput screening of structural rearrangements.</p></div

Fluorescent in situ hybridization of the pIRS-PCR probe, derived from hybrid Hb1.8, on normal pig metaphase chromosomes.

<p>pIRS-PCR products label the entire length of chromosomes 1, 8, 10, 16, X, and partially chromosomes 3, 4 and 5. Chromosomes 3 and 4 are indicated with yellow and blue arrows respectively (a), and are magnified in (b). On these two chromosomes, signals are observed on the p arm for SSC3 and in the p15-q22 region for SSC4, confirming the presence of the translocated derivative chromosome der(4) in Hb1.8 hybrid clone.</p

Distribution of discordant mate-pairs linking SSC3 to SSC4.

<p>Each circle represents a single mate-pair with one read mapped to SSC3 (x-axis) and the other read mapped to SSC4 (y-axis). In total, more than 12,000 mate-pairs link the two chromosomes, underlining the high background noise due to false-positive discordant mate-pairs. The dark green circle, where the two dashed lines cross, highlights a concentration of mate-pairs linking the two regions involved in the translocation.</p

Results obtained by PCR-screening of each hybrid clone with microsatellite markers.

<p>Results obtained by PCR-screening of each hybrid clone with microsatellite markers.</p

Karyotype of the balanced reciprocal translocation t(3;4)(p1.3;q1.5).

<p>(a): G-band karyotype obtained from a LW female, heterozygous for the translocation. The breakpoint positions on SSC3 and SSC4 are indicated with arrows. (b): Schematic representation of SSC3 and SSC4 pairs including one normal SSC3 or SSC4 chromosome (N) and one translocated chromosome der(3) or der(4).</p