Detection of DNA methylation signatures through the lens of genomic imprinting

Genomic imprinting represents an original model of epigenetic regulation resulting in a parent-of-origin expression. Despite the critical role of imprinted genes in mammalian growth, metabolism and neuronal function, there is no molecular tool specifically targeting them for a systematic evaluation. We show here that enzymatic methyl-seq consistently outperforms the bisulfite-based standard in capturing 165 candidate regions for genomic imprinting in the pig. This highlights the potential for a turnkey, fully customizable and reliable capture tool of genomic regions regulated by cytosine methylation in any population of interest. For the field of genomic imprinting, it opens up the possibility of detecting multilocus imprinting variations across the genome, with implications for basic research, agrigenomics and clinical practice.

Genome structures resolve the early diversification of teleost fishes

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

Rabbit targeted genomic sequences after heterologous hybridization using human exome

Causal mutations for major genes that underlie a broad range of morphological traits are often located within exons of genes that then affect protein functions. Non-model organism genetic studies are not easy to perform due to the lack of genome-wide molecular tools such as SNP genotyping array. Genotyping-By-Sequencing (GBS) methods offer an alternative. Consequently, we used this approach that is focused on the exome to target and identify major genes in rabbit populations.

Data description

We used a heterologous enrichment method before sequencing, allowing us to capture the rabbit exome using the marketed human panel since mammal protein coding genes are well conserved across the phylogenic tree of species. This targeted strategy was performed on 52 French rabbits from 5 different French strains (Californian, New-Zealand, Castor, Chinchilla and Laghmere). We generated 3.4 billion sequencing reads and approximately 29–140 million of reads per DNA sample. The expected exome coverage per sample ranged between 118 and 566X.

The present dataset could be useful for the scientific community working on rabbit species in order to (i) improve the annotation of the rabbit reference genome Oryctolagus cuniculus (OryCun2.0), (ii) enrich the characterization of polymorphisms segregating in rabbits and (iii) evaluate the genetic biodiversity in different rabbit strains.

Raw sequences were deposited in the European Nucleotide Archive (ENA) at the European Molecular Biology Laboratory- European Bioinformatics Institute (EMBL-EBI) data portal under bioproject accession number PRJEB37917.

Generation of a chromosome‐level genome assembly for Pacific halibut ( Hippoglossus stenolepis ) and characterization of its sex‐determining genomic region

The Pacific halibut (Hippoglossus stenolepis) is a key species in the North Pacific Ocean and Bering Sea ecosystems, where it also supports important fisheries. However, the lack of genomic resources limits our understanding of evolutionary, environmental and anthropogenic forces affecting key life history characteristics of Pacific halibut and prevents the application of genomic tools in fisheries management and conservation efforts. In the present study, we report on the first generation of a high‐quality chromosome‐level assembly of the Pacific halibut genome, with an estimated size of 602 Mb, 24 chromosome‐length scaffolds that contain 99.8% of the assembly and a N₅₀ scaffold length of 27.3 Mb. In the first application of this important resource, we conducted genome‐wide analyses of sex‐specific genetic variation by pool sequencing and characterized a potential sex‐determining region in chromosome 9 with a high density of female‐specific SNPs. Within this region, we identified the bmpr1ba gene as a potential candidate for master sex‐determining (MSD) gene. bmpr1ba is a member of the TGF‐β family that in teleosts has provided the largest number of MSD genes, including a paralogue of this gene in Atlantic herring. The genome assembly constitutes an essential resource for future studies on Pacific halibut population structure and dynamics, evolutionary history and responses to environmental and anthropogenic influences. Furthermore, the genomic location of the sex‐determining region in Pacific halibut has been identified and a putative candidate MSD gene has been proposed, providing further support for the rapid evolution of sex‐determining mechanisms in teleost fish.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

A duplicated copy of id2b is an unusual sex-determining candidate gene on the Y chromosome of arapaima (Arapaima gigas)

Arapaima gigas is one of the largest freshwater fish species of high ecological and economic importance. Overfishing and habitat destruction are severe threats to the remaining wild populations. By incorporating a chromosomal Hi-C contact map, we improved the arapaima genome assembly to chromosome-level, revealing an unexpected high degree of chromosome rearrangements during evolution of the bonytongues (Osteoglossiformes). Combining this new assembly with pool-sequencing of male and female genomes, we identified id2bbY, a duplicated copy of the inhibitor of DNA binding 2b (id2b) gene on the Y chromosome as candidate male sex-determining gene. A PCR-test for id2bbY was developed, demonstrating that this gene is a reliable male-specific marker for genotyping. Expression analyses showed that this gene is expressed in juvenile male gonads. Its paralog, id2ba, exhibits a male-biased expression in immature gonads. Transcriptome analyses and protein structure predictions confirm id2bbY as a prime candidate for the master sex-determiner. Acting through the TGFβ signaling pathway, id2bbY from arapaima would provide the first evidence for a link of this family of transcriptional regulators to sex determination. Our study broadens our current understanding about the evolution of sex determination genetic networks and provide a tool for improving arapaima aquaculture for commercial and conservation purposes.

Siberian sturgeon multi-tissue reference transcriptome database

Motivation: Siberian sturgeon is a long lived and late maturing fish farmed for caviar production in 50 countries. Functional genomics enable to find genes of interest for fish farming. In the absence of a reference genome, a reference transcriptome is very useful for sequencing based functional studies.

Results: We present here a high-quality transcriptome assembly database built using RNA-seq reads coming from brain, pituitary, gonadal, liver, stomach, kidney, anterior kidney, heart, embryonic and pre-larval tissues. It will facilitate crucial research on topics such as puberty, reproduction, growth, food intake and immunology. This database represents a major contribution to the publicly available sturgeon transcriptome reference datasets.

Availability: The database is publicly available at http://siberiansturgeontissuedb.sigenae.org

Supplementary information:  Supplementary data are available at Database online.

Synteny-Guided Resolution of Gene Trees Clarifies the Functional Impact of Whole-Genome Duplications

Whole-genome duplications (WGDs) have major impacts on the evolution of species, as they produce new gene copies contributing substantially to adaptation, isolation, phenotypic robustness, and evolvability. They result in large, complex gene families with recurrent gene losses in descendant species that sequence-based phylogenetic methods fail to reconstruct accurately. As a result, orthologs and paralogs are difficult to identify reliably in WGD-descended species, which hinders the exploration of functional consequences of WGDs. Here, we present Synteny-guided CORrection of Paralogies and Orthologies (SCORPiOs), a novel method to reconstruct gene phylogenies in the context of a known WGD event. WGDs generate large duplicated syntenic regions, which SCORPiOs systematically leverages as a complement to sequence evolution to infer the evolutionary history of genes. We applied SCORPiOs to the 320-My-old WGD at the origin of teleost fish. We find that almost one in four teleost gene phylogenies in the Ensembl database (3,394) are inconsistent with their syntenic contexts. For 70% of these gene families (2,387), we were able to propose an improved phylogenetic tree consistent with both the molecular substitution distances and the local syntenic information. We show that these synteny-guided phylogenies are more congruent with the species tree, with sequence evolution and with expected expression conservation patterns than those produced by state-of-the-art methods. Finally, we show that synteny-guided gene trees emphasize contributions of WGD paralogs to evolutionary innovations in the teleost clade.

Multi-tissue Siberian sturgeon RNA sequencing data

Siberian sturgeon, Acipenser baerii, is a commercially valuable fish for flesh and caviar production and a threatened species. We produced transcriptomic data for ten tissues with relevance to puberty, reproduction, early development, growth and food intake. The data includes RNA-Seq read sets of brain, pituitary, anterior-kidney, kidney, stomach, liver, heart, embryonic, pre-larval, and immature gonad sequences. Tissues were collected from sex differentiated fish (17 to 42 months of age, 66 to 85 cm) RNA was extracted and sequenced. Our purpose is to facilitate fundamental studies of sturgeon physiology to wild and aquaculture populations management.

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.

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.

Functional genomics of the digestive tract in broilers

Background

The sustainability of poultry farming relies on the development of more efficient and autonomous production systems in terms of feed supply. This implies a better integration of adaptive traits in breeding programs, including digestive efficiency, in order to favor the use of a wider variety of feedstuffs. The aim of the project was to improve the understanding of genes involved in digestive functions by characterizing the transcriptome of different sections of the digestive tract: the junction between the proventriculus and the gizzard, the gizzard, the gastroduodenal junction, and the jejunum.

Results

Total RNA from the four tissues were sequenced on a HiSeq2500 for six 23-day-old chickens from a second generation (F2) cross between two lines that were divergent for their digestive efficiency (D+/D-). Bioinformatics and biostatistics analyses of the RNA-seq data showed a total of 11,040 differentially expressed transcripts between the four tissues. In total, seven clusters of genes with markedly different expression profiles were identified. Functional analysis on gene groups was performed using “Gene Ontology” and semantic similarity. It showed a significant enrichment of body immune defenses in the jejunum, and an enrichment of transcriptional activity in the gizzard. Moreover, an interesting enrichment for neurohormonal control of muscle contraction was found for the two gizzard’s junctions.

Conclusion

This analysis allows us to draw the first molecular portrait of the different sections of the digestive tract, which will serve as a basis for future studies on the genetic and physiological control of the response of the animal to feed variations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5344-z) contains supplementary material, which is available to authorized users.

Compacting and correcting Trinity and Oases RNA-Seq de novo assemblies

Background

De novo transcriptome assembly of short reads is now a common step in expression analysis of organisms lacking a reference genome sequence. Several software packages are available to perform this task. Even if their results are of good quality it is still possible to improve them in several ways including redundancy reduction or error correction. Trinity and Oases are two commonly used de novo transcriptome assemblers. The contig sets they produce are of good quality. Still, their compaction (number of contigs needed to represent the transcriptome) and their quality (chimera and nucleotide error rates) can be improved.

Results

We built a de novo RNA-Seq Assembly Pipeline (DRAP) which wraps these two assemblers (Trinity and Oases) in order to improve their results regarding the above-mentioned criteria. DRAP reduces from 1.3 to 15 fold the number of resulting contigs of the assemblies depending on the read set and the assembler used. This article presents seven assembly comparisons showing in some cases drastic improvements when using DRAP. DRAP does not significantly impair assembly quality metrics such are read realignment rate or protein reconstruction counts.

Conclusion

Transcriptome assembly is a challenging computational task even if good solutions are already available to end-users, these solutions can still be improved while conserving the overall representation and quality of the assembly. The de novo RNA-Seq Assembly Pipeline (DRAP) is an easy to use software package to produce compact and corrected transcript set. DRAP is free, open-source and available under GPL V3 license at http://www.sigenae.org/drap.

Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database

With more than 30,000 species, ray-finned fish represent approximately half of vertebrates. The evolution of ray-finned fish was impacted by several whole genome duplication (WGD) events including a teleost-specific WGD event (TGD) that occurred at the root of the teleost lineage about 350 million years ago (Mya) and more recent WGD events in salmonids, carps, suckers and others. In plants and animals, WGD events are associated with adaptive radiations and evolutionary innovations. WGD-spurred innovation may be especially relevant in the case of teleost fish, which colonized a wide diversity of habitats on earth, including many extreme environments. Fish biodiversity, the use of fish models for human medicine and ecological studies, and the importance of fish in human nutrition, fuel an important need for the characterization of gene expression repertoires and corresponding evolutionary histories of ray-finned fish genes. To this aim, we performed transcriptome analyses and developed the PhyloFish database to provide (i) de novo assembled gene repertoires in 23 different ray-finned fish species including two holosteans (i.e. a group that diverged from teleosts before TGD) and 21 teleosts (including six salmonids), and (ii) gene expression levels in ten different tissues and organs (and embryos for many) in the same species. This resource was generated using a common deep RNA sequencing protocol to obtain the most exhaustive gene repertoire possible in each species that allows between-species comparisons to study the evolution of gene expression in different lineages. The PhyloFish database described here can be accessed and searched using RNAbrowse, a simple and efficient solution to give access to RNA-seq de novo assembled transcripts.

The BioMart community portal: an innovative alternative to large, centralized data repositories

The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

High-throughput sequencing analyses of XX genital ridges lacking FOXL2 reveal DMRT1 up-regulation before SOX9 expression during the sex-reversal process in goats

FOXL2 loss of function in goats leads to the early transdifferentiation of ovaries into testes, then to the full sex reversal of XX homozygous mutants. By contrast, Foxl2 loss of function in mice induces an arrest of follicle formation after birth, followed by complete female sterility. In order to understand the molecular role of FOXL2 during ovarian differentiation in the goat species, putative FOXL2 target genes were determined at the earliest stage of gonadal sex-specific differentiation by comparing the mRNA profiles of XX gonads expressing the FOXL2 protein or not. Of these 163 deregulated genes, around two-thirds corresponded to testicular genes that were up-regulated when FOXL2 was absent, and only 19 represented female-associated genes, down-regulated in the absence of FOXL2. FOXL2 should therefore be viewed as an antitestis gene rather than as a female-promoting gene. In particular, the key testis-determining gene DMRT1 was found to be up-regulated ahead of SOX9, thus suggesting in goats that SOX9 primary up-regulation may require DMRT1. Overall, our results equated to FOXL2 being an antitestis gene, allowing us to propose an alternative model for the sex-determination process in goats that differs slightly from that demonstrated in mice.

RNAbrowse: RNA-Seq de novo assembly results browser

Transcriptome analysis based on a de novo assembly of next generation RNA sequences is now performed routinely in many laboratories. The generated results, including contig sequences, quantification figures, functional annotations and variation discovery outputs are usually bulky and quite diverse. This article presents a user oriented storage and visualisation environment permitting to explore the data in a top-down manner, going from general graphical views to all possible details. The software package is based on biomart, easy to install and populate with local data. The software package is available under the GNU General Public License (GPL) at http://bioinfo.genotoul.fr/RNAbrowse.

Induced pluripotent stem cells derived from rabbits exhibit some characteristics of naïve pluripotency

Not much is known about the molecular and functional features of pluripotent stem cells (PSCs) in rabbits. To address this, we derived and characterized 2 types of rabbit PSCs from the same breed of New Zealand White rabbits: 4 lines of embryonic stem cells (rbESCs), and 3 lines of induced PSCs (rbiPSCs) that were obtained by reprogramming adult skin fibroblasts. All cell lines required fibroblast growth factor 2 for their growth and proliferation. All rbESC lines showed molecular and functional properties typically associated with primed pluripotency. The cell cycle of rbESCs had a prolonged G1 phase and a DNA damage checkpoint before entry into the S phase, which are the 2 features typically associated with the somatic cell cycle. In contrast, the rbiPSC lines exhibited some characteristics of naïve pluripotency, including resistance to single-cell dissociation by trypsin, robust activity of the distal enhancer of the mouse Oct4 gene, and expression of naïve pluripotency-specific genes, as defined in rodents. According to gene expression profiles, rbiPSCs were closer to the rabbit inner cell mass (ICM) than rbESCs. Furthermore, rbiPSCs were capable of colonizing the ICM after aggregation with morulas. Therefore, we propose that rbiPSCs self-renew in an intermediate state between naïve and primed pluripotency, which represents a key step toward the generation of bona fide naïve PSC lines in rabbits.

Induced pluripotent stem cells derived from rabbits exhibit some characteristics of naïve pluripotency

Not much is known about the molecular and functional features of pluripotent stem cells (PSCs) in rabbits. To address this, we derived and characterized 2 types of rabbit PSCs from the same breed of New Zealand White rabbits: 4 lines of embryonic stem cells (rbESCs), and 3 lines of induced PSCs (rbiPSCs) that were obtained by reprogramming adult skin fibroblasts. All cell lines required fibroblast growth factor 2 for their growth and proliferation. All rbESC lines showed molecular and functional properties typically associated with primed pluripotency. The cell cycle of rbESCs had a prolonged G1 phase and a DNA damage checkpoint before entry into the S phase, which are the 2 features typically associated with the somatic cell cycle. In contrast, the rbiPSC lines exhibited some characteristics of naïve pluripotency, including resistance to single-cell dissociation by trypsin, robust activity of the distal enhancer of the mouse Oct4 gene, and expression of naïve pluripotency-specific genes, as defined in rodents. According to gene expression profiles, rbiPSCs were closer to the rabbit inner cell mass (ICM) than rbESCs. Furthermore, rbiPSCs were capable of colonizing the ICM after aggregation with morulas. Therefore, we propose that rbiPSCs self-renew in an intermediate state between naïve and primed pluripotency, which represents a key step toward the generation of bona fide naïve PSC lines in rabbits.

Fine mapping of complex traits in non-model species: using next generation sequencing and advanced intercross lines in Japanese quail

Background

As for other non-model species, genetic analyses in quail will benefit greatly from a higher marker density, now attainable thanks to the evolution of sequencing and genotyping technologies. Our objective was to obtain the first genome wide panel of Japanese quail SNP (Single Nucleotide Polymorphism) and to use it for the fine mapping of a QTL for a fear-related behaviour, namely tonic immobility, previously localized on Coturnix japonica chromosome 1. To this aim, two reduced representations of the genome were analysed through high-throughput 454 sequencing: AFLP (Amplified Fragment Length Polymorphism) fragments as representatives of genomic DNA, and EST (Expressed Sequence Tag) as representatives of the transcriptome.

Results

The sequencing runs produced 399,189 and 1,106,762 sequence reads from cDNA and genomic fragments, respectively. They covered over 434 Mb of sequence in total and allowed us to detect 17,433 putative SNP. Among them, 384 were used to genotype two Advanced Intercross Lines (AIL) obtained from three quail lines differing for duration of tonic immobility. Despite the absence of genotyping for founder individuals in the analysis, the previously identified candidate region on chromosome 1 was refined and led to the identification of a candidate gene.

Conclusions

These data confirm the efficiency of transcript and AFLP-sequencing for SNP discovery in a non-model species, and its application to the fine mapping of a complex trait. Our results reveal a significant association of duration of tonic immobility with a genomic region comprising the DMD (dystrophin) gene. Further characterization of this candidate gene is needed to decipher its putative role in tonic immobility in Coturnix.

Transcriptomic profiling of proteases and antiproteases in the liver of sexually mature hens in relation to vitellogenesis

Background

Most egg yolk precursors are synthesized by the liver, secreted into the blood and transferred into oocytes, to provide nutrients and bioactive molecules for the avian embryo. Three hundred and sixteen distinct proteins have been identified in egg yolk. These include 37 proteases and antiproteases, which are likely to play a role in the formation of the yolk (vitellogenesis), as regulators of protein metabolism. We used a transcriptomic approach to define the protease and antiprotease genes specifically expressed in the hen liver in relation to vitellogenesis by comparing sexually mature and pre-laying chickens showing different steroid milieu.

Results

Using a 20 K chicken oligoarray, a total of 582 genes were shown to be over-expressed in the liver of sexually mature hens (1.2 to 67 fold-differences). Eight of the top ten over-expressed genes are known components of the egg yolk or perivitelline membrane. This list of 582 genes contains 12 proteases and 3 antiproteases. We found that “uncharacterized protein LOC419301/similar to porin” (GeneID:419301), an antiprotease and “cathepsin E-A-like/similar to nothepsin” (GeneID:417848), a protease, were the only over-expressed candidates (21-fold and 35-fold difference, respectively) that are present in the egg yolk. Additionally, we showed the 4-fold over-expression of “ovochymase-2/similar to oviductin” (GeneID:769290), a vitelline membrane-specific protease.

Conclusions

Our approach revealed that three proteases and antiproteases are likely to participate in the formation of the yolk. The role of the other 12 proteases and antiproteases which are over-expressed in our model remains unclear. At least 1/3 of proteases and antiproteases identified in egg yolk and vitelline membrane proteomes are expressed similarly in the liver regardless of the maturity of hens, and have been initially identified as regulators of haemostasis and inflammatory events. The lack of effect of sex steroids on these genes expressed in the liver but the products of which are found in the yolk suggests that these may be passively incorporated into the yolk rather than actively produced for that purpose. These results raise the question of the biological significance of egg yolk proteases and antiproteases, and more generally of all minor proteins that have been identified in egg yolk.

Development and characterisation of an expressed sequence tags (EST)-derived single nucleotide polymorphisms (SNPs) resource in rainbow trout

Background

There is considerable interest in developing high-throughput genotyping with single nucleotide polymorphisms (SNPs) for the identification of genes affecting important ecological or economical traits. SNPs are evenly distributed throughout the genome and are likely to be functionally relevant. In rainbow trout, in silico screening of EST databases represents an attractive approach for de novo SNP identification. Nevertheless, EST sequencing errors and assembly of EST paralogous sequences can lead to the identification of false positive SNPs which renders the reliability of EST-derived SNPs relatively low. Further validation of EST-derived SNPs is therefore required. The objective of this work was to assess the quality of and to validate a large number of rainbow trout EST-derived SNPs.

Results

A panel of 1,152 EST-derived SNPs was selected from the INRA Sigenae SNP database and was genotyped in standard and double haploid individuals from several populations using the Illumina GoldenGate BeadXpress assay. High-quality genotyping data were obtained for 958 SNPs representing a genotyping success rate of 83.2 %, out of which, 350 SNPs (36.5 %) were polymorphic in at least one population and were designated as true SNPs. They also proved to be a potential tool to investigate genetic diversity of the species, as the set of SNP successfully sorted individuals into three main groups using STRUCTURE software. Functional annotations revealed 28 non-synonymous SNPs, out of which four substitutions were predicted to affect protein functions. A subset of 223 true SNPs were polymorphic in the two INRA mapping reference families and were integrated into the INRA microsatellite-based linkage map.

Conclusions

Our results represent the first study of EST-derived SNPs validation in rainbow trout, a species whose genome sequences is not yet available. We designed several specific filters in order to improve the genotyping yield. Nevertheless, our selection criteria should be further improved in order to reduce the observed high rate of false positive SNPs which results from the occurrence of whole genome duplications.

Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources

Background

Gene expression profiling studies of mastitis in ruminants have provided key but fragmented knowledge for the understanding of the disease. A systematic combination of different expression profiling studies via meta-analysis techniques has the potential to test the extensibility of conclusions based on single studies. Using the program Pointillist, we performed meta-analysis of transcription-profiling data from six independent studies of infections with mammary gland pathogens, including samples from cattle challenged in vivo with S. aureus, E. coli, and S. uberis, samples from goats challenged in vivo with S. aureus, as well as cattle macrophages and ovine dendritic cells infected in vitro with S. aureus. We combined different time points from those studies, testing different responses to mastitis infection: overall (common signature), early stage, late stage, and cattle-specific.

Results

Ingenuity Pathway Analysis of affected genes showed that the four meta-analysis combinations share biological functions and pathways (e.g. protein ubiquitination and polyamine regulation) which are intrinsic to the general disease response. In the overall response, pathways related to immune response and inflammation, as well as biological functions related to lipid metabolism were altered. This latter observation is consistent with the milk fat content depression commonly observed during mastitis infection. Complementarities between early and late stage responses were found, with a prominence of metabolic and stress signals in the early stage and of the immune response related to the lipid metabolism in the late stage; both mechanisms apparently modulated by few genes, including XBP1 and SREBF1.

The cattle-specific response was characterized by alteration of the immune response and by modification of lipid metabolism. Comparison of E. coli and S. aureus infections in cattle in vivo revealed that affected genes showing opposite regulation had the same altered biological functions and provided evidence that E. coli caused a stronger host response.

Conclusions

This meta-analysis approach reinforces previous findings but also reveals several novel themes, including the involvement of genes, biological functions, and pathways that were not identified in individual studies. As such, it provides an interesting proof of principle for future studies combining information from diverse heterogeneous sources.

PhyleasProg: a user-oriented web server for wide evolutionary analyses

Evolutionary analyses of biological data are becoming a prerequisite in many fields of biology. At a time of high-throughput data analysis, phylogenetics is often a necessary complementary tool for biologists to understand, compare and identify the functions of sequences. But available bioinformatics tools are frequently not easy for non-specialists to use. We developed PhyleasProg (http://phyleasprog.inra.fr), a user-friendly web server as a turnkey tool dedicated to evolutionary analyses. PhyleasProg can help biologists with little experience in evolutionary methodologies by analysing their data in a simple and robust way, using methods corresponding to robust standards. Via a very intuitive web interface, users only need to enter a list of Ensembl protein IDs and a list of species as inputs. After dynamic computations, users have access to phylogenetic trees, positive/purifying selection data (on site and branch-site models), with a display of these results on the protein sequence and on a 3D structure model, and the synteny environment of related genes. This connection between different domains of phylogenetics opens the way to new biological analyses for the discovery of the function and structure of proteins.

Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg

BACKGROUND

As uricoletic animals, chickens produce cleidoic eggs, which are self-contained bacteria-resistant biological packages for extra-uterine development of the chick embryo. The eggshell constitutes a natural physical barrier against bacterial penetration if it forms correctly and remains intact. The eggshell's remarkable mechanical properties are due to interactions among mineral components and the organic matrix proteins. The purpose of our study was to identify novel eggshell proteins by examining the transcriptome of the uterus during calcification of the eggshell. An extensive bioinformatic analysis on genes over-expressed in the uterus allowed us to identify novel eggshell proteins that contribute to the egg's natural defenses.

RESULTS

Our 14 K Del-Mar Chicken Integrated Systems microarray was used for transcriptional profiling in the hen's uterus during eggshell deposition. A total of 605 transcripts were over-expressed in the uterus compared with the magnum or white isthmus across a wide range of abundance (1.1- to 79.4-fold difference). The 605 highly-expressed uterine transcripts correspond to 469 unique genes, which encode 437 different proteins. Gene Ontology (GO) analysis was used for interpretation of protein function. The most over-represented GO terms are related to genes encoding ion transport proteins, which provide eggshell mineral precursors. Signal peptide sequence was found for 54 putative proteins secreted by the uterus during eggshell formation. Many functional proteins are involved in calcium binding or biomineralization--prerequisites for interacting with the mineral phase during eggshell fabrication. While another large group of proteins could be involved in proper folding of the eggshell matrix. Many secreted uterine proteins possess antibacterial properties, which would protect the egg against microbial invasion. A final group includes proteases and protease inhibitors that regulate protein activity in the acellular uterine fluid where eggshell formation takes place.

CONCLUSIONS

Our original study provides the first detailed description of the chicken uterus transcriptome during formation of the eggshell. We have discovered a cache of about 600 functional genes and identified a large number of encoded proteins secreted into uterine fluid for fabrication of the eggshell and chemical protection of the egg. Some of these uterine genes could prove useful as biological markers for genetic improvement of phenotypic traits (i.e., egg and eggshell quality).

Revealing the dynamics of gene expression during embryonic genome activation and first differentiation in the rabbit embryo with a dedicated array screening

Early mammalian development is characterized by extensive changes in nuclear functions that result from epigenetic modifications of the newly formed embryonic genome. While the first embryonic cells are totipotent, this status spans only a few cell cycles. At the blastocyst stage, the embryo already contains differentiated trophectoderm cells and pluripotent inner cell mass cells. Concomitantly, the embryonic genome becomes progressively transcriptionally active. During this unique period of development, the gene expression pattern has been mainly characterized in the mouse, in which embryonic genome activation (EGA) spans a single cell cycle after abrupt epigenetic modifications. To further characterize this period, we chose to analyze it in the rabbit, in which, as in most mammals, EGA is more progressive and occurs closer to the first cell differentiation events. In this species, for which no transcriptomic arrays were available, we focused on genes expressed at EGA and first differentiation and established a 2,000-gene dedicated cDNA array. Screening this with pre-EGA, early post-EGA, and blastocyst embryos divided genes into seven clusters of expression according to their regulation during this period and revealed their dynamics of expression during EGA and first differentiation. Our results point to transient properties of embryo transcriptome at EGA, due not only to the transition between maternal and embryonic transcripts but also to the transient expression of a subset of embryonic genes whose functions remained largely uncharacterized. They also provide a first view of the functional consequences of the changes in gene expression program.

Identification of transcripts involved in meiosis and follicle formation during ovine ovary development

Background

The key steps in germ cell survival during ovarian development are the entry into meiosis of oogonies and the formation of primordial follicles, which then determine the reproductive lifespan of the ovary. In sheep, these steps occur during fetal life, between 55 and 80 days of gestation, respectively. The aim of this study was to identify differentially expressed ovarian genes during prophase I meiosis and early folliculogenesis in sheep.

Results

In order to elucidate the molecular events associated with early ovarian differentiation, we generated two ovary stage-specific subtracted cDNA libraries using SSH. Large-scale sequencing of these SSH libraries identified 6,080 ESTs representing 2,535 contigs. Clustering and assembly of these ESTs resulted in a total of 2,101 unique sequences depicted in 1,305 singleton (62.11%) and 796 contigs (37.9%) ESTs (clusters). BLASTX evaluation indicated that 99% of the ESTs were homologous to various known genes/proteins in a broad range of organisms, especially ovine, bovine and human species. The remaining 1% which exhibited any homology to known gene sequences was considered as novel. Detailed study of the expression patterns of some of these genes using RT-PCR revealed new promising candidates for ovary differentiation genes in sheep.

Conclusion

We showed that the SSH approach was relevant to determining new mammalian genes which might be involved in oogenesis and early follicle development, and enabled the discovery of new potential oocyte and granulosa cell markers for future studies. These genes may have significant implications regarding our understanding of ovarian function in molecular terms, and for the development of innovative strategies to both promote and control fertility.

Molecular evidence for a critical period in mural trophoblast development in bovine blastocysts

Embryonic and extra-embryonic lineages are separated at the blastocyst stage in the mouse at the onset of implantation but well ahead of implantation in most mammals. To provide information on the development of the trophoblast lineage in late-implanting bovine embryos, we combined the use of molecular markers defining embryonic and extra-embryonic lineages in the mouse with a transcriptomic approach dedicated to the early steps of the elongation process, a characteristic feature of blastocyst development in ruminants. In this study, we present molecular evidence for differences between the cow and the mouse in the programming of trophoblast differentiation. This different programming encompasses: (i) the expression of epiblast specifying genes (Oct-4, Nanog) in bovine trophoblast cells at the onset of elongation, (ii) the transcription of proliferation markers in early elongating blastocysts, (iii) the early detection of trophoblast-specific transcripts related to extra-embryonic tissue's differentiation (Hand1, Ets2, IFN-tau) and (iv) the identification of a new transcript (c12) which displays a reciprocal pattern to that of Oct-4 and Nanog genes in the embryonic cells and for which no equivalent has thus far been found in the mouse. Altogether, these results tended to show that early elongation is a critical transition in bovine trophoblast development.

Specialized microbial databases for inductive exploration of microbial genome sequences

Background

The enormous amount of genome sequence data asks for user-oriented databases to manage sequences and annotations. Queries must include search tools permitting function identification through exploration of related objects.

Methods

The GenoList package for collecting and mining microbial genome databases has been rewritten using MySQL as the database management system. Functions that were not available in MySQL, such as nested subquery, have been implemented.

Results

Inductive reasoning in the study of genomes starts from "islands of knowledge", centered around genes with some known background. With this concept of "neighborhood" in mind, a modified version of the GenoList structure has been used for organizing sequence data from prokaryotic genomes of particular interest in China. GenoChore , a set of 17 specialized end-user-oriented microbial databases (including one instance of Microsporidia, Encephalitozoon cuniculi, a member of Eukarya) has been made publicly available. These databases allow the user to browse genome sequence and annotation data using standard queries. In addition they provide a weekly update of searches against the world-wide protein sequences data libraries, allowing one to monitor annotation updates on genes of interest. Finally, they allow users to search for patterns in DNA or protein sequences, taking into account a clustering of genes into formal operons, as well as providing extra facilities to query sequences using predefined sequence patterns.

Conclusion

This growing set of specialized microbial databases organize data created by the first Chinese bacterial genome programs (ThermaList, Thermoanaerobacter tencongensis, LeptoList, with two different genomes of Leptospira interrogans and SepiList, Staphylococcus epidermidis) associated to related organisms for comparison.

Contribution to high-resolution mapping in pigs with 101 type I markers and progress in comparative map between humans and pigs

In the frame of the European program GenetPig, we localized on the Pig map 105 coding sequences (type I markers) from different origins, using INRA-University of Minnesota porcine Radiation Hybrid Panel (IMpRH, 101 markers) and somatic cell hybrid panel (SCHP, 93 markers, of which only four were not also mapped using IMpRH). Thus, we contributed to the improvement of the porcine high-resolution map, and we complemented the integration between the RH and cytogenetic maps. IMpRH tools allowed us to map 101 new markers relatively to reference markers of the first generation radiation hybrid map. Ninety out of 101 markers are linked to an already mapped marker with a LOD score greater than 4.8. Seventy-eight markers were informative for comparative mapping. Comparison of marker positions on the RH map with those obtained on the cytogenetic map or those expected by Human-Pig comparative map data suggested to us to be cautious with markers linked with a LOD lower than 6. These results allowed us to specify chromosomal fragments well conserved between humans and pigs and also to suggest new correspondences (Sscr1-Hsap3, Sscr9-Hsap9, Sscr13-Hsap11, Sscr15-Hsap6) confirmed by FISH on pig chromosomes. We examined in more detail the comparative map between Hsap12 and Sscr5 considering gene order, which suggests that rearrangements have occurred within the conserved synteny.

In the QTL region surrounding porcine MHC, gene order is conserved with human genome

On the porcine genome, the region surrounding the Major Histocompatibility Complex, also called Swine Leukocyte Antigens (SLA), is of particular interest not only owing to itq role in the control of immune response, but also because of its influence on many traits such as growth, fatness, and meat quality. To help in the identification of responsible genes, detailed comparative maps of the MHC region in mammalian species and powerful mapping tools allowing accurate ordering of genes and markers in this region are needed. In this report, we describe the use of the recently developed IMpRH radiation hybrid panel, to construct a higher density radiation hybrid map of swine Sscr 7p-q12, containing 23 additional loci. Our results show that the gene order is conserved between the two MHC-containing regions, even if an inversion is observed above the QTL region in the region containing DEK, SCA1, and EDN1 genes. The framework map produced shows that the IMpRH panel permits the ordering of genes and markers in the three MHC classes and would thus allow accurate localization of ESTs and candidate genes.

IMpRH server: an RH mapping server available on the Web

SUMMARY

The INRA-Minnesota Porcine Radiation Hybrid (IMpRH) Server provides both a mapping tool (IMpRH mapping tool) and a database (IMpRH database) of officially submitted results. The mapping tool permits the mapping of a new marker relatively to markers previously mapped on the IMpRH panel. The IMpRH database is the official database for submission of new results and queries. The database not only permits the sharing of public data but also semi-private and private data.