- Open Access
Characterization of the endogenous retrovirus insertion in CYP19A1 associated with henny feathering in chicken
Mobile DNAvolume 10, Article number: 38 (2019)
Henny feathering in chickens is determined by a dominant mutation that transforms male-specific plumage to female-like plumage. Previous studies indicated that this phenotype is caused by ectopic expression in skin of CYP19A1 encoding aromatase that converts androgens to estrogen and thereby inhibits the development of male-specific plumage. A long terminal repeat (LTR) from an uncharacterized endogenous retrovirus (ERV) insertion was found in an isoform of the CYP19A1 transcript from henny feathering chicken. However, the complete sequence and the genomic position of the insertion were not determined.
We used publicly available whole genome sequence data to determine the flanking sequences of the ERV, and then PCR amplified the entire insertion and sequenced it using Nanopore long reads and Sanger sequencing. The 7524 bp insertion contains an intact endogenous retrovirus that was not found in chickens representing 31 different breeds not showing henny feathering or in samples of the ancestral red junglefowl. The sequence shows over 99% sequence identity to the avian leukosis virus ev-1 and ev-21 strains, suggesting a recent integration. The ERV 3’LTR, containing a powerful transcriptional enhancer and core promoter with TATA box together with binding sites for EFIII and Ig/EBP inside the CYP19A1 5′ untranslated region, was detected partially in an aromatase transcript, which present a plausible explanation for ectopic expression of aromatase in non-ovarian tissues underlying the henny feathering phenotype.
We demonstrate that the henny feathering allele harbors an insertion of an intact avian leukosis virus at the 5’end of CYP19A1. The presence of this ERV showed complete concordance with the henny feathering phenotype both within a pedigree segregating for this phenotype and across breeds.
Plumage color is a striking and variable aspect of sexual dimorphism in avian species. Males often exhibit a showy plumage during mating season, which is an example of a compromise between sexual selection and avoiding predation by camouflage color outside the breeding season [1, 2]. Understanding the molecular mechanisms of male feathering therefore involves basic principles in developmental and evolutionary biology. The wild ancestor of the domestic chicken, the red junglefowl, shows a spectacular sexual dimorphism like many other pheasant species. However, some domestic chickens carry the dominant henny feathering allele that makes the plumage of males resemble females . This phenotype offers a model to study the molecular basis for sexual dimorphism in avian species and has therefore been extensively studied . The henny feathering trait was first observed by Sir John Sebright in 1800 , and became a distinguishing feature of the Sebright Bantam breed (Fig. 1). Other henny feathered chicken breeds such as the Golden Campine likely received this trait from Sebright due to interbreeding by chicken fanciers . The henny feathering (Hf) mutation is desirable when breeding fancy chickens because it allows males to show intrafeather patterns which are normally interrupted in males homozygous for the wild-type allele (hf+). In many other show chicken breeds, only females can meet the criteria for exhibition because of their uniformity of intrafeather patterns. However, in Sebright and Golden Campine, males are preferred because Hf allows them to express uniform intrafeather patterning while retaining other male characters, like larger combs and elongated tail feathers .
The plumage in hf+ male chickens has a shiny appearance with pointy tips in saddle and hackle feathers. Hens have a more ‘dun’ plumage with duller, rounder feathers. Testis removal does not affect the appearance of these feathers, however ovary removal in females induces a male-shiny phenotype . The typical feathering of hens is produced by the action of estrogen, predominantly produced in the ovaries through aromatase conversion from androgens . In early studies of the henny feathering phenotype, castration of gonads [10, 11], castration followed by treatment with different androgens [12, 13], as well as transplantation of gonads  and skin  suggested that the changes induced by Hf must reside in the skin, decreasing testosterone and increasing estrogen, resulting in a feminized plumage . This is consistent with current knowledge that Hf is caused by ectopic expression of aromatase in the skin of Sebright and Campine chickens [4, 16,17,18,19,20,21]. In adult wild-type chickens, aromatase is detectable only in ovaries and hypothalamus . In males carrying the dominant Hf allele, such conversion also happens in non-ovarian tissues such as skin . Hf is co-dominant with respect to aromatase activity, as heterozygotes have an enzyme activity intermediate between the two homozygotes. It is dominant with respect to feathering, as a single dose of the Hf allele leads to the production of adequate estrogen in skin to induce henny feathering . Furthermore, a study on transgenic cocks overexpressing aromatase also resulted in female feathering resembling the Hf phenotype . In other studies, after the injection of an aromatase inhibitor during embryonic stage, female chickens developed testes and had the behavior and physical appearance of males [9, 23].
A long terminal repeat (LTR) from an uncharacterized inherited endogenous retrovirus (ERV) was found in the cDNA of CYP19A1, which encodes aromatase in chicken. This LTR was hypothesized to act as a promoter inducing ubiquitous upregulation of CYP19A1 and thus inducing aromatase activity in non-ovarian tissues . However, the ERV was not characterized in detail, and the mechanism causing ectopic aromatase expression is not fully understood. Here we used publicly available whole genome sequence (WGS) data from chickens with and without henny feathering to identify the flanking sequence of the insertion, and long-read Nanopore as well as Sanger sequencing to characterize the locus. We show that the insertion is indeed an intact 7524 bp recently integrated ERV, which is related to avian leukosis viruses and missing in the ancestral red junglefowl as well as in related chicken breeds.
Insertion site identified using WGS data
Guided by the localization of the putative ERV adjacent to CYP19A1 , we used publicly available whole genome sequence data from individually sequenced chickens and from pooled sequencing in total representing 32 breeds of domestic chicken, and in addition five individually sequenced and two pools of red junglefowl (Additional file 1: Table S1). We focused the analysis on the 200 kb region surrounding CYP19A1 to identify structural variants present in Sebright chickens (Hf/−) and not present in chickens from other breeds (hf+/hf+). Only one structural variant met this criterion, and was present at position 9,683,879 bp on chromosome 10 of GalGal6, which is in the 5′ untranslated region (5’UTR) of CYP19A1 (Fig. 2a). The sequence of soft-clipped reads surrounding the insertion position confirmed the presence of an LTR (Additional file 5: Figure S1).
Sequencing of the Hf associated ERV insertion
PCR amplification of the Hf insertion produced an amplicon of about 7.5 kb from all Hf/− chickens whereas only short amplicons (163 bp) were obtained from wild-type chickens (Fig. 3). The polished de novo contig from long-read sequencing (GenBank accession number MK937054) was parsed for retroviral similarity, revealing high sequence identity to avian leukosis virus. A sequence comparison of the henny feathering associated ERV (Hf_ERV), described in this study, with 140 published avian leukosis viruses identified 50 retroviruses with over 90% identity (Additional file 2: Table S2). The highest sequence identity (99.4%) was obtained against the ev-21 strain, which shares an origin with the Rous sarcoma virus (Additional file 6: Figure S2). Examination of the Hf_ERV sequence motifs revealed intact gag, pro, pol and env genes flanked by two identical 274 bp LTRs (Fig. 2b). Each LTR contains a TATA box promoter sequence , as well as binding sites for the avian serum response factor EFIII and Ig/EBP  (Fig. 2c), a ubiquitously expressed immunoglobulin enhancer binding protein. Immediately flanking the LTRs, we found identical 6 bp sequences, GAGGAC, identified as the chromosomal target site duplications (TSD) formed during integration (Fig. 2b).
5’RACE of CYP19A1 reveals start of transcription in Hf_ERV
5’RACE experiments using CYP19A1-specific primers (Additional file 3: Table S3) and mRNA isolated from the skin of a one-day-old Sebright chick (genotyped as Hf/Hf) revealed only one transcript isoform (GenBank accession number MN124451). The sequence data showed that the transcription start is 24 bp downstream of the 3’LTR TATA box. It includes the remaining 99 bp of the 3’LTR together with the TSD, as well as partial 5’UTR (4 bp) of the wild-type CYP19A1 transcript (NM_001364699.1), and the CYP19A1 coding sequence (CDS) (Fig. 2a). These sequence features suggest that the 3’LTR TATA box of Hf_ERV acts as the promoter that induce the transcription of the Hf_ERV and CYP19A1 fusion transcript, which uses the same splicing sites as the wild-type transcript. The cDNA sequence was fully consistent with the genomic sequence of the Hf_ERV insertion.
Hf_ERV-induced transcripts are detected in Sebright but not in wild-type chicken
We explored the CYP19A1 expression pattern across six tissues (skin, lung, intestine, hypothalamus, muscle, and liver) and using two sets of primer combinations, one encompassing only coding sequences and expected to amplify all CYP19A1 transcripts and the other only detecting HF_ERV-induced transcripts (Fig. 4). In wild-type (hf+/hf+) one-day-old chicks, the expression of the CYP19A1 coding region was only detected in hypothalamus (Fig. 4a). For Sebright chicks a very similar expression pattern was observed with the two sets of primers, with high ectopic expression in skin, lung, and intestine, and low expression in hypothalamus, muscle, and liver (Fig. 4a, b). The results strongly suggest that the Hf_ERV-induced transcript, starting at the 3’LTR TATA box, is the only major isoform responsible for ectopic expression of CYP19A1 in Hf chicks. The level of ectopic expression was higher in the homozygote (Hf/Hf) than in the heterozygote (Hf/hf+) as expected.
Segregation of the henny feathering phenotype in a backcross population
To further explore the genotype-phenotype association for henny feathering we analyzed segregation data in a backcross pedigree. Four F0 Sebright females homozygous for Hf_ERV (Hf/Hf) were mated with three Silver Spangled Hamburg males (hf +/hf +), which produced 17 F1 females. These were backcrossed to F0 males to produce 60 backcross males. The henny feathering phenotype was observed in 27 progenies, while 33 were wild-type, consistent with an expected 1:1 segregation (P = 0.44). A fifth F0 Sebright female was heterozygous for Hf_ERV and produced 4 F1 females, one (ID: 544) was heterozygous (Hf/hf+) while the other three (ID: 541, 545, 549) were homozygous wild-type (Fig. 5). The segregation of the henny feathering phenotype among their back-cross males was in perfect agreement with the genotype determined based on the Hf_ERV insertion. The result from this pedigree is fully consistent with the dominant inheritance of henny feathering and shows that Hf is not fixed in this Sebright population. We estimated the allele frequency of Hf to 0.80 in a sample of Sebright chickens (n = 20). One of the two Campine chickens that we genotyped was homozygous mutant (Hf/Hf) while the other was homozygous wild-type (hf+/hf+, Additional file 4: Table S4), indicating that Hf is not fixed in this breed either. None of these samples had phenotypic information because they were females or this specific phenotype was not recorded.
Hf was previously mapped to the linkage group E29C09W09 (cited by Carefoot ), which is now known to reside on chicken chromosome 10. Previous studies indicated that an LTR from a, putative ERV upstream of CYP19A1 might cause henny feathering by acting as a cis-regulatory mutation driving ectopic expression of aromatase [21, 27]. In our study, a 7524 bp ERV within the 5’UTR of CYP19A1 was completely associated with Hf. The genomic sequence documented in our study matches 83 bp of the cDNA sequence reported by Matsumine et al.  whereas 45 bp of the cDNA shows multiple sequence differences (Additional file 7: Figure S3), possibly due to sequencing errors in the previous sequence.
Our study demonstrates that ERVs in the chicken genome can have profound phenotypic effects. Other examples include a retroviral insertion in the tyrosinase gene causing the recessive white plumage color  and an EAV-HP insertion in the promoter region of SLCO1B3, which upregulates gene expression causing the blue eggshell phenotype [29, 30]. We found two identical 274 bp LTRs in the Hf_ERV, showing high sequence identity to avian leukosis virus strain ev-21, each containing a powerful transcriptional enhancer and core promoter with a TATA box (Fig. 2c). The 5’RACE experiment strongly suggested that the proximity of the promoter in the 3’LTR, and its binding sites for EFIII and Ig/EBP, to CYP19A1 (Fig. 2) are the drivers of ectopic expression of aromatase and thereby the conversion of androgen to estrogen in the skin and other non-ovarian tissues. Given the identical LTRs, it is conceivable that expression of ERV genes could be induced from the 5’LTR. Hf_ERV has intact gag, pro, pol, and env genes, as well as intact TSDs immediately flanking the two identical LTRs, indicating a recent insertion (Fig. 2b). The Hf insertion most likely occurred in one of the close ancestors of the “hen-tailed Bantam cock” found by Sir John Sebright in 1800 . It also suggests that the Hf_ERV could still be active and potentially replicate further. We observed that the Hf allele was not fixed neither in the Sebright nor in the Campine. This is probably due to reduced fertility of males that are homozygous for Hf [17, 31]. Male heterozygotes produce more offspring, which thus reduces selection against the hf+ allele. The reduced male fertility is likely caused by the ectopic expression of aromatase in testis [17, 32], or possibly expression of Hf_ERV which could reduce fertility and hatchability of the host chicken . Therefore, functional studies are needed to determine the direct effects of the Hf_ERV insertion.
Two isoforms of CYP19A1 transcripts have been documented in wild-type chickens, NM_001364699.1 and NM_001001761.3. Their coding regions are identical, however NM_001364699.1 is longer because of an extended 5’UTR (Fig. 2a). Although both transcripts are expressed in cerebrum, only the shorter transcript is detected in ovary , which is consistent with previous reports that aromatase is detectable only in ovary and hypothalamus in wild-type chickens . In both Hf/− and hf+/hf+ chickens, the shorter (NM_001001761.3) isoform can be detected in the ovary but not in fibroblasts . However, because a transcript (MN124451) corresponding to NM_001364699.1, with an extended UTR and partial LTR sequence (Fig. 2a), can be detected in both ovary and fibroblasts from Hf/− chickens , it is most likely responsible for the aromatase activity in all extragonadal tissues. This is supported by our qPCR result showing that Hf_ERV-induced transcripts are detected at high levels in skin, lung, and intestine, but also at low level in liver, muscle, and hypothalamus in Sebright chickens. The observed CYP19A1 expression pattern in Sebright chicks (skin > lung > intestine > > muscle > liver; Fig. 4) matches previous reports about estrogen formation  and aromatase activity  in Sebright chickens. We propose that the effect of the Hf mutation is changing the expression pattern of the long isoform, represented by NM_001364699.1, from hypothalamus-specific to more broadly expressed and possibly without affecting the ovary-specific transcript.
This study shows that the henny feathering trait in chicken is associated with the insertion of an intact avian leukosis virus in the 5’UTR of CYP19A1. The insertion of strong viral promotors in this region appears as a plausible causal mutation for the ectopic expression of CYP19A1 underlying henny feathering.
A trait mapping population was initiated from 8 chickens purchased from Murray McMurray Hatchery (www.mcmurrayhatchery.com, Webster City, Iowa, USA) to examine the segregation of Hf. They consisted of 3 Silver Spangled Hamburg males and 5 Silver Sebright females crossed to generate the F1 generation. Matings between 21 F1 females and 3 F0 Hamburg males produced 86 backcross males. Photos for phenotyping the backcrosses were taken (focusing on the hackle and saddle feathers which show the most distinct sexual dimorphism) at 12 weeks. The putative ERV insertion associated with Hf was characterized using 18 DNA samples, 16 Sebright and 2 Campine chickens, and 17 of these were sequenced using Nanopore long reads (Additional file 4: Table S4).
Liver tissues for DNA preparation and tissue samples (liver, lung, hypothalamus, muscle, intestine, and dorsal skin) for RNA isolation were collected from two one-day-old Silver Sebright chicks purchased from Ideal Poultry (http://www.idealpoultry.com/, Cameron, Texas, USA). The same tissues for RNA samples were collected from two one-day-old red junglefowls from a colony of red junglefowl kept in Texas A&M University’s Poultry Research Center.
Whole genome sequence (WGS) analyses
Illumina paired-end FASTQ data for 82 individuals or pooled samples from public databases, including one Sebright pool (Additional file 1: Table S1), aligned to the red junglefowl genome assembly version GalGal6 using BWA, sorted with SAMtools, and variants were called with GATK HaplotypeCaller 3.8 according to Broad best practices . Structural variants were called with Lumpy in single-sample mode .
Individual DNA samples were isolated from blood or liver using Qiagen, Puregene Tissue Core Kit B, DNA Isolation Protocol for Avian Blood with minor modifications. Each of the backcrossed males, F1 females, F0 and 18 chickens from Sebright or Campine lines (Additional file 4: Table S4) were genotyped by PCR (standard protocol for TAKARA PrimeSTAR GXL DNA Polymerase) with the forward primer HFEV_F and reverse primer HFEV_R2 (Additional file 3: Table S3).
Amplicon sequencing and assembly
Purified amplicons, encompassing the entire Hf insertion from 17 chickens, were barcoded using the Oxford Nanopore (ONT) Rapid Barcoding Kit (#SQK-RBK004) and sequenced on a single R9.5.1 flowcell. Raw nanopore FAST5 reads were converted to FASTQ using Albacore v2.3.4 (ONT), and assembled de novo using Canu 1.8 . An 85% majority consensus of the Hf insertion was derived from all 17 individuals, and the amplicon from one individual was manually polished by Sanger sequencing. Sanger primers are listed in Additional file 3: Table S3. PCR products for Sanger sequencing were generated via two-step nested PCR, which used purified amplicon from the genotyping PCR as the template to avoid amplification of homologous regions in the chicken genome. The polished sequence was used in BLAST searches to identify the most similar sequences, avian leukosis viruses. The consensus was aligned with 140 published avian leukosis viruses (Additional file 2: Table S2), and a maximum likelihood phylogeny was generated using RAxML 8.2.12 using the GTR + gamma nucleotide substitution model with 100 bootstraps . Annotation of the sequence was performed with RepeatMasker for LTRs and the EBI-EMBL Pfam database for retroviral protein coding sequences. Retroviral sequence features were characterized by RetroTector .
5′ rapid amplification of cDNA ends (5’RACE)
To determine the 5′ ends of the ectopically expressed CYP19A1 transcript in the skin of henny feathering chicken, 5’RACE experiments were performed with the kit Rapid Amplification of cDNA Ends (Invitrogen), according to the manufacturer’s protocol. The CYP19A1 gene-specific primers were designed for cDNA synthesis (Hf_RACE_R1, Additional file 3: Table S3) and for subsequent PCR reactions and Sanger sequencing (Hf_RACE_R2, Additional file 3: Table S3).
Quantitative real-time RT-PCR
Total RNA was extracted using Quick-RNA Miniprep Plus Kit (Zymo Research). First-strand cDNA was synthesized using SuperScript™ IV VILO™ Master Mix (Invitrogen). qPCRs were conducted with PowerUp™ SYBR™ Green Master Mix (Applied Biosystems) according to the manufacturer’s protocol. The products detected with the Roche LightCycler®480 using the standard protocol. Sequences of primers for the Hf_ERV induced transcript (Hf_qF and Hf_qR), the coding region of CYP19A1 (CYP_qF and CYP_qR) and for house-keeping gene (GAPDH_qF and GAPDH_qR) are listed in Additional file 3: Table S3. Each PCR reaction has three technical replicates.
Availability of data and materials
The annotated consensus sequence for Hf_ERV has been deposited in GenBank under the accession number MK937054.
The 5′ partial sequence for Hf_ERV-induced CYP19A1 transcripts, based on 5’RACE, has been deposited in GenBank under the accession number MN124451.
The accession numbers of the 140 published avian leukosis virus sequences used in this article are included within Additional file 2: Table S2.
- CYP19A1 :
Cytochrome P450 family 19 subfamily A member 1
Enhancer factor III
- Env :
- Gag :
Group specific antigen
Immunoglobulin/enhancer binding protein
Long terminal repeat
- Pol :
- Pro :
Rapid amplification of cDNA ends
Target site duplications
Whole genome sequence
Kimball RT, Ligon JD. Evolution of avian plumage dichromatism from a proximate perspective. Am Nat. 1999;154:182–93.
Dunn PO, Armenta JK, Whittingham LA. Natural and sexual selection act on different axes of variation in avian plumage color. Sci Adv. 2015;1:e1400155.
Morgan TH. The genetic factor for hen-feathering in the Sebright bantam. Biol Bull. 1920;39:257–9.
Wilson JD, George FW, Leshin M. Genetic control of extraglandular aromatase activity in the chicken. Steroids. 1987;50:235–44.
Tegetmeier WB. The poultry book: comprising the breeding and management of profitable and ornamental poultry; to which is added “the standard of excellence in exhibition birds.” London: G. Routledge; 1873. p. 284.
George FW, Matsumine H, Mcphaul MJ, Somes RG, Wilson JD. Inheritance of the henny feathering trait in the Golden Campine chicken: evidence for allelism with the gene that causes henny feathering in the sebright bantam. J Hered. 1990;81:107–10.
Eliot TS. The influence of the gonads on the plumage of Sebright bantams. Physiol Zool. 1928;1:286–324.
Owens IPF, Short RV. Hormonal basis of sexual dimorphism in birds: implications for new theories of sexual selection. Trends Ecol Evol. 1995;10:44–7.
Elbrecht A, Smith RG. Aromatase enzyme activity and sex determination in chickens. Science. 1992;255:467–70.
Morgan TH. Demonstration of the effects of castration on Seabright cockerels. Exp Biol Med. 1917;15:3–4.
Morgan TH. The effects of castration of hen-feathered Campines. Biol Bull. 1920;39:231–47.
Gallagher TF, Domm LV, Koch FC. The problem of hen-feathering in Sebright cocks. J Biol Chem. 1933;100:47.
George FW, Noble JF, Wilson JD. Female feathering in Sebright cocks is due to conversion of testosterone to estradiol in skin. Science. 1981;213:557–9.
Roxas HA. Gonad cross-transplantation in Sebright and Leghorn fowls. J Exp Zool. 1926;46:63–119.
Danforth CH. Gene H and testosterone in the fowl. In: Essays in biology, in honor of Herbert M Evans. Los Angeles: University of California Press; 1943. p. 157–167.
Leshin M, Baron J, George FW, Wilson JD. Increased estrogen formation and aromatase activity in fibroblasts cultured from the skin of chickens with the henny feathering trait. J Biol Chem. 1981;256:4341–4.
George FW, Wilson JD. Pathogenesis of the henny feathering trait in the sebright bantam chicken. J Clin Invest. 1980;66:57–65.
Somes RG, George FW, Baron J, Noble JF, Wilson JD. Inheritance of the henny-feathering trait of the Sebright bantam chicken. J Hered. 1984;75:99–102.
Wilson JD, Leshin M, George FW. The sebright bantam chicken and the genetic control of extraglandular aromatase. Endocr Rev. 1987;8:363–76.
Matsumine H, Wilson JD, McPhaul MJ. Sebright and Campine chickens express aromatase P-450 messenger RNA inappropriately in extraglandular tissues and in skin fibroblasts. Mol Endocrinol. 1990;4:905–11.
Matsumine H, Herbst MA, Ou SHI, Wilson JD, McPhaul MJ. Aromatase mRNA in the extragonadal tissues of chickens with the henny-feathering trait is derived from a distinctive promoter structure that contains a segment of a retroviral long terminal repeat. J Biol Chem. 1991;266:19900–7.
Lambeth LS, Morris KR, Wise TG, Cummins DM, O’Neil TE, Cao Y, et al. Transgenic chickens overexpressing aromatase have high estrogen levels but maintain a predominantly male phenotype. Endocrinology. 2016;157:83–90.
Abinawanto, Shimada K, Yoshida K, Saito N. Effects of aromatase inhibitor on sex differentiation and levels of P450 (17 alpha) and P450 arom messenger ribonucleic acid of gonads in chicken embryos. Gen Comp Endocrinol. 1996;102:241–6.
Ju G, Cullen BR. The role of avian retroviral LTRs in the regulation of gene expression and viral replication. Adv Virus Res. 1985;30:179–223.
Ruddell A. Transcription regulatory elements of the avian retroviral long terminal repeat. Virology. 1995;206:1–7.
Carefoot WC. Hen-feathering mutation HF*H may act as a eumelanising factor and modify the expression of autosomal barring. Br Poult Sci. 2002;43:391–4.
Leshin M. 5-Azacytidine and sodium butyrate induce expression of aromatase in fibroblasts from chickens carrying the henny feathering trait but not from wild-type chickens. Proc Natl Acad Sci. 1985;82:3005–9.
Chang CM, Coville JL, Coquerelle G, Gourichon D, Oulmouden A, Tixier-Boichard M. Complete association between a retroviral insertion in the tyrosinase gene and the recessive white mutation in chickens. BMC Genomics. 2006;7:19.
Wang Z, Qu L, Yao J, Yang X, Li G, Zhang Y, et al. An EAV-HP insertion in 5′ flanking region of SLCO1B3 causes blue eggshell in the chicken. PLoS Genet. 2013;9:1003183.
Hanotte O, Han J-L, Alcalde JA, Wragg D, Mwacharo JM, Gourichon D, et al. Endogenous retrovirus EAV-HP linked to blue egg phenotype in Mapuche fowl. Kozak CA, editor. PLoS One. 2013;8:e71393.
Hadden DR, Mayes FJ, Bayne BEF, Finlay L, Leslie H, Sheridan B. Elevation of serum estradiol in hen-feathered male sebright bantams. J Endocrinol. 1986;108:84.
Fowler KA, Gill K, Kirma N, Dillehay DL, Tekmal RR. Overexpression of aromatase leads to development of testicular Leydig cell tumors: an in vivo model for hormone-mediated testicular cancer. Am J Pathol. 2000;156:347–53.
Bacon LD, Fulton JE, Kulkarni GB. Methods for evaluating and developing commercial chicken strains free of endogenous subgroup E avian leukosis virus. Avian Pathol. 2004;33:233–43.
Schmid M, Smith J, Burt DW, Aken BL, Antin PB, Archibald AL, et al. Third report on chicken genes and chromosomes 2015. Cytogenet Genome Res. 2015;145:78–179.
George FW, Wilson JD. Developmental pattern of increased aromatase activity in the Sebright bantam chicken. Endocrinology. 1982;110:1203–7.
Poplin R, Ruano-Rubio V, DePristo MA, Fennell TJ, Carneiro MO, Auwera GA Van der, et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv. 2017;201178.
Layer RM, Chiang C, Quinlan AR, Hall IM. LUMPY: a probabilistic framework for structural variant discovery. Genome Biol. 2014;15:R84.
Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017;27:722–36.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
Sperber GO, Airola T, Jern P, Blomberg J. Automated recognition of retroviral sequences in genomic data - RetroTector©. Nucleic Acids Res. 2007;35:4964–76.
We thank members in Dr. Paul Siegel’s laboratory for husbandry and sampling of the chickens, and for DNA preparations. We thank Dr. Giri Athrey for providing the red junglefowls for RNA samples.
The work was supported by The Hagler Institute for Advanced Study at Texas A&M University.
Procedures for the mapping population were approved by the Institutional Animal Care and Use Committee at Virginia Polytechnic Institute and State University. Procedures for the tissue sample collections followed by the guidelines of Institutional Animal Care and Use Committee at Texas A&M University College of Veterinary Medicine & Biomedical Sciences.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Public whole genome sequence data used in this study. (XLSX 11 kb)
Table S2. Sequence comparison of the henny feathering endogenous retrovirus sequence (Hf_ERV) with 140 published avian leukosis viruses. (XLSX 16 kb)
Table S3. PCR primers used in this study. (XLSX 10 kb)
Table S4. List of DNA samples. (XLSX 9 kb)
Figure S1. Aligned reads showing the insertion site of Hf_ERV on chicken chromosome 10. Aligned reads are shown in gray; multicolored “rainbow” reads are soft-clipped. Left clip and right clip sequences match the 5′ and 3′ of the LTR sequence, respectively. Discordant reads are in cyan which are well aligned in this figure (Chr10), while their mate pairs are aligned to Chr1, either to LOC770705 or LOC107052718, both retrovirus related genes, because there is no retrovirus sequence at this target site on chromosome 10 in the reference genome. Altogether, this alignment suggests a retrovirus insertion located between the left and right clips, which contains the 6 bp target site duplication (TSD) sequence. Reads are obtained from BioSample Accession SAMEA104432213 (Pool of Silver Sebright Bantam) and aligned to GalGal6. (DOCX 214 kb)
Figure S2. Maximum likelihood phylogeny for the endogenous retrovirus associated with henny feathering (Hf_ERV). Hf_ERV and the 50 other retrovirus sequences with above 90% sequence identity to Hf_ERV. Bootstrap values are reported as percentages. Grey box indicates the clade containing Hf_ERV, ev-1 and ev21. (DOCX 475 kb)
Figure S3. Comparison of sequences reported in this and previous studies. Alignment between a previously reported partial cDNA sequence of CYP19A1 in Sebright chickens  and the retrovirus insertion detected in the present study (LTR). The base pair positions of the LTR sequence correspond to the complete LTR (5′ to 3′), which is 274 bp in length. (DOCX 140 kb)