Weighted gene co-expression network analysis reveals potential candidate genes affecting drip loss in pork

Drip loss is an essential evaluation indicator for pork quality. It is closely related to other meat quality indicators, including water-holding capacity, water loss rate and pH value at 45 min (pH₁) and 24 h post-mortem (pH₂), and is influenced by environmental and genetic factors and their interactions. We previously conducted differentially expressed gene analysis to identify candidate genes affecting drip loss using eight individuals with extremely high- and low-drip loss selected from 28 purebred Duroc pigs. Using 28 identical samples, in the present study, we performed weighted gene co-expression network analysis with drip loss and drip loss-related traits, including water-holding capacity, water loss rate, pH₁ and pH₂. A total of 25 modules were identified, and five of them correlated with at least two drip loss or drip loss-related traits. After functional enrichment analysis of genes in the five modules, three modules were found to be critical, as their genes were significantly involved in amino acid metabolism, immune response and apoptosis, which have potential relationships with drip loss. Furthermore, we identified five candidate genes affecting drip loss in one critical module, AASS, BCKDHB, ALDH6A1, MUT and MCCC1, as they overlapped with differentially expressed genes detected in our previous study, exhibited protein–protein interactions and had potential biological functions in affecting drip loss according to the literature. The outcomes of the present study enhance our understanding of the molecular mechanisms underlying drip loss and will aid in improving the pork quality.     

RNA‐seq analysis reveals new candidate genes for drip loss in a Pietrain × Duroc × Landrace × Yorkshire population

Drip loss, one of the most important meat quality traits, is characterized by low heritability. To date, the genetic factors affecting the drip loss trait have not been clearly elucidated. The objective of this study was to identify critical candidate genes affecting drip loss. First, we generated a Pietrain × Duroc × Landrace × Yorkshire commercial pig population and obtained phenotypic values for the drip loss trait. Furthermore, we constructed two RNA libraries from pooled samples of longissimus dorsi muscles with the highest (H group) and lowest (L group) drip loss and identified the differentially expressed genes (DEGs) between these extreme phenotypes using RNA‐seq technology. In total, 25 883 genes were detected in the H and L group libraries, and none was specifically expressed in only one library. Comparative analysis of gene expression levels found that 150 genes were differentially expressed, of which 127 were upregulated and 23 were downregulated in the H group relative to the L group. In addition, 68 drip loss quantitative trait loci (QTL) overlapping with 63 DEGs were identified, and these QTL were distributed mainly on chromosomes 1, 2, 5 and 6. Interestingly, the triadin (TRDN) gene, which is involved in muscle contraction and fat deposition, and the myostatin (MSTN) gene, which has a role in muscle growth, were localized to more than two drip loss QTL, suggesting that both are critical candidate genes responsible for drip loss.     

A functional genomics approach to evaluate candidate genes located in a QTL interval for milk production traits on BTA6

The potential genetic and economic advantage of marker-assisted selection for enhanced production in dairy cattle has provided an impetus to conduct numerous genome scans in order to identify associations between DNA markers and future productive potential. One area of focus has been a quantitative trait locus on bovine chromosome 6 (BTA6) found to be associated with milk yield, milk protein and fat percentage, which has been subsequently fine-mapped to six positional candidate genes. Subsequent investigations have yet to resolve which of the potential positional candidate genes is responsible for the observed associations with productive performance. In this study, we analysed candidate gene expression and the effects of gene knockdown on expression of β- and κ-casein mRNA in a small interfering RNA transfected bovine in vitro mammosphere model. From our expression studies in vivo , we observed that four of the six candidates ( ABCG2 , SPP1 , PKD2 and LAP3 ) exhibited differential expression in bovine mammary tissue over the lactation cycle, but in vitro functional studies indicate that inhibition of only one gene, SPP1 , had a significant impact on milk protein gene expression. These data suggest that the gene product of SPP1 (also known as osteopontin) has a significant role in the modulation of milk protein gene expression. While these findings do not exclude other positional candidates from influencing lactation, they support the hypothesis that the gene product of SPP1 is a significant lactational regulatory molecule.     

A statistical multiprobe model for analyzing cis and trans genes in genetical genomics experiments with short-oligonucleotide arrays

Short-oligonucleotide arrays typically contain multiple probes per gene. In genetical genomics applications a statistical model for the individual probe signals can help in separating "true" differential mRNA expression from "ghost" effects caused by polymorphisms, misdesigned probes, and batch effects. It can also help in detecting alternative splicing, start, or termination.     

A novel microarray approach reveals new tissue-specific signatures of known and predicted mammalian microRNAs

Microarrays to examine the global expression levels of microRNAs (miRNAs) in a systematic in-parallel manner have become important tools to help unravel the functions of miRNAs and to understand their roles in RNA-based regulation and their implications in human diseases. We have established a novel miRNA-specific microarray platform that enables the simultaneous expression analysis of both known and predicted miRNAs obtained from human or mouse origin. Chemically modified 2′-O-(2-methoxyethyl)-(MOE) oligoribonucleotide probes were arrayed onto Evanescent Resonance (ER) microchips by robotic spotting. Supplementing the complementary probes against miRNAs with carefully designed mismatch controls allowed for accurate sequence-specific determination of miRNA expression profiles obtained from a panel of mouse tissues. This revealed new expression signatures of known miRNAs as well as of novel miRNAs previously predicted using bioinformatic methods. Systematic confirmation of the array data with northern blotting and, in particular, real-time PCR suggests that the described microarray platform is a powerful tool to analyze miRNA expression patterns with rapid throughput and high fidelity.     

Genome-wide association study revealed a promising region and candidate genes for eggshell quality in an F2 resource population

Background

Eggshell is subject to quality loss with aging process of laying hens, and damaged eggshells result in economic losses of eggs. However, the genetic architecture underlying the dynamic eggshell quality remains elusive. Here, we measured eggshell quality traits, including eggshell weight (ESW), eggshell thickness (EST) and eggshell strength (ESS) at 11 time points from onset of laying to 72 weeks of age and conducted comprehensive genome-wide association studies (GWAS) in 1534 F₂ hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX).

Results

ESWs at all ages exhibited moderate SNP-based heritability estimates (0.30 ~ 0.46), while the estimates for EST (0.21 ~ 0.31) and ESS (0.20 ~ 0.27) were relatively low. Eleven independent univariate genome-wide screens for each trait totally identified 1059, 1026 and 1356 significant associations with ESW, EST and ESS, respectively. Most significant loci were in a region spanning from 57.3 to 71.4 Mb of chromosome 1 (GGA1), which together account for 8.4 ~ 16.5 % of the phenotypic variance for ESW from 32 to 72 weeks of age, 4.1 ~ 6.9 % and 2.95 ~ 16.1 % for EST and ESS from 40 to 72 weeks of age. According to linkage disequilibrium (LD) and conditional analysis, the significant SNPs in this region were in extremely strong linkage disequilibrium status. Ultimately, two missense SNPs in GGA1 and one in GGA4 were considered as promising loci on three independent genes including ITPR2, PIK3C2G, and NCAPG. The homozygotes of advantageously effective alleles on PIK3C2G and ITPR2 possessed the best eggshell quality and could partly counteract the negative effect of aging process. NCAPG had certain effect on eggshell quality for young hens.

Conclusions

Identification of the promising region as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic eggshell quality and has the practical significance in breeding program for the improvement of eggshell quality, especially at the later part of laying cycle.

Electronic supplementary material

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

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long‐term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome‐wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue‐specific expression patterns indicated that the up‐regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up‐regulated selection genes were identified in leaves and seeds than in the other organs. Genome‐wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large‐scale characterization of key agronomic traits in this crop species.     

Genetic analysis of 103 candidate genes for coronary artery disease and associated phenotypes in a founder population reveals a new association between endothelin-1 and high-density lipoprotein cholesterol

Coronary artery disease (CAD) is a major health concern in both developed and developing countries. With a heritability estimated at ~50%, there is a strong rationale to better define the genetic contribution to CAD. This project involves the analysis of 884 individuals from 142 families (with average sibships of 5.7) as well as 558 case and control subjects from the Saguenay Lac St-Jean region of northeastern Quebec, with the use of 1,536 single-nucleotide polymorphisms (SNPs) in 103 candidate genes for CAD. By use of clusters of SNPs to generate multiallelic haplotypes at candidate loci for segregation studies within families, suggestive linkage for high-density lipoprotein (HDL) cholesterol is observed on chromosome 1p36.22. Furthermore, several associations that remain significant after Bonferroni correction are observed with lipoprotein-related traits as well as plasma concentrations of adiponectin. Of note, HDL cholesterol levels are associated with an amino acid substitution (lysine/asparagine) at codon 198 (rs5370) of endothelin-1 (EDN1) in a sex-specific manner, as well as with a SNP (rs2292318) located 7.7 kb upstream of lecithin cholesterol acyl-transferase (LCAT). Whereas the other observed associations are described in the current literature, these two are new. Using an independent validation sample of 806 individuals, we confirm the EDN1 association (P<.005), whereas the LCAT association was nonsignificant (P=.12).     

A note on a conditional-likelihood approach for family-based association studies of candidate genes

The family-based association study design is a variation of the case-control study design, where unaffected family members instead of unrelated subjects are sampled as controls. This variation is useful in assessing the effects of candidate genes on disease, because it avoids false associations caused by admixture of populations. A complication of this design is that because of an inherited genotypic correlation among family members, the genotypic distributions between cases and relative controls may be distorted by the ascertainment criteria of families, which could involve not only cases and relative controls, but also other relatives. Analyzing such data naively may lead to biased estimates of relative risk. In this note, we will discuss the consistency of a conditional-likelihood approach. We show analytically that maximum conditional-likelihood estimators are consistent for the true relative risks, if genotypes for family members are exchangeable under the sampling process, for example, sibling clusters. Besides being straightforward conceptually and computationally, this approach is robust to ascertainment bias and naturally accommodates genetic heterogeneity across families.     

Systems genetics approach reveals candidate genes for parasite resistance from quantitative trait loci studies in agricultural species

A systems genetics approach combining pathway analysis of quantitative trait loci (QTL) and gene expression information has provided strong evidence for common pathways associated with genetic resistance to internal parasites. Gene data, collected from published QTL regions in sheep, cattle, mice, rats and humans, and microarray data from sheep, were converted to human Entrez Gene IDs and compared to the KEGG pathway database. Selection of pathways from QTL data was based on a selection index that ensured that the selected pathways were in all species and the majority of the projects overall and within species. Pathways with either up- and down-regulated genes, primarily up-regulated genes or primarily down-regulated genes, were selected from gene expression data. After comparing the data sets independently, the pathways from each data set were compared and the common set of pathways and genes was identified. Comparisons within data sets identified 21 pathways from QTL data and 66 pathways from gene expression data. Both selected sets were enriched with pathways involved in immune functions, disease and cell responses to signals. The analysis identified 14 pathways that were common between QTL and gene expression data, and four directly associated with IFNγ or MHCII, with 31 common genes, including three MHCII genes. In conclusion, a systems genetics approach combining data from multiple QTL and gene expression projects led to the discovery of common pathways associated with genetic resistance to internal parasites. This systems genetics approach may prove significant for the discovery of candidate genes for many other multifactorial, economically important traits.     

Fun in facets: A flexible new tool for population genomics in R

The landscape of analytical tools for population genomics continues to evolve. However, these tools are scattered across programming languages, making them largely inaccessible for many biologists. In this issue of Molecular Ecology Resources, Hemstrom and Jones, 2022 (Mol Ecol Resour; 962) introduce a new R package, snpR. This package combines a large number of existing analyses, to provide a one-stop shop for population genomics. F-statistics, admixture analyses, effective population size inferences, genome-wide association studies (GWAS), and parentage analyses are all implemented natively within the package. A variety of third-party software can also be run without leaving the R environment. The authors pay particular attention to data structure – avoiding redundancy – and allowing analyses to be run across multiple sample or single-nucleotide polymorphism (SNP) groupings. Because of its great accessibility and wide range of analyses, snpR has the potential to become a favourite within the Molecular Ecology community.     

A novel approach for identifying candidate imprinted genes through sequence analysis of imprinted and control genes

Through the sequence analysis of 27 imprinted human genes and a set of 100 control genes we have developed a novel approach for identifying candidate imprinted genes based on the differences in sequence composition observed. The imprinted genes were found to be associated with significantly reduced numbers of short interspersed transposable element (SINE) Alus and mammalian-wide interspersed repeat (MIR) repeat elements, as previously reported. In addition, a significant association between imprinted genes and increased numbers of low-complexity repeats was also evident. Numbers of the Alu classes AluJ and AluS were found to be significantly depleted in some parts of the flanking regions of imprinted genes. A recent study has proposed that there is active selection against SINE elements in imprinted regions. Alternatively, there may be differences in the rates of insertion of Alu elements. Our study indicates that this difference extends both upstream and downstream of the coding region. This and other consistent differences between the sequence characteristics of imprinted and control genes has enabled us to develop discriminant analysis, which can be used to screen the genome for candidate imprinted genes. We have applied this function to a number of genes whose imprinting status is disputed or uncertain.     

A genome scan reveals QTL for growth, fatness, leanness and meat quality in a Duroc-Pietrain resource population

We performed a genome-wide QTL scan for production traits in a line cross between Duroc and Pietrain breeds of pigs, which included 585 F(2) progeny produced from 31 full-sib families genotyped with 106 informative microsatellites. A linkage map covering all 18 autosomes and spanning 1987 Kosambi cM was constructed. Thirty-five phenotypic traits including body weight, growth, carcass composition and meat quality traits were analysed using least square regression interval mapping. Twenty-four QTL exceeded the genome-wide significance threshold, while 47 QTL reached the suggestive threshold. These QTL were located at 28 genomic regions on 16 autosomal chromosomes and QTL in 11 regions were significant at the genome-wide level. A QTL affecting pH value in loin was detected on SSC1 between marker-interval S0312-S0113 with strong statistical support (P < 3.0 x 10(-14)); this QTL was also associated with meat colour and conductivity. QTL for carcass composition and average daily gain was also found on SSC1, suggesting multiple QTL. Seventeen genomic segments had only a single QTL that reached at least suggestive significance. Forty QTL exhibited additive inheritance whereas 31 QTL showed (over-) dominance effects. Two QTL for trait backfat thickness were detected on SSC2; a significant paternal effect was found for a QTL in the IGF2 region while another QTL in the middle of SSC2 showed Mendelian expression.     

A novel bioinformatics approach identifies candidate genes for the synthesis and feruloylation of arabinoxylan

Arabinoxylans (AXs) are major components of graminaceous plant cell walls, including those in the grain and straw of economically important cereals. Despite some recent advances in identifying the genes encoding biosynthetic enzymes for a number of other plant cell wall polysaccharides, the genes encoding enzymes of the final stages of AX synthesis have not been identified. We have therefore adopted a novel bioinformatics approach based on estimation of differential expression of orthologous genes between taxonomic divisions of species. Over 3 million public domain cereal and dicot expressed sequence tags were mapped onto the complete sets of rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) genes, respectively. It was assumed that genes in cereals involved in AX biosynthesis would be expressed at high levels and that their orthologs in dicotyledonous plants would be expressed at much lower levels. Considering all rice genes encoding putative glycosyl transferases (GTs) predicted to be integral membrane proteins, genes in the GT43, GT47, and GT61 families emerged as much the strongest candidates. When the search was widened to all other rice or Arabidopsis genes predicted to encode integral membrane proteins, cereal genes in Pfam family PF02458 emerged as candidates for the feruloylation of AX. Our analysis, known activities, and recent findings elsewhere are most consistent with genes in the GT43 families encoding beta-1,4-xylan synthases, genes in the GT47 family encoding xylan alpha-1,2- or alpha-1,3-arabinosyl transferases, and genes in the GT61 family encoding feruloyl-AX beta-1,2-xylosyl transferases.     

Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula

Bottom‐up evolutionary approaches, including geographically explicit population genomic analyses, have the power to reveal the mechanistic basis of adaptation. Here, we conduct a population genomic analysis in the model legume, Medicago truncatula, to characterize population genetic structure and identify symbiosis‐related genes showing evidence of spatially variable selection. Using RAD‐seq, we generated over 26,000 SNPs from 191 accessions from within three regions of the native range in Europe. Results from STRUCTURE analysis identify five distinct genetic clusters with divisions that separate east and west regions in the Mediterranean basin. Much of the genetic variation is maintained within sampling sites, and there is evidence for isolation by distance. Extensive linkage disequilibrium was identified, particularly within populations. We conducted genetic outlier analysis with F_ST‐based genome scans and a Bayesian modeling approach (PCAdapt). There were 70 core outlier loci shared between these distinct methods with one clear candidate symbiosis related gene, DMI1. This work sets that stage for functional experiments to determine the important phenotypes that selection has acted upon and complementary efforts in rhizobium populations.     

Plasticity genes and plasticity costs: a new approach using an Arabidopsis recombinant inbred population

Earlier flowering is triggered by vernalization in some but not all Arabidopsis ecotypes, often reflecting allelic variation at the FRIGIDA (FRI) locus. Using a recombinant inbred (RI) population polymorphic at FRI, we examined fitness consequences of variation for plasticity. Flowering and fitness were scored for 68 RI genotypes following full and partial vernalization treatments. Within-environment and mixed-model anovas estimated variance components for a genotype effect and a G x E term, respectively. Selection analyses examined whether delayed bolting increases fitness; a plasticity costs analysis asked whether increased plasticity lowers fitness. We also explored whether trait QTL had environment-specific effects, colocated in the immediate vicinity of FRI, or overlapped with fitness QTL. Selection may favor fri alleles and constitutive early flowering, especially in conditions that only partially vernalize plants. Plasticity costs, detected only after partial vernalization and only marginally significant, were nonetheless consistent with FRI-FLC function. We discuss how information about QTL with environment-specific effects, fitness QTL, and knowledge about plasticity genes can improve interpretation of selection or plasticity cost analyses.