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

C. Meslin, Colette Desert, I. Callebaut, A. Djari, C. Kropp, Frederique Pitel, S. Leroux, P. Martin, P. Froment, E. Guilbert, et al.

► To cite this version:

C. Meslin, Colette Desert, I. Callebaut, A. Djari, C. Kropp, et al. Expanding Duplication of Free Fatty Acid Receptor-2 (GPR43) Genes in the Chicken Genome. Genome Biology and Evolution, 2015, 7 (5), pp.1332-1348. 10.1093/gbe/evv072. hal-01157691

HAL Id: hal-01157691

https://hal.science/hal-01157691

Submitted on 4 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Expanding Duplication of Free Fatty Acid Receptor-2 (GPR43) Genes in the Chicken Genome

Camille Meslin<sup>1,2,3,4,†</sup>, Colette Desert<sup>5,6,†</sup>, Isabelle Callebaut<sup>7</sup>, Anis Djari<sup>8</sup>, Christophe Klopp<sup>8</sup>, Frédérique Pitel<sup>9</sup>, Sophie Leroux<sup>9</sup>, Pascal Martin<sup>10</sup>, Pascal Froment<sup>1,2,3,4</sup>, Edith Guilbert<sup>1,2,3,4</sup>, Florence Gondret<sup>5,6</sup>, Sandrine Lagarrigue<sup>5,6</sup>, and Philippe Monget<sup>1,2,3,4,*</sup>

Accepted: April 22, 2015

## Abstract

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 non-synonymous/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.

**Key words:** FFAR, evolution, gene conversion, positive selection, chicken, duplication.

## Introduction

Free fatty acid receptors (FFAR) are members of the G-protein-coupled receptor (GPCR) superfamily and are activated by free fatty acids (FFAs). Five receptors of this subfamily have been yet identified: FFAR1 (GPR40), FFAR2 (GPR43), FFAR3

(GPR41), FFAR4 (GPR120), and GPR84. They are characterized by their respective ligands, their pattern of expression, and their biological functions. GPR84 and FFAR4 are activated by medium-chain and unsaturated long-chain FFAs, respectively, whereas FFAR1 is activated by both medium- and long-chain

© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals permissions@oup.com

<footnote n="1"><sup>1</sup>UMR85 Physiologie de la Reproduction et des Comportements, INRA, Nouzilly, France</footnote>

<footnote n="2"><sup>2</sup>UMR7247, CNRS, Nouzilly, France</footnote>

<footnote n="3"><sup>3</sup>Université François Rabelais de Tours, France</footnote>

<footnote n="4"><sup>4</sup>Institut Français du Cheval et de l'Equitation, Nouzilly, France</footnote>

<footnote n="5"><sup>5</sup>INRA, UMR1348 Physiologie, Environnement et Génétique pour l'Animal et les Systèmes d'élevage, Saint-Gilles, France</footnote>

<footnote n="6"><sup>6</sup>Agrocampus-Ouest, UMR1348, Rennes, France</footnote>

<footnote n="7"><sup>7</sup>IMPMC, UMR CNRS 7590, Museum National d'Histoire Naturelle, IRD UMR 206, Sorbonne Universités—UPMC Université Paris 06, France</footnote>

<footnote n="8"><sup>8</sup>INRA, BIA, CS 52627, Castanet-Tolosan, France</footnote>

<footnote n="9"><sup>9</sup>UMR INRA/<sup>h2</sup>INPT ENSAT/<sup>h3</sup>INPT ENVT, UMR1388 Génétique, Physiologie et Systèmes d'élevage, INRA, Castanet Tolosan, France</footnote>

<footnote n="10"><sup>10</sup>INRA, UR 0066 Pharmacologie-Toxicologie, Toulouse, France</footnote>

<footnote>*Corresponding author: E-mail: pmonget@tours.inra.fr.</footnote>

<footnote>†These authors contributed equally to this work.</footnote>
FFAs. In contrast, both FFAR2 and FFAR3 are selectively activated by short-chain FFAs (SCFAs) from one to six carbon chain length (Ulven 2012; Hudson et al. 2013). The pattern of expression of these receptors also differs, as previously shown in human and rodent species (Stoddart et al. 2008; Ulven 2012). FFAR1 expression has been mainly reported in pancreas (notably in β-cells), FFAR3 expression has been observed in many tissues with the highest level in white adipose tissue, and GPR84 is predominantly expressed in hematopoietic tissues and bone marrow (Wang et al. 2006). FFAR4 is widely expressed in various tissues and cell types including intestine, macrophages, adipose tissue, taste buds, brain, pancreas, lung, thymus, and pituitary (Ichimura et al. 2009). Finally, FFAR2 is highly expressed in immune cells such as neutrophils, monocytes and peripheral blood mononuclear cells (PBMC), but has been also detected in bone marrow, spleen, skeletal muscle, heart, adipose tissue, and intestine (Stoddart et al. 2008). Bovine FFAR2 was found in almost all tissues tested (Wang et al. 2009), with variations in adipose tissue according to the age of animals (Smith et al. 2012). Pig FFAR2 was also detected in adipose tissue (Hou and Sun 2008) and intestine (Colombo et al. 2012). To our knowledge, only one study reported the expression of a chicken FFAR, claimed to be FFAR1 in chicken hepatocytes in vitro (Suh et al. 2008).

Functional studies have highlighted different roles of FFAR family in human health. In particular, human loss-of-function variants of mouse FFAR2 and human FFAR4 have been shown to increase the risk to develop obesity (Ichimura et al. 2012; Kimura et al. 2013). FFAR2 deficiency protects from high-fat diet-induced obesity and dyslipidemia, at least partly through increased energy expenditure (Bjursell et al. 2011). Mice over-expressing FFAR2 specifically in adipose tissue remain lean even when fed a high-fat diet. Furthermore, SCFAs have been described as key molecules produced by gut microbial fermentation of soluble fibers and the activation of FFAR2 by SCFAs has been involved in the regulation of energy balance (Kimura et al. 2013). Taken together, it is suggested that FFAR2 may play a key role in lipid metabolism, glucose tolerance, immune regulation, and may be involved in the crosstalk between gut microbiota and whole-body homeostasis.

Beside rodents, other animal organisms are now recognized for their potential interest in a better understanding of developmental biology, physiology, and human diseases. The chicken was the first avian species and domestic animal selected for complete genome sequencing and assembly. Chicken exhibit "natural" hyperglycemia but no signs of insulin resistance (Simon 1989), making them a valuable model to understand the regulation of energy homeostasis (Burt 2007). Like humans, de novo synthesis of lipids occurs mainly in the liver, whereas fat is deposited mainly at the visceral (abdominal) location. As stated above, only one study has reported the expression in vitro of a chicken FFAR to date,

which was claimed to be FFAR1 in chicken hepatocytes (Suh et al. 2008). Twenty-six chicken genes encoding paralogs of FFAR2 were accessible in old versions of Ensembl database (Ensembl release 70) but have been removed in the current version (Ensembl release 77).

In this article, we experimentally confirmed the existence of FFAR2 paralogs, examined their patterns of expression in different tissues, and studied their evolution by gene conversion and positive selection. Comparisons were made with pigs, where FFAR2 has been previously detected in adipose tissue (Hou and Sun 2008) and intestine (Colombo et al. 2012).

# Materials and Methods

## Animals and Experimental Procedures

### Chicken

Broilers were reared together in a closed building at the experimental unit PEAT (Pôle d'Expérimentation Avicole de Tours, INRA) under standard condition. They were fed ad libitum using conventional diet for a minimum of 4 h after overnight fasting and then weighed and sacrificed by electrical stunning in the experimental processing plant. Following sacrifice, blood was collected from all animals for DNA extraction and tissues were collected, quickly frozen in liquid nitrogen and stored at -80 °C until RNA extraction. Peripheral blood mononuclear cells (PBMC) were isolated from the blood using Ficoll Histopaque (Sigma, #10771) density gradient centrifugation according to standard protocols, before liquid nitrogen congelation.

In order to check the link between FFAR2 expression and adiposity, another animal protocol was used: Chickens from two experimental broiler lines divergently selected for abdominal fat (Leclercq et al. 1980) were fed with a conventional diet and slaughtered at 9 weeks of age for the reverse transcriptase (RT)-qPCR measure and at 14 weeks of age for the RNA-Seq study.

Concerning experiments on testis, seminiferous tubules (containing germ cells and Sertoli cells) and Sertoli cells preparations were obtained from 6-week-old (immature animals) chickens (ISA Brown, Institut de Sélection Animale, Saint Brieuc, France) and were adapted and modified from other methodologies published in rodents species (Ellingson and Yao 1970; Foster et al. 1984; Toebosch et al. 1989; Zwain and Cheng 1994; Staub et al. 2000; Guibert et al. 2011). At this age, the population of germ cells detectable is spermatogonia and type I spermatocytes. Testes were decapsulated, slightly minced, and incubated in a shaking water bath during 15 min at 37°C with DNase I and collagenase (Sigma, l'Isle d'Abeau Chesnes, France). Cells were centrifuged to remove collagenase and then cells were allowed to sediment by gravity in order to separate seminiferous tubules, in the pellet, and Leydig cells, in the supernatant. The fragments of seminiferous tubules were slightly digested by a last
collagenase bath. Chicken tubules were frozen at -80°C. Then, Sertoli cells preparations were obtained by two successive digestion of seminiferous tubules (two collagenase bathes 0.6 and 0.8 mg/ml for 15 min at 37°C each digestion) added with DNase (20 µg/ml) and followed by 0.1% hyaluronidase treatment (Sigma, l'Isle d'Abeau Chesnes, France) for 10 min at 37°C, to reduce peritubular cell contamination. The purity of the Sertoli cell preparation was averaged 80-85%. Indeed, contamination of Sertoli cell preparations with germ cells was less than 10% and the percentage of peritubular myoid cells, evaluated by alkaline phosphatase staining, was close to 8% of the total cell population (Guibert et al. 2011). Sertoli cells were frozen at -80°C.

DNA samples from wild jungle fowls from Thailand (*Gallus gallus* and *Gallus spadiceus*) (one individual per breed) were also kindly provided by Michèle Tixier-Boichard from UMR GABI, INRA, AgroParisTech, 78352 Jouy-en-Josas.

All experiments were conducted under Licence No. 37-123 from the Veterinary Services (Indre et Loire, France) and in accordance with guidelines for care and use of animals in Agricultural Research and Teaching (French Agricultural Agency and Scientific Research Agency).

### Pig

Male (n = 4) and female (n = 4) pigs of a crossbred genotype (Pietrain × (Large White × Landrace)) were considered. From 40 kg body weight to slaughter, pigs were reared in isolated cages in a temperature-controlled room, and fed ad libitum a standard cereal-based diet. At 98.2 ± 0.9 kg body weight (i.e., 151 ± 12 days of age), pigs were killed 2 h after their last meal intake, by exsanguination following electronarcosis. The liver and perirenal adipose tissue (visceral location) were taken off and weighed. Samples were immediately prepared, frozen in liquid nitrogen, and stored at -76°C until analysis.

### Gene Dosage Analysis by qPCR

Variable amounts of chicken genomic DNA (6.25, 3.125, 1.563, 0.781, and 0.391 ng/µl) were used to quantify FFAR2 copy number in a 96-well plate. PCR primers were determined to amplify all the chicken FFAR2 sequences and were forward (F)-GCCCCATAGCAAACTTCT and reverse (R)-GGGCAGCCAT AAAGAGAG. As a control, four single-copy genes were amplify: FADS1, ACOX1, CPT1A, and ETS1. The primers were FADS1_F-CAGCACCACGCGAAACC, FADS1_R-TCTACAGAG AGCTTCTTTCCCAAAG, ACOX1_F-TCATCCGGTCTCTGATTG TAGGA, ACOX1_R-GCACTATAGCGGATGGCAATG, CPT1A_F-CCCTGAAAATGCTGCTTTCCTA, CPT1A_R-TGGTG CCTGCAGAAAGTTTG, ETS1_F-CAGCATCAGCACAAAGCA G, and ETS1_R-CAGCCAACCCAACCAAAG. Quantitative real-time PCR was performed using SsoFast EvaGreen Supermix on a CFX96 Touch Real-Time PCR system (Bio-

Rad) in 15 µl final volume (comprising 1 × SsoFast EvaGreen supermix, 0.3 µM forward and reverse primers, and 5 µl of DNA). PCR cycling program consisted of 5 min at 98°C, followed by 44 cycles of 5 s at 98°C and 10 s at 59°C. An additional step was used (from 65 to 95°C during 10 s, with 0.5°C increment) for dissociation curve analysis. A same fluorescence threshold was fixed for all the plate. For each PCR, post-PCR melting curves confirmed the specificity of single-target amplification. The amplification efficiencies were similar for the five reactions, so the $2^{-\Delta\Delta Cq}$ method (Livak and Schmittgen 2001) could be used with the average Cq of the four control genes as calibrator. Cq is defined as the number of cycles required for the fluorescent signal to cross the threshold, which is when the signal exceeds the background level. $\Delta Cq$ then shows the difference of expression between two genes.

### RNA Isolation

For all tissue samples except the pancreas, total RNA was extracted with TRIzol reagent (Invitrogen, Cergy Pontoise, France) according to the manufacturer's instructions. Pancreatic total RNA was prepared by the guanidinium thio-cyanate extraction procedure (Chirgwin et al. 1979). The quantity and quality of RNA extracted were evaluated by 1% agarose gel electrophoresis and spectrophotometrically using a nanodrop measuring system. No sign of RNA degradation was observed. The TURBO DNA-free kit (Ambion) was used for the DNase treatment.

### Real-Time RT PCR

One microgram of total RNA was reverse-transcribed using iScript cDNA Synthesis Kit (Bio-Rad). cDNA templates (diluted in a 1:10 ratio) were used to perform PCR quantification (utilizing conditions described above). For the measurement of tissue expression, three biological replicates were used. Data were normalized to levels of 18 S rRNA (primers sequences: F-TTAAGTCCCTGCCCTTTGTACAC; R-CGATCCGAGGACCTCA CTAAAC). Average Cq of the samples with weak expression was used as calibrator. Specific primers able to target five FFAR2 paralogs were F-CTCTTTATGGCTGCCCTCAG and R-GTAGCCCAGGCTTGGTTGG. They amplified the genes ENSGALG00000022595, ENSGALG00000022164, ENSGALG00000022608, ENSGALG00000022422, and ENSGALG0000022399 (supplementary data S1, Supplementary Material online).

In tissue-specific expression measure, reference genes were either GAPDH for FFAR2 quantification in adipose tissue of fat and lean lines (primers sequences: F-GCTAAGGCTGTGGGGA AAGT; R-TCAGCAGCAGCCTTCACTAC) or RPL15 (Ribosomal Protein L15) for expression measure in Sertoli cells, seminiferous tubules, and testis (F-TGTGATGCGTTTCCTCCTTGG; R-CC ATAGGTTGCACCTTTTGGG).
## RH Mapping

PCR amplifications were carried out for each marker (FFAR2: primers used above, MAG: F-GGCAGCCC ATACCCTAAAAG, R-GCTCGCTGAAGCTGTACTGG, USF2: F-GCGCAGGAGGGATAAAATC, R-CCCCGTCTTG CTGTTGTC) in 15 µl reactions containing 25 ng DNA from the chickRH6 panel (Morisson et al. 2002), 0.4 µM of each primer, 0.25 units Taq polymerase (GoTaq, Promega), 1.5 mM MgCl2, and 0.2 mM dNTP on a GeneAmp PCR System 9700 thermocycler (Applied Biosystems). The first 5-min denaturation step was followed by 35 cycles, each consisting of denaturation at 94°C for 30 s, annealing at Tm for 30 s and elongation at 72°C for 30 s. PCR products were analyzed on 2% agarose gels, electrophoresed in 1× TBE (Tris/Borate/EDTA) buffer, and visualized by staining with ethidium bromide. Mapping of the markers on the RH panel was performed through the ChickRH server (http://chickrh.toulouse.inra.fr/, last accessed May 7, 2015). Distances and two-point LODs (logarithm of odds) were calculated through the Carthagene software (De Givry et al. 2005). Maps were drawn with MapChart 2.0 (Voorrips 2002). FFAR2 fragment from each positive hybrid was sequenced on an ABI 3730 sequencer (Applied Biosystem).

## Genome DNA Resequencing Data

Genome DNA resequencing was performed from 18 birds belonging to European broiler lines. The 18 DNA-Seq libraries were prepared using the TruSeq DNA Sample Preparation Kit (Illumina, San Diego, CA) according to the manufacturer's instructions. Briefly, paired-end libraries with a 250 bp insert size were generated. The libraries were quantified using QPCR Library Quantification Kit (Agilent), checked on a High Sensitivity DNA Chip (Agilent), and sequenced in pair end 2 × 100 bp on a HiSeq 2000 (Illumina) using a TruSeq v3 Kit. Sequencing produced per animal an average 20.4 Gb corresponding to a sequencing depth of 19.7 ×. For each animal, reads were aligned against the 26 gene sequences encoding paralogs of FFAR2 identified in 2011 by Ensembl, using BWA software (Li and Durbin 2009). Only reads with a unique mapping hit and a minimal quality score of 30 were kept.

## RNA-Seq Data

RNA-Seq was performed on liver and adipose tissue mRNA of 14-week-old prepubertal chicken (n = 8: four males and four females; two issued from fat lines and two from lean lines per sex) and pig (4.5 months) (four males and four females) (SRA accession number SRP042257). RNA-Seq data from 4.5-day chick embryos (n = 20) were obtained from a previous study (SRA accession number SRP033603; Frésard et al. 2014).

The 32 mRNA-Seq libraries were prepared using the TruSeq RNA Sample Preparation Kit (Illumina) according to the manufacturer's instructions. Briefly, Poly-A RNA was purified from 4 µg of total RNA using oligo(dT) magnetic beads, fragmented into 150-400 bp pieces using divalent cations at 94°C for 8 min, and retrotranscribed to double-stranded cDNA using random primers. The resulting cDNA was purified using Agencourt AMPure XP beads (Beckman Coulter, Villepinte, France). Then, cDNA was subjected to end-repair and phosphorylation and subsequent purification was performed using Agencourt AMPure XP beads (Beckman Coulter). These repaired cDNA fragments were 3'-adenylated producing cDNA fragments with a single "A" base overhang at their 3'-ends for subsequent adapter-ligation. Illumina adapters containing indexing tags were ligated to the ends of these 3'-adenylated cDNA fragments followed by two purification steps using Agencourt AMPure XP beads (Beckman Coulter). Ten rounds of PCR amplification were performed to enrich the adapter-modified cDNA library using primers complementary to the ends of the adapters. The PCR products were purified using Agencourt AMPure XP beads (Beckman Coulter) and size-selected (200 ± 25 bp) on a 2% agarose Invitrogen EGel (Thermo Scientific). Libraries were then checked on an Agilent Technologies 2100 Bioanalyzer using the Agilent High Sensitivity DNA Kit and quantified by qPCR with the QPCR NGS Library Quantification kit (Agilent Technologies). After quantification, tagged cDNA libraries were pooled in equal ratios and a final qPCR check was performed postpooling. The pooled libraries were used for 2 × 100 bp paired-end sequencing on one lane of the Illumina HiSeq2000 with a TruSeq SBS v3-HS Kit (Illumina). After sequencing, the samples were demultiplexed and the indexed adapter sequences were trimmed using the CASAVA v1.8.2 software (Illumina). Sequencing produced an average of 6.67 Gb (i.e., 33 million of pair end reads) per animal. Reads were then mapped to the chicken reference transcript database (Ensembl cDNA release 64) using BWA software (bwa-0.5.9rc1; Li and Durbin 2009). The number of reads mapped to each transcript was determined using samtools (0.1.17) and an in house script (Li et al. 2009).

## Gene Conversion

The 26 genes of the FFAR2 cluster used for the positive selection calculation were tested for intralocus gene conversion by using GENECONV version 1.81 (Sawyer 1989), which is a well-established method for detecting partial gene conversion (Posada 2002). The analysis was performed on the codon-based alignment performed by the Muscle algorithm (Edgar 2004) in the PhyleasProg pipeline (Busset et al. 2011) and permitted to the program to look for pair of sequences which are sufficiently similar to be suggestive of nonreciprocal transfer of genetic information between the different sequences. Three *P* values
were calculated and compared to assess the significance of the results. Strong evidence of gene conversion was noted when a fragment had a *P* value that was less than 0.05 for at least two different types of statistical test. In each alignment, indels and missing data were ignored. All polymorphic sites were tested for evidence of gene conversion using adjusted mismatch penalties of 0 (g0) to detect recent gene conversion events. The figure for gene conversion events detected with the g0 option was made using Circos (Krzywinski et al. 2009).

## Positive Selection Calculation

The inference of positive selection was performed on the tree of paralogous genes by PhyleasProg with branch-site and site models of codeml of the PAML package, with the fast option (Busset et al. 2011). The multiple sequence alignment was carefully examined to avoid all false positive results. In particular, amino acids predicted to be under positive selection that were at the boundary of the alignments were not considered because they are doubtful. Both branch-site and site models are designed to identify amino acids under positive selection. However, the site model allows the ω ratio to vary among sites (among amino acids in the protein). The branch-site model on the other side allows ω to vary both among sites in the protein and across branches on the tree and therefore aims to detect positive selection affecting a few sites along particular lineages, sites that would not be detected by using the site model. For the use of branch-site models, each branch of the phylogenetic tree was tested for positive selection. So we performed multiple test corrections by controlling for the false discovery rate using the R package QVALUE (Storey and Tibshirani 2003). Results are considered significant with a threshold of q=10% of false positives. Sites with posterior probabilities of Bayes Empirical Bayes analyses superior to 95% or 99% were considered as positively selected. No overlap was found between the two models because the branch-site model is detecting positive selection on a selected branch whereas the site model is detecting positive selection affecting the whole phylogenetic tree.

## 3D Structure Modeling

A model of the 3D structure of FFAR2 has been constructed using Modeller (v9.10) (Martí-Renom et al. 2000), the alignment shown in *supplementary data S7*, *Supplementary Material online*, deduced from a HHpred search (Söding et al. 2005) and, as template, the 3D structure of proteinase-activated receptor 1 (pdb 3vw; Zhang et al. 2012). The local reliability was evaluated using VERIFY3D (http://services.mbi.ucla.edu/Verify_3D/, last accessed May 6, 2015; Bowie et al. 1991), whereas the stereochemical quality was evaluated with PROCHECK (Laskowski et al. 1993), 3D structures were manipulated

using Chimera (Pettersen et al. 2004), and the sequence alignment was presented using ESPript (Robert and Gouet 2014).

## Results

### FFAR2 Gene Was Highly Duplicated in the Chicken Genome

In the version 70 of Ensembl database, 26 chicken genes encoding paralogs of FFAR2 were accessible (Ensembl accession numbers are provided in *supplementary data S1*, *Supplementary Material online*). To check the existence of these multiple FFAR2 paralogs in the chicken genome, complementary analyses were performed in this study.

### Validation by Quantitative Polymerase Chain Reaction

We first performed real-time polymerase chain reaction (PCR) using series of genomic DNA dilutions obtained from two French commercial broiler chickens. Primer pairs were designed to amplify all chicken FFAR2 sequences, so called "universal primers" of FFAR2. Four single copy genes were used as controls. The differences of quantitative cycle (ΔCq values) between the single-copy genes and FFAR2 genes were found between 4.27 and 4.70 (fig. 1). These results confirm the massive duplication of FFAR2 genes in the chicken genome, and estimate a relative copy number of 23.7 ± 1.2 paralogs (mean ± SEM; *n* = 5). A similar experiment was also performed using genomic DNA dilutions from two ancestral chicken lines from Thailand. For these populations, we also found a massive duplication of FFAR2 (fig. 1) with a copy number close-similar to that found in the European domestic line (relative copy number of 23.0 ± 1.5 (mean ± SEM; *n* = 4).

### Validation by Genomic DNA Resequencing

Whole-genome resequencing data were generated from 18 chickens with a 20× sequencing depth per individual. For each animal, reads were then mapped against the 26 gene sequences encoding paralogs of FFAR2 identified by Ensembl, v70 (*supplementary data S1*, *Supplementary Material online*), providing that only reads which had a unique hit after mapping were kept. For each animal, we were able to map reads to all 26 gene sequences, thus confirming existence of FFAR2 paralogs in the chicken genome. The percentages of identity between FFAR2 nucleotide sequences were found to range from 87.2% to 99.0% (*supplementary data S3*, *Supplementary Material online*). Examination of pairwise alignments of nucleotide sequences of FFAR2 paralogs that shared an identity rate above 98.5% (*supplementary data S3*, *Supplementary Material online*, in red) allowed detecting from 9.8 to 11.9 single nucleotide polymorphisms (SNP)/kb. Importantly, Basic Local Alignment Search Tool (BLAST)
<figure><img src="image_2.png" /><figcaption>**Fig. 1.**—FFAR2 has numerous copies in European broiler and in ancestral chicken genome. qPCR on genomic DNA shows clear difference in ΔCq between FFAR2 genes and the genes carrying only one copy per genome (four genes in European lines and three genes in Ancestral lines). qPCR was performed using "universal" primers able to amplify the 26 sequences of FFAR2 (see supplementary data S1, Supplementary Material online). FFAR2 error bars indicate standard deviation between two individual chickens. For the single copy genes, error bars represent standard deviation between the Cq measures for the three or four genes from two individual chickens.</figcaption></figure>

searches against turkey, finch, and quail genome and expressed sequence tags (ESTs) database did not detect any duplication of FFAR2 genes.

### Validation by RNA Sequencing Data in Whole Embryos and Prepupertal Tissues

The availability of RNA sequencing (RNA-Seq) data from whole body or specific tissues can be valuable to prove expression of particular genes. We produced on average 47 million reads per chicken whole embryos and 20 and 16 million reads for liver and adipose tissue, respectively, per prepupertal chickens. Reads were then mapped to the chicken reference transcript database (Ensembl cDNA release 64). Among the 26 Ensembl FFAR2 paralogs, 11 were observed with at least 10 reads in each tissue (supplementary data S2, Supplementary Material online). These results confirm existence of many expressed FFAR2 paralogs in chicken. It remains however to determine whether other paralogs could be detected when other tissues or other physiological conditions are considered.

### Chicken FFAR2 Paralogs Have Highly Similar Sequences and Are Located in the Same Genomic Region

#### Amino Acid and Protein Sequences

The percentages of identity shared by the chicken FFAR2 amino acid sequences reached 83.1–98.6% (table 1), these identities being much greater than for other more distant chicken paralogs of the same (F2R, F2RI1, F2RL2, F2RL3, P2RY8, from 21.2% to 49.3%). FFAR1 and FFAR3 genes

were not found in the chicken genome. In order to compare with other species, percentages of identity analysis were determined on nucleic and protein sequences of FFAR paralogs (FFAR1, FFAR2, FFAR3, FFAR4, and GPR84) in human, mouse, and pig species (supplementary data S4–S6, Supplementary Material online). The percentages of identity between FFAR1, -2, -3, and -4 were below 58% and 48% for nucleic and protein sequences, respectively (in green in supplementary data S5 and S6, Supplementary Material online). Moreover, these different FFAR genes are shared by most of mammals, suggesting that they have duplicated before divergence between these species, and that the duplications that led to FFAR2 chicken paralogs are more recent.

The construction of the phylogenetic tree of the 26 FFAR2 amino acid sequences in chicken confirmed that they belonged to the FFAR2 subclass of FFAR encoding genes (fig. 2).

Alignment of the chicken protein sequences of FFAR2 (supplementary data S7, Supplementary Material online) showed that the two arginine and two histidine, previously described as central for the binding of fatty acid ligands (Sum et al. 2007; Stoddart et al. 2008; Hudson et al. 2013), are well conserved in chicken FFAR2, as in human FFAR2 but also human FFAR1 and FFAR3 sequences (fig. 3). The 3D fold is well conserved between FFAR1, FFAR2 and FFAR3, as highlighted by the overall conservation of amino acids over the whole alignment (fig. 3).

Hydrophobic amino acids, which make up the core of the fold, are especially well conserved. The first arginine residue (Arg¹⁸⁰ in human FFAR2) was found in the 26 sequences
<table><caption>**Table 1**
Percentage of Identity of Amino Acids Sequences between FFAR-2 Chicken Paralogs</caption><thead><tr><th>Ensembl Protein ID</th><th>Ensembl Gene name</th><th>Ensembl Protein ID (ENSGALP000000)</th></tr></thead><tbody><tr><td>ENSGALP00000035053</td><td>Novel Ensembl prediction</td><td>35053 34827 15122 34808 34780 15906 28197 35396 34731 34729 24919 28377 35394 19013 35229 35368 34818 28102 24962 35207 35169 35315 38108 35408 34929 34941 25875 37768 04047 34061 24124 24125 37783</td></tr><tr><td>ENSGALP00000034837</td><td>Novel Ensembl prediction</td><td>91.3 100.0</td></tr><tr><td>ENSGALP00000015212</td><td>Novel Ensembl prediction</td><td>96.9 94.2 100.0</td></tr><tr><td>ENSGALP00000034808</td><td>Novel Ensembl prediction</td><td>94.6 92.4 95.8 100.0</td></tr><tr><td>ENSGALP00000034780</td><td>Novel Ensembl prediction</td><td>88.1 88.8 90.1 89.1 100.0</td></tr><tr><td>ENSGALP00000019806</td><td>Novel Ensembl prediction</td><td>89.1 88.9 88.8 89.2 98.2 100.0</td></tr><tr><td>ENSGALP00000028197</td><td>Novel Ensembl prediction</td><td>89.9 89.1 90.5 90.2 82.6 97.5 100.0</td></tr><tr><td>ENSGALP00000035396</td><td>Novel Ensembl prediction</td><td>86.9 85.1 84.4 85.7 91.3 91.3 92.0 100.0</td></tr><tr><td>ENSGALP00000034731</td><td>Novel Ensembl prediction</td><td>89.8 85.6 92.3 90.0 92.9 92.8 93.1 90.9 100.0</td></tr><tr><td>ENSGALP00000034729</td><td>Novel Ensembl prediction</td><td>85.3 85.1 85.7 86.1 92.4 94.2 95.7 84.2 92.3 100.0</td></tr><tr><td>ENSGALP00000034919</td><td>Novel Ensembl prediction</td><td>86.3 86.8 87.2 91.7 93.5 95.0 90.6 92.3 88.8 100.0</td></tr><tr><td>ENSGALP00000028377</td><td>Novel Ensembl prediction</td><td>89.3 88.5 90.7 90.7 94.4 95.1 96.0 94.1 94.9 96.7 96.7 100.0</td></tr><tr><td>ENSGALP00000035394</td><td>Novel Ensembl prediction</td><td>87.8 90.0 89.4 89.3 93.6 94.4 94.5 93.4 94.1 96.0 97.4 88.5 100.0</td></tr><tr><td>ENSGALP00000019013</td><td>Novel Ensembl prediction</td><td>84.0 85.5 85.9 86.9 88.8 89.9 92.1 91.4 89.2 86.5 87.8 94.7 93.9 100.0</td></tr><tr><td>ENSGALP00000035229</td><td>Novel Ensembl prediction</td><td>83.8 88.8 88.1 85.6 86.2 85.6 86.7 91.9 87.5 88.3 89.4 90.5 90.5 90.4 100.0</td></tr><tr><td>ENSGALP00000035368</td><td>Novel Ensembl prediction</td><td>83.3 85.1 83.3 85.4 87.2 88.9 91.0 88.4 88.5 84.9 90.6 93.0 93.0 88.2 89.0 100.0</td></tr><tr><td>ENSGALP00000034818</td><td>Novel Ensembl prediction</td><td>87.4 88.1 87.4 86.4 89.6 90.8 91.4 92.2 90.4 90.6 93.0 92.1 88.6 87.8 92.6 100.0</td></tr><tr><td>ENSGALP00000028102</td><td>Novel Ensembl prediction</td><td>86.7 87.7 85.2 87.0 92.0 92.2 93.1 90.1 91.1 90.4 91.8 94.5 94.5 92.2 90.8 92.8 95.9 100.0</td></tr><tr><td>ENSGALP00000034962</td><td>Novel Ensembl prediction</td><td>86.8 87.0 87.5 88.3 91.6 93.8 96.0 85.7 92.3 88.1 89.1 96.7 95.2 89.5 88.9 88.2 92.6 95.8 100.0</td></tr><tr><td>ENSGALP00000035207</td><td>Novel Ensembl prediction</td><td>86.1 85.8 85.3 87.5 91.7 91.1 92.0 91.6 90.8 91.3 92.6 94.1 92.7 93.4 93.6 94.3 93.6 94.3 96.0 88.8 93.8 93.9 95.3 95.3 96.7 96.7 100.0</td></tr><tr><td>ENSGALP00000035169</td><td>Novel Ensembl prediction</td><td>85.9 85.0 84.8 86.0 89.9 90.9 91.7 91.4 90.6 88.7 91.7 93.8 87.5 88.6 92.1 92.9 92.1 88.2 95.6 100.0</td></tr><tr><td>ENSGALP00000035315</td><td>Novel Ensembl prediction</td><td>87.6 88.7 87.4 88.7 90.8 91.1 92.0 91.7 90.2 93.1 94.6 93.4 93.4 90.6 88.8 93.8 93.9 95.3 95.3 96.7 96.7 100.0</td></tr><tr><td>ENSGALP0000002808</td><td>Novel Ensembl prediction</td><td>89.5 87.7 89.1 89.7 92.8 93.0 93.1 89.1 92.5 91.4 92.1 93.8 93.0 92.2 88.6 90.8 93.9 93.8 94.1 94.0 90.8 93.8 100.0</td></tr><tr><td>ENSGALP00000035408</td><td>Novel Ensembl prediction</td><td>87.2 86.6 87.0 87.9 92.8 92.3 93.2 92.9 91.8 92.0 94.3 93.8 93.8 93.1 91.6 92.0 93.1 94.6 95.0 94.3 92.3 94.2 95.0 100.0</td></tr><tr><td>ENSGALP00000034929</td><td>Novel Ensembl prediction</td><td>87.8 87.0 86.7 87.7 92.4 91.9 92.0 90.1 91.5 90.1 92.7 93.4 94.1 92.5 91.9 89.5 92.7 94.7 92.1 94.3 90.1 94.2 94.4 94.9 98.0 100.0</td></tr><tr><td>ENSGALP00000034941</td><td>Novel Ensembl prediction</td><td>87.4 86.0 85.9 88.0 92.0 91.5 92.4 90.4 91.1 89.8 92.4 93.8 93.8 92.5 92.3 89.8 93.1 95.1 92.4 94.6 90.5 94.6 94.1 98.3 98.4 100.0</td></tr><tr><td>ENSGALP00000026875</td><td>P2RY8</td><td>26.5 26.0 28.4 26.7 28.2 28.1 29.0 26.1 28.7 26.0 25.9 28.5 26.1 29.5 24.4 28.0 27.2 25.2 25.3 24.9 27.7 26.6 26.8 26.6 26.9 100.0</td></tr><tr><td>ENSGALP00000037768</td><td>F2RL2</td><td>26.5 26.4 25.9 26.7 26.0 27.8 27.9 26.3 27.7 26.5 25.7 27.0 26.6 25.8 27.7 24.3 27.2 26.3 25.9 24.9 27.0 27.4 26.8 25.9 26.3 35.2 100.0</td></tr><tr><td>ENSGALP0000004047</td><td>Novel Ensembl prediction</td><td>24.4 24.3 23.5 25.4 24.9 25.5 25.3 23.2 25.4 22.8 23.2 24.3 23.9 22.1 25.8 23.2 24.5 25.0 24.4 23.8 24.3 24.5 24.3 24.5 24.3 24.6 32.0 38.1 100.0</td></tr><tr><td>ENSGALP00000034861</td><td>F2RL3</td><td>24.2 23.7 23.6 24.0 23.0 23.0 22.4 23.3 21.9 22.5 23.2 22.3 23.1 21.1 23.5 23.8 22.2 22.5 21.7 23.3 23.8 23.4 23.1 23.5 30.0 40.4 41.5 100.0</td></tr><tr><td>ENSGALP00000024124</td><td>F2RL4</td><td>25.4 25.0 26.4 25.0 27.0 26.6 27.1 25.9 27.2 23.9 25.7 25.8 24.8 26.4 24.2 28.4 26.0 24.6 24.7 25.1 26.1 25.6 25.9 26.0 35.5 31.6 35.4 32.4 100.0</td></tr><tr><td>ENSGALP00000034861</td><td>F2RL1</td><td>27.2 27.0 28.3 27.0 28.6 29.0 29.4 27.0 28.8 27.0 26.9 28.6 28.6 29.2 25.0 27.7 27.2 25.7 25.8 25.3 27.4 26.8 27.4 26.5 26.8 37.5 35.2 40.1 39.3 35.1 100.0</td></tr><tr><td>ENSGALP00000037783</td><td>F2RL3</td><td>25.6 25.5 27.7 26.1 24.4 25.0 25.5 25.3 24.7 24.8 24.5 24.3 26.2 23.6 26.5 25.2 23.9 22.8 23.9 25.8 25.5 25.2 25.8 39.1 34.6 49.3 44.0 37.8 100</td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.**—Phylogenetic tree of FFAR2 paralogous genes. The phylogenetic tree was reconstructed using PhyML (Guindon and Gascuel 2003). Bootstrap values are given when nodes are strongly supported (>70%). The scale represents the substitution rate. The sequence with an * has a specific insertion of eight amino acids (DNGSEADG) at the following positions: ENSGALP00000034780 (163-170), ENSGALP00000028197 (152-159), and ENSGALP00000019806 (158-165).</figcaption></figure>

(located in position 172 of the alignment as indicated by a blue star, supplementary data S7, Supplementary Material online, and fig. 3), whereas the second one (Arg<sup>255</sup> in human FFAR2) was present in 25 out of 26 sequences (position 247 of the alignment indicated by a light blue star). The protein ENSGALP00000034818 is truncated before the arginine residue, which strongly suggests that the receptor function of this gene is altered. One histidine residue (His<sup>242</sup> in human FFAR2), considered as playing an important role in

fatty acid length selectivity, was conserved in chicken FFAR2 (position 234 of the alignment as indicated by a red star). This histidine was present only in human FFAR2 and human FFAR3, but was replaced by an asparagine in human FFAR1. The other histidine residue that contributes to recognition of the carboxylate (His<sup>140</sup> in human FFAR2) was found in 24 out of 26 sequences (position 133 of the alignment, green star); noticeable exceptions concerned proteins ENSGALP00000035229 and ENSGALP00000034731, which could exhibit differences
<figure><img src="image_4.png" /><figcaption>**Fig. 3.**—Amino acids that are critical for ligand binding are conserved within the FFAR gene family. Amino acid sequences of one of the chicken FFAR2 paralogs, and FFAR1, -2 and -3 human sequences were aligned. The well-conserved amino acids that are critical of ligand binding are indicated with colored stars. Sequence identities are reported white on a black background, whereas similarities are boxed. Secondary structures, as deduced from our modeling (supplementary data S7, Supplementary Material online, and fig. 8), are reported above the sequences.</figcaption></figure>

in their specificity/affinity for ligands. Analysis of the sequences obtained after DNA resequencing of chickens confirmed the lack of this His¹³³ in the ENSGALP00000035229 sequence. Of note, a SNP was found at this position in the ENSGALP00000034731 with either a proline or a histidine (50% of both variants in the studied populations).

## Chromosomal Localization

"Universal" primers were used to amplify a consensus fragment of FFAR2, giving rise to a unique band on agarose gels (supplementary data S8, Supplementary Material online). Sequencing of this band confirmed that we amplified several FFAR2 fragments. These fragments were mapped near the SEQ0067 marker (accession number AJ862640) and the

ROS0264 marker (31.7 cR, LOD = 10.8, accession number XM_003643130; fig. 4), which is located inside the predicted LIG1 gene (ligase I, DNA, ATP-dependent). These two markers are located on linkage group E64 (ROS0264: chrLGE64:721094 + 721330, SEQ0067: chrLGE64:709088 + 709284, UCSC Galgal4).

Two other markers in the vicinity of FFAR2 in the human genome were analyzed on the RH panel: MAG (myelin-associated glycoprotein) and USF2 (upstream transcription factor 2, c-fos interacting). Both map to the same linkage group, not linked to E64 (fig. 4). This suggests that, contrary to the human genome, MAG and USF2 are not located near FFAR2 genes in the chicken genome, suggesting a genomic rearrangement.
<figure><img src="image_5.png" /><figcaption>**Fig. 4.**—Localization on human chromosome 19 of genes from the FFAR2 region. (A) Fragment of human chromosome 19 from 35.9 to 48.6 Mb and (B) FFAR2 region on human chromosome 19. Maps are given in Mb, from assembly version GRCh37/hg19. At the left of each map are shown blocks of conserved synteny between HSA19 and the chicken genome, obtained through RH mapping. Genes with an asterisk are not found in the current version of the chicken genome (Galgal4).</figcaption></figure>

## Tissue Expression of FFAR2 Paralogs

### Chicken FFAR2 Genes Are Highly Expressed in the Testis

The pattern of expression of FFAR2 mRNA was analyzed by quantitative PCR (qPCR), using "universal" primers able to amplify all the 26 mRNA paralogs on 18 tissues sampled from chickens. FFAR2 was expressed at low levels in all these tissues (Cq between 27 and 36). The highest expression level was observed in testis, spleen, PBMC, abdominal adipose tissue, intestine, and lung (fig. 5A). Despite the high percentage of sequence identity between paralogs, we were also able to design primer pairs able to specifically amplify five different paralogs. The specificity in amplification was confirmed by sequencing. No significant difference in expression levels between the subgroup of these five paralogs was detected, and the higher expression levels of these paralogs were confirmed in testis, spleen, PBMC, and adipose tissue.

A particular attention was then paid to testis, which was divided in seminiferous tubules and Sertoli cells. A similar expression level of FFAR2 was observed in Sertoli cells and seminiferous tubules, which are composed of Sertoli cells and

germ cells (spermatogonia and type I spermatocyte stage). Expression level is higher in the whole testis (fig. 5B).

### Expression of FFAR2 Genes Is Inversely Related to Adiposity

To explore whether the expression of FFAR2 paralogs can be dependent of body adiposity, FFAR2 expression levels were measured using qPCR, in abdominal (visceral) adipose tissue from two experimental chicken lines divergently selected for abdominal fatness (Leclercq et al. 1980). Expression of FFAR2 genes was 2.7-fold higher (P = 0.008) in adipose tissue of the lean line compared with the fat line (fig. 5C). Similar results were obtained using adipose chicken RNA-Seq data, for which total FFAR2 read numbers were markedly higher (3.78 vs. 1.34; P = 0.034) in the lean line than in the fat line.

### Species-Specific Expression of FFAR Genes in Tissues Involved in Lipid Metabolism

We compared the expression of chicken FFAR genes with porcine FFAR genes. Expression level of the genes of the
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5.**—FFAR2 mRNA expression in different tissues in chicken (A), in chicken Sertoli cells, seminiferous tubules, and whole testis (B), in adipose tissue of lean and fat divergent lines (C). Relative tissue abundances of FFAR2 genes were examined by RT-qPCR using "universal" primers able to amplify the 26 sequences of FFAR2. The ribosomal subunit 18 S, RPL15, and GAPDH were used as reference genes. Y axis indicates relative level of mRNA expressed in comparison to the lowest level observed (Uropygial gland [A], seminiferous tubules [B], and fat chickens [C]). Error bars indicate standard deviation between three biological replicates (A), three independent experiments performed in triplicate (B), and ten animals (C). Different letters indicate significant expression difference (*P* < 0.05).</figcaption></figure>

FFAR family in pig is interesting on several aspects. First, pig is another domesticated species that show divergence along the phylogenetic tree with chicken, and second, although de novo lipid synthesis occurs in the liver in chicken, it takes place in the adipose tissue in the pig. Expression levels of the FFAR family genes were compared between chicken and pig in both liver and adipose tissues (fig. 6 and *supplementary data S4, Supplementary Material online*). Although we observed a higher expression of FFAR2 paralogs in adipose tissue (39 reads, 2.56 fragments per kilobase of exon per million fragments mapped [FPKM]) compare with the liver (12 reads, 0.7 FPKM) in the chicken, the expression of the porcine FFAR2 was low in both tissues (<10 reads). Expression levels for other members of the FFAR family revealed that FFAR4 was highly

expressed in the pig (632 reads, 39 FPKM). Finally, very low expression levels were observed for porcine FFAR1 and FFAR3 in both adipose tissue (0 and 0.1 reads, respectively) and liver (0 and 0.4 reads, respectively).

## Estimation of Gene Conversion and Positive Selection in FFAR2 Gene

To determine the selective pressures that may have triggered the evolution of the FFAR2 gene clusters in the chicken genome, we conducted two different analyses, the gene conversion events detection and the positive selection calculation. Both analyses were conducted using the multiple amino acid sequence alignment and the phylogenetic tree (fig. 2 and *supplementary data S7, Supplementary Material online*).
<figure><img src="image_8.png" /><figcaption>**Fig. 6.**—FFAR expression in chicken and pig using RNA-Seq in adipose tissue (A) and liver (B). Values correspond to the mean of eight individual expression counts (four males and four females). Chicken FFAR2 reads are the sum of all the paralog expression counts. Errors bars indicate standard deviation between individual reads.</figcaption></figure>

Gene conversion can occur between sequences if they present enough sequence identity. By using GENECONV tool, we searched for statistical evidence of this phenomenon. We used the option g0 that allows no mismatch between sequences involved in a gene conversion event, so fragments resulting from gene conversion events contain no polymorphism. The control of "randomize sites," which randomizes the order of polymorphic sites before analysis, detected no gene conversion events in the cluster, so that the results of subsequent GENECONV analyses for FFAR2 are reliable. Each ribbon on figure 7 represents a unique conversion event between FFAR2 paralogs in chicken (supplementary data S9, Supplementary Material online). We found that gene conversion is a pervasive evolutionary mechanism in the FFAR2 cluster in the chicken genome with 24 out to 26 paralogs involved in at least one recent gene conversion (fig. 7 and supplementary data S9, Supplementary Material online). Moreover, by using the site models implemented in PhyleasProg, which allows the ω ratio to vary among sites, we detected eight different positions under positive selection scattered along the sequence of the FFAR2 paralogs. By using the branch-site models, which estimate different ω values among branches and among sites, we identified three amino acids specifically under positive selection for one of the paralogs (Table 2).

Taken together, this analysis of gene conversion events on FFAR2 cluster showed that extensive gene conversion led to sequence homogenization. However, the amino acids under positive selection may be involved in the diversification of function of the different paralogs.

### Mapping of Amino Acids under Positive Selection on the 3D Structure of Chicken FFAR2

A model of the 3D structure of chicken FFAR2 (ENSGALP00000034780) was made, using the experimental

3D structure of the proteinase-activated receptor 1 as a template (pdb 3vw7, 25% identity, HH-Pred probability 100, E value $4.5 \times 10^{-45}$). Despite the low levels of sequence identity, the protein sequences were well aligned, with clear anchor points within the transmembrane helices and insertions occurring within loops (see the corresponding alignment in supplementary data S10, Supplementary Material online). The four basic amino acids mentioned above to be crucial for receptor functionality were clearly located in the ligand-binding groove, at the top of the transmembrane assembly (fig. 8). One of the amino acids under positive selection (Y246) participates in this groove, whereas the other amino acids described here did not appear to be in critical positions. However, two of them (M80 and H182) lied in the extracellular loops, in the vicinity of the ligand-binding groove.

### Discussion

The chicken was the first avian species and domestic animal selected for complete genome sequencing and assembly (International Chicken Genome Sequencing Consortium [ICGSC] Hillier et al. 2004). The chicken genome (as other bird genomes) however displays a reduced size compared with mammals. Moreover, phylogenetic analyses of chicken, mammalian, and fish paralogs support the hypothesis that loss of paralogs occurred much more frequently in chicken than in mammals (Hughes and Friedman 2008). However, several gene duplications in chicken have also been described, such as genes encoding Toll-receptors (ICGSC 2004; Temperley et al. 2008). Specific epimutations have been also recently associated with genes related to the bone morphogenic protein, toll receptor, and melanogenesis signaling pathways in birds (Skinner et al. 2014). Beside these pathways, genes
<figure><img src="image_9.png" /><figcaption>**Fig. 7.**—Gene conversion events in the FFAR2 cluster. Each line represents a gene and each ribbon represents a gene conversion event between two genes of the FFAR2 cluster. The thickness of the ribbon represents the significance of each identified fragment. The biggest, medium and smallest ribbons represent the fragments identified with three, two and one significant P values, respectively. For clarity, Ensembl gene IDs are shortened. The correspondence between reduced gene ID, gene ID, and protein ID can be found in the *supplementary data S1*, *Supplementary Material* online. The coordinates of each fragment can be found in the *supplementary data S9*, *Supplementary Material* online. Note that only the 24 genes involved in at least one gene conversion event are represented.</figcaption></figure>

related to energy metabolism could have been affected dif-
ferently in birds compared with mammals, because chicken
exhibit normal hyperglycemia without signs for insulin resis-
tance. Recently, the subfamily of FFAR has received a consid-
erable interest in literature related to energy metabolism,
obesity and diabetes in human and rodents (Hara et al.
2014; Ichimura et al. 2012), because of their potential interest
as key therapeutic targets (Hara et al. 2014). A computational
analysis of the chicken genome notably reported an expansion
of eight homologs to human FFAR2 gene (GPR43; Lagerström
et al. 2006). Here, we prove the existence of more than 20
genes encoding FFAR2 paralogs in the chicken genome. First,

qPCR experiments on genomic DNA using different chicken
populations estimated a relative copy number of 23 ± 1.5 of
these paralogs in the chicken genome. Second, analysis of
genome DNA resequencing and RNA-Seq data confirmed
the massive duplication of FFAR2. Although having strong
identity (above 98.5%), 9.8–11.9 substitutions between
FFAR2 nucleotide sequences (each containing approximately
1,000 nucleotides) were detected. This is more than the 5
SNP/kb previously estimated for the entire chicken genome
(Wong et al. 2004). This is also much more than the approx-
imately 2.7 SNP/kb observed in coding regions between broi-
lers and red jungle fowl (Wong et al. 2004). This suggests that
<table><caption>**Table 2**
Parameter Estimates and Likelihood Scores for Branch-Site and Site Models</caption><thead><tr><th>Gene</th><th>Model</th><th>P<sup>a</sup></th><th>Estimates of Parameters<sup>b</sup></th><th>2Δ/</th><th>Positively Selected Sites (BEB)<sup>c</sup></th></tr></thead><tbody><tr><td>Branch-site models:</td><td>Null</td><td>2,170.978924</td><td>ρ0 = 0.42803, (ρ1 = 0.17627),<br>ω0 = 0.03667, (ω1 = 1)</td><td>7.97**</td><td>Not allowed</td></tr><tr><td>ENSGALG 00000022004</td><td>Alternative</td><td>2,166.991593</td><td>ρ0 = 0.66474, ρ1 = 0.27047, (ρ2 = 0.0648),<br>ω0 = 0.03763, (ω1 = 1), ω2 = ∞</td><td></td><td>Three sites P > 99%: 113E,<br>115H, and 201T</td></tr><tr><td rowspan="2">Site models</td><td>M8a</td><td>2,173.482314</td><td>ρ0 = 0.70692, ρ1 = 0.29308,<br>p = 4.58561 q = 99</td><td rowspan="2">20.45***</td><td rowspan="2">Not allowed</td></tr><tr><td>M8a</td><td>2,163.258223</td><td>ρ0 = 0.85604, ρs = 0.14396, p = 0.24218,<br>q = 1.07143, ωs = 2.58370</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Fig. 8.**—Ribbon representation of the 3D model of FFAR2. The model was made by homology with the 3D structure of the proteinase-activated receptor 1 (pdb 3vw7). This last structure was obtained in complex with the antagonist vorapaxar, whose binding site is shown here light gray shaded. The four critical amino acids (corresponding to H140, R180, H242, and R255 in fig. 3 after correction of the alignment) are presented in yellow, whereas amino acids under positive selection, as reported in table 2, are indicated in orange, red, and purple.</figcaption></figure>

true paralogs rather than allelic variants of FFAR2 genes exist in the chicken genome. Furthermore, sequence comparison of these FFAR2 chicken paralogs further showed that the four basic amino acids responsible of fatty acid binding (Hudson et al. 2013) were conserved among paralogs, suggesting that most of these paralogs are functional.

In mammalian genome, FFAR2 gene is located in a locus that also contains genes encoding two other types of the FFAR family (FFAR1, FFAR3), HAMP coding for hepcidin, involved in the regulation of iron metabolism, and USF2, the upstream transcription factor involved in lipid homeostasis. By exploring the chicken genome by reciprocal tBLASTn analysis on genome and EST database, we did not find traces of FFAR1 and FFAR3 genes in chicken. This suggests that the only one study that has reported the expression of a chicken FFAR to date, which was claimed to be FFAR1 (Suh et al. 2008), would refer rather to FFAR2. The hepcidin gene is absent from chicken, fish, and reptile genomes. Although chicken USF2 sequence is available for chicken (http://www.ncbi.nlm.nih.gov/protein/47420475, last accessed May 6, 2015), USF2 has not been still assigned to a chicken chromosome. The amplification pattern obtained here on the RH panel however indicates that the FFAR2 locus is not mapped near the USF2 linkage group: Very few hybrid clones were both positive for USF2 and the FFAR2 locus, indicating a different location for the FFAR2 locus and the portion of the genome orthologous to human HSA19 bearing USF2. Altogether, this suggests that this particular locus encompassing FFAR2 has been subjected to a specific rearrangement in chicken. We cannot discriminate between the 26 FFAR2 fragments, but sequencing of each hybrid clone amplicon proved that there are several different FFAR2 fragments in each clone, leading to the hypothesis that the ancestral FFAR2 gene has been duplicated in the same locus. Assuming a unique position for the FFAR2 locus, it

<footnote><sup>a</sup>Log-likelihood values.</footnote>

<footnote><sup>b</sup>Null and alternative models: ρ0, ρ1, and ρ2 are the proportions of codons subject to purifying selection, neutral evolution, and positive selection, respectively. ω0, ω1, and ω2 represented dN/dS for each class (purifying, neutral, and positive selection, respectively). M8 and M8a models: ωs, Average dN/dS ratio for sites subject to positive selection (models M2a and M8); p and q, shape parameters for the beta distribution of ω (models M7, M8, and MEC); ρ0, ρ1, and ρs are the proportions of codons subject to purifying selection, neutral evolution, and positive selection, respectively.</footnote>

<footnote><sup>c</sup>Site numbers and amino acids refer to the ENSGALP00000034780 sequence (gene ID: ENSGALG00000022004).</footnote>

<footnote>**Significant at P < 0.01 and ***P < 0.001.</footnote>
belongs to chicken linkage group E64, a small group (799,899 bp) probably corresponding to chicken microchro-mosome 31 (Masabanda et al. 2004). Altogether, this con-firms that there is no overall synteny conservation between human chromosome 19, where the FFAR2 gene is mapped, and the chicken genome, with only several small groups of synteny mainly localized on chicken microchromosomes or unknown linkage groups (Morisson et al. 2007). The duplica-tion of the FFAR2 locus may have occurred during the rear-rangements of these genomic fragments.

Interestingly, this expansion among FFAR2 genes seems to be chicken-specific, because tBLASTn analysis on ESTs from finch, turkey, or quail species found only one significant match in this study. We cannot totally rule out the possibility that further sequencing of other bird genomes may reveal FFAR2 duplication similar to that in chicken. However, as we found no trace of a massive duplication of FFAR2 in neither turkey nor quail but found FFAR2 duplicates in the ancestral line from Thailand, we can date the massive duplication event of FFAR2 between 30 Ma (estimated divergence between chicken and turkey) and 8,000 years ago (chicken domestication) (Dimcheff et al. 2002; Jarvis et al. 2014).

Although all FFAR2 paralogs seem to encode full-length proteins, the functionality of those genes remains unknown. Further investigations are needed to understand physiological significance of such a massive duplication in chicken. In human and rodents, FFAR2 has been described as a receptor specific of small-chain fatty acids (1-6 carbon chain length), so that it is activated by physiological concentrations of volatile FFAs acetate, propionate, and butyrate which are gut-micro-biote-derived bacterial fermentation products (for review, Hara et al. 2013). Whether chicken produce concentrations of small-chain fatty acids that differ from mammals would be interesting to investigate. The FFAR2 has been shown to mediate the beneficial effects associated with high soluble fiber diet and protect from high-fat diet-induced obesity and dyslipidemia. In particular, in vivo in the mouse, acetate is able to reduce levels of nonesterified fatty acids in serum, these effects being abolished in FFAR2-/- mice (Ge et al. 2008). Because chicken in modern breeding are usually fed cereal-based diets without major sources of fibers, the physiological significance of the massive FFAR2 duplication remains questionable. Activation of FFAR2 by short-chain fatty acids may lead to the inhibition of lipolysis and the decrease of FFA levels in serum, a phenomenon which may be amplified by the presence of these numerous FFAR2 paralogs in chicken.

Overall, whether the duplication of FFAR2 gene is involved in a particular regulation of energy metabolism in chicken compared with mammals has to be investigated. In human, specific duplications of amylase gene were described, and copy number of this gene was shown to be positively correlated with salivary amylase protein level (Perry et al. 2007). Moreover, individuals from populations with high-starch

diets have more copies than those with traditionally low-starch diets. Thus, the number of copies of amylase genes may be selected to adapt populations to their traditional diet. An adaptation to diet is then an interesting hypothesis to explain FFAR2 gene numerous duplications.

To gain insights into the functions of FFAR2 in chicken, we decided to investigate the tissue-expression of FFAR2 and its different paralogs. Considering the very high percentage of identity between nucleotide sequences of the FFAR2 paralogs, it is however impossible to produce primers that are specific of each of them. We demonstrated for the first time the expression of FFAR2 gene in testis, spleen, PBMC, adipose tissue, intestine, and lung of the chicken, with testis expressing the highest levels of these genes. Therefore, we also evaluated FFAR2 expression in different compartments of the chicken testis. We found equivalent mRNA levels in seminiferous tubules and Sertoli cell. Whether these genes are predominantly expressed in Sertoli cells or other somatic or germ cells remains questionable. Testis is known for permissive transcription (transcription of DNA caused by biochemical happenstance rather than usefulness) (Soumillon et al. 2013), so the physiological significance of expression of FFAR2 genes still remains to be experimentally explored. Further investigations of the role of FFAR2 paralogs in testis require specific experimental designs modulating male fertility after meal supplementation, together with in vitro modulation of Sertoli and/or Leydig cells steroidogenesis.

Importantly, FFAR2 was expressed in chicken and pig adipose tissue, as already described in mouse in this tissue (Hong et al. 2005; Kimura et al. 2013). However, because there is an evidence of an implication of FFAR2 on obesity (Dewulf et al. 2011; Kimura et al. 2013), we looked at the expression levels of FFAR2 in visceral adipose tissue between chicken lines divergently selected for abdominal fatness. We showed that FFAR2 was significantly overexpressed in the lean line compared with the fat line. This result is consistent with previous study in mouse (Kimura et al. 2013), showing that FFAR2 knockout mice are obese, whereas transgenic mice overexpressing adipose tissue-specific FFAR2 are lean.

Suh et al. (2008) reported the expression of a chicken FFAR claimed to be FFAR1 in primary cultured chicken hepatocytes. In their study, FFAR1 was activated with linoleic acid (C18:2), and FFAR1 protein expression was blocked by "FFAR1-specific siRNA" in vitro. However the corresponding sequence is not available, so it is not possible to check this experiment. Whether FFAR1 does exist in the chicken remains to be investigated.

Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has a homogenized part of the length of the chicken FFAR2 paralogs, suggesting that extensive gene conversion is responsible for the particularly high degree of sequence similarity between these genes, in particular the high conservation of the four amino acids known to be critical of ligand binding (Hudson
et al. 2013). Moreover, elevated nonsynonymous/synonymous substitution ratios on some amino acids within (Y248) or in the vicinity (Met80, H182) of the ligand-binding groove suggest that positive selection may have reduced the effective rate of gene conversion in this region, and may contribute to diversify/specialize the function of some FFAR2 paralogs in chicken.

In conclusion, this work brings significant new data on the evolution in the FFAR2 subfamily in chicken, by massive duplication, gene conversion, and positive selection. These dominant processes of their sequence evolution might have led to a sort of mixing of homogenization and diversification of protein functions. Further investigations are needed to study whether these recent chicken-specific events are the consequence of domestication. It is important to determine which might be the functions of these FFAR2 paralogs in chicken. Namely, it would be interesting to test whether those paralogs, if actually functional GPCRs, activate specific pathways and by which agonists they can be activated. Present data suggest roles of FFAR2 in testis and adipose tissue expansion in chicken.

## Supplementary Material

Supplementary data S1–S10 are available at *Genome Biology and Evolution* online (http://www.gbe.oxfordjournals.org/).

## Acknowledgments

This work was supported by INRA (Chick-seq and SOSrnaSEQ programs) and by French ANR (Programs FATINTEGER and EpiBird). The authors thank Mireille Morisson for providing the ChickRH panel.

## Literature Cited

Bjursell M, et al. 2011. Improved glucose control and reduced body fat mass in free fatty acid receptor 2-deficient mice fed a high-fat diet. *J Physiol Endocrinol Metab* 300(1):E211–E220.

Bowie JU, Luthy R, Eisenberg D. 1991. A method to identify protein sequences that fold into a known three-dimensional structure. *Science* 253:164–170.

Burt D. 2007. Emergence of the chicken as a model organism: implications for agriculture and biology. *Poult Sci.* 86:1460–1471.

Busset J, Cabau C, Meslin C, Pascal G. 2011. PhyleasProg: a user-oriented web server for wide evolutionary analyses. *Nucleic Acids Res.* 39: W479–W485.

Chirgwin JM, Przybyla AE, MacDonald RJ, Rutter WJ. 1979. Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease. *Biochemistry* 18:5294–5299.

Colombo M, Trevisi P, Gandolfi G, Bosi P. 2012. Assessment of the presence of chemosensing receptors based on bitter and fat taste in the gastrointestinal tract of young pig. *J Anim Sci.* 90:128–130.

De Givry S, Bouchez M, Chabrier P, Milan D, Schiex T. 2005. Carh ta Gene: multipopulation integrated genetic and radiation hybrid mapping. *Bioinformatics* 21:1703–1704.

Dewulf EM, et al. 2011. Inulin-type fructans with prebiotic properties counteract GPR43 overexpression and PPARγ-related adipogenesis in

the white adipose tissue of high-fat diet-fed mice. *J Nutr Biochem.* 22: 712–722.

Dimcheff DE, Drovetski SV, Mindell DP. 2002. Phylogeny of Tetraoninae and other galliform birds using mitochondrial 12S and ND2 genes. *Mol Phylogenet Evol.* 24:203–215.

Edgar RC. 2004. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. *BMC Bioinformatics* 5:113.

Ellingson D, Yao K. 1970. Growth and observations of Chinese hamster seminiferous epithelium in vitro. *J Cell Sci.* 6:195–205.

Foster P, Creasy DM, Foster JR, Gray T. 1984. Testicular toxicity produced by ethylene glycol monomethyl and monoethyl ethers in the rat. *Environ Health Perspect.* 57:207.

Frésard L, et al. 2014. Transcriptome-wide investigation of genomic imprinting in chicken. *Nucleic Acids Res.* 42:3768–3782.

Ge H, et al. 2008. Activation of G protein-coupled receptor 43 in adipocytes leads to inhibition of lipolysis and suppression of plasma free fatty acids. *Endocrinology* 149:4519–4526.

Guibert E, Brière S, Pelletier R, Brillard J, Froment P. 2011. Characterization of chicken Sertoli cells in vitro. *Poult Sci.* 90:1276–1286.

Guindon S, Gascuel O. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. *Syst Biol.* 52: 696–704.

Hara T, et al. 2014. Role of free fatty acid receptors in the regulation of energy metabolism. *Biochim Biophys Acta Mol Cell Biol Lipids.* 1841: 1292–1300.

Hara T, Kimura I, Inoue D, Ichimura A, Hirasawa A. 2013. Free fatty acid receptors and their role in regulation of energy metabolism. *Rev Physiol Biochem Pharmacol.* 164:77–116.

Hillier LW, et al. 2004. International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. *Nature* 432(7018):69–716.

Hong Y-H, et al. 2005. Acetate and propionate short chain fatty acids stimulate adipogenesis via GPCR43. *Endocrinology* 146:5092–5099.

Hou Z, Sun C. 2008. Transcriptional expression of GPR43 gene in adipose tissue and primary cultured adipocytes of pig. *Chin J Biotechnol* [Sheng wu gong cheng xue bao]. 24:1361–1366.

Hudson BD, Murdoch H, Milligan G. 2013. Minireview: the effects of species ortholog and SNP variation on receptors for free fatty acids. *Mol Endocrinol.* 27:1177–1187.

Hughes AL, Friedman R. 2008. Genome size reduction in the chicken has involved massive loss of ancestral protein-coding genes. *Mol Biol Evol.* 25:2681–2688.

Ichimura A, et al. 2012. Dysfunction of lipid sensor GPR120 leads to obesity in both mouse and human. *Nature* 483(7389):350–354.

Ichimura A, Hirasawa A, Hara T, Tsujimoto G. 2009. Free fatty acid receptors act as nutrient sensors to regulate energy homeostasis. *Prostaglandins other Lipid Mediat.* 89:82–88.

Jarvis ED, et al. 2014. Whole-genome analyses resolve early branches in the tree of life of modern birds. *Science* 346:1320–1331.

Kimura I, et al. 2013. The gut microbiota suppresses insulin-mediated fat accumulation via the short-chain fatty acid receptor GPR43. *Nat Commun.* 4:1829.

Krzywinski M, et al. 2009. Circos: an information aesthetic for comparative genomics. *Genome Res.* 19:1639–1645.

Lagerström MC, et al. 2006. The G protein-coupled receptor subset of the chicken genome. *PLoS Comput Biol.* 2:e54.

Laskowski RA, MacArthur MW, Moss DS, Thornton JM. 1993. PROCHECK: a program to check the stereochemical quality of protein structures. *J Appl Crystallogr.* 26:283–291.

Leclercq B, Blum J, Boyer J. 1980. Selecting broilers for low or high abdominal fat: initial observations. *Br Poult Sci.* 21:107–113.

Li H, Durbin R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* 25:1754–1760.
Li H, et al. 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25:2078-2079.

Livak KJ, Schmittgen TD. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2⁻ΔΔCT method. Methods 25:402-408.

Martí-Renom MA, et al. 2000. Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 29: 291-325.

Masabanda JS, et al. 2004. Molecular cytogenetic definition of the chicken genome: the first complete avian karyotype. Genetics 166:1367-1373.

Morisson M, et al. 2002. ChickRH6: a chicken whole-genome radiation hybrid panel. Genet Sel Evol. 34:521.

Morisson M, et al. 2007. The chicken RH map: current state of progress and microchromosome mapping. Cytogenet Genome Res. 117: 14-21.

Perry GH, et al. 2007. Diet and the evolution of human amylase gene copy number variation. Nat Genet. 39:1256-1260.

Pettersen EF, et al. 2004. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 25: 1605-1612.

Posada D. 2002. Evaluation of methods for detecting recombination from DNA sequences: empirical data. Mol Biol Evol. 19:708-717.

Robert X, Gouet P. 2014. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 42(W1): W320-W324.

Sawyer S. 1989. Statistical tests for detecting gene conversion. Mol Biol Evol. 6:526-538.

Simon J. 1989. Chicken as a useful species for the comprehension of insulin action. Crit Rev Poult Biol. 2(2):121-148.

Skinner MK, et al. 2014. Epigenetics and the Evolution of Darwin's Finches. Genome Biol Evol. 6(8):1972-1989.

Smith S, et al. 2012. Adipogenic gene expression and fatty acid composition in subcutaneous adipose tissue depots of Angus steers between 9 and 16 months of age. J Anim Sci. 90:2505-2514.

Söding J, Biegert A, Lupas AN. 2005. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 33:W244-W248.

Soumillon M, et al. 2013. Cellular source and mechanisms of high transcriptome complexity in the mammalian testis. Cell Rep. 3: 2179-2190.

Staub C, et al. 2000. The whole meiotic process can occur in vitro in untransformed rat spermatogenic cells. Exp Cell Res. 260:85-95.

Stoddart LA, Smith NJ, Jenkins L, Brown AJ, Milligan G. 2008. Conserved polar residues in transmembrane domains V, VI, and VII of free fatty acid receptor 2 and free fatty acid receptor 3 are required for the binding and function of short chain fatty acids. J Biol Chem. 283: 32913-32924.

Storey JD, Tibshirani R. 2003. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 100:9440-9445.

Suh HN, Huong HT, Song CH, Lee JH, Han HJ. 2008. Linoleic acid stimulates gluconeogenesis via Ca2+/PLC, cPLA2, and PPAR pathways through GPR40 in primary cultured chicken hepatocytes. Am J Physiol Cell Physiol. 295:C1518-C1527.

Sum CS, et al. 2007. Identification of residues important for agonist recognition and activation in GPR40. J Biol Chem. 282(40): 29248-29255.

Temperley ND, Berlin S, Paton IR, Griffin DK, Burt DW. 2008. Evolution of the chicken Toll-like receptor gene family: a story of gene gain and gene loss. BMC Genomics 9:62.

Toebosch A, Brussée R, Verkerk A, Grootegoed J. 1989. Quantitative evaluation of the maintenance and development of spermatocytes and round spermatids in cultured tubule fragments from immature rat testis. Int J Androl. 12:360-374.

Ulven T. 2012. Short-chain free fatty acid receptors FFA2/GPR43 and FFA3/GPR41 as new potential therapeutic targets. Front Endocrinol. 3:111; doi: 10.3389/fendo.2012.00111.

Voorrips R. 2002. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 93:77-78.

Wang A, Gu Z, Heid B, Akers R, Jiang H. 2009. Identification and characterization of the bovine G protein-coupled receptor GPR41 and GPR43 genes. J Dairy Sci. 92:2696-2705.

Wang J, Wu X, Simonavicius N, Tian H, Ling L. 2006. Medium-chain fatty acids as ligands for orphan G protein-coupled receptor GPR84. J Biol Chem. 281:34457-34464.

Wong GK, et al. 2004. A genetic variation map for chicken with 2.8 million single-nucleotide polymorphisms. Nature 432:717-722.

Zhang C, et al. 2012. High-resolution crystal structure of human protease-activated receptor 1. Nature 492:387-392.

Zwain IH, Cheng CY. 1994. Rat seminiferous tubular culture medium contains a biological factor that inhibits Leydig cell steroidogenesis: its purification and mechanism of action. Mol Cell Endocrinol. 104: 213-227.

Associate editor: Dan Graur