Integration of GWAS and eGWAS to screen candidate genes underlying green head traits in male ducks

Sexually dimorphic plumage coloration is widespread in birds. The male possesses more brightly colored feathers than the female. Dark green head feathers comprise one of the most typical appearance characteristics of the male Ma duck compared with the female. However, there are noticeable individual differences observed in these characteristics. Herein, genome-wide association studies (GWAS) were employed to investigate the genetic basis of individual differences in male duck green head-related traits. Our results showed that 165 significant SNPs were associated with green head traits. Meanwhile, 71 candidate genes were detected near the significant SNPs, including four genes (CACNA1I, WDR59, GNAO1 and CACNA2D4) related to the individual differences in the green head traits of male ducks. Additionally, the eGWAS identified three SNPs located within two candidate genes (LOC101800026 and SYNPO2) associated with TYRP1 gene expression, and might be important regulators affecting the expression level of TYRP1 in the head skin of male ducks. Our data also suggested that transcription factor MXI1 might regulate the expression of TYRP1, thereby causing differences in the green head traits among male ducks. This study provided primary data for further analysis of the genetic regulation of duck feather color.

Gene expression networks and functionally enriched pathways involved in the response of domestic chicken to acute heat stress

Heat stress in poultry houses, especially in warm areas, is one of the main environmental factors that restrict the growth of broilers or laying performance of layers, suppresses the immune system, and deteriorates egg quality and feed conversion ratio. The molecular mechanisms underlying the response of chicken to acute heat stress (AHS) have not been comprehensively elucidated. Therefore, the main object of the current work was to investigate the liver gene expression profile of chickens under AHS in comparison with their corresponding control groups, using four RNA-seq datasets. The meta-analysis, GO and KEGG pathway enrichment, WGCNA, machine-learning, and eGWAS analyses were performed. The results revealed 77 meta-genes that were mainly related to protein biosynthesis, protein folding, and protein transport between cellular organelles. In other words, under AHS, the expression of genes involving in the structure of rough reticulum membrane and in the process of protein folding was adversely influenced. In addition, genes related to biological processes such as "response to unfolded proteins," "response to reticulum stress" and "ERAD pathway" were differentially regulated. We introduce here a couple of genes such as HSPA5, SSR1, SDF2L1, and SEC23B, as the most significantly differentiated under AHS, which could be used as bio-signatures of AHS. Besides the mentioned genes, the main findings of the current work may shed light to the identification of the effects of AHS on gene expression profiling of domestic chicken as well as the adaptive response of chicken to environmental stresses.

Expression-based Genome-wide Association Study Links OPN and IL1-RA With Newly Diagnosed Type 1 Diabetes in Children

CONTEXT

Islet autoantibodies (IAbs) are currently the most reliable indicators of islet autoimmunity. However, IAbs do not fully meet the need for the prediction and intervention of type 1 diabetes (T1D). Serological proteins should be great sources for biomarkers.

OBJECTIVE

This work aimed to identify new proteomic biomarkers with the technology of an expression-based genome-wide association study (eGWAS) in children newly diagnosed with T1D.

METHODS

In an attempt to identify additional biomarkers, we performed an eGWAS using microarray data from 169 arrays of the pancreatic islets of T1D rodents (78 T1D cases and 91 controls). We ranked all 16 099 protein-coding genes by the likelihood of differential expression in the pancreatic islets. Our top 20 secreted proteins were screened in 170 children including 100 newly diagnosed T1D, and 50 type 2 diabetes (T2D) and 20 age-matched healthy children. With 6 proteins showing significance, we further conducted a validation study using the second independent set of 400 samples from children including 200 newly diagnosed with T1D, 100 T2D, and 100 age-matched controls.

RESULTS

We identified 2 serum proteins that were significantly changed in T1D vs both control and T2D, and 5 serum proteins were significantly changed both in T1D and T2D vs control. Serum osteopontin (OPN) levels were uniquely higher in T1D (T1D vs controls, P = 1.29E-13 ~ 9.38E-7, T1D vs T2D, P = 2.65E-8 ~ 1.58E-7) with no difference between T2D and healthy control individuals. Serum interleukin 1 receptor antagonist (IL-1RA) levels were lower in T1D compared both with T2D (P = 3.36E-9~0.0236) and healthy participants (P = 1.09E-79 ~ 2.00E-12).

CONCLUSION

Our results suggest that OPN and IL1-RA could be candidates for useful biomarkers for T1D in children.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the Allele Frequency Spectrum (AFS) of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a super-abundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of co-expressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BACKGROUND

Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.

RESULTS

The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.

CONCLUSIONS

This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.