Differentially expressed genes and gene networks involved in pig ovarian follicular atresia.

Ovarian folliculogenesis corresponds to the development of follicles leading to either ovulation or degeneration, this latter process being called atresia. Even if atresia involves apoptosis, its mechanism is not well understood. The objective of this study was to analyze global gene expression in pig granulosa cells of ovarian follicles during atresia. The transcriptome analysis was performed on a 9,216 cDNA microarray to identify gene networks and candidate genes involved in pig ovarian follicular atresia. We found 1,684 significantly regulated genes to be differentially regulated between small healthy follicles and small atretic follicles. Among them, 287 genes had a fold-change higher than two between the two follicle groups. Eleven genes (DKK3, GADD45A, CAMTA2, CCDC80, DAPK2, ECSIT, MSMB, NUPR1, RUNX2, SAMD4A, and ZNF628) having a fold-change higher than five between groups could likely serve as markers of follicular atresia. Moreover, automatic confrontation of deregulated genes with literature data highlighted 93 genes as regulatory candidates of pig granulosa cell atresia. Among these genes known to be inhibitors of apoptosis, stimulators of apoptosis, or tumor suppressors INHBB, HNF4, CLU, different interleukins (IL5, IL24), TNF-associated receptor (TNFR1), and cytochrome-c oxidase (COX) were suggested as playing an important role in porcine atresia. The present study also enlists key upstream regulators in follicle atresia based on our results and on a literature review. The novel gene candidates and gene networks identified in the current study lead to a better understanding of the molecular regulation of ovarian follicular atresia.

Transcriptome-wide investigation of genomic imprinting in chicken.

Genomic imprinting is an epigenetic mechanism by which alleles of some specific genes are expressed in a parent-of-origin manner. It has been observed in mammals and marsupials, but not in birds. Until now, only a few genes orthologous to mammalian imprinted ones have been analyzed in chicken and did not demonstrate any evidence of imprinting in this species. However, several published observations such as imprinted-like QTL in poultry or reciprocal effects keep the question open. Our main objective was thus to screen the entire chicken genome for parental-allele-specific differential expression on whole embryonic transcriptomes, using high-throughput sequencing. To identify the parental origin of each observed haplotype, two chicken experimental populations were used, as inbred and as genetically distant as possible. Two families were produced from two reciprocal crosses. Transcripts from 20 embryos were sequenced using NGS technology, producing ∼200 Gb of sequences. This allowed the detection of 79 potentially imprinted SNPs, through an analysis method that we validated by detecting imprinting from mouse data already published. However, out of 23 candidates tested by pyrosequencing, none could be confirmed. These results come together, without a priori, with previous statements and phylogenetic considerations assessing the absence of genomic imprinting in chicken.

Muscle transcriptomic investigation of late fetal development identifies candidate genes for piglet maturity.

BACKGROUND: In pigs, the perinatal period is the most critical time for survival. Piglet maturation, which occurs at the end of gestation, leads to a state of full development after birth. Therefore, maturity is an important determinant of early survival. Skeletal muscle plays a key role in adaptation to extra-uterine life, e.g. glycogen storage and thermoregulation. In this study, we performed microarray analysis to identify the genes and biological processes involved in piglet muscle maturity. Progeny from two breeds with extreme muscle maturity phenotypes were analyzed at two time points during gestation (gestational days 90 and 110). The Large White (LW) breed is a selected breed with an increased rate of mortality at birth, whereas the Meishan (MS) breed produces piglets with extremely low mortality at birth. The impact of the parental genome was analyzed with reciprocal crossed fetuses. RESULTS: Microarray analysis identified 12,326 differentially expressed probes for gestational age and genotype. Such a high number reflects an important transcriptomic change that occurs between 90 and 110 days of gestation. 2,000 probes, corresponding to 1,120 unique annotated genes, involved more particularly in the maturation process were further studied. Functional enrichment and graph inference studies underlined genes involved in muscular development around 90 days of gestation, and genes involved in metabolic functions, such as gluconeogenesis, around 110 days of gestation. Moreover, a difference in the expression of key genes, e.g. PCK2, LDHA or PGK1, was detected between MS and LW just before birth. Reciprocal crossing analysis resulted in the identification of 472 genes with an expression preferentially regulated by one parental genome. Most of these genes (366) were regulated by the paternal genome. Among these paternally regulated genes, some known imprinted genes, such as MAGEL2 or IGF2, were identified and could have a key role in the maturation process. CONCLUSION: These results reveal the biological mechanisms that regulate muscle maturity in piglets. Maturity is also under the conflicting regulation of the parental genomes. Crucial genes, which could explain the biological differences in maturity observed between LW and MS breeds, were identified. These genes could be excellent candidates for a key role in the maturity.

Sperm DNA methylation analysis in swine reveals conserved and species-specific methylation patterns and highlights an altered methylation at the GNAS locus in infertile boars.

Male infertility is an increasing health issue in today's society for both human and livestock populations. In livestock, male infertility slows the improvement of animal selection programs and agricultural productivity. There is increasing evidence that epigenetic marks play an important role in the production of good-quality sperm. We therefore screened for specific or common epigenetic signatures of livestock infertility. To do so, we compared DNA methylation level in sperm DNA from fertile and infertile boars. We evaluated first the global level of sperm DNA methylation and found no difference between the two groups of boars. We then selected 42 loci of interest, most of them known to be imprinted in human or mice, and assessed the imprinting status of five of them not previously described in swine tissues: WT1, CNTN3, IMPACT, QPCT, and GRB10. DNA methylation level was then quantified in fertile and infertile boars at these 42 loci. Results from fertile boars indicated that the methylation level of the selected loci is highly conserved between pig, human, and mice, with a few exceptions, including the POU5F1 (OCT4) promoter and RTL1. Comparison between fertile and infertile boars revealed that one imprinted region, the GNAS locus, shows an increase in sperm DNA methylation in three out of eight infertile boars with low semen quality. This increase in DNA methylation is associated with an altered expression of the genes belonging to the GNAS locus, suggesting a new role for GNAS in the proper formation of functional gametes.

Cryptic patterning of avian skin confers a developmental facility for loss of neck feathering.

Vertebrate skin is characterized by its patterned array of appendages, whether feathers, hairs, or scales. In avian skin the distribution of feathers occurs on two distinct spatial levels. Grouping of feathers within discrete tracts, with bare skin lying between the tracts, is termed the macropattern, while the smaller scale periodic spacing between individual feathers is referred to as the micropattern. The degree of integration between the patterning mechanisms that operate on these two scales during development and the mechanisms underlying the remarkable evolvability of skin macropatterns are unknown. A striking example of macropattern variation is the convergent loss of neck feathering in multiple species, a trait associated with heat tolerance in both wild and domestic birds. In chicken, a mutation called Naked neck is characterized by a reduction of body feathering and completely bare neck. Here we perform genetic fine mapping of the causative region and identify a large insertion associated with the Naked neck trait. A strong candidate gene in the critical interval, BMP12/GDF7, displays markedly elevated expression in Naked neck embryonic skin due to a cis-regulatory effect of the causative mutation. BMP family members inhibit embryonic feather formation by acting in a reaction-diffusion mechanism, and we find that selective production of retinoic acid by neck skin potentiates BMP signaling, making neck skin more sensitive than body skin to suppression of feather development. This selective production of retinoic acid by neck skin constitutes a cryptic pattern as its effects on feathering are not revealed until gross BMP levels are altered. This developmental modularity of neck and body skin allows simple quantitative changes in BMP levels to produce a sparsely feathered or bare neck while maintaining robust feather patterning on the body.

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.

Integrative mapping analysis of chicken microchromosome 16 organization.

BACKGROUND: The chicken karyotype is composed of 39 chromosome pairs, of which 9 still remain totally absent from the current genome sequence assembly, despite international efforts towards complete coverage. Some others are only very partially sequenced, amongst which microchromosome 16 (GGA16), particularly under-represented, with only 433 kb assembled for a full estimated size of 9 to 11 Mb. Besides the obvious need of full genome coverage with genetic markers for QTL (Quantitative Trait Loci) mapping and major genes identification studies, there is a major interest in the detailed study of this chromosome because it carries the two genetically independent MHC complexes B and Y. In addition, GGA16 carries the ribosomal RNA (rRNA) genes cluster, also known as the NOR (nucleolus organizer region). The purpose of the present study is to construct and present high resolution integrated maps of GGA16 to refine its organization and improve its coverage with genetic markers. RESULTS: We developed 79 STS (Sequence Tagged Site) markers to build a physical RH (radiation hybrid) map and 34 genetic markers to extend the genetic map of GGA16. We screened a BAC (Bacterial Artificial Chromosome) library with markers for the MHC-B, MHC-Y and rRNA complexes. Selected clones were used to perform high resolution FISH (Fluorescent In Situ Hybridization) mapping on giant meiotic lampbrush chromosomes, allowing meiotic mapping in addition to the confirmation of the order of the three clusters along the chromosome. A region with high recombination rates and containing PO41 repeated elements separates the two MHC complexes. CONCLUSIONS: The three complementary mapping strategies used refine greatly our knowledge of chicken microchromosome 16 organisation. The characterisation of the recombination hotspots separating the two MHC complexes demonstrates the presence of PO41 repetitive sequences both in tandem and inverted orientation. However, this region still needs to be studied in more detail.

QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines.

Quantitative trait loci (QTL) for metabolic and body composition traits were mapped at 7 and 9 wk, respectively, in an F(2) intercross between high-growth and low-growth chicken lines. These lines also diverged for abdominal fat percentage (AFP) and plasma insulin-like growth factor-I (IGF-I), insulin, and glucose levels. Genotypings were performed with 129 microsatellite markers covering 21 chromosomes. A total of 21 QTL with genomewide level of significance were detected by single-trait analyses for body weight (BW), breast muscle weight (BMW) and percentage (BMP), AF weight (AFW) and percentage (AFP), shank length (ShL) and diameter (ShD), fasting plasma glucose level (Gluc), and body temperature (T(b)). Other suggestive QTL were identified for these parameters and for plasma IGF-I and nonesterified fatty acid levels. QTL controlling adiposity and Gluc were colocalized on GGA3 and GGA5 and QTL for BW, ShL and ShD, adiposity, and T(b) on GGA4. Multitrait analyses revealed two QTL controlling Gluc and AFP on GGA5 and Gluc and T(b) on GGA26. Significant effects of the reciprocal cross were observed on BW, ShD, BMW, and Gluc, which may result from mtDNA and/or maternal effects. Most QTL regions for Gluc and adiposity harbor genes for which alleles have been associated with increased susceptibility to diabetes and/or obesity in humans. Identification of genes responsible for these metabolic QTL will increase our understanding of the constitutive "hyperglycemia" found in chickens. Furthermore, a comparative approach could provide new information on the genetic causes of diabetes and obesity in humans.

Addition of the microchromosome GGA25 to the chicken genome sequence assembly through radiation hybrid and genetic mapping.

BACKGROUND: The publication of the first draft chicken sequence assembly became available in 2004 and was updated in 2006. However, this does not constitute a definitive and complete sequence of the chicken genome, since the microchromosomes are notably under-represented. In an effort to develop maps for the microchromosomes absent from the chicken genome assembly, we developed radiation hybrid (RH) and genetic maps with markers isolated from sequence currently assigned to "chromosome Unknown" (chrUn). The chrUn is composed of sequence contigs not assigned to named chromosomes. To identify and map sequence belonging to the microchromosomes we used a comparative mapping strategy, and we focused on the small linkage group E26C13. RESULTS: In total, 139 markers were analysed with the chickRH6 panel, of which 120 were effectively assigned to the E26C13 linkage group, the remainder mapping elsewhere in the genome. The final RH map is composed of 22 framework markers extending over a 245.6 cR distance. A corresponding genetic map was developed, whose length is 103 cM in the East Lansing reference population. The E26C13 group was assigned to GGA25 (Gallus gallus chromosome 25) by FISH (fluorescence in situ hybridisation) mapping. CONCLUSION: The high-resolution RH framework map obtained here covers the entire chicken chromosome 25 and reveals the existence of a high number of intrachromosomal rearrangements when compared to the human genome. The strategy used here for the characterization of GGA25 could be used to improve knowledge on the other uncharacterized small, yet gene-rich microchromosomes.

Construction of a radiation hybrid map of chicken chromosome 2 and alignment to the chicken draft sequence.

BACKGROUND: The ChickRH6 whole chicken genome radiation hybrid (RH) panel recently produced has already been used to build radiation hybrid maps for several chromosomes, generating comparative maps with the human and mouse genomes and suggesting improvements to the chicken draft sequence assembly. Here we present the construction of a RH map of chicken chromosome 2. Markers from the genetic map were used for alignment to the existing GGA2 (Gallus gallus chromosome 2) linkage group and EST were used to provide valuable comparative mapping information. Finally, all markers from the RH map were localised on the chicken draft sequence assembly to check for eventual discordances. RESULTS: Eighty eight microsatellite markers, 10 genes and 219 EST were selected from the genetic map or on the basis of available comparative mapping information. Out of these 317 markers, 270 gave reliable amplifications on the radiation hybrid panel and 198 were effectively assigned to GGA2. The final RH map is 2794 cR6000 long and is composed of 86 framework markers distributed in 5 groups. Conservation of synteny was found between GGA2 and eight human chromosomes, with segments of conserved gene order of varying lengths. CONCLUSION: We obtained a radiation hybrid map of chicken chromosome 2. Comparison to the human genome indicated that most of the 8 groups of conserved synteny studied underwent internal rearrangements. The alignment of our RH map to the first draft of the chicken genome sequence assembly revealed a good agreement between both sets of data, indicative of a low error rate.

A high-resolution radiation hybrid map of chicken chromosome 5 and comparison with human chromosomes.

BACKGROUND: The resolution of radiation hybrid (RH) maps is intermediate between that of the genetic and BAC (Bacterial Artificial Chromosome) contig maps. Moreover, once framework RH maps of a genome have been constructed, a quick location of markers by simple PCR on the RH panel is possible. The chicken ChickRH6 panel recently produced was used here to construct a high resolution RH map of chicken GGA5. To confirm the validity of the map and to provide valuable comparative mapping information, both markers from the genetic map and a high number of ESTs (Expressed Sequence Tags) were used. Finally, this RH map was used for testing the accuracy of the chicken genome assembly for chromosome 5. RESULTS: A total of 169 markers (21 microsatellites and 148 ESTs) were typed on the ChickRH6 RH panel, of which 134 were assigned to GGA5. The final map is composed of 73 framework markers extending over a 1315.6 cR distance. The remaining 61 markers were placed alongside the framework markers within confidence intervals. CONCLUSION: The high resolution framework map obtained in this study has markers covering the entire chicken chromosome 5 and reveals the existence of a high number of rearrangements when compared to the human genome. Only two discrepancies were observed in relation to the sequence assembly recently reported for this chromosome.

Non PCR-amplified Transcripts and AFLP fragments as reduced representations of the quail genome for 454 Titanium sequencing.

BACKGROUND: SNP (Single Nucleotide Polymorphism) discovery is now routinely performed using high-throughput sequencing of reduced representation libraries. Our objective was to adapt 454 GS FLX based sequencing methodologies in order to obtain the largest possible dataset from two reduced representations libraries, produced by AFLP (Amplified Fragment Length Polymorphism) for genomic DNA, and EST (Expressed Sequence Tag) for the transcribed fraction of the genome. FINDINGS: The expressed fraction was obtained by preparing cDNA libraries without PCR amplification from quail embryo and brain. To optimize the information content for SNP analyses, libraries were prepared from individuals selected in three quail lines and each individual in the AFLP library was tagged. Sequencing runs produced 399,189 sequence reads from cDNA and 373,484 from genomic fragments, covering close to 250 Mb of sequence in total. CONCLUSIONS: Both methods used to obtain reduced representations for high-throughput sequencing were successful after several improvements.The protocols may be used for several sequencing applications, such as de novo sequencing, tagged PCR fragments or long fragment sequencing of cDNA.

Detection of novel quantitative trait loci for cutaneous melanoma by genome-wide scan in the MeLiM swine model.

Human cutaneous melanoma is a complex trait inherited in about 10% of cases. Although 2 high-risk genes, CDKN2A and CDK4, and 1 low risk gene, MC1R, have been identified, susceptibility genes remain to be discovered. Here, we attempted to determine new genomic regions linked to melanoma using the pig MeLiM strain, which develops hereditary cutaneous melanomas. We applied quantitative trait loci (QTL) mapping method to a significant genome-wide scan performed on 331 backcross pigs derived from this strain. QTLs were detected at chromosome-wide level for a melanoma synthetic trait corresponding to the development of melanoma. The peak positions on Sus scrofa chromosomes (SSC) were at 49.4 and 88.0 cM (SSC1), 56.0 cM (SSC13), 86.5 cM (SSC15) and 39.8 cM (SSC17), and, on SSC2, at 16.9 cM, in families derived from F1 males only (p < 0.05, except for SSC13, p < 0.01). Analysis of 7 precise specific traits revealed highly significant QTLs on SSC10 (ulceration), on SSC12 (presence of melanoma at birth), on SSC13 (lesion type), and on SSC16 and SSC17 (number of aggressive melanomas) at the respective positions 42.0, 95.6, 81.0, 45.3 and 44.8 cM (p < 0.001 and p < 0.05 respectively at the chromosome- and genome-wide levels). We also showed that MeLiM MC1R*2 allele, which determines black coat colour in pigs, predisposes significantly to melanoma. Interactions were observed between MC1R and markers located on SSC1 (p < 0.05). Taken together, these results indicate that MeLiM swine is a model for human multigenic diseases. Comparative mapping revealed human regions of interest to search for new melanoma susceptibility candidates.

Fatness QTL on chicken chromosome 5 and interaction with sex.

Quantitative trait loci (QTL) affecting fatness in male chickens were previously identified on chromosome 5 (GGA5) in a three-generation design derived from two experimental chicken lines divergently selected for abdominal fat weight. A new design, established from the same pure lines, produced 407 F2 progenies (males and females) from 4 F1-sire families. Body weight and abdominal fat were measured on the F2 at 9 wk of age. In each sire family, selective genotyping was carried out for 48 extreme individuals for abdominal fat using seven microsatellite markers from GGA5. QTL analyses confirmed the presence of QTL for fatness on GGA5 and identified a QTL by sex interaction. By crossing one F1 sire heterozygous at the QTL with lean line dams, three recombinant backcross 1 (BC1) males were produced and their QTL genotypes were assessed in backcross 2 (BC2) progenies. These results confirmed the QTL by sex interaction identified in the F2 generation and they allow mapping of the female QTL to less than 8 Mb at the distal part of the GGA5. They also indicate that fat QTL alleles were segregating in both fat and lean lines.