The chicken RH map: current state of progress and microchromosome mapping

Mapping of leptin and its syntenic genes to chicken chromosome 1p

Background

Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly.

Results

A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = − 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes.

Conclusion

Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0543-1) contains supplementary material, which is available to authorized users.

Expanding Duplication of Free Fatty Acid Receptor-2 (GPR43) Genes in the Chicken Genome

Free fatty acid receptors (FFAR) belong to a family of five G-protein coupled receptors that are involved in the regulation of lipid metabolism, so that their loss of function increases the risk of obesity. The aim of this study was to determine the expansion of genes encoding paralogs of FFAR2 in the chicken, considered as a model organism for developmental biology and biomedical research. By estimating the gene copy number using quantitative polymerase chain reaction, genomic DNA resequencing, and RNA sequencing data, we showed the existence of 23 ± 1.5 genes encoding FFAR2 paralogs in the chicken genome. The FFAR2 paralogs shared an identity from 87.2% up to 99%. Extensive gene conversion was responsible for this high degree of sequence similarities between these genes, and this concerned especially the four amino acids known to be critical for ligand binding. Moreover, elevated nonsynonymous/synonymous substitution ratios on some amino acids within or in close-vicinity of the ligand-binding groove suggest that positive selection may have reduced the effective rate of gene conversion in this region, thus contributing to diversify the function of some FFAR2 paralogs. All the FFAR2 paralogs were located on a microchromosome in a same linkage group. FFAR2 genes were expressed in different tissues and cells such as spleen, peripheral blood mononuclear cells, abdominal adipose tissue, intestine, and lung, with the highest rate of expression in testis. Further investigations are needed to determine whether these chicken-specific events along evolution are the consequence of domestication and may play a role in regulating lipid metabolism in this species.

Evolutionary genomics and adaptive evolution of the Hedgehog gene family (Shh, Ihh and Dhh) in vertebrates

The Hedgehog (Hh) gene family codes for a class of secreted proteins composed of two active domains that act as signalling molecules during embryo development, namely for the development of the nervous and skeletal systems and the formation of the testis cord. While only one Hh gene is found typically in invertebrate genomes, most vertebrates species have three (Sonic hedgehog – Shh; Indian hedgehog – Ihh; and Desert hedgehog – Dhh), each with different expression patterns and functions, which likely helped promote the increasing complexity of vertebrates and their successful diversification. In this study, we used comparative genomic and adaptive evolutionary analyses to characterize the evolution of the Hh genes in vertebrates following the two major whole genome duplication (WGD) events. To overcome the lack of Hh-coding sequences on avian publicly available databases, we used an extensive dataset of 45 avian and three non-avian reptilian genomes to show that birds have all three Hh paralogs. We find suggestions that following the WGD events, vertebrate Hh paralogous genes evolved independently within similar linkage groups and under different evolutionary rates, especially within the catalytic domain. The structural regions around the ion-binding site were identified to be under positive selection in the signaling domain. These findings contrast with those observed in invertebrates, where different lineages that experienced gene duplication retained similar selective constraints in the Hh orthologs. Our results provide new insights on the evolutionary history of the Hh gene family, the functional roles of these paralogs in vertebrate species, and on the location of mutational hotspots.

Mapping and genotypic analysis of the NK-lysin gene in chicken

Background

Antimicrobial peptides (AMP) are important elements of the first line of defence against pathogens in animals. NK-lysin is a cationic AMP that plays a critical role in innate immunity. The chicken NK-lysin gene has been cloned and its antimicrobial and anticancer activity has been described but its location in the chicken genome remains unknown. Here, we mapped the NK-lysin gene and examined the distribution of a functionally significant single nucleotide polymorphism (SNP) among different chicken inbred lines and heritage breeds.

Results

A 6000 rad radiation hybrid panel (ChickRH6) was used to map the NK-lysin gene to the distal end of chromosome 22. Two additional genes, the adipocyte enhancer-binding protein 1-like gene (AEBP1) and the DNA polymerase delta subunit 2-like (POLD2) gene, are located in the same NW_003779909 contig as NK-lysin, and were thus indirectly mapped to chromosome 22 as well. Previously, we reported a functionally significant SNP at position 271 of the NK-lysin coding sequence in two different chicken breeds. Here, we examined this SNP and found that the A allele appears to be more common than the G allele in these heritage breeds and inbred lines.

Conclusions

The chicken NK-lysin gene mapped to the distal end of chromosome 22. Two additional genes, AEBP1 and POLD2, were indirectly mapped to chromosome 22 also. SNP analyses revealed that the A allele, which encodes a peptide with a higher antimicrobial activity, is more common than the G allele in our tested inbred lines and heritage breeds.

Mapping genes to chicken microchromosome 16 and discovery of olfactory and scavenger receptor genes near the major histocompatibility complex

Trisomy mapping is a powerful method for assigning genes to chicken microchromosome 16 (GGA 16). The single chicken nucleolar organizer region (NOR), the 2 major histocompatibility complex regions (MHC-Y and MHC-B), and CD1 genes were all previously assigned to GGA 16 using trisomy mapping. Here, we combined array comparative genomic hybridization with trisomy mapping to screen unassigned genomic scaffolds (consigned temporarily to chrUn_random) for sequences originating from GGA 16. A number of scaffolds mapped to GGA 16. Among these were scaffolds that contain genes for olfactory (OR) and cysteine-rich domain scavenger (SRCR) receptors, along with a number of genes that encode putative immunoglobulin-like receptors and other molecules. We used high-resolution cytogenomic analyses to confirm assignment of OR and SRCR genes to GGA 16 and to pinpoint members of these gene families to the q-arm in partially overlapping regions between the centromere and the NOR. Southern blots revealed sequence polymorphism within the OR/SRCR region and linkage with the MHC-Y region, thereby providing evidence for conserved linkage between OR genes and the MHC within birds. This work localizes OR genes to the vicinity of the chicken MHC and assigns additional genes, including immune defense genes, to GGA 16.

A duck RH panel and its potential for assisting NGS genome assembly

Background

Owing to the low cost of the high throughput Next Generation Sequencing (NGS) technology, more and more species have been and will be sequenced. However, de novo assemblies of large eukaryotic genomes thus produced are composed of a large number of contigs and scaffolds of medium to small size, having no chromosomal assignment. Radiation hybrid (RH) mapping is a powerful tool for building whole genome maps and has been used for several animal species, to help assign sequence scaffolds to chromosomes and determining their order.

Results

We report here a duck whole genome RH panel obtained by fusing female duck embryonic fibroblasts irradiated at a dose of 6,000 rads, with HPRT-deficient Wg3hCl₂ hamster cells. The ninety best hybrids, having an average retention of 23.6% of the duck genome, were selected for the final panel. To allow the genotyping of large numbers of markers, as required for whole genome mapping, without having to cultivate the hybrid clones on a large scale, three different methods involving Whole Genome Amplification (WGA) and/or scaling down PCR volumes by using the Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM for genotyping were tested. RH maps of APL12 and APL22 were built, allowing the detection of intrachromosomal rearrangements when compared to chicken. Finally, the panel proved useful for checking the assembly of sequence scaffolds and for mapping EST located on one of the smallest microchromosomes.

Conclusion

The Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM genotyping by quantitative PCR provides a rapid and cost-effective method for building RH linkage groups. Although the vast majority of genotyped markers exhibited a picture coherent with their associated scaffolds, a few of them were discordant, pinpointing potential assembly errors. Comparative mapping with chicken chromosomes GGA21 and GGA11 allowed the detection of the first chromosome rearrangements on microchromosomes between duck and chicken. As in chicken, the smallest duck microchromosomes appear missing in the assembly and more EST data will be needed for mapping them. Altogether, this underlines the added value of RH mapping to improve genome assemblies.

An unbiased approach to identify genes involved in development in a turtle with temperature-dependent sex determination

BACKGROUND

Many reptiles exhibit temperature-dependent sex determination (TSD). The initial cue in TSD is incubation temperature, unlike genotypic sex determination (GSD) where it is determined by the presence of specific alleles (or genetic loci). We used patterns of gene expression to identify candidates for genes with a role in TSD and other developmental processes without making a priori assumptions about the identity of these genes (ortholog-based approach). We identified genes with sexually dimorphic mRNA accumulation during the temperature sensitive period of development in the Red-eared slider turtle (Trachemys scripta), a turtle with TSD. Genes with differential mRNA accumulation in response to estrogen (estradiol-17β; E(2)) exposure and developmental stages were also identified.

RESULTS

Sequencing 767 clones from three suppression-subtractive hybridization libraries yielded a total of 581 unique sequences. Screening a macroarray with a subset of those sequences revealed a total of 26 genes that exhibited differential mRNA accumulation: 16 female biased and 10 male biased. Additional analyses revealed that C16ORF62 (an unknown gene) and MALAT1 (a long noncoding RNA) exhibited increased mRNA accumulation at the male producing temperature relative to the female producing temperature during embryonic sexual development. Finally, we identified four genes (C16ORF62, CCT3, MMP2, and NFIB) that exhibited a stage effect and five genes (C16ORF62, CCT3, MMP2, NFIB and NOTCH2) showed a response to E(2) exposure.

CONCLUSIONS

Here we report a survey of genes identified using patterns of mRNA accumulation during embryonic development in a turtle with TSD. Many previous studies have focused on examining the turtle orthologs of genes involved in mammalian development. Although valuable, the limitations of this approach are exemplified by our identification of two genes (MALAT1 and C16ORF62) that are sexually dimorphic during embryonic development. MALAT1 is a noncoding RNA that has not been implicated in sexual differentiation in other vertebrates and C16ORF62 has an unknown function. Our results revealed genes that are candidates for having roles in turtle embryonic development, including TSD, and highlight the need to expand our search parameters beyond protein-coding genes.

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.

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.