A first complete catalog of highly expressed genes in eight chicken tissues reveals uncharacterized gene families specific for the chicken testis.

Based on next-generation sequencing, we established a repertoire of differentially overexpressed genes (DoEG) in eight adult chicken tissues: testis, brain, lung, liver, kidney, muscle, heart, and intestine. With 4499 DoEG, the testis had the highest number and proportion of DoEG compared with the seven somatic tissues. The testis DoEG set included the highest proportion of long noncoding RNAs (lncRNAs; 1851, representing 32% of the lncRNA genes in the whole genome) and the highest proportion of protein-coding genes (2648, representing 14.7% of the protein-coding genes in the whole genome). The main significantly enriched Gene Ontology terms related to the protein-coding genes are "reproductive process," "tubulin binding," "microtubule cytoskeleton." By using real-time quantitative reverse transcription polymerase chain reaction, we confirmed the overexpression of genes that encode proteins already described in chicken sperm (such as CABYR, SPAT18 and CDK5RAP2) but whose testis origin had not been confirmed previously. Moreover, we demonstrated the overexpression of vertebrate orthologs of testis genes not yet described in the adult chicken testis (such as NEK2, AK7 and CCNE2). Using clustering according to primary sequence homology, we found that 67% (1737) of the 2648 testis protein-coding genes were unique genes. This proportion was significantly higher than the somatic tissues except muscle. We clustered the other 911 testis protein-coding genes into 495 families, from which 47 had all paralogs overexpressed in the testis. Among these 47 testis-specific families, eight contain uncharacterized members without orthologs in other metazoans except birds: these families are then specific for chickens (probably more broadly birds).

CoNIC Challenge: Pushing the frontiers of nuclear detection, segmentation, classification and counting.

Nuclear detection, segmentation and morphometric profiling are essential in helping us further understand the relationship between histology and patient outcome. To drive innovation in this area, we setup a community-wide challenge using the largest available dataset of its kind to assess nuclear segmentation and cellular composition. Our challenge, named CoNIC, stimulated the development of reproducible algorithms for cellular recognition with real-time result inspection on public leaderboards. We conducted an extensive post-challenge analysis based on the top-performing models using 1,658 whole-slide images of colon tissue. With around 700 million detected nuclei per model, associated features were used for dysplasia grading and survival analysis, where we demonstrated that the challenge's improvement over the previous state-of-the-art led to significant boosts in downstream performance. Our findings also suggest that eosinophils and neutrophils play an important role in the tumour microevironment. We release challenge models and WSI-level results to foster the development of further methods for biomarker discovery.

Combined Multiplexed Phage Display, High-Throughput Sequencing, and Functional Assays as a Platform for Identifying Modulatory VHHs Targeting the FSHR.

Developing modulatory antibodies against G protein-coupled receptors is challenging. In this study, we targeted the follicle-stimulating hormone receptor (FSHR), a significant regulator of reproduction, with variable domains of heavy chain-only antibodies (VHHs). We built two immune VHH libraries and submitted them to multiplexed phage display approaches. We used next-generation sequencing to identify 34 clusters of specifically enriched sequences that were functionally assessed in a primary screen based on a cAMP response element (CRE)-dependent reporter gene assay. In this assay, 23 VHHs displayed negative or positive modulation of FSH-induced responses, suggesting a high success rate of the multiplexed strategy. We then focused on the largest cluster identified (i.e., PRC1) that displayed positive modulation of FSH action. We demonstrated that PRC1 specifically binds to the human FSHR and human FSHR/FSH complex while potentiating FSH-induced cAMP production and Gs recruitment. We conclude that the improved selection strategy reported here is effective for rapidly identifying functionally active VHHs and could be adapted to target other challenging membrane receptors. This study also led to the identification of PRC1, the first potential positive modulator VHH reported for the human FSHR.

Genes encoding teleost orthologues of human signal transduction proteins remain duplicated or triplicated more frequently than the whole genome.

Cell signalling involves a myriad of proteins, many of which belong to families of related proteins, and these proteins display a huge number of interactions. One of the events that has led to the creation of new genes is whole genome duplication (WGD), a phenomenon that has made some major innovations possible. In addition to the two WGDs that happened before gnathostome radiation, teleost genomes underwent one (the 3WGD group) or two (the 4WGD group) extra WGD after separation from the lineage leading to holostei. In the present work, we studied in 63 teleost species whether the orthologues of human genes involved in 47 signalling pathways (HGSP) remain more frequently duplicated, triplicated or in the singleton state compared with the whole genome. We found that these genes have remained duplicated and triplicated more frequently in teleost of the 3WGD and 4WGD groups, respectively. Moreover, by examining pairs of interacting gene products in terms of conserved copy numbers, we found a majority of the 1:1 and 1:2 proportions in the 3WGD group (between 54% and 60%) and of the 2:2 and 2:4 proportions in the 4WGD group (30%). In both groups, we observed the 0:n proportion at a mean of approximately 10%, and we found some pseudogenes in the concerned genomes. Finally, the proportions were very different between the studied pathways. The n:n (i.e. same) proportion concerned 20%-65% of the interactions, depending on the pathways, and the n:m (i.e. different) proportion concerned 34%-70% of the interactions. Among the n:n proportion, the 1:1 ratio is most represented (25.8%) and among the n:m ratios, the 1:2 is most represented (25.0%). We noted the absence of gene loss for the JAK-STAT, FoxO and glucagon pathways. Overall, these results show that the teleost gene orthologues of HGSP remain duplicated (3WGD) and triplicated (4WGD) more frequently than the whole genome, although some genes have been lost, and the proportions have not always been maintained.

Seasonal remodeling of the progenitor pool and its distribution in the ewe mediobasal hypothalamus.

Recent studies have reported the presence of adult neurogenesis in the arcuate nucleus periventricular space (pvARH) and in the median eminence (ME), two structures involved in reproductive function. In sheep, a seasonal mammal, decreasing daylight in autumn induces a higher neurogenic activity in these two structures. However, the different types of neural stem and progenitor cells (NSCs/NPCs) that populate the arcuate nucleus and median eminence, as well as their location, have not been evaluated. Here, using semi-automatic image analyzing processes, we identified and quantified the different populations of NSCs/NPCs, showing that, during short days, higher densities of [SOX2 +] cells are found in pvARH and ME. In the pvARH, higher densities of astrocytic and oligodendrocitic progenitors mainly contribute to these variations. The different populations of NSCs/NPCs were mapped according to their position relative to the third ventricle and their proximity to the vasculature. We showed that [SOX2 +] cells extended deeper into the hypothalamic parenchyma during short days. Similarly, [SOX2 +] cells were found further from the vasculature in the pvARH and the ME, at this time of year, indicating the existence of migratory signals. The expression levels of neuregulin transcripts (NRGs), whose proteins are known to stimulate proliferation and adult neurogenesis and to regulate progenitor migration, as well as the expression levels of ERBB mRNAs, cognate receptors for NRGs, were assessed. We showed that mRNA expression changed seasonally in pvARH and ME, suggesting that the ErbB-NRG system is potentially involved in the photoperiodic regulation of neurogenesis in seasonal adult mammals.

Expanding duplication of the testis PHD Finger Protein 7 (PHF7) gene in the chicken genome.

Gene duplications increase genetic and phenotypic diversity and occur in complex genomic regions that are still difficult to sequence and assemble. PHD Finger Protein 7 (PHF7) acts during spermiogenesis for histone-to-histone protamine exchange and is a determinant of male fertility in Drosophila and the mouse. We aimed to explore and characterise in the chicken genome the expanding family of the numerous orthologues of the unique mouse Phf7 gene (highly expressed in the testis), observing the fact that this information is unclear and/or variable according to the versions of databases. We validated nine primer pairs by in silico PCR for their use in screening the chicken bacterial artificial chromosome (BAC) library to produce BAC-derived probes to detect and localise PHF7-like loci by fluorescence in situ hybridisation (FISH). We selected nine BAC that highlighted nine chromosomal regions for a total of 10 distinct PHF7-like loci on five Gallus gallus chromosomes: Chr1 (three loci), Chr2 (two loci), Chr12 (one locus), Chr19 (one locus) and ChrZ (three loci). We sequenced the corresponding BAC by using high-performance PacBio technology. After assembly, we performed annotation with the FGENESH program: there were a total of 116 peptides, including 39 PHF7-like proteins identified by BLASTP. These proteins share a common exon-intron core structure of 8-11 exons. Phylogeny revealed that the duplications occurred first between chromosomal regions and then inside each region. There are other duplicated genes in the identified BAC sequences, suggesting that these genomic regions exhibit a high rate of tandem duplication. We showed that the PHF7 gene, which is highly expressed in the rooster testis, is a highly duplicated gene family in the chicken genome, and this phenomenon probably concerns other bird species.

Seasonal vascular plasticity in the mediobasal hypothalamus of the adult ewe.

Sheep, like most seasonal mammals, exhibit a cyclic adaptive reproductive physiology that allows ewes to give birth to their progeny during the spring when environmental conditions are favorable to their survival. This process relies on the detection of day length (or photoperiod) and is associated with profound changes in cellular plasticity and gene expression in the hypothalamic-pituitary-gonadal axis, mechanisms that are suggested to participate in the seasonal adaptation of neuroendocrine circuits. Recently, pituitary vascular growth has been proposed as a seasonally regulated process in which the vascular endothelial growth factor A (VEGFA), a well-known angiogenic cytokine, is suspected to play a crucial role. However, whether this mechanism is restricted to the pituitary gland or also occurs in the mediobasal hypothalamus (MBH), a crucial contributor to the control of the reproductive function, remains unexplored. Using newly developed image analysis tools, we showed that the arcuate nucleus (ARH) of the MBH exhibits an enhanced vascular density during the long photoperiod or non-breeding season, associated with higher expression of VEGFA. In the median eminence (ME), a structure connecting the MBH to the pituitary gland, higher VEGFA, kinase insert domain receptor (KDR/VEGFR2) and plasmalemma vesicle-associated protein (PLVAP) gene expressions were detected during the long photoperiod. We also found that VEGFA and its receptor, VEGFR2, are expressed by neurons and tanycytes in both the ARH and ME. Altogether, these data show variations in the MBH vasculature according to seasons potentially through a VEGFA-dependent pathway, paving the way for future studies aiming to decipher the role of these changes in the hypothalamic control of seasonal reproduction.

Genes Encoding Teleost Orthologs of Human Haploinsufficient and Monoallelically Expressed Genes Remain in Duplicate More Frequently Than the Whole Genome.

Gene dosage is an important issue both in cell and evolutionary biology. Most genes are present in two copies or alleles in diploid eukariotic cells. The most outstanding exception is monoallelic gene expression (MA) that concerns genes localized on the X chromosome or in regions undergoing parental imprinting in eutherians, and many other genes scattered throughout the genome. In diploids, haploinsufficiency (HI) implies that a single functional copy of a gene in a diploid organism is insufficient to ensure a normal biological function. One of the most important mechanisms ensuring functional innovation during evolution is whole genome duplication (WGD). In addition to the two WGDs that have occurred in vertebrate genomes, the teleost genomes underwent an additional WGD, after their divergence from tetrapods. In the present work, we have studied on 57 teleost species whether the orthologs of human MA or HI genes remain more frequently in duplicates or returned more frequently in singleton than the rest of the genome. Our results show that the teleost orthologs of HI human genes remained more frequently in duplicate than the rest of the genome in all of the teleost species studied. No signal was observed for the orthologs of genes mapping to the human X chromosome or subjected to parental imprinting. Surprisingly, the teleost orthologs of the other human MA genes remained in duplicate more frequently than the rest of the genome for most teleost species. These results suggest that the teleost orthologs of MA and HI human genes also undergo selective pressures either related to absolute protein amounts and/or of dosage balance issues. However, these constraints seem to be different for MA genes in teleost in comparison with human genomes.

The piggyBac-derived protein 5 (PGBD5) transposes both the closely and the distantly related piggyBac-like elements Tcr-pble and Ifp2.

The vertebrate piggyBac derived transposase 5 (PGBD5) encodes a domesticated transposase, which is active and able to transpose its distantly related piggyBac-like element (pble), Ifp2. This raised the question whether PGBD5 would be more effective at mobilizing a phylogenetically closely related pble element. We aimed to identify the pble most closely related to the pgbd5 gene. We updated the landscape of vertebrate pgbd genes to develop efficient filters and identify the most closely related pble to each of these genes. We found that Tcr-pble is phylogenetically the closest pble to the pgbd5 gene. Furthermore, we evaluated the capacity of two murine and human PGBD5 isoforms, Mm523 and Hs524, to transpose both Tcr-pble and Ifp2 elements. We found that both pbles could be transposed by Mm523 with similar efficiency. However, integrations of both pbles occurred through both proper transposition and improper PGBD5-dependent recombination. This suggested that the ability of PGBD5 to bind both pbles may not be based on the primary sequence of element ends, but may involve recognition of inner DNA motifs, possibly related to palindromic repeats. In agreement with this hypothesis, we identified internal palindromic repeats near the end of 24 pble sequences, which display distinct sequences.

Variations in genome size between wild and domesticated lineages of fowls belonging to the Gallus gallus species.

Efforts to elucidate the causes of biological differences between wild fowls and their domesticated relatives, the chicken, have to date mainly focused on the identification of single nucleotide mutations. Other types of genomic variations have however been demonstrated to be important in avian evolution and associated to variations in phenotype. They include several types of sequences duplicated in tandem that can vary in their repetition number. Here we report on genome size differences between the red jungle fowl and several domestic chicken breeds and selected lines. Sequences duplicated in tandem such as rDNA, telomere repeats, satellite DNA and segmental duplications were found to have been significantly re-shaped during domestication and subsequently by human-mediated selection. We discuss the extent to which changes in genome organization that occurred during domestication agree with the hypothesis that domesticated animal genomes have been shaped by evolutionary forces aiming to adapt them to anthropized environments.

Thermal Manipulation During Embryogenesis Impacts H3K4me3 and H3K27me3 Histone Marks in Chicken Hypothalamus.

Changes in gene activity through epigenetic alterations induced by early environmental challenges during embryogenesis are known to impact the phenotype, health, and disease risk of animals. Learning how environmental cues translate into persisting epigenetic memory may open new doors to improve robustness and resilience of developing animals. It has previously been shown that the heat tolerance of male broiler chickens was improved by cyclically elevating egg incubation temperature. The embryonic thermal manipulation enhanced gene expression response in muscle (P. major) when animals were heat challenged at slaughter age, 35 days post-hatch. However, the molecular mechanisms underlying this phenomenon remain unknown. Here, we investigated the genome-wide distribution, in hypothalamus and muscle tissues, of two histone post-translational modifications, H3K4me3 and H3K27me3, known to contribute to environmental memory in eukaryotes. We found 785 H3K4me3 and 148 H3K27me3 differential peaks in the hypothalamus, encompassing genes involved in neurodevelopmental, metabolic, and gene regulation functions. Interestingly, few differences were identified in the muscle tissue for which differential gene expression was previously described. These results demonstrate that the response to embryonic thermal manipulation (TM) in chicken is mediated, at least in part, by epigenetic changes in the hypothalamus that may contribute to the later-life thermal acclimation.

Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies.

BACKGROUND: More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced. RESULTS: The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected. CONCLUSIONS: High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.

Genes Encoding Teleost Fish Ligands and Associated Receptors Remained in Duplicate More Frequently than the Rest of the Genome.

Signaling through ligand/receptor interactions is a widespread mechanism across all living taxa. During evolution, however, there has been a diversification in multigene families and changes in their interaction patterns. Among the events that led to the creation of new genes is the whole-genome duplication, which made possible some major innovations. Teleost fishes descended from a common ancestor which underwent one such whole-genome duplication. In our study, we investigated the effect of complete genome duplication on the evolution of ligand-receptor pairs in teleosts. We selected ten teleost species and used bioinformatics programs and phylogenetic tools in order to study the evolution of the human ligands and receptors that have orthologous genes in fishes, as well as the rest of the fish genomes. We established that since the complete duplication of the fish genomes, the conservation in duplicate copy of ligand and receptor genes is higher than expected. However, the ligand/receptor pair partners did not necessarily evolve in the same way, and a lot of situations occurred in which one of the partners returned in singleton copy when the other one was maintained in duplicate. This suggests that changes in interaction partners may have taken place during the evolution of teleosts. Moreover, the fate of the ligands and receptor coding genes is partly congruent with the phylogeny of teleosts. However, some incongruences can be observed. We suggest that these incongruences are correlated to the environment.

Perturbation of mRNP biogenesis reveals a dynamic landscape of the Rrp6-dependent surveillance machinery trafficking along the yeast genome.

Eukaryotic cells have evolved a nuclear quality control (QC) system to monitor the co-transcriptional mRNA processing and packaging reactions that lead to the formation of export-competent ribonucleoprotein particles (mRNPs). Aberrant mRNPs that fail to pass the QC steps are retained in the nucleus and eliminated by the exonuclease activity of Rrp6. It is still unclear how the surveillance system is precisely coordinated both physically and functionally with the transcription machinery to detect the faulty events that may arise at each step of transcript elongation and mRNP formation. To dissect the QC mechanism, we previously implemented a powerful assay based on global perturbation of mRNP biogenesis in yeast by the bacterial Rho helicase. By monitoring model genes, we have shown that the QC process is coordinated by Nrd1, a component of the NNS complex (Nrd1-Nab3-Sen1) involved in termination, processing and decay of ncRNAs which is recruited by the CTD of RNAP II. Here, we have extended our investigations by analyzing the QC behaviour over the whole yeast genome. We performed high-throughput RNA sequencing (RNA-seq) to survey a large collection of mRNPs whose biogenesis is affected by Rho action and which can be rescued upon Rrp6 depletion. This genome-wide perspective was extended by generating high-resolution binding landscapes (ChIP-seq) of QC components along the yeast chromosomes before and after perturbation of mRNP biogenesis. Our results show that perturbation of mRNP biogenesis redistributes the QC components over the genome with a significant hijacking of Nrd1 and Nab3 from genomic loci producing ncRNAs to Rho-affected protein-coding genes, triggering termination and processing defects of ncRNAs.

Relative effects of location relative to the corpus luteum and lactation on the transcriptome of the bovine oviduct epithelium.

BACKGROUND: Lactation and associated metabolic stresses during the post-partum period have been shown to impair fertility in dairy cows. The oviduct plays key roles in embryo development and the establishment of pregnancy in cattle. The aim of this study was to investigate the effects of lactation and location relative to the corpus luteum (CL) on the transcriptome of the bovine oviduct epithelium. RESULTS: An original animal model was used. At 60 days post-partum, Holstein lactating (n = 4) and non-lactating (i.e. never milked after calving; n = 5) cows, as well as control nulliparous heifers (n = 5), were slaughtered on Day 3 following induced estrus, and epithelial samples from the oviductal ampulla and isthmus ipsilateral and contralateral to the corpus luteum (CL) were recovered for RNA sequencing. In the oviduct ipsilateral to the CL, differentially expressed genes (DEGs) were identified between heifers compared with both postpartum cow groups. However, only 15 DEGs were identified between post-partum lactating and non-lactating cows in the ipsilateral isthmus and none were identified in the ipsilateral ampulla. In contrast, 192 and 2583 DEGs were identified between ipsilateral and contralateral ampulla and isthmus, respectively. In both regions, more DEGs were identified between ipsilateral and contralateral oviducts in non-lactating cows and heifers than in lactating cows. Functional annotation of the DEGs associated with comparisons between metabolic groups highlighted a number of over-represented biological functions and cell pathways including immune response and cholesterol/steroid biosynthesis. CONCLUSIONS: Gene expression in the oviduct epithelium, particularly in the isthmus, was more affected by the location relative to the CL than by lactation at Day 3 post-estrus. Furthermore, the effect of the proximity to the CL was modulated by the metabolic status of the cow.

But where did the centromeres go in the chicken genome models?

The chicken genome was the third vertebrate to be sequenced. To date, its sequence and feature annotations are used as the reference for avian models in genome sequencing projects developed on birds and other Sauropsida species, and in genetic studies of domesticated birds of economic and evolutionary biology interest. Therefore, an accurate description of this genome model is important to a wide number of scientists. Here, we review the location and features of a very basic element, the centromeres of chromosomes in the galGal5 genome model. Centromeres are elements that are not determined by their DNA sequence but by their epigenetic status, in particular by the accumulation of the histone-like protein CENP-A. Comparison of data from several public sources (primarily marker probes flanking centromeres using fluorescent in situ hybridization done on giant lampbrush chromosomes and CENP-A ChIP-seq datasets) with galGal5 annotations revealed that centromeres are likely inappropriately mapped in 9 of the 16 galGal5 chromosome models in which they are described. Analysis of karyology data confirmed that the location of the main CENP-A peaks in chromosomes is the best means of locating the centromeres in 25 galGal5 chromosome models, the majority of which (16) are fully sequenced and assembled. This data re-analysis reaffirms that several sources of information should be examined to produce accurate genome annotations, particularly for basic structures such as centromeres that are epigenetically determined.

Anti-angiogenic VEGFAxxxb transcripts are not expressed in the medio-basal hypothalamus of the seasonal sheep.

This study investigated Vegfa expression in the pars tuberalis (PT) of the pituitary and medio-basal hypothalamus (MBH) of sheep, across seasons and reproductive states. It has recently been proposed that season impacts alternative splicing of Vegfa mRNA in the PT, which shifts the balance between angiogenic VEGFAxxx and anti-angiogenic VEGFAxxxb isoforms (with xxx the number of amino acids of the mature VEGFA proteins) to modulate seasonal breeding. Here, we used various RT-PCR methodologies and analysis of RNAseq datasets to investigate seasonal variation in expression and splicing of the ovine Vegfa gene. Collectively, we identify 5 different transcripts for Vegfa within the ewe PT/MBH, which correspond to splicing events previously described in mouse and human. All identified transcripts encode angiogenic VEGFAxxx isoforms, with no evidence for alternative splicing within exon 8. These findings led us to investigate in detail how "Vegfaxxxb-like" PCR products could be generated by RT-PCR and misidentified as endogenous transcripts, in sheep and human HEK293 cells. In conclusion, our findings do not support the existence of anti-angiogenic VEGFAxxxb isoforms in the ovine PT/MBH and shed new light on the interpretation of prior studies, which claimed to identify Vegfaxxxb isoforms by RT-PCR.

An Assessment of Fixed and Native Chromatin Preparation Methods to Study Histone Post-Translational Modifications at a Whole Genome Scale in Skeletal Muscle Tissue.

BACKGROUND: Genomic loci associated with histone marks are typically analyzed by immunoprecipitation of the chromatin followed by quantitative-PCR (ChIP-qPCR) or high throughput sequencing (ChIP-seq). Chromatin can be either cross-linked (X-ChIP) or used in the native state (N-ChIP). Cross-linking of DNA and proteins helps stabilizing their interactions before analysis. Despite X-ChIP is the most commonly used method, muscle tissue fixation is known to be relatively inefficient. Moreover, no protocol described a simple and reliable preparation of skeletal muscle chromatin of sufficient quality for subsequent high-throughput sequencing. Here we aimed to set-up and compare both chromatin preparation methods for a genome-wide analysis of H3K27me3, a broad-peak histone mark, using chicken P. major muscle tissue. RESULTS: Fixed and unfixed chromatin were prepared from chicken muscle tissues (Pectoralis major). Chromatin fixation, shearing by sonication or digestion and immunoprecipitation performed equivalently. High-quality Illumina reads were obtained (q30 > 93%). The bioinformatic analysis of the data was performed using epic, a tool based on SICER, and MACS2. Forty millions of reads were analyzed for both X-ChIP-seq and N-ChIP-seq experiments. Surprisingly, H3K27me3 X-ChIP-seq analysis led to the identification of only 2000 enriched regions compared to about 15,000 regions identified in the case of N-ChIP-seq. N-ChIP-seq peaks were more consistent between replicates compared to X-ChIP-seq. Higher N-ChIP-seq enrichments were confirmed by ChIP-qPCR at the PAX5 and SOX2 loci known to be enriched for H3K27me3 in myotubes and at the loci of common regions of enrichment identified in this study. CONCLUSIONS: Our findings suggest that the preparation of muscle chromatin for ChIP-seq in cross-linked conditions can compromise the systematic analysis of broad histone marks. Therefore, native chromatin preparation should be preferred to cross-linking when a ChIP experiment has to be performed on skeletal muscle tissue, particularly when a broad source signal is considered.

Deep landscape update of dispersed and tandem repeats in the genome model of the red jungle fowl, Gallus gallus, using a series of de novo investigating tools.

BACKGROUND: The program RepeatMasker and the database Repbase-ISB are part of the most widely used strategy for annotating repeats in animal genomes. They have been used to show that avian genomes have a lower repeat content (8-12 %) than the sequenced genomes of many vertebrate species (30-55 %). However, the efficiency of such a library-based strategies is dependent on the quality and completeness of the sequences in the database that is used. An alternative to these library based methods are methods that identify repeats de novo. These alternative methods have existed for a least a decade and may be more powerful than the library based methods. We have used an annotation strategy involving several complementary de novo tools to determine the repeat content of the model genome galGal4 (1.04 Gbp), including identifying simple sequence repeats (SSRs), tandem repeats and transposable elements (TEs). RESULTS: We annotated over one Gbp. of the galGal4 genome and showed that it is composed of approximately 19 % SSRs and TEs repeats. Furthermore, we estimate that the actual genome of the red jungle fowl contains about 31-35 % repeats. We find that library-based methods tend to overestimate TE diversity. These results have a major impact on the current understanding of repeats distributions throughout chromosomes in the red jungle fowl. CONCLUSIONS: Our results are a proof of concept of the reliability of using de novo tools to annotate repeats in large animal genomes. They have also revealed issues that will need to be resolved in order to develop gold-standard methodologies for annotating repeats in eukaryote genomes.

DensityMap: a genome viewer for illustrating the densities of features.

BACKGROUND: Several tools are available for visualizing genomic data. Some, such as Gbrowse and Jbrowse, are very efficient for small genomic regions, but they are not suitable for entire genomes. Others, like Phenogram and CViT, can be used to visualise whole genomes, but are not designed to display very dense genomic features (eg: interspersed repeats). We have therefore developed DensityMap, a lightweight Perl program that can display the densities of several features (genes, ncRNA, cpg, etc.) along chromosomes on the scale of the whole genome. A critical advantage of DensityMap is that it uses GFF annotation files directly to compute the densities of features without needing additional information from the user. The resulting picture is readily configurable, and the colour scales used can be customized for a best fit to the data plotted. RESULTS: DensityMap runs on Linux architecture with few requirements so that users can easily and quickly visualize the distributions and densities of genomic features for an entire genome. The input is GFF3-formated data representing chromosomes (linkage groups or pseudomolecules) and sets of features which are used to calculate representations in density maps. In practise, DensityMap uses a tilling window to compute the density of one or more features and the number of bases covered by these features along chromosomes. The densities are represented by colour scales that can be customized to highlight critical points. DensityMap can compare the distributions of features; it calculates several chromosomal density maps in a single image, each of which describes a different genomic feature. It can also use the genome nucleotide sequence to compute and plot a density map of the GC content along chromosomes. CONCLUSIONS: DensityMap is a compact, easily-used tool for displaying the distribution and density of all types of genomic features within a genome. It is flexible enough to visualize the densities of several types of features in a single representation. The images produced are readily configurable and their SVG format ensures that they can be edited.