Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs

Background

The recent completion of the swine genome sequencing project and development of a high density porcine SNP array has made genome-wide association (GWA) studies feasible in pigs.

Methodology/Principal Findings

Using Illumina''s PorcineSNP60 BeadChip, we performed a pilot GWA study in 820 commercial female pigs phenotyped for backfat, loin muscle area, body conformation in addition to feet and leg (FL) structural soundness traits. A total of 51,385 SNPs were jointly fitted using Bayesian techniques as random effects in a mixture model that assumed a known large proportion (99.5%) of SNPs had zero effect. SNP annotations were implemented through the Sus scrofa Build 9 available from pig Ensembl. We discovered a number of candidate chromosomal regions, and some of them corresponded to QTL regions previously reported. We not only have identified some well-known candidate genes for the traits of interest, such as MC4R (for backfat) and IGF2 (for loin muscle area), but also obtained novel promising genes, including CHCHD3 (for backfat), BMP2 (for loin muscle area, body size and several FL structure traits), and some HOXA family genes (for overall leg action). The candidate regions responsible for body conformation and FL structure soundness did not overlap greatly which implied that these traits were controlled by different genes. Functional clustering analyses classified the genes into categories related to bone and cartilage development, muscle growth and development or the insulin pathway suggesting the traits are regulated by common pathways or gene networks that exert roles at different spatial and temporal stages.

Conclusions/Significance

This study is one of the earliest GWA reports on important quantitative traits in pigs, and the findings will contribute to the further biological function analysis of the identified candidate genes and potential utilization of them in marker assisted selection.     

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.     

Identification of microRNAs associated with allergic airway disease using a genetically diverse mouse population

Background

Allergic airway diseases (AADs) such as asthma are characterized in part by granulocytic airway inflammation. The gene regulatory networks that govern granulocyte recruitment are poorly understood, but evidence is accruing that microRNAs (miRNAs) play an important role. To identify miRNAs that may underlie AADs, we used two complementary approaches that leveraged the genotypic and phenotypic diversity of the Collaborative Cross (CC) mouse population. In the first approach, we sought to identify miRNA expression quantitative trait loci (eQTL) that overlap QTL for AAD-related phenotypes. Specifically, CC founder strains and incipient lines of the CC were sensitized and challenged with house dust mite allergen followed by measurement of granulocyte recruitment to the lung. Total lung RNA was isolated and miRNA was measured using arrays for CC founders and qRT-PCR for incipient CC lines.

Results

Among CC founders, 92 miRNAs were differentially expressed. We measured the expression of 40 of the most highly expressed of these 92 miRNAs in the incipient lines of the CC and identified 18 eQTL corresponding to 14 different miRNAs. Surprisingly, half of these eQTL were distal to the corresponding miRNAs, and even on different chromosomes. One of the largest-effect local miRNA eQTL was for miR-342-3p, for which we identified putative causal variants by bioinformatic analysis of the effects of single nucleotide polymorphisms on RNA structure. None of the miRNA eQTL co-localized with QTL for eosinophil or neutrophil recruitment. In the second approach, we constructed putative miRNA/mRNA regulatory networks and identified three miRNAs (miR-497, miR-351 and miR-31) as candidate master regulators of genes associated with neutrophil recruitment. Analysis of a dataset from human keratinocytes transfected with a miR-31 inhibitor revealed two target genes in common with miR-31 targets correlated with neutrophils, namely Oxsr1 and Nsf.

Conclusions

miRNA expression in the allergically inflamed murine lung is regulated by genetic loci that are smaller in effect size compared to mRNA eQTL and often act in trans. Thus our results indicate that the genetic architecture of miRNA expression is different from mRNA expression. We identified three miRNAs, miR-497, miR-351 and miR-31, that are candidate master regulators of genes associated with neutrophil recruitment. Because miR-31 is expressed in airway epithelia and is predicted to target genes with known links to neutrophilic inflammation, we suggest that miR-31 is a potentially novel regulator of airway inflammation.

Electronic supplementary material

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

Prevalence of the Cryptosporidium pig genotype II in pigs from the Yangtze River Delta, China

Background

Cryptosporidium spp. is prevalent globally, pigs are an important Cryptosporidium reservoir. In China, little data regarding rates of Cryptosporidium infections in pigs are available. The present study was therefore aimed at characterizing the distribution of Cryptosporidium species in pigs from two different cities, Shaoxing and Shanghai, from the Yangtze River delta.

Methodology/Principal Findings

Nested PCR to amplify the 18S rRNA locus on DNA extracted from fecal samples (n = 94) revealed the positive rate of Cryptosporidium in pigs from two cities was approximately 17.0%. The positive rates in Shanghai and Shaoxing were 14.3% and 25.0% respectively. Amplified sequences were verified by sequencing. The identified strain belonged to the C. pig genotype II using BLAST analysis in the NCBI database.

Conclusion/Significance

Our finding of Cryptosporidium pig genotype II in pigs in the Yangtze River delta area suggests that pig farms in this region must be considered a public health threat and proper control measures be introduced.     

Methylomic analysis identifies frequent DNA methylation of zinc finger protein 582 (ZNF582) in cervical neoplasms

Background

Despite of the trend that the application of DNA methylation as a biomarker for cancer detection is promising, clinically applicable genes are few. Therefore, we looked for novel hypermethylated genes for cervical cancer screening.

Methods and Findings

At the discovery phase, we analyzed the methylation profiles of human cervical carcinomas and normal cervixes by methylated DNA immunoprecipitation coupled to promoter tiling arrays (MeDIP-on-chip). Methylation-specific PCR (MSP), quantitative MSP and bisulfite sequencing were used to verify the methylation status in cancer tissues and cervical scrapings from patients with different severities. Immunohistochemical staining of a cervical tissue microarray was used to confirm protein expression. We narrowed to three candidate genes: DBC1, PDE8B, and ZNF582; their methylation frequencies in tumors were 93%, 29%, and 100%, respectively. At the pre-validation phase, the methylation frequency of DBC1 and ZNF582 in cervical scraping correlated significantly with disease severity in an independent cohort (n = 330, both P<0.001). For the detection of cervical intraepithelial neoplasia 3 (CIN3) and worse, the area under the receiver operating characteristic curve (AUC) of ZNF582 was 0.82 (95% confidence interval  = 0.76–0.87).

Conclusions

Our study shows ZNF582 is frequently methylated in CIN3 and worse lesions, and it is demonstrated as a potential biomarker for the molecular screening of cervical cancer.     

Genome-wide linkage in a highly consanguineous pedigree reveals two novel loci on chromosome 7 for non-syndromic familial Premature Ovarian Failure

Background

The human condition known as Premature Ovarian Failure (POF) is characterized by loss of ovarian function before the age of 40. A majority of POF cases are sporadic, but 10–15% are familial, suggesting a genetic origin of the disease. Although several causal mutations have been identified, the etiology of POF is still unknown for about 90% of the patients.

Methodology/Principal Findings

We report a genome-wide linkage and homozygosity analysis in one large consanguineous Middle-Eastern POF-affected family presenting an autosomal recessive pattern of inheritance. We identified two regions with a LOD_max of 3.26 on chromosome 7p21.1-15.3 and 7q21.3-22.2, which are supported as candidate regions by homozygosity mapping. Sequencing of the coding exons and known regulatory sequences of three candidate genes (DLX5, DLX6 and DSS1) included within the largest region did not reveal any causal mutations.

Conclusions/Significance

We detect two novel POF-associated loci on human chromosome 7, opening the way to the identification of new genes involved in the control of ovarian development and function.     

Joint QTL linkage mapping for multiple-cross mating design sharing one common parent

Background

Nested association mapping (NAM) is a novel genetic mating design that combines the advantages of linkage analysis and association mapping. This design provides opportunities to study the inheritance of complex traits, but also requires more advanced statistical methods. In this paper, we present the detailed algorithm of a QTL linkage mapping method suitable for genetic populations derived from NAM designs. This method is called joint inclusive composite interval mapping (JICIM). Simulations were designed on the detected QTL in a maize NAM population and an Arabidopsis NAM population so as to evaluate the efficiency of the NAM design and the JICIM method.

Principal Findings

Fifty-two QTL were identified in the maize population, explaining 89% of the phenotypic variance of days to silking, and nine QTL were identified in the Arabidopsis population, explaining 83% of the phenotypic variance of flowering time. Simulations indicated that the detection power of these identified QTL was consistently high, especially for large-effect QTL. For rare QTL having significant effects in only one family, the power of correct detection within the 5 cM support interval was around 80% for 1-day effect QTL in the maize population, and for 3-day effect QTL in the Arabidopsis population. For smaller-effect QTL, the power diminished, e.g., it was around 50% for maize QTL with an effect of 0.5 day. When QTL were linked at a distance of 5 cM, the likelihood of mapping them as two distinct QTL was about 70% in the maize population. When the linkage distance was 1 cM, they were more likely mapped as one single QTL at an intermediary position.

Conclusions

Because it takes advantage of the large genetic variation among parental lines and the large population size, NAM is a powerful multiple-cross design for complex trait dissection. JICIM is an efficient and specialty method for the joint QTL linkage mapping of genetic populations derived from the NAM design.     

Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 (HMGA2) as a positional candidate gene

Background

Ear size and shape are distinct conformation characteristics of pig breeds. Previously, we identified a significant quantitative trait locus (QTL) influencing ear surface on pig chromosome 5 in a White Duroc × Erhualian F₂ resource population. This QTL explained more than 17% of the phenotypic variance.

Methods

Four new markers on pig chromosome 5 were genotyped across this F₂ population. RT-PCR was performed to obtain expression profiles of different candidate genes in ear tissue. Standard association test, marker-assisted association test and F-drop test were applied to determine the effects of single nucleotide polymorphisms (SNP) on ear size. Three synthetic commercial lines were also used for the association test.

Results

We refined the QTL to an 8.7-cM interval and identified three positional candidate genes i.e. HMGA2, SOX5 and PTHLH that are expressed in ear tissue. Seven SNP within these three candidate genes were selected and genotyped in the F₂ population. Of the seven SNP, HMGA2 SNP ({"type":"entrez-nucleotide","attrs":{"text":"JF748727","term_id":"347311053","term_text":"JF748727"}}JF748727: g.2836 A > G) showed the strongest association with ear size in the standard association test and marker-assisted association test. With the F-drop test, F value decreased by more than 97% only when the genotypes of HMGA2 g.2836 A > G were included as a fixed effect. Furthermore, the significant association between g.2836 A > G and ear size was also demonstrated in the synthetic commercial Sutai pig line. The haplotype-based association test showed that the phenotypic variance explained by HMGA2 was similar to that explained by the QTL and at a much higher level than by SOX5. More interestingly, HMGA2 is also located within the dog orthologous chromosome region, which has been shown to be associated with ear type and size.

Conclusions

HMGA2 was the closest gene with a potential functional effect to the QTL or marker for ear size on chromosome 5. This study will contribute to identify the causative gene and mutation underlying this QTL.     

The effectiveness of online cognitive behavioral treatment in routine clinical practice

Context

Randomized controlled trails have identified online cognitive behavioral therapy as an efficacious intervention in the management of common mental health disorders.

Objective

To assess the effectiveness of online CBT for different mental disorders in routine clinical practice.

Design

An uncontrolled before-after study, with measurements at baseline, posttest, 6-week follow-up, and 1-year follow-up.

Participants & Setting

1500 adult patients (female: 67%; mean age: 40 years) with a GP referral for psychotherapy were treated at a Dutch online mental health clinic for symptoms of depression (n = 413), panic disorder (n = 139), posttraumatic stress (n = 478), or burnout (n = 470).

Interventions

Manualized, web-based, therapist-assisted CBT, of which the efficacy was previously demonstrated in a series of controlled trials. Standardized duration of treatment varied from 5 weeks (online CBT for Posttraumatic stress) to 16 weeks (online CBT for Depression).

Main Outcome Measures

Validated self-report questionnaires of specific and general psychopathology, including the Beck Depression Inventory, the Impact of Event Scale, the Panic Disorder Severity Scale-Self Report, the Oldenburg Burnout Inventory, and the Depression Anxiety Stress Scales.

Results

Treatment adherence was 71% (n = 1071). Study attrition was 21% at posttest, 33% at 6-week FU and 65% at 1-year FU. Mixed-model repeated measures regression identified large short-term reductions in all measures of primary symptoms (d = 1.9±0.2 to d = 1.2±0.2; P<.001), which sustained up to one year after treatment. At posttest, rates of reliable improvement and recovery were 71% and 52% in the completer sample (full sample: 55%/40%). Patient satisfaction was high.

Conclusions

Results suggest that online therapist-assisted CBT may be as effective in routine practice as it is in clinical trials. Although pre-treatment withdrawal and long-term outcomes require further study, results warrant continued implementation of online CBT.     

High throughput interrogation of somatic mutations in high grade serous cancer of the ovary

Background

Epithelial ovarian cancer is the most lethal of all gynecologic malignancies, and high grade serous ovarian cancer (HGSC) is the most common subtype of ovarian cancer. The objective of this study was to determine the frequency and types of point somatic mutations in HGSC using a mutation detection protocol called OncoMap that employs mass spectrometric-based genotyping technology.

Methodology/Principal Findings

The Center for Cancer Genome Discovery (CCGD) Program at the Dana-Farber Cancer Institute (DFCI) has adapted a high-throughput genotyping platform to determine the mutation status of a large panel of known cancer genes. The mutation detection protocol, termed OncoMap has been expanded to detect more than 1000 mutations in 112 oncogenes in formalin-fixed paraffin-embedded (FFPE) tissue samples. We performed OncoMap on a set of 203 FFPE advanced staged HGSC specimens. We isolated genomic DNA from these samples, and after a battery of quality assurance tests, ran each of these samples on the OncoMap v3 platform. 56% (113/203) tumor samples harbored candidate mutations. Sixty-five samples had single mutations (32%) while the remaining samples had ≥2 mutations (24%). 196 candidate mutation calls were made in 50 genes. The most common somatic oncogene mutations were found in EGFR, KRAS, PDGRFα, KIT, and PIK3CA. Other mutations found in additional genes were found at lower frequencies (<3%).

Conclusions/Significance

Sequenom analysis using OncoMap on DNA extracted from FFPE ovarian cancer samples is feasible and leads to the detection of potentially druggable mutations. Screening HGSC for somatic mutations in oncogenes may lead to additional therapies for this patient population.     

Genomic library screens for genes involved in n-butanol tolerance in Escherichia coli

Background

n-Butanol is a promising emerging biofuel, and recent metabolic engineering efforts have demonstrated the use of several microbial hosts for its production. However, most organisms have very low tolerance to n-butanol (up to 2% (v/v)), limiting the economic viability of this biofuel. The rational engineering of more robust n-butanol production hosts relies upon understanding the mechanisms involved in tolerance. However, the existing knowledge of genes involved in n-butanol tolerance is limited. The goal of this study is therefore to identify E. coli genes that are involved in n-butanol tolerance.

Methodology/Principal Findings

Using a genomic library enrichment strategy, we identified approximately 270 genes that were enriched or depleted in n-butanol challenge. The effects of these candidate genes on n-butanol tolerance were experimentally determined using overexpression or deletion libraries. Among the 55 enriched genes tested, 11 were experimentally shown to confer enhanced tolerance to n-butanol when overexpressed compared to the wild-type. Among the 84 depleted genes tested, three conferred increased n-butanol resistance when deleted. The overexpressed genes that conferred the largest increase in n-butanol tolerance were related to iron transport and metabolism, entC and feoA, which increased the n-butanol tolerance by 32.8±4.0% and 49.1±3.3%, respectively. The deleted gene that resulted in the largest increase in resistance to n-butanol was astE, which enhanced n-butanol tolerance by 48.7±6.3%.

Conclusions/Significance

We identified and experimentally verified 14 genes that decreased the inhibitory effect of n-butanol tolerance on E. coli. From the data, we were able to expand the current knowledge on the genes involved in n-butanol tolerance; the results suggest that an increased iron transport and metabolism and decreased acid resistance may enhance n-butanol tolerance. The genes and mechanisms identified in this study will be helpful in the rational engineering of more robust biofuel producers.     

Real-time PCR improves Helicobacter pylori detection in patients with peptic ulcer bleeding

Background and Aims

Histological and rapid urease tests to detect H. pylori in biopsy specimens obtained during peptic ulcer bleeding episodes (PUB) often produce false-negative results. We aimed to examine whether immunohistochemistry and real-time PCR can improve the sensitivity of these biopsies.

Patients and Methods

We selected 52 histology-negative formalin-fixed paraffin-embedded biopsy specimens obtained during PUB episodes. Additional tests showed 10 were true negatives and 42 were false negatives. We also selected 17 histology-positive biopsy specimens obtained during PUB to use as controls. We performed immunohistochemistry staining and real-time PCR for 16S rRNA, ureA, and 23S rRNA for H. pylori genes on all specimens.

Results

All controls were positive for H. pylori on all PCR assays and immunohistochemical staining. Regarding the 52 initially negative biopsies, all PCR tests were significantly more sensitive than immunohistochemical staining (p<0.01). Sensitivity and specificity were 55% and 80% for 16S rRNA PCR, 43% and 90% for ureA PCR, 41% and 80% for 23S rRNA PCR, and 7% and 100% for immunohistochemical staining, respectively. Combined analysis of PCR assays for two genes were significantly more sensitive than ureA or 23S rRNA PCR tests alone (p<0.05) and marginally better than 16S rRNA PCR alone. The best combination was 16S rRNA+ureA, with a sensitivity of 64% and a specificity of 80%.

Conclusions

Real-time PCR improves the detection of H. pylori infection in histology-negative formalin-fixed paraffin-embedded biopsy samples obtained during PUB episodes. The low reported prevalence of H. pylori in PUB may be due to the failure of conventional tests to detect infection.     

Oncogenic CagA promotes gastric cancer risk via activating ERK signaling pathways: a nested case-control study

Background

CagA cellular interaction via activation of the ERK signaling pathway may be a starting point in the development of gastric cancer. This study aimed to evaluate whether genes involved in ERK downstream signaling pathways activated by CagA are susceptible genetic markers for gastric cancer.

Methods

In the discovery phase, a total of 580 SNPs within +/−5 kbp of 30 candidate genes were genotyped to examine an association with gastric cancer risk in the Korean Multi-center Cancer Cohort (100 incident gastric cancer case-control sets). The most significant SNPs (raw or permutated p value<0.02) identified in the discovery analysis were re-evaluated in the extension phase using unconditional logistic regression model (400 gastric cancer case-control sets). Combined analyses including pooled- and meta-analysis were conducted to summarize all the results.

Results

24 SNPs in eight genes (ERK, Dock180, C3G, Rap1, Src, CrkL, Mek and Crk) were significantly associated with gastric cancer risk in the individual SNP analyses in the discovery phase (p<0.05). In the extension analyses, ERK rs5999749, Dock180 rs4635002 and C3G rs7853122 showed marginally significant gene-dose effects for gastric cancer. Consistently, final combined analysis presented the SNPs as significantly associated with gastric cancer risk (OR = 1.56, [95% CI: 1.19–2.06], OR = 0.61, [95% CI: 0.43–0.87], OR = 0.59, [95% CI: 0.54–0.76], respectively).

Conclusions

Our findings suggest that ERK rs5999749, Dock180 rs4635002 and C3G rs7853122 are genetic determinants in gastric carcinogenesis.     

Genetic variation in cell death genes and risk of non-Hodgkin lymphoma

Background

Non-Hodgkin lymphomas are a heterogeneous group of solid tumours that constitute the 5^th highest cause of cancer mortality in the United States and Canada. Poor control of cell death in lymphocytes can lead to autoimmune disease or cancer, making genes involved in programmed cell death of lymphocytes logical candidate genes for lymphoma susceptibility.

Materials and Methods

We tested for genetic association with NHL and NHL subtypes, of SNPs in lymphocyte cell death genes using an established population-based study. 17 candidate genes were chosen based on biological function, with 123 SNPs tested. These included tagSNPs from HapMap and novel SNPs discovered by re-sequencing 47 cases in genes for which SNP representation was judged to be low. The main analysis, which estimated odds ratios by fitting data to an additive logistic regression model, used European ancestry samples that passed quality control measures (569 cases and 547 controls). A two-tiered approach for multiple testing correction was used: correction for number of tests within each gene by permutation-based methodology, followed by correction for the number of genes tested using the false discovery rate.

Results

Variant rs928883, near miR-155, showed an association (OR per A-allele: 2.80 [95% CI: 1.63–4.82]; p_F = 0.027) with marginal zone lymphoma that is significant after correction for multiple testing.

Conclusions

This is the first reported association between a germline polymorphism at a miRNA locus and lymphoma.