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

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.

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.

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

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.

Favoring Expression of Yak Alleles in Interspecies F1 Hybrids of Cattle and Yak Under High-Altitude Environments

Both cis- and trans-regulation could cause differential expression between the parental alleles in diploid species that might have broad biological implications. Due to the relatively distant genetic divergence between cattle and yak, as well as their differential adaptation to high-altitude environments, we investigated genome-wide allelic differential expression (ADE) in their F1 hybrids using Nanopore long-read RNA-seq technology. From adult F1 hybrids raised in high-altitude, ten lung and liver tissues were individually sequenced for producing 31.6 M full-length transcript sequences. Mapping against autosomal homologous regions between cattle and yak, we detected 17,744 and 14,542 protein-encoding genes expressed in lung and liver tissues, respectively. According to the parental assignments of transcript sequences, a total of 3,381 genes were detected to show ADE in at least one sample. There were 186 genes showing ubiquitous ADE in all the studied animals, and among them 135 and 37 genes had consistent higher expression of yak and cattle alleles, respectively. Functional analyses revealed that the genes with favoring expression of yak alleles have been involved in the biological progresses related with hypoxia adaptation and immune response. In contrast, the genes with favoring expression of cattle alleles have been enriched into different biological progresses, such as secretion of endocrine hormones and lipid metabolism. Our results would support unequal contribution of parental genes to environmental adaptation in the F1 hybrids of cattle and yak.

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.

Sex identification in embryos and adults of Darwin's finches

Darwin’s finches are an iconic example of adaptive radiation and evolution under natural selection. Comparative genetic studies using embryos of Darwin’s finches have shed light on the possible evolutionary processes underlying the speciation of this clade. Molecular identification of the sex of embryonic samples is important for such studies, where this information often cannot be inferred otherwise. We tested a fast and simple chicken embryo protocol to extract DNA from Darwin’s finch embryos. In addition, we applied minor modifications to two of the previously reported PCR primer sets for CHD1, a gene used for sexing adult passerine birds. The sex of all 29 tested embryos of six species of Darwin’s finches was determined successfully by PCR, using both primer sets. Next to embryos, hatchlings and fledglings are also impossible to distinguish visually. This extends to juveniles of sexually dimorphic species which are yet to moult in adult-like plumage and beak colouration. Furthermore, four species of Darwin’s finches are monomorphic, males and females looking alike. Therefore, sex assessment in the field can be a source of error, especially with respect to juveniles and mature monomorphic birds outside of the mating season. We caught 567 juveniles and adults belonging to six species of Darwin’s finches and only 44% had unambiguous sex-specific morphology. We sexed 363 birds by PCR: individuals sexed based on marginal sex specific morphological traits; and birds which were impossible to classify in the field. PCR revealed that for birds with marginal sex specific traits, sexing in the field produced a 13% error rate. This demonstrates that PCR based sexing can improve field studies on Darwin’s finches, especially when individuals with unclear sex-related morphology are involved. The protocols used here provide an easy and reliable way to sex Darwin’s finches throughout ontogeny, from embryos to adults.

Investigation of allele specific expression in various tissues of broiler chickens using the detection tool VADT

Differential abundance of allelic transcripts in a diploid organism, commonly referred to as allele specific expression (ASE), is a biologically significant phenomenon and can be examined using single nucleotide polymorphisms (SNPs) from RNA-seq. Quantifying ASE aids in our ability to identify and understand cis-regulatory mechanisms that influence gene expression, and thereby assist in identifying causal mutations. This study examines ASE in breast muscle, abdominal fat, and liver of commercial broiler chickens using variants called from a large sub-set of the samples (n = 68). ASE analysis was performed using a custom software called VCF ASE Detection Tool (VADT), which detects ASE of biallelic SNPs using a binomial test. On average ~ 174,000 SNPs in each tissue passed our filtering criteria and were considered informative, of which ~ 24,000 (~ 14%) showed ASE. Of all ASE SNPs, only 3.7% exhibited ASE in all three tissues, with ~ 83% showing ASE specific to a single tissue. When ASE genes (genes containing ASE SNPs) were compared between tissues, the overlap among all three tissues increased to 20.1%. Our results indicate that ASE genes show tissue-specific enrichment patterns, but all three tissues showed enrichment for pathways involved in translation.

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.

Allele-specific expression and gene regulation help explain transgressive thermal tolerance in non-native hybrids of the endangered California tiger salamander (Ambystoma californiense)

Hybridization between native and non-native species is an ongoing global conservation threat. Hybrids that exhibit traits and tolerances that surpass parental values are of particular concern, given their potential to outperform native species. Effective management of hybrid populations requires an understanding of both physiological performance and the underlying mechanisms that drive transgressive hybrid traits. Here, we explore several aspects of the hybridization between the endangered California tiger salamander (Ambystoma californiense; CTS) and the introduced barred tiger salamander (Ambystoma mavortium; BTS). We assayed critical thermal maximum (CTMax) to compare the ability of CTS, BTS and F1 hybrids to tolerate acute thermal stress, and found that hybrids exhibit a wide range of CTMax values, with 33% (4/12) able to tolerate temperatures greater than either parent. We then quantified the genomic response, measured at the RNA transcript level, of each salamander, to explore the mechanisms underlying thermal tolerance strategies. We found that CTS and BTS have strikingly different values and tissue-specific patterns of overall gene expression, with hybrids expressing intermediate values. F1 hybrids display abundant and variable degrees of allele-specific expression (ASE), likely arising from extensive compensatory evolution in gene regulatory mechanisms between CTS and BTS. We found evidence that the proportion of genes with allelic imbalance in individual hybrids correlates with their CTMax, suggesting a link between ASE and expanded thermal tolerance that may contribute to the success of hybrid salamanders in California. Future climate change may further complicate management of CTS if hybrid salamanders are better equipped to deal with rising temperatures.

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.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

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.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.

Identification and Comparison of Imprinted Genes Across Plant Species

Genomic imprinting is a phenomenon that occurs in flowering plants and mammals, whereby a gene is expressed in a parent-of-origin-specific manner. Although imprinting has now been examined genome-wide in a number of species using RNA-seq, the analyses used to assess imprinting vary between studies, making consistent comparisons between species difficult. Here we present a simple, easy-to-use bioinformatic pipeline for imprinting analyses suitable for any tissue, including plant endosperm. All relevant scripts can be downloaded. As an illustrative example, we reanalyze published data from A. thaliana and Z. mays endosperm using the pipeline and then demonstrate how to use the results to assess the conservation of imprinting between these species. We also introduce the Plant Imprinting Database, a repository for published imprinting datasets in plants that can be used to view, compare, and download data.

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.

Embryo microinjection of the lecithotrophic sea urchin Heliocidaris erythrogramma

Microinjection is a common embryological technique used for many types of experiments, including lineage tracing, manipulating gene expression, or genome editing. Injectable reagents include mRNA overexpression, mis-expression, or dominant-negative experiments to examine a gene of interest, a morpholino antisense oligo to prevent translation of an mRNA or spliceoform of interest and CRISPR-Cas9 reagents. Thus, the technique is broadly useful for basic embryological studies, constructing gene regulatory networks, and directly testing hypotheses about cis-regulatory and coding sequence changes underlying the evolution of development. However, the methods for microinjection in typical planktotrophic marine invertebrates may not work well in the highly modified eggs and embryos of lecithotrophic species. This protocol is optimized for the lecithotrophic sea urchin Heliocidaris erythrogramma.

Analysis of the Progeny of Sibling Matings Reveals Regulatory Variation Impacting the Transcriptome of Immune Cells in Commercial Chickens

There is increasing recognition that the underlying genetic variation contributing to complex traits influences transcriptional regulation and can be detected at a population level as expression quantitative trait loci. At the level of an individual, allelic variation in transcriptional regulation of individual genes can be detected by measuring allele-specific expression in RNAseq data. We reasoned that extreme variants in gene expression could be identified by analysis of inbred progeny with shared grandparents. Commercial chickens have been intensively selected for production traits. Selection is associated with large blocks of linkage disequilibrium with considerable potential for co-selection of closely linked "hitch-hiker alleles" affecting traits unrelated to the feature being selected, such as immune function, with potential impact on the productivity and welfare of the animals. To test this hypothesis that there is extreme allelic variation in immune-associated genes we sequenced a founder population of commercial broiler and layer birds. These birds clearly segregated genetically based upon breed type. Each genome contained numerous candidate null mutations, protein-coding variants predicted to be deleterious and extensive non-coding polymorphism. We mated selected broiler-layer pairs then generated cohorts of F2 birds by sibling mating of the F1 generation. Despite the predicted prevalence of deleterious coding variation in the genomic sequence of the founders, clear detrimental impacts of inbreeding on survival and post-hatch development were detected in only one F2 sibship of 15. There was no effect on circulating leukocyte populations in hatchlings. In selected F2 sibships we performed RNAseq analysis of the spleen and isolated bone marrow-derived macrophages (with and without lipopolysaccharide stimulation). The results confirm the predicted emergence of very large differences in expression of individual genes and sets of genes. Network analysis of the results identified clusters of co-expressed genes that vary between individuals and suggested the existence of trans-acting variation in the expression in macrophages of the interferon response factor family that distinguishes the parental broiler and layer birds and influences the global response to lipopolysaccharide. This study shows that the impact of inbreeding on immune cell gene expression can be substantial at the transcriptional level, and potentially opens a route to accelerate selection using specific alleles known to be associated with desirable expression levels.

Identification of Single Nucleotide Polymorphism in Red Clover (Trifolium pratense L.) Using Targeted Genomic Amplicon Sequencing and RNA-seq

Red clover (Trifolium pratense L.) is a diploid, naturally cross-pollinated, cool-season species. As a perennial forage legume, red clover is mostly cultivated in temperate regions worldwide. Being a non-model crop species, genomic resources for red clover have been underdeveloped. Thus far, genomic analysis used in red clover has mainly relied on simple sequence repeat (SSR) markers. However, SSR markers are sparse in the genome and it is often difficult to unambiguously map them using short reads generated by next generation sequencing technology. Single nucleotide polymorphisms (SNPs) have been successfully applied in genomics assisted breeding in several agriculturally important species. Due to increasing importance of legumes in forage production, there is a clear need to develop SNP based markers for red clover that can be applied in breeding applications. In this study, we first developed an analytical pipeline that can confidently identify SNPs in a set of 72 different red clover genotypes using sequences generated by targeted amplicon sequencing. Then, with the same filtering stringency used in this pipeline, we used sequences from publicly available RNA-seq data to identify confident SNPs in different red clover varieties. Using this strategy, we have identified a total of 69,975 SNPs across red clover varieties. Among these, 28% (19,116) of them are missense mutations. Using Medicago truncatula as the reference, we annotated the regions affected by these missense mutations. We identified 2,909 protein coding regions with missense mutations. Pathway analysis of these coding regions indicated several biological processes impacted by these mutations. Specifically, three domains (homeobox domain, pentatricopeptide repeat containing plant-like, and regulator of Vps4 activity) were identified with five or more missense SNPs. These domain might also be a functional contributor in the molecular mechanisms of self-incompatibility in red clover. Future in-depth sequence diversity analysis of these three genes may yield valuable insights into the molecular mechanism involved in self-incompatibility in red clover.

Variant analysis pipeline for accurate detection of genomic variants from transcriptome sequencing data

The wealth of information deliverable from transcriptome sequencing (RNA-seq) is significant, however current applications for variant detection still remain a challenge due to the complexity of the transcriptome. Given the ability of RNA-seq to reveal active regions of the genome, detection of RNA-seq SNPs can prove valuable in understanding the phenotypic diversity between populations. Thus, we present a novel computational workflow named VAP (Variant Analysis Pipeline) that takes advantage of multiple RNA-seq splice aware aligners to call SNPs in non-human models using RNA-seq data only. We applied VAP to RNA-seq from a highly inbred chicken line and achieved high accuracy when compared with the matching whole genome sequencing (WGS) data. Over 65% of WGS coding variants were identified from RNA-seq. Further, our results discovered SNPs resulting from post transcriptional modifications, such as RNA editing, which may reveal potentially functional variation that would have otherwise been missed in genomic data. Even with the limitation in detecting variants in expressed regions only, our method proves to be a reliable alternative for SNP identification using RNA-seq data. The source code and user manuals are available at https://modupeore.github.io/VAP/.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

RNA-Seq Analysis Reveals Spatial and Sex Differences in Pectoralis Major Muscle of Broiler Chickens Contributing to Difference in Susceptibility to Wooden Breast Disease

Wooden Breast Disease (WBD) is a novel myopathy affecting the pectoralis major muscle of modern broiler chickens. The etiology of WBD is not currently known, but has been linked to increased feed efficiency, growth rate, and muscle yield in broiler chickens. Differential effect of WBD has been detected between regions of the P. major and between sexes of broilers-male birds and the cranial aspect of the muscle tend to be more severely affected by the disease than females and the caudal aspect. This study aimed to characterize biological differences in the P. major between regions of the muscle and sexes of birds. Samples were taken from the cranial and caudal aspects of P. major muscles of 3-week-old, unaffected male and female birds for RNA sequencing. RNA was extracted and used for preparation of cDNA libraries, which were sequenced by the Delaware Biotechnology Institute (DBI) using HiSeq2500. Sequence reads were aligned to the chicken reference genome with HISAT, and genes were analyzed for differential expression between regions of the breast muscle and sexes of birds using CuffDiff. Functional analysis was performed on differentially expressed genes (DEGs) between sex groups using DAVID and Ingenuity Pathway Analysis (IPA). There were 12 DEGs between cranial and caudal samples, and 260 between male and female birds. Out of the 260 genes differentially expressed between sexes, 189 were upregulated in males. Of this subset, 103 genes (55%) were located on the Z-chromosome. There was increased expression of genes involved in fat metabolism and oxidative stress responses in the cranial region of the P. major muscle, as well as increased expression of fat metabolism, oxidative stress response, antiangiogenesis, and connective tissue proliferation genes in male broilers. These results support the hypothesis that there are biological characteristics in male broilers and the cranial region of the breast muscle that may make them more susceptible to WBD, as well as raising the possibility of a metabolic switch in modern broiler chickens that may be more prominent in males.

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.

Inheritance patterns of the transcriptome in hybrid chickens and their parents revealed by expression analysis

Although many phenotypic traits of chickens have been well documented, the genetic patterns of gene expression levels in chickens remain to be determined. In the present study, we crossed two chicken breeds, White Leghorn (WL) and Cornish (Cor), which have been selected for egg and meat production, respectively, for a few hundred years. We evaluated transcriptome abundance in the brain, muscle, and liver from the day-old progenies of pure-bred WL and Cor, and the hybrids of these two breeds, by RNA-Seq in order to determine the inheritance patterns of gene expression. Comparison among expression levels in the different groups revealed that most of the genes showed conserved expression patterns in all three examined tissues and that brain had the highest number of conserved genes, which indicates that conserved genes are predominantly important compared to others. On the basis of allelic expression analysis, in addition to the conserved genes, we identified the extensive presence of additive, dominant (Cor dominant and WL dominant), over-dominant, and under-dominant genes in all three tissues in hybrids. Our study is the first to provide an overview of inheritance patterns of the transcriptome in layers and broilers, and we also provide insights into the genetics of chickens at the gene expression level.

Survey of allele specific expression in bovine muscle

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3’UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

RNA-Seq Analyses Identify Additivity as the Predominant Gene Expression Pattern in F1 Chicken Embryonic Brain and Liver

The superior performance of hybrids to parents, termed heterosis, has been widely utilized in animal and plant breeding programs, but the molecular mechanism underlying heterosis remains an enigma. RNA-Seq provides a novel way to investigate heterosis at the transcriptome-wide level, because gene expression functions as an intermediate phenotype that contributes to observable traits. Here we compared embryonic gene expression between chicken hybrids and their inbred parental lines to identify inheritance patterns of gene expression. Inbred Fayoumi and Leghorn were crossed reciprocally to obtain F1 fertile eggs. RNA-Seq was carried out using 24 brain and liver samples taken from day 12 embryos, and the differentially expressed (DE) genes were identified by pairwise comparison among the hybrids, parental lines, and mid-parent expression values. Our results indicated the expression levels of the majority of the genes in the F1 cross are not significantly different from the mid-parental values, suggesting additivity as the predominant gene expression pattern in the F1. The second and third prevalent gene expression patterns are dominance and over-dominance. Additionally, we found only 7⁻20% of the DE genes exhibit allele-specific expression in the F1, suggesting that trans regulation is the main driver for differential gene expression and thus contributes to heterosis effect in the F1 crosses.

Homeologue-specific expression divergence in the recently formed tetraploid Capsella bursa-pastoris (Brassicaceae)

Following allopolyploid formation, extensive genome evolution occurs, with the eventual loss of many homeologous gene copies. Although this process of diploidization has occurred many times independently, the evolutionary forces determining the probability and rate of gene loss remain poorly understood. Here, we conduct genome and transcriptome sequencing in a broad sample of Chinese accessions of Capsella bursa-pastoris, a recently formed allotetraploid. Our whole genome data reveal three groups of these accessions: an Eastern group from low-altitude regions, a Western group from high-altitude regions, and a much more differentiated Northwestern group. Population differentiation in total expression was limited among closely related populations; by contrast, the relative expression of the two homeologous copies closely mirrors the genome-wide SNP divergence. Consistent with this, we observe a negative correlation between expression changes in the two homeologues. However, genes showing population genomic evidence for adaptive evolution do not show an enrichment for expression divergence between homeologues, providing no clear evidence for adaptive shifts in relative gene expression. Overall, these patterns suggest that neutral drift may contribute to the population differentiation in the expression of the homeologues, and drive eventual gene loss over longer periods of time.

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.