Transcriptome-wide investigation of genomic imprinting in chicken

United by conflict: Convergent signatures of parental conflict in angiosperms and placental mammals

Endosperm in angiosperms and placenta in eutherians are convergent innovations for efficient embryonic nutrient transfer. Despite advantages, this reproductive strategy incurs metabolic costs that maternal parents disproportionately shoulder, leading to potential inter-parental conflict over optimal offspring investment. Genomic imprinting-parent-of-origin-biased gene expression-is fundamental for endosperm and placenta development and has convergently evolved in angiosperms and mammals, in part, to resolve parental conflict. Here, we review the mechanisms of genomic imprinting in these taxa. Despite differences in the timing and spatial extent of imprinting, these taxa exhibit remarkable convergence in the molecular machinery and genes governing imprinting. We then assess the role of parental conflict in shaping evolution within angiosperms and eutherians using four criteria: 1) Do differences in the extent of sibling relatedness cause differences in the inferred strength of parental conflict? 2) Do reciprocal crosses between taxa with different inferred histories of parental conflict exhibit parent-of-origin growth effects? 3) Are these parent-of-origin growth effects caused by dosage-sensitive mechanisms and do these loci exhibit signals of positive selection? 4) Can normal development be restored by genomic perturbations that restore stoichiometric balance in the endosperm/placenta? Although we find evidence for all criteria in angiosperms and eutherians, suggesting that parental conflict may help shape their evolution, many questions remain. Additionally, myriad differences between the two taxa suggest that their respective biologies may shape how/when/where/to what extent parental conflict manifests. Lastly, we discuss outstanding questions, highlighting the power of comparative work in quantifying the role of parental conflict in evolution.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Reinforcement of repressive marks in the chicken primordial germ cell epigenetic signature: divergence from basal state resetting in mammals

BACKGROUND

In mammals, primordial germ cells (PGCs), the embryonic precursors of the germline, arise from embryonic or extra-embryonic cells upon induction by the surrounding tissues during gastrulation, according to mechanisms which are elucidated in mice but remain controversial in primates. They undergo genome-wide epigenetic reprogramming, consisting of extensive DNA demethylation and histone post-translational modification (PTM) changes, toward a basal, euchromatinized state. In contrast, chicken PGCs are specified by preformation before gastrulation based on maternally-inherited factors. They can be isolated from the bloodstream during their migration to the genital ridges. Our prior research highlighted differences in the global epigenetic profile of cultured chicken PGCs compared with chicken somatic cells and mammalian PGCs. This study investigates the acquisition and evolution of this profile during development.

RESULTS

Quantitative analysis of global DNA methylation and histone PTMs, including their distribution, during key stages of chicken early development revealed divergent PGC epigenetic changes compared with mammals. Unlike mammalian PGCs, chicken PGCs do not undergo genome-wide DNA demethylation or exhibit a decrease in histone H3 lysine 9 dimethylation. However, chicken PGCs show 5‑hydroxymethylcytosine loss, macroH2A redistribution, and chromatin decompaction, mirroring mammalian processes. Chicken PGCs initiate their epigenetic signature during migration, progressively accumulating high global levels of H3K9me3, with preferential enrichment in inactive genome regions. Despite apparent global chromatin decompaction, abundant heterochromatin marks, including repressive histone PTMs, HP1 variants, and DNA methylation, persists in chicken PGCs, contrasting with mammalian PGCs.

CONCLUSIONS

Chicken PGCs' epigenetic signature does not align with the basal chromatin state observed in mammals, suggesting a departure from extensive epigenetic reprogramming. Despite disparities in early PGC development, the persistence of several epigenetic features shared with mammals implies their involvement in chromatin-regulated germ cell properties, with the distinctive elevation of chicken-specific H3K9me3 potentially participating in these processes.

Livestock species as emerging models for genomic imprinting

Genomic imprinting is an epigenetically-regulated process of central importance in mammalian development and evolution. It involves multiple levels of regulation, with spatio-temporal heterogeneity, leading to the context-dependent and parent-of-origin specific expression of a small fraction of the genome. Genomic imprinting studies have therefore been essential to increase basic knowledge in functional genomics, evolution biology and developmental biology, as well as with regard to potential clinical and agrigenomic perspectives. Here we offer an overview on the contribution of livestock research, which features attractive resources in several respects, for better understanding genomic imprinting and its functional impacts. Given the related broad implications and complexity, we promote the use of such resources for studying genomic imprinting in a holistic and integrative view. We hope this mini-review will draw attention to the relevance of livestock genomic imprinting studies and stimulate research in this area.

Avian sperm-borne RNAs: optimisation of a new isolation protocol

1. Sperm-borne RNAs are involved in sperm and embryonic protein translation, the regulation of early development and the epigenetic inheritance of the paternal phenotype. Sperm-borne RNA purification protocols generally include a cell purification stage to discard contamination by somatic cells. In avian species, no protocol is currently available to isolate all the populations composing sperm-borne RNAs.2. This study evaluated the presence of somatic cells in semen samples of chickens and quails using visual examination after fluorescent nuclei staining. The efficiency of somatic cell lysis buffer (SCLB) on chicken liver cells and its impacts on chicken sperm cell integrity was explored. Three different approaches were tested to isolate RNA: two developed for mammalian sperm cells and a commercial kit for somatic cells. The efficiency and reliability of each approach was determined based on RNA quality and purity. Eventually, the presence of miRNA and mRNA in purified avian sperm-borne RNAs was investigated by RT-(q)PCR.3. No somatic cells were found in chicken and quail semen. The SCLB totally lysed chicken liver cells but also induced sperm cell necrosis. Consequently, this treatment wasn't performed on samples prior to RNA isolation. Among the tested RNA purification protocols, the commercial one was the least variable and isolated RNA with the highest purity levels. No DNA contamination was observed. Furthermore, the samples contained miRNA and mRNA already known as present in mammalian sperm cells (gga-miR-100-5p, gga-miR-191-5p, GAPDH and PLCZ1), but mRNAs associated with leucocytes (CD4) and Sertoli cells (SOX4, CLDN11) were not detected. This protocol was successfully applied to quail sperm cells.4. Altogether, the study reveals that it is unnecessary to pre-treat samples to remove somatic cell contamination before RNA purification and successfully describes an isolation protocol for sperm-borne RNAs, including small non-coding and long coding RNAs, in two distinct avian species highly valuable as biological models.

Placental imprinting of SLC22A3 in the IGF2R imprinted domain is conserved in therian mammals

Background

The eutherian IGF2R imprinted domain is regulated by an antisense long non-coding RNA, Airn, which is expressed from a differentially methylated region (DMR) in mice. Airn silences two neighbouring genes, Solute carrier family 22 member 2 (Slc22a2) and Slc22a3, to establish the Igf2r imprinted domain in the mouse placenta. Marsupials also have an antisense non-coding RNA, ALID, expressed from a DMR, although the exact function of ALID is currently unknown. The eutherian IGF2R DMR is located in intron 2, while the marsupial IGF2R DMR is located in intron 12, but it is not yet known whether the adjacent genes SLC22A2 and/or SLC22A3 are also imprinted in the marsupial lineage. In this study, the imprinting status of marsupial SLC22A2 and SLC22A3 in the IGF2R imprinted domain in the chorio-vitelline placenta was examined in a marsupial, the tammar wallaby.

Results

In the tammar placenta, SLC22A3 but not SLC22A2 was imprinted. Tammar SLC22A3 imprinting was evident in placental tissues but not in the other tissues examined in this study. A putative promoter of SLC22A3 lacked DNA methylation, suggesting that this gene is not directly silenced by a DMR on its promoter as seen in the mouse. Based on immunofluorescence, we confirmed that the tammar SLC22A3 is localised in the endodermal cell layer of the tammar placenta where nutrient trafficking occurs.

Conclusions

Since SLC22A3 is imprinted in the tammar placenta, we conclude that this placental imprinting of SLC22A3 has been positively selected after the marsupial and eutherian split because of the differences in the DMR location. Since SLC22A3 is known to act as a transporter molecule for nutrient transfer in the eutherian placenta, we suggest it was strongly selected to control the balance between supply and demand of nutrients in marsupial as it does in eutherian placentas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00465-4.

Evolution of the short form of DNMT3A, DNMT3A2, occurred in the common ancestor of mammals

Genomic imprinting is found in marsupial and eutherian mammals, but not in monotremes. While the primary regulator of genomic imprinting in eutherians is differential DNA methylation between parental alleles, conserved imprinted genes in marsupials tend to lack DNA methylation at their promoters. DNA methylation at eutherian imprinted genes is mainly catalysed by a DNA methyltransferase (DNMT) enzyme, DNMT3A. There are two isoforms of eutherian DNMT3A: DNMT3A and DNMT3A2. DNMT3A2 is the primary isoform for establishing DNA methylation at eutherian imprinted genes and is essential for eutherian genomic imprinting. In this study, we investigated whether DNMT3A2 is also present in the two other mammalian lineages, marsupials and monotremes. We identified DNMT3A2 in both marsupials and monotremes, although imprinting has not been identified in monotremes. By analysing genomic sequences and transcriptome data across vertebrates, we concluded that the evolution of DNMT3A2 occurred in the common ancestor of mammals. In addition, DNMT3A/3A2 gene and protein expression during gametogenesis showed distinct sexual dimorphisms in a marsupial, the tammar wallaby, and this pattern coincided with the sex-specific DNA methylation reprogramming in this species as it does in mice. Our results show that DNMT3A2 is present in all mammalian groups and suggests that the basic DNMT3A/3A2-based DNA methylation mechanism is conserved at least in therian mammals.

Genomic Imprinting in the New Omics Era: A Model for Systems-Level Approaches

Genomic imprinting represents a noteworthy inheritance mechanism leading to allele-specific regulations dependent of the parental origin. Imprinted loci are especially involved in essential mammalian functions related to growth, development and behavior. In this mini-review, we first offer a summary of current representations associated with genomic imprinting through key results of the three last decades. We then outline new perspectives allowed by the spread of new omics technologies tackling various interacting levels of imprinting regulations, including genomics, transcriptomics and epigenomics. We finally discuss the expected contribution of new omics data to unresolved big questions in the field.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

Whole-genome resequencing reveals loci with allelic transmission ratio distortion in F1 chicken population

Allelic transmission ratio distortion (TRD) is the significant deviation from the expected ratio under Mendelian inheritance theory, which may be resulted from multiple disrupted biological processes, including germline selection, meiotic drive, gametic competition, imprint error, and embryo lethality. However, it is less known that whether or what extent the allelic TRD is present in farm animals. In this study, whole-genome resequencing technology was applied to reveal TRD loci in chicken by constructing a full-sib F₁ hybrid population. Through the whole-genome resequencing data of two parents (30 ×) and 38 offspring (5 ×), we detected a total of 2850 TRD SNPs (p-adj < 0.05) located within 400 genes showing TRD, and all of them were unevenly distributed on macrochromosomes and microchromosomes. Our findings suggested that TRD in the chicken chromosome 16 might play an important role in chicken immunity and disease resistance and the MYH1F with significant TRD and allele-specific expression could play a key role in the fast muscle development. In addition, functional enrichment analyses revealed that many genes (e.g., TGFBR2, TGFBR3, NOTCH1, and NCOA1) with TRD were found in the significantly enriched biological process and InterPro terms in relation to embryonic lethality and germline selection. Our results suggested that TRD is considerably prevalent in the chicken genome and has functional implications.

Epigenetics of Animal Personality: DNA Methylation Cannot Explain the Heritability of Exploratory Behavior in a Songbird

The search for the hereditary mechanisms underlying quantitative traits traditionally focused on the identification of underlying genomic polymorphisms such as single-nucleotide polymorphisms. It has now become clear that epigenetic mechanisms, such as DNA methylation, can consistently alter gene expression over multiple generations. It is unclear, however, if and how DNA methylation can stably be transferred from one generation to the next and can thereby be a component of the heritable variation of a trait. In this study, we explore whether DNA methylation responds to phenotypic selection using whole-genome and genome-wide bisulfite approaches. We assessed differential erythrocyte DNA methylation patterns between extreme personality types in the Great Tit (Parus major). For this, we used individuals from a four-generation artificial bi-directional selection experiment and siblings from eight F2 inter-cross families. We find no differentially methylated sites when comparing the selected personality lines, providing no evidence for the so-called epialleles associated with exploratory behavior. Using a pair-wise sibling design in the F2 intercrosses, we show that the genome-wide DNA methylation profiles of individuals are mainly explained by family structure, indicating that the majority of variation in DNA methylation in CpG sites between individuals can be explained by genetic differences. Although we found some candidates explaining behavioral differences between F2 siblings, we could not confirm this with a whole-genome approach, thereby confirming the absence of epialleles in these F2 intercrosses. We conclude that while epigenetic variation may underlie phenotypic variation in behavioral traits, we were not able to find evidence that DNA methylation can explain heritable variation in personality traits in Great Tits.

Tissue-specific patterns of regulatory changes underlying gene expression differences among Ficedula flycatchers and their naturally occurring F1 hybrids

Changes in interacting cis- and trans-regulatory elements are important candidates for Dobzhansky-Muller hybrid incompatibilities and may contribute to hybrid dysfunction by giving rise to misexpression in hybrids. To gain insight into the molecular mechanisms and determinants of gene expression evolution in natural populations, we analyzed the transcriptome from multiple tissues of two recently diverged Ficedula flycatcher species and their naturally occurring F₁ hybrids. Differential gene expression analysis revealed that the extent of differentiation between species and the set of differentially expressed genes varied across tissues. Common to all tissues, a higher proportion of Z-linked genes than autosomal genes showed differential expression, providing evidence for a fast-Z effect. We further found clear signatures of hybrid misexpression in brain, heart, kidney, and liver. However, while testis showed the highest divergence of gene expression among tissues, it showed no clear signature of misexpression in F₁ hybrids, even though these hybrids were found to be sterile. It is therefore unlikely that incompatibilities between cis-trans regulatory changes explain the observed sterility. Instead, we found evidence that cis-regulatory changes play a significant role in the evolution of gene expression in testis, which illustrates the tissue-specific nature of cis-regulatory evolution bypassing constraints associated with pleiotropic effects of genes.

The Landscape of Genomic Imprinting at the Porcine SGCE/PEG10 Locus from Methylome and Transcriptome of Parthenogenetic Embryos

In mammals, imprinted genes often exist in the form of clusters in specific chromosome regions. However, in pigs, genomic imprinting of a relatively few genes and clusters has been identified, and genes within or adjacent to putative imprinted clusters need to be investigated including those at the SGCE/PEG10 locus. The objective of this study was to, using porcine parthenogenetic embryos, investigate imprinting status of genes within the genomic region spans between the COL1A2 and ASB4 genes in chromosome 9. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) were conducted with normal and parthenogenetic embryos, and methylome and transcriptome were analyzed. As a result, differentially methylated regions (DMRs) between the embryos were identified, and parental allele-specific expressions of the SGCE and PEG10 genes were verified. The pig imprinted interval was limited between SGCE and PEG10, since both the COL1A2 and CASD1 genes at the centromere-proximal region and the genes between PPP1R9A and ASB4 toward the telomere were non-imprinted and biallelically expressed. Consequently, our combining analyses of methylome, transcriptome, and informative polymorphisms revealed the boundary of imprinting cluster at the SGCE/PEG10 locus in pig chromosome 9 and consolidated the landscape of genomic imprinting in pigs.

Transcriptome analysis revealed misregulated gene expression in blastoderms of interspecific chicken and Japanese quail F1 hybrids

Hybrid incompatibility, such as sterility and inviability, prevents gene flow between closely-related populations as a reproductive isolation barrier. F₁ hybrids between chickens and Japanese quail (hereafter, referred to as quail), exhibit a high frequency of developmental arrest at the preprimitive streak stage. To investigate the molecular basis of the developmental arrest at the preprimitive streak stage in chicken–quail F₁ hybrid embryos, we investigated chromosomal abnormalities in the hybrid embryos using molecular cytogenetic analysis. In addition, we quantified gene expression in parental species and chicken- and quail-derived allele-specific expression in the hybrids at the early blastoderm and preprimitive streak stages by mRNA sequencing. Subsequently, we compared the directions of change in gene expression, including upregulation, downregulation, or no change, from the early blastoderm stage to the preprimitive streak stage between parental species and their hybrids. Chromosome analysis revealed that the cells of the hybrid embryos contained a fifty-fifty mixture of parental chromosomes, and numerical chromosomal abnormalities were hardly observed in the hybrid cells. Gene expression analysis revealed that a part of the genes that were upregulated from the early blastoderm stage to the preprimitive streak stage in both parental species exhibited no upregulation of both chicken- and quail-derived alleles in the hybrids. GO term enrichment analysis revealed that these misregulated genes are involved in various biological processes, including ribosome-mediated protein synthesis and cell proliferation. Furthermore, the misregulated genes included genes involved in early embryonic development, such as primitive streak formation and gastrulation. These results suggest that numerical chromosomal abnormalities due to a segregation failure does not cause the lethality of chicken–quail hybrid embryos, and that the downregulated expression of the genes that are involved in various biological processes, including translation and primitive streak formation, mainly causes the developmental arrest at the preprimitive streak stage in the hybrids.

Imprinting analysis by droplet digital PCR coupled with locked nucleic acid TaqMan probes

Imprinted genes are differentially expressed in a parent-of-origin-specific manner. Parental origin of the alleles is discriminated by intragenic DNA polymorphisms. Comparisons of parental allelic expression have been analysed by semiquantitative RT-PCR. Here, we developed a novel quantitative method for allelic expression of the imprinted gene Ube3a, which inactivation and mutations cause Angelman syndrome and predominantly expressed by the maternal allele in neuronal tissues. In this method, cDNA was amplified by droplet digital PCR (ddPCR) coupled with allele-specific locked nucleic acid (LNA) TaqMan probes, which labelled by FAM and HEX were designed to detect the SNPs in the target regions. ddPCR assay demonstrated that the sense transcript of Ube3a was equally expressed from both parental alleles in adult tissues except neuronal tissues, where Ube3a expression from the paternal allele was about 10 to 14% of total Ube3a expression in adult brain, and 20% in spinal cord. The antisense transcript of Ube3a was expressed at 60% to 70% of the sense transcript of Ube3a in adult brain. Changes in the Ube3a transcripts during postnatal brain development were also evaluated by ddPCR. The ddPCR method is far more reliable and simpler to use than semiquantitative PCR to analyse skewed or faint allelic expression of imprinted genes.

Zygotic genome activation in the chicken: a comparative review

Maternal RNAs and proteins in the oocyte contribute to early embryonic development. After fertilization, these maternal factors are cleared and embryonic development is determined by an individual's own RNAs and proteins, in a process called the maternal-to-zygotic transition. Zygotic transcription is initially inactive, but is eventually activated by maternal transcription factors. The timing and molecular mechanisms involved in zygotic genome activation (ZGA) have been well-described in many species. Among birds, a transcriptome-based understanding of ZGA has only been explored in chickens by RNA sequencing of intrauterine embryos. RNA sequencing of chicken intrauterine embryos, including oocytes, zygotes, and Eyal-Giladi and Kochav (EGK) stages I-X has enabled the identification of differentially expressed genes between consecutive stages. These studies have revealed that there are two waves of ZGA: a minor wave at the one-cell stage (shortly after fertilization) and a major wave between EGK.III and EGK.VI (during cellularization). In the chicken, the maternal genome is activated during minor ZGA and the paternal genome is quiescent until major ZGA to avoid transcription from supernumerary sperm nuclei. In this review, we provide a detailed overview of events in intrauterine embryonic development in birds (and particularly in chickens), as well as a transcriptome-based analysis of ZGA.

Genome-wide analysis of spatiotemporal allele-specific expression in F1 hybrids of meat- and egg-type chickens

In diploid organisms, each gene locus is composed of two parental alleles, which would interact with each other for determining the phenotypic variation. Better understanding of the allele-specific expression (ASE) in farm animals is much important to explore the genetic basis underlying economically important traits, which have been poorly understood yet. In this study, genome-wide analysis was applied to explore the spatiotemporal pattern of ASE in the F1 hybrids of chicken. First, meat- and egg-type chickens were selected for producing a full-sib F1 hybrid population (n = 57). Then, genome resequencing of two parents and 38 offspring were performed and liver and breast muscle samples (n = 38) were subjected to strand-specific RNA sequencing (ssRNA-seq) for ASE detection at 1, 28, and 56 days of age, respectively. The results accurately identified a total of 465 informative genes that could be distinguished with respect to their parental origins. There were 0.4% - 4.1% of informative genes showing ASE, and 57 of them were found across different tissues and time points. Besides, most ASE genes in chickens were tissue-specific, and no matter what the time-point pattern of one ASE gene, the same parental allele of this gene almost showed consistently higher or lower expression across all time points in the same type tissue. In conclusion, this study indicated that most of ASE genes were tissue-specific and time-dependent.

Evolution of cis- and trans-regulatory divergence in the chicken genome between two contrasting breeds analyzed using three tissue types at one-day-old

Background

Gene expression variation is a key underlying factor influencing phenotypic variation, and can occur via cis- or trans-regulation. To understand the role of cis- and trans-regulatory variation on population divergence in chicken, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibit major differences in body size and reproductive traits, and used them to determine the degree of cis versus trans variation in the brain, liver, and muscle tissue of male and female 1-day-old specimens.

Results

We provided an overview of how transcriptomes are regulated in hybrid progenies of two contrasting breeds based on allele specific expression analysis. Compared with cis-regulatory divergence, trans-acting genes were more extensive in the chicken genome. In addition, considerable compensatory cis- and trans-regulatory changes exist in the chicken genome. Most importantly, stronger purifying selection was observed on genes regulated by trans-variations than in genes regulated by the cis elements.

Conclusions

We present a pipeline to explore allele-specific expression in hybrid progenies of inbred lines without a specific reference genome. Our research is the first study to describe the regulatory divergence between two contrasting breeds. The results suggest that artificial selection associated with domestication in chicken could have acted more on trans-regulatory divergence than on cis-regulatory divergence.

Transcriptional regulatory divergence underpinning species-specific learned vocalization in songbirds

Learning of most motor skills is constrained in a species-specific manner. However, the proximate mechanisms underlying species-specific learned behaviors remain poorly understood. Songbirds acquire species-specific songs through learning, which is hypothesized to depend on species-specific patterns of gene expression in functionally specialized brain regions for vocal learning and production, called song nuclei. Here, we leveraged two closely related songbird species, zebra finch, owl finch, and their interspecific first-generation (F1) hybrids, to relate transcriptional regulatory divergence between species with the production of species-specific songs. We quantified genome-wide gene expression in both species and compared this with allele-specific expression in F1 hybrids to identify genes whose expression in song nuclei is regulated by species divergence in either cis- or trans-regulation. We found that divergence in transcriptional regulation altered the expression of approximately 10% of total transcribed genes and was linked to differential gene expression between the two species. Furthermore, trans-regulatory changes were more prevalent than cis-regulatory and were associated with synaptic formation and transmission in song nucleus RA, the avian analog of the mammalian laryngeal motor cortex. We identified brain-derived neurotrophic factor (BDNF) as an upstream mediator of trans-regulated genes in RA, with a significant correlation between individual variation in BDNF expression level and species-specific song phenotypes in F1 hybrids. This was supported by the fact that the pharmacological overactivation of BDNF receptors altered the expression of its trans-regulated genes in the RA, thus disrupting the learned song structures of adult zebra finch songs at the acoustic and sequence levels. These results demonstrate functional neurogenetic associations between divergence in region-specific transcriptional regulation and species-specific learned behaviors.

Genome-Wide Analysis of Allele-Specific Expression Patterns in Seventeen Tissues of Korean Cattle (Hanwoo)

The functional hemizygosity could be caused by the MAE of a given gene and it can be one of the sources to affect the phenotypic variation in cattle. We aimed to identify MAE genes across the transcriptome in Korean cattle (Hanwoo). For three Hanwoo family trios, the transcriptome data of 17 tissues were generated in three offspring. Sixty-two MAE genes had a monoallelic expression in at least one tissue. Comparing genotypes among each family trio, the preferred alleles of 18 genes were identified (maternal expression, n = 9; paternal expression, n = 9). The MAE genes are involved in gene regulation, metabolic processes, and immune responses, and in particular, six genes encode transcription factors (FOXD2, FOXM1, HTATSF1, SCRT1, NKX6-2, and UBN1) with tissue-specific expression. In this study, we report genome-wide MAE genes in seventeen tissues of adult cattle. These results could help to elucidate epigenetic effects on phenotypic variation in Hanwoo.

The Role of Epigenetics in Placental Development and the Etiology of Preeclampsia

In this review, we comprehensively present the function of epigenetic regulations in normal placental development as well as in a prominent disease of placental origin, preeclampsia (PE). We describe current progress concerning the impact of DNA methylation, non-coding RNA (with a special emphasis on long non-coding RNA (lncRNA) and microRNA (miRNA)) and more marginally histone post-translational modifications, in the processes leading to normal and abnormal placental function. We also explore the potential use of epigenetic marks circulating in the maternal blood flow as putative biomarkers able to prognosticate the onset of PE, as well as classifying it according to its severity. The correlation between epigenetic marks and impacts on gene expression is systematically evaluated for the different epigenetic marks analyzed.

Assessment of maternal and parent of origin effects in genetic variation of economic traits in Iranian native fowl

1. The objective of the study was to investigate the influence of maternal and parent of origin effects (POE) on genetic variation of Iranian native fowl on economic traits. 2. Studied traits were body weights at birth (BW0), at eight (BW8) and 12 weeks of age (BW12), age (ASM) and weight at sexual maturity (WSM), egg number (EN) and average egg weight (AEW). 3. Several models, including additive, maternal additive genetics, permanent environmental effects and POE were compared using Wombat software. Bayesian Information Criterion (BIC) was used to identify the best model for each trait. The chance of reranking of birds between models was investigated using Spearman correlation and Wilcoxon rank test. 4. Based on the best model, direct heritability estimates for BW0, BW8, BW12, ASM, WSM, EN and AEW traits were 0.05, 0.21, 0.23, 0.30, 0.39, 0.22 and 0.38, respectively. Proportion of variance due to paternal POE for BW8 was 4% and proportion of variance due to maternal POE for BW12 was 5%. 5. Estimated maternal heritability for BW0 was 0.30 and for BW8 and BW12 were 0.00 and 0.01, respectively, which shows that maternal heritability was reduced by age. 6. Based on the results, considering POE for BW8 and BW12 and maternal genetic effects for BW0 improved the accuracy of estimations and avoid reranking of birds for these traits.

Revisiting mechanisms and functions of prenatal hormone-mediated maternal effects using avian species as a model

Maternal effects can adaptively modulate offspring developmental trajectories in variable but predictable environments. Hormone synthesis is sensitive to environmental factors, and maternal hormones are thus a powerful mechanism to transfer environmental cues to the next generation. Birds have become a key model for the study of hormone-mediated maternal effects because the embryo develops outside the mother's body, facilitating the measurement and manipulation of prenatal hormone exposure. At the same time, birds are excellent models for the integration of both proximate and ultimate approaches, which is key to a better understanding of the evolution of hormone-mediated maternal effects. Over the past two decades, a surge of studies on hormone-mediated maternal effects has revealed an increasing number of discrepancies. In this review, we discuss the role of the environment, genetic factors and social interactions in causing these discrepancies and provide a framework to resolve them. We also explore the largely neglected role of the embryo in modulating the maternal signal, as well as costs and benefits of hormone transfer and expression for the different family members. We conclude by highlighting fruitful avenues for future research that have opened up thanks to new theoretical insights and technical advances in the field. This article is part of the theme issue 'Developing differences: early-life effects and evolutionary medicine'.

Scanning the genomes of parents for imprinted loci acting in their un-genotyped progeny

Depending on their parental origin, alleles at imprinted loci are fully or partially inactivated through epigenetic mechanisms. Their effects contribute to the broader class of parent-of-origin effects. Standard methodology for mapping imprinted quantitative trait loci in association studies requires phenotypes and parental origin of marker alleles (ordered genotypes) to be simultaneously known for each individual. As such, many phenotypes are known from un-genotyped offspring in ongoing breeding programmes (e.g. meat animals), while their parents have known genotypes but no phenotypes. By theoretical considerations and simulations, we showed that the limitations of standard methodology can be overcome in such situations. This is achieved by first estimating parent-of-origin effects, which then serve as dependent variables in association analyses, in which only imprinted loci give a signal. As a theoretical foundation, the regression of parent-of-origin effects on the number of B-alleles at a biallelic locus — representing the un-ordered genotype — equals the imprinting effect. The applicability to real data was demonstrated for about 1800 genotyped Brown Swiss bulls and their un-genotyped fattening progeny. Thus, this approach unlocks vast data resources in various species for imprinting analyses and offers valuable clues as to what extent imprinted loci contribute to genetic variability.

Investigation of allele-specific expression of genes involved in adipogenesis and lipid metabolism suggests complex regulatory mechanisms of PPARGC1A expression in porcine fat tissues

Background

The expression of genes involved in regulating adipogenesis and lipid metabolism may affect economically important fatness traits in pigs. Allele-specific expression (ASE) reflects imbalance between allelic transcript levels and can be used to identify underlying cis-regulatory elements. ASE has not yet been intensively studied in pigs. The aim of this investigation was to analyze the differential allelic expression of four genes, PPARA, PPARG, SREBF1, and PPARGC1A, which are involved in the regulation of fat deposition in porcine subcutaneous and visceral fat and longissimus dorsi muscle.

Results

Quantification of allelic proportions by pyrosequencing revealed that both alleles of PPARG and SREBF1 are expressed at similar levels. PPARGC1A showed the greatest ASE imbalance in fat deposits in Polish Large White (PLW), Polish Landrace and Pietrain pigs; and PPARA in PLW pigs. Significant deviations of mean PPARGC1A allelic transcript ratio between cDNA and genomic DNA were detected in all tissues, with the most pronounced difference (p < 0.001) in visceral fat of PLW pigs. To search for potential cis-regulatory elements affecting ASE in the PPARGC1A gene we analyzed the effects of four SNPs (rs337351686, rs340650517, rs336405906 and rs345224049) in the promoter region, but none were associated with ASE in the breeds studied. DNA methylation analysis revealed significant CpG methylation differences between samples showing balanced (allelic transcript ratio ≈1) and imbalanced allelic expression for CpG site at the genomic position in chromosome 8 (SSC8): 18527678 in visceral fat (p = 0.017) and two CpG sites (SSC8:18525215, p = 0.030; SSC8:18525237, p = 0.031) in subcutaneous fat.

Conclusions

Our analysis of differential allelic expression suggests that PPARGC1A is subjected to cis-regulation in porcine fat tissues. Further studies are necessary to identify other regulatory elements localized outside the PPARGC1A proximal promoter region.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0696-6) contains supplementary material, which is available to authorized users.

Zygotic gene activation in the chicken occurs in two waves, the first involving only maternally derived genes

The first wave of transcriptional activation occurs after fertilisation in a species-specific pattern. Despite its importance to initial embryonic development, the characteristics of transcription following fertilisation are poorly understood in Aves. Here, we report detailed insights into the onset of genome activation in chickens. We established that two waves of transcriptional activation occurred, one shortly after fertilisation and another at Eyal-Giladi and Kochav Stage V. We found 1544 single nucleotide polymorphisms across 424 transcripts derived from parents that were expressed in offspring during the early embryonic stages. Surprisingly, only the maternal genome was activated in the zygote, and the paternal genome remained silent until the second-wave, regardless of the presence of a paternal pronucleus or supernumerary sperm in the egg. The identified maternal genes involved in cleavage that were replaced by bi-allelic expression. The results demonstrate that only maternal alleles are activated in the chicken zygote upon fertilisation, which could be essential for early embryogenesis and evolutionary outcomes in birds.

The early stages of animal development involve a handover of genetic control. Initially, the egg cell is maintained by genetic information inherited from the mother, but soon after fertilization it starts to depend on its own genes instead. Activating genes inside the fertilized egg cell (zygote) so that they can take control of development is known as zygotic genome activation.

Despite the fact that birds are often used to study how embryos develop, zygotic genome activation in birds is not well understood. Fertilization in birds, including chickens, is different to mammals in that it requires multiple sperm to fertilize an egg cell. As such, zygotic genome activation in birds is likely to differ from that in mammals.

By examining gene expression in embryos from mixed-breed chickens, Hwang, Seo et al. showed that there are two stages of zygotic genome activation in chickens. The genes derived from the mother become active in the first stage, while genes from the father become active in the second stage. Genome activation in birds is therefore very different to the same process in mammals, which involves genome activation of both parents from the first stage. This extra level of control may help to prevent genetic complications resulting from the presence of multiple sperm, each of which carries a different set of genes from the father.

Canine MAS1: monoallelic expression is suggestive of an imprinted gene

Imprinted genes are epigenetically modified in a parent-of-origin dependent manner and as a consequence are differentially expressed, with one allele typically expressed while the other is repressed. In canine, the insulin like growth factor 2 receptor gene (IGF2R) is imprinted with predominant expression of the maternally inherited allele. Because imprinted genes usually occur in clusters, we examined the allelic expression pattern of the gene encoding the canine Mas receptor (MAS1), which is located upstream of IGF2R on canine chromosome 1 and is highly conserved in mammals. In this report we describe monoallelic expression of canine MAS1 in the neonatal umbilical cord of several individuals and we identify the expressed allele as maternally inherited. These data suggest that canine MAS1 is an imprinted gene.

DNA methylation and expression of proopiomelanocortin (POMC) gene in the hypothalamus of three-week-old chickens show sex-specific differences

Increased availability and improved sequence annotation of the chicken (Gallus gallus f. domestica) genome have sparked interest in the bird as a model system to investigate translational embryonic development and health/disease outcomes. However, the epigenetics of this bird genome remain unclear. The aim of this study was to determine the levels of gene expression and DNA methylation at the proopiomelanocortin (POMC) gene in the hypothalamus of 3-week-old chickens. POMC is a key player in the control of the stress response, food intake, and metabolism. DNA methylation of the promoter, CpG island, and gene body regions of POMC were measured. Our data illustrate the pattern, variability, and functionality of DNA methylation for POMC expression in the chicken. Our findings show correlation of methylation pattern and gene expression along with sex-specific differences in POMC. Overall, these novel data highlight the promising potential of the chicken as a model and also the need for breeders and researchers to consider sex ratios in their studies.

Transgenerational epigenetic inheritance in birds

While it has been shown that epigenetics accounts for a portion of the variability of complex traits linked to interactions with the environment, the real contribution of epigenetics to phenotypic variation remains to be assessed. In recent years, a growing number of studies have revealed that epigenetic modifications can be transmitted across generations in several animal species. Numerous studies have demonstrated inter- or multi-generational effects of changing environment in birds, but very few studies have been published showing epigenetic transgenerational inheritance in these species. In this review, we mention work conducted in parent-to-offspring transmission analyses in bird species, with a focus on the impact of early stressors on behaviour. We then present recent advances in transgenerational epigenetics in birds, which involve germline linked non-Mendelian inheritance, underline the advantages and drawbacks of working on birds in this field and comment on future directions of transgenerational studies in bird species.

Vertical transmission of positive-sense single-stranded RNA viruses in plants as a model for arboviral induced teratogenesis

Teratogenic viruses have increased public health importance with the emergence of Zika virus and a recent decline in rubella virus vaccination. Of the seven viruses known to cause birth defects in humans, three are mosquito-borne pathogens. Ethical oversight, compliance, rising costs, and the need for specialized training slow the pace of study of these human pathogens compared to study of similar teratogenic viruses in plants. Plant viruses have served as models for human viruses which can be applied to animal systems. This review describes the similar features of plant and animal teratogenic arboviruses and the common systems and barriers that are encountered during vertical transmission in the host.

Allele-Specific Expression Analysis Does Not Support Sex Chromosome Inactivation on the Chicken Z Chromosome

Heterogametic sex chromosomes have evolved many times independently, and in many cases, the loss of functional genes from the sex-limited Y or W chromosome leaves only one functional gene copy on the corresponding X or Z chromosome in the heterogametic sex. Because gene dose often correlates with gene expression level, this difference in gene dose between males and females for X- or Z-linked genes in some cases has selected for chromosome-wide transcriptional dosage compensation mechanisms to counteract any reduction in expression in the heterogametic sex. These mechanisms are thought to restore the balance between sex-linked loci and the autosomal genes they interact with, and this also typically results in equal expression between the sexes. However, dosage compensation in many other species is incomplete, and in the case of birds average expression from males (ZZ) remains higher than in females (ZW). Interestingly, recent reports in chickens and related species have shown that the Z chromosome is expressed less in males than would be expected from two copies of the chromosome, and recent data from cell-based approaches on 11 loci in chicken have suggested that one Z chromosome is partially inactivated in males, in a mechanism thought to be homologous to X inactivation in therian mammals. In the present study, we use controlled crosses in three tissues to test for the presence of Z inactivation in males, which would be expected to bias transcription to the active gene copy (allele-specific expression). We show that for the vast majority of genes on the chicken Z chromosome, males express both parental alleles at statistically similar levels, indicating no Z chromosome inactivation. For those Z chromosome loci with detectable ASE in males, we show that the most likely cause is cis-regulatory variation, rather than Z chromosome inactivation. Taken together, our results indicate that unlike the X chromosome in mammals, Z inactivation does not affect an appreciable number of loci in chicken.

RNA-Seq Analyses Identify Frequent Allele Specific Expression and No Evidence of Genomic Imprinting in Specific Embryonic Tissues of Chicken

Epigenetic and genetic cis-regulatory elements in diploid organisms may cause allele specific expression (ASE) – unequal expression of the two chromosomal gene copies. Genomic imprinting is an intriguing type of ASE in which some genes are expressed monoallelically from either the paternal allele or maternal allele as a result of epigenetic modifications. Imprinted genes have been identified in several animal species and are frequently associated with embryonic development and growth. Whether genomic imprinting exists in chickens remains debatable, as previous studies have reported conflicting evidence. Albeit no genomic imprinting has been reported in the chicken embryo as a whole, we interrogated the existence or absence of genomic imprinting in the 12-day-old chicken embryonic brain and liver by examining ASE in F1 reciprocal crosses of two highly inbred chicken lines (Fayoumi and Leghorn). We identified 5197 and 4638 ASE SNPs, corresponding to 18.3% and 17.3% of the genes with a detectable expression in the embryonic brain and liver, respectively. There was no evidence detected of genomic imprinting in 12-day-old embryonic brain and liver. While ruling out the possibility of imprinted Z-chromosome inactivation, our results indicated that Z-linked gene expression is partially compensated between sexes in chickens.

Accurate characterization of the IFITM locus using MiSeq and PacBio sequencing shows genetic variation in Galliformes

Background

Interferon inducible transmembrane (IFITM) proteins are effectors of the immune system widely characterized for their role in restricting infection by diverse enveloped and non-enveloped viruses. The chicken IFITM (chIFITM) genes are clustered on chromosome 5 and to date four genes have been annotated, namely chIFITM1, chIFITM3, chIFITM5 and chIFITM10. However, due to poor assembly of this locus in the Gallus Gallus v4 genome, accurate characterization has so far proven problematic. Recently, a new chicken reference genome assembly Gallus Gallus v5 was generated using Sanger, 454, Illumina and PacBio sequencing technologies identifying considerable differences in the chIFITM locus over the previous genome releases.

Methods

We re-sequenced the locus using both Illumina MiSeq and PacBio RS II sequencing technologies and we mapped RNA-seq data from the European Nucleotide Archive (ENA) to this finalized chIFITM locus. Using SureSelect probes capture probes designed to the finalized chIFITM locus, we sequenced the locus of a different chicken breed, namely a White Leghorn, and a turkey.

Results

We confirmed the Gallus Gallus v5 consensus except for two insertions of 5 and 1 base pair within the chIFITM3 and B4GALNT4 genes, respectively, and a single base pair deletion within the B4GALNT4 gene. The pull down revealed a single amino acid substitution of A63V in the CIL domain of IFITM2 compared to Red Jungle fowl and 13, 13 and 11 differences between IFITM1, 2 and 3 of chickens and turkeys, respectively. RNA-seq shows chIFITM2 and chIFITM3 expression in numerous tissue types of different chicken breeds and avian cell lines, while the expression of the putative chIFITM1 is limited to the testis, caecum and ileum tissues.

Conclusions

Locus resequencing using these capture probes and RNA-seq based expression analysis will allow the further characterization of genetic diversity within Galliformes.

Electronic supplementary material

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

Detection of QTL for traits related to adaptation to sub-optimal climatic conditions in chickens

BACKGROUND

Growth traits can be used as indicators of adaptation to sub-optimal conditions. The current study aimed at identifying quantitative trait loci (QTL) that control performance under variable temperature conditions in chickens.

METHODS

An F2 population was produced by crossing the Taiwan Country chicken L2 line (selected for body weight, comb area, and egg production) with an experimental line of Rhode Island Red layer R- (selected for low residual feed consumption). A total of 844 animals were genotyped with the 60 K Illumina single nucleotide polymorphism (SNP) chip. Whole-genome interval linkage mapping and a genome-wide association study (GWAS) were performed for body weight at 0, 4, 8, 12, and 16 weeks of age, shank length at 8 weeks of age, size of comb area at 16 weeks of age, and antibody response to sheep red blood cells at 11 weeks of age (7 and 14 days after primary immunization). Relevant genes were identified based on functional annotation of candidate genes and potentially relevant SNPs were detected by comparing whole-genome sequences of several birds between the parental lines.

RESULTS

Whole-genome QTL analysis revealed 47 QTL and 714 effects associated with 178 SNPs were identified by GWAS with 5% Bonferroni genome-wide significance. Little overlap was observed between the QTL and GWAS results, with only two chromosomal regions detected by both approaches, i.e. one on GGA24 (GGA for Gallus gallus chromosome) for BW04 and one on GGAZ for six growth-related traits. Based on whole-genome sequence, differences between the parental lines based on several birds were screened in the genome-wide QTL regions and in a region detected by both methods, resulting in the identification of 106 putative candidate genes with a total of 15,443 SNPs, of which 41 were missense and 1698 were not described in the dbSNP archive.

CONCLUSIONS

The QTL detected in this study for growth and morphological traits likely influence adaptation of chickens to sub-tropical climate. Using whole-genome sequence data, we identified candidate SNPs for further confirmation of QTL in the F2 design. A strong QTL effect found on GGAZ underlines the importance of sex-linked inheritance for growth traits in chickens.

Snail2 and Zeb2 repress P-cadherin to define embryonic territories in the chick embryo

Snail and Zeb transcription factors induce epithelial-to-mesenchymal transition (EMT) in embryonic and adult tissues by direct repression of E-cadherin transcription. The repression of E-cadherin transcription by the EMT inducers Snail1 and Zeb2 plays a fundamental role in defining embryonic territories in the mouse, as E-cadherin needs to be downregulated in the primitive streak and in the epiblast, concomitant with the formation of mesendodermal precursors and the neural plate, respectively. Here, we show that in the chick embryo, E-cadherin is weakly expressed in the epiblast at pre-primitive streak stages where it is substituted for by P-cadherin We also show that Snail2 and Zeb2 repress P-cadherin transcription in the primitive streak and the neural plate, respectively. This indicates that E- and P-cadherin expression patterns evolved differently between chick and mouse. As such, the Snail1/E-cadherin axis described in the early mouse embryo corresponds to Snail2/P-cadherin in the chick, but both Snail factors and Zeb2 fulfil a similar role in chick and mouse in directly repressing ectodermal cadherin genes to contribute to the delamination of mesendodermal precursors at gastrulation and the proper specification of the neural ectoderm during neural induction.

Disentangling Genetic and Prenatal Maternal Effects on Offspring Size and Survival

Organizational processes during prenatal development can have long-term effects on an individual's phenotype. Because these early developmental stages are sensitive to environmental influences, mothers are in a unique position to alter their offspring's phenotype by differentially allocating resources to their developing young. However, such prenatal maternal effects are difficult to disentangle from other forms of parental care, additive genetic effects, and/or other forms of maternal inheritance, hampering our understanding of their evolutionary consequences. Here we used divergent selection lines for high and low prenatal maternal investment and their reciprocal line crosses in a precocial bird-the Japanese quail (Coturnix japonica)-to quantify the relative importance of genes and prenatal maternal effects in shaping offspring phenotype. Maternal but not paternal origin strongly affected offspring body size and survival throughout development. Although the effects of maternal egg investment faded over time, they were large at key life stages. Additionally, there was evidence for other forms of maternal inheritance affecting offspring phenotype at later stages of development. Our study is among the first to successfully disentangle prenatal maternal effects from all other sources of confounding variation and highlights the important role of prenatal maternal provisioning in shaping offspring traits closely linked to fitness.

Matrilineal inheritance of a key mediator of prenatal maternal effects

Sex-linkage is predicted to evolve in response to sex-specific or sexually antagonistic selection. In line with this prediction, most sex-linked genes are associated with reproduction in the respective sex. In addition to traits directly involved in fertility and fecundity, mediators of maternal effects may be predisposed to evolve sex-linkage, because they indirectly affect female fitness through their effect on offspring phenotype. Here, we test for sex-linked inheritance of a key mediator of prenatal maternal effects in oviparous species, the transfer of maternally derived testosterone to the eggs. Consistent with maternal inheritance, we found that in Japanese quail (Coturnix japonica) granddaughters resemble their maternal (but not their paternal) grandmother in yolk testosterone deposition. This pattern of resemblance was not due to non-genetic priming effects of testosterone exposure during prenatal development, as an experimental manipulation of yolk testosterone levels did not affect the females' testosterone transfer to their own eggs later in life. Instead, W chromosome and/or mitochondrial variation may underlie the observed matrilineal inheritance pattern. Ultimately, the inheritance of mediators of maternal effects along the maternal line will allow for a fast and direct response to female-specific selection, thereby affecting the dynamics of evolutionary processes mediated by maternal effects.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

Embryonic development and inviability phenotype of chicken-Japanese quail F1 hybrids

Interspecific hybrid incompatibility, including inviability and sterility, is important in speciation; however, its genetic basis remains largely unknown in vertebrates. Crosses between male chickens and female Japanese quails using artificial insemination can generate intergeneric hybrids; however, the hatching rate is low, and hatched hybrids are only sterile males. Hybrid development is arrested frequently during the early embryonic stages, and the sex ratio of living embryos is male-biased. However, the development and sex ratio of hybrid embryos have not been comprehensively analyzed. In the present study, we observed delayed embryonic development of chicken-quail hybrids during the early stage, compared with that of chickens and quails. The survival rate of hybrids decreased markedly during the blastoderm-to-pre-circulation stage and then decreased gradually through the subsequent stages. Hybrid females were observed at more than 10 d of incubation; however, the sex ratio of hybrids became male-biased from 10 d of incubation. Severely malformed embryos were observed frequently in hybrids. These results suggest that developmental arrest occurs at various stages in hybrid embryos, including a sexually non-biased arrest during the early stage and a female-biased arrest during the late stage. We discuss the genetic basis for hybrid inviability and its sex bias.

The Extent of mRNA Editing Is Limited in Chicken Liver and Adipose, but Impacted by Tissular Context, Genotype, Age, and Feeding as Exemplified with a Conserved Edited Site in COG3

RNA editing is a posttranscriptional process leading to differences between genomic DNA and transcript sequences, potentially enhancing transcriptome diversity. With recent advances in high-throughput sequencing, many efforts have been made to describe mRNA editing at the transcriptome scale, especially in mammals, yielding contradictory conclusions regarding the extent of this phenomenon. We show, by detailed description of the 25 studies focusing so far on mRNA editing at the whole-transcriptome scale, that systematic sequencing artifacts are considered in most studies whereas biological replication is often neglected and multi-alignment not properly evaluated, which ultimately impairs the legitimacy of results. We recently developed a rigorous strategy to identify mRNA editing using mRNA and genomic DNA sequencing, taking into account sequencing and mapping artifacts, and biological replicates. We applied this method to screen for mRNA editing in liver and white adipose tissue from eight chickens and confirm the small extent of mRNA recoding in this species. Among the 25 unique edited sites identified, three events were previously described in mammals, attesting that this phenomenon is conserved throughout evolution. Deeper investigations on five sites revealed the impact of tissular context, genotype, age, feeding conditions, and sex on mRNA editing levels. More specifically, this analysis highlighted that the editing level at the site located on COG3 was strongly regulated by four of these factors. By comprehensively characterizing the mRNA editing landscape in chickens, our results highlight how this phenomenon is limited and suggest regulation of editing levels by various genetic and environmental factors.

Next-Generation Sequencing Techniques Reveal that Genomic Imprinting Is Absent in Day-Old Gallus gallus domesticus Brains

Genomic imprinting is a phenomenon characterized by parent-of-origin-specific gene expression. While widely documented in viviparous mammals and plants, imprinting in oviparous birds remains controversial. Because genomic imprinting is temporal- and tissue-specific, we investigated this phenomenon only in the brain tissues of 1-day-old chickens (Gallus gallus). We used next-generation sequencing technology to compare four transcriptomes pooled from 11 chickens, generated from reciprocally crossed families, to the DNA sequences of their parents. Candidate imprinted genes were then selected from these sequence alignments and subjected to verification experiments that excluded all but one SNP. Subsequent experiments performed with two new sets of reciprocally crossed families resulted in the exclusion of that candidate SNP as well. Attempts to find evidence of genomic imprinting from long non-coding RNAs yielded negative results. We therefore conclude that genomic imprinting is absent in the brains of 1-day-old chickens. However, due to the temporal and tissue specificity of imprinting, our results cannot be extended to all growth stages and tissue types.

Genome-Wide Characterization of RNA Editing in Chicken Embryos Reveals Common Features among Vertebrates

RNA editing results in a post-transcriptional nucleotide change in the RNA sequence that creates an alternative nucleotide not present in the DNA sequence. This leads to a diversification of transcription products with potential functional consequences. Two nucleotide substitutions are mainly described in animals, from adenosine to inosine (A-to-I) and from cytidine to uridine (C-to-U). This phenomenon is described in more details in mammals, notably since the availability of next generation sequencing technologies allowing whole genome screening of RNA-DNA differences. The number of studies recording RNA editing in other vertebrates like chicken is still limited. We chose to use high throughput sequencing technologies to search for RNA editing in chicken, and to extend the knowledge of its conservation among vertebrates. We performed sequencing of RNA and DNA from 8 embryos. Being aware of common pitfalls inherent to sequence analyses that lead to false positive discovery, we stringently filtered our datasets and found fewer than 40 reliable candidates. Conservation of particular sites of RNA editing was attested by the presence of 3 edited sites previously detected in mammals. We then characterized editing levels for selected candidates in several tissues and at different time points, from 4.5 days of embryonic development to adults, and observed a clear tissue-specificity and a gradual increase of editing level with time. By characterizing the RNA editing landscape in chicken, our results highlight the extent of evolutionary conservation of this phenomenon within vertebrates, attest to its tissue and stage specificity and provide support of the absence of non A-to-I events from the chicken transcriptome.

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

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

Acquired alterations of hypothalamic gene expression of insulin and leptin receptors and glucose transporters in prenatally high-glucose exposed three-week old chickens do not coincide with aberrant promoter DNA methylation

Background

Prenatal exposures may have a distinct impact for long-term health, one example being exposure to maternal ‘diabesity’ during pregnancy increasing offspring ‘diabesity’ risk. Malprogramming of the central nervous regulation of body weight, food intake and metabolism has been identified as a critical mechanism. While concrete disrupting factors still remain unclear, growing focus on acquired epigenomic alterations have been proposed. Due to the independent development from the mother, the chicken embryo provides a valuable model to distinctively establish causal factors and mechanisms.

Aim

The aim of this study was to determine the effects of prenatal hyperglycemia on postnatal hypothalamic gene expression and promoter DNA methylation in the chicken.

Methods and Findings

To temporarily induce high-glucose exposure in chicken embryos, 0.5 ml glucose solution (30 mmol/l) were administered daily via catheter into a vessel of the chorioallantoic egg membrane from days 14 to 17 of incubation. At three weeks of postnatal age, body weight, total body fat, blood glucose, mRNA expression (INSR, LEPR, GLUT1, GLUT3) as well as corresponding promoter DNA methylation were determined in mediobasal hypothalamic brain slices (Nucleus infundibuli hypothalami). Although no significant changes in morphometric and metabolic parameters were detected, strongly decreased mRNA expression occurred in all candidate genes. Surprisingly, however, no relevant alterations were observed in respective promoter methylation.

Conclusion

Prenatal hyperglycemia induces strong changes in later hypothalamic expression of INSR, LEPR, GLUT1, and GLUT3 mRNA. While the chicken provides an interesting approach for developmental malprogramming, the classical expression regulation via promoter methylation was not observed here. This may be due to alternative/interacting brain mechanisms or the thus far under-explored bird epigenome.

Genome-wide identification of allele-specific expression in response to Streptococcus suis 2 infection in two differentially susceptible pig breeds

Although allele expression imbalance has been recognized in many species, and strongly linked to diseases, no whole transcriptome allele imbalance has been detected in pigs during pathogen infections. The pathogen Streptococcus suis 2 (SS2) causes serious zoonotic disease. Different pig breeds show differential susceptibility/resistance to pathogen infection, but the biological insight is little known. Here we analyzed allele-specific expression (ASE) using the spleen transcriptome of four pigs belonging to two phenotypically different breeds after SS2 infection. The comparative analysis of allele specific SNPs between control and infected animals revealed 882 and 1096 statistically significant differentially expressed allele SNPs (criteria: ratio ≧ 2 or ≦ 0.5) in Landrace and Enshi black pig, respectively. Twenty nine allelically imbalanced SNPs were further verified by Sanger sequencing, and later six SNPs were quantified by pyrosequencing assay. The pyrosequencing results are in agreement with the RNA-seq results, except two SNPs. Looking at the role of ASE in predisposition to diseases, the discovery of causative variants by ASE analysis might help the pig industry in long term to design breeding programs for improving SS2 resistance.