Genome-Wide Association Studies for Growth and Meat Production Traits in Sheep

Background

Growth and meat production traits are significant economic traits in sheep. The aim of the study is to identify candidate genes affecting growth and meat production traits at genome level with high throughput single nucleotide polymorphisms (SNP) genotyping technologies.

Methodology and Results

Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep for 11 growth and meat production traits (birth weight, weaning weight, 6-month weight, eye muscle area, fat thickness, pre-weaning gain, post-weaning gain, daily weight gain, height at withers, chest girth, and shin circumference). After quality control, 319 sheep and 48,198 SNPs were analyzed by TASSEL program in a mixed linear model (MLM). 36 significant SNPs were identified for 7 traits, and 10 of them reached genome-wise significance level for post-weaning gain. Gene annotation was implemented with the latest sheep genome Ovis_aries_v3.1 (released October 2012). More than one-third SNPs (14 out of 36) were located within ovine genes, others were located close to ovine genes (878bp-398,165bp apart). The strongest new finding is 5 genes were thought to be the most crucial candidate genes associated with post-weaning gain: s58995.1 was located within the ovine genes MEF2B and RFXANK, OAR3_84073899.1, OAR3_115712045.1 and OAR9_91721507.1 were located within CAMKMT, TRHDE, and RIPK2 respectively. GRM1, POL, MBD5, UBR2, RPL7 and SMC2 were thought to be the important candidate genes affecting post-weaning gain too. Additionally, 25 genes at chromosome-wise significance level were also forecasted to be the promising genes that influencing sheep growth and meat production traits.

Conclusions

The results will contribute to the similar studies and facilitate the potential utilization of genes involved in growth and meat production traits in sheep in future.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

Accuracy of genotype imputation in Nelore cattle

Background

Genotype imputation from low-density (LD) to high-density single nucleotide polymorphism (SNP) chips is an important step before applying genomic selection, since denser chips tend to provide more reliable genomic predictions. Imputation methods rely partially on linkage disequilibrium between markers to infer unobserved genotypes. Bos indicus cattle (e.g. Nelore breed) are characterized, in general, by lower levels of linkage disequilibrium between genetic markers at short distances, compared to taurine breeds. Thus, it is important to evaluate the accuracy of imputation to better define which imputation method and chip are most appropriate for genomic applications in indicine breeds.

Methods

Accuracy of genotype imputation in Nelore cattle was evaluated using different LD chips, imputation software and sets of animals. Twelve commercial and customized LD chips with densities ranging from 7 K to 75 K were tested. Customized LD chips were virtually designed taking into account minor allele frequency, linkage disequilibrium and distance between markers. Software programs FImpute and BEAGLE were applied to impute genotypes. From 995 bulls and 1247 cows that were genotyped with the Illumina® BovineHD chip (HD), 793 sires composed the reference set, and the remaining 202 younger sires and all the cows composed two separate validation sets for which genotypes were masked except for the SNPs of the LD chip that were to be tested.

Results

Imputation accuracy increased with the SNP density of the LD chip. However, the gain in accuracy with LD chips with more than 15 K SNPs was relatively small because accuracy was already high at this density. Commercial and customized LD chips with equivalent densities presented similar results. FImpute outperformed BEAGLE for all LD chips and validation sets. Regardless of the imputation software used, accuracy tended to increase as the relatedness between imputed and reference animals increased, especially for the 7 K chip.

Conclusions

If the Illumina® BovineHD is considered as the target chip for genomic applications in the Nelore breed, cost-effectiveness can be improved by genotyping part of the animals with a chip containing around 15 K useful SNPs and imputing their high-density missing genotypes with FImpute.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0069-1) contains supplementary material, which is available to authorized users.     

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

Background

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

Results

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

Conclusions

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Electronic supplementary material

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

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

Preliminary Investigations on Intradiscal Pressures during Daily Activities: An In Vivo Study Using the Merino Sheep

Purpose

Currently, no studies exist, which attest the suitability of the ovine intervertebral disc as a biomechanical in vivo model for preclinical tests of new therapeutic strategies of the human disc. By measuring the intradiscal pressure in vivo, the current study attempts to characterize an essential biomechanical parameter to provide a more comprehensive physiological understanding of the ovine intervertebral disc.

Methods

Intradiscal pressure (IDP) was measured for 24 hours within the discs L2-L3 and L4-L5 via a piezo-resistive pressure sensor in one merino sheep. The data were divided into an activity and a recovery phase and the corresponding average pressures for both phases were determined. Additionally, IDPs for different static and dynamic activities were analyzed and juxtaposed to human data published previously. After sacrificing the sheep, the forces corresponding to the measured IDPs were examined ex vivo in an axial compression test.

Results

The temporal patterns of IDP where pressure decreased during activity and increased during rest were comparable between humans and sheep. However, large differences were observed for different dynamic activities such as standing up or walking. Here, IDPs averaged 3.73 MPa and 1.60 MPa respectively, approximately two to four times higher in the ovine disc compared to human. These IDPs correspond to lower ex vivo derived axial compressive forces for the ovine disc in comparison to the human disc. For activity and rest, average ovine forces were 130 N and 58 N, compared to human forces of 400-600 N and 100 N, respectively.

Conclusions

In vivo IDPs were found to be higher in the ovine than in the human disc. In contrast, axial forces derived ex vivo were markedly lower in comparison to humans. Both should be considered in future preclinical tests of intradiscal therapies using the sheep. The techniques used in the current study may serve as a protocol for measuring IDP in a variety of large animal models.     

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array

Background

High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far.

Results

We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel.

Conclusions

The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-823) contains supplementary material, which is available to authorized users.     

A Neuropeptide Y Variant (rs16139) Associated with Major Depressive Disorder in Replicate Samples from Chinese Han Population

Objective

This study aimed to investigate the single nucleotide polymorphisms (SNPs) of neuropeptide Y (NPY) and major depressive disorder (MDD) in Chinese Han population.

Design

Prospective and randomized studies were carried out.

Patients

A total of 700 patients (324 male and 376 female; mean age = 40±14.9 years) with depression who met the diagnostic criteria of Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) and 673 healthy controls (313 male and 360 female; mean age = 41.9±17.2 years) were used to investigate the relationship between SNPs of NPY and the pathogenesis of MDD. A total of 417 patients (195 male and 202 female; mean age = 36±14.2 years) diagnosed with MDD and 314 healthy controls (153 male and 161 female; mean age = 37.9±14.2 years) from Chinese Han population were used to verify the relationship between SNPs of NPY and the pathogenesis of MDD.

Intervention and outcome

Ligase detection reactions were performed to detect the SNP sites of NPY. A series of statistical methods was carried out to investigate the correlation between the NPY gene SNP and MDD.

Results

Statistical analysis showed a significant correlation between the SNP sites rs16139 in NPY and the morbidity of depression. Patients with MDD have a lower frequency of A-allele in rs16139 in replicate samples from Chinese Han population. However, the frequency varied between male and female patients.

Conclusion

The gene polymorphism loci rs16139 was closely related to MDD in Chinese Han population.     

The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes

Background

The extent of linkage disequilibrium (LD) between molecular markers impacts genome-wide association studies and implementation of genomic selection. The availability of high-density single nucleotide polymorphism (SNP) genotyping platforms makes it possible to investigate LD at an unprecedented resolution. In this work, we characterised LD decay in breeds of beef cattle of taurine, indicine and composite origins and explored its variation across autosomes and the X chromosome.

Findings

In each breed, LD decayed rapidly and r² was less than 0.2 for marker pairs separated by 50 kb. The LD decay curves clustered into three groups of similar LD decay that distinguished the three main cattle types. At short distances between markers (< 10 kb), taurine breeds showed higher LD (r² = 0.45) than their indicine (r² = 0.25) and composite (r² = 0.32) counterparts. This higher LD in taurine breeds was attributed to a smaller effective population size and a stronger bottleneck during breed formation. Using all SNPs on only the X chromosome, the three cattle types could still be distinguished. However for taurine breeds, the LD decay on the X chromosome was much faster and the background level much lower than for indicine breeds and composite populations. When using only SNPs that were polymorphic in all breeds, the analysis of the X chromosome mimicked that of the autosomes.

Conclusions

The pattern of LD mirrored some aspects of the history of breed populations and showed a sharp decay with increasing physical distance between markers. We conclude that the availability of the HD chip can be used to detect association signals that remained hidden when using lower density genotyping platforms, since LD dropped below 0.2 at distances of 50 kb.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.     

Linkage disequilibrium in crossbred and pure line chickens

Background

Both genome-wide association (GWA) studies and genomic selection depend on the level of non-random association of alleles at different loci, i.e. linkage disequilibrium (LD), across the genome. Therefore, characterizing LD is of fundamental importance to implement both approaches. In this study, using a 60K single nucleotide polymorphism (SNP) panel, we estimated LD and haplotype structure in crossbred broiler chickens and their component pure lines (one male and two female lines) and calculated the consistency of LD between these populations.

Results

The average level of LD (measured by r²) between adjacent SNPs across the chicken autosomes studied here ranged from 0.34 to 0.40 in the pure lines but was only 0.24 in the crossbred populations, with 28.4% of adjacent SNP pairs having an r² higher than 0.3. Compared with the pure lines, the crossbred populations consistently showed a lower level of LD, smaller haploblock sizes and lower haplotype homozygosity on macro-, intermediate and micro-chromosomes. Furthermore, correlations of LD between markers at short distances (0 to 10 kb) were high between crossbred and pure lines (0.83 to 0.94).

Conclusions

Our results suggest that using crossbred populations instead of pure lines can be advantageous for high-resolution QTL (quantitative trait loci) mapping in GWA studies and to achieve good persistence of accuracy of genomic breeding values over generations in genomic selection. These results also provide useful information for the design and implementation of GWA studies and genomic selection using crossbred populations.

Electronic supplementary material

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

Single nucleotide polymorphism and haplotype effects associated with somatic cell score in German Holstein cattle

Background

To better understand the genetic determination of udder health, we performed a genome-wide association study (GWAS) on a population of 2354 German Holstein bulls for which daughter yield deviations (DYD) for somatic cell score (SCS) were available. For this study, we used genetic information of 44 576 informative single nucleotide polymorphisms (SNPs) and 11 725 inferred haplotype blocks.

Results

When accounting for the sub-structure of the analyzed population, 16 SNPs and 10 haplotypes in six genomic regions were significant at the Bonferroni threshold of P ≤ 1.14 × 10^-6. The size of the identified regions ranged from 0.05 to 5.62 Mb. Genomic regions on chromosomes 5, 6, 18 and 19 coincided with known QTL affecting SCS, while additional genomic regions were found on chromosomes 13 and X. Of particular interest is the region on chromosome 6 between 85 and 88 Mb, where QTL for mastitis traits and significant SNPs for SCS in different Holstein populations coincide with our results. In all identified regions, except for the region on chromosome X, significant SNPs were present in significant haplotypes. The minor alleles of identified SNPs on chromosomes 18 and 19, and the major alleles of SNPs on chromosomes 6 and X were favorable for a lower SCS. Differences in somatic cell count (SCC) between alternative SNP alleles reached 14 000 cells/mL.

Conclusions

The results support the polygenic nature of the genetic determination of SCS, confirm the importance of previously reported QTL, and provide evidence for the segregation of additional QTL for SCS in Holstein cattle. The small size of the regions identified here will facilitate the search for causal genetic variations that affect gene functions.     

Quantitative trait loci for resistance to trichostrongylid infection in Spanish Churra sheep

Background

For ruminants reared on grazing systems, gastrointestinal nematode (GIN) parasite infections represent the class of diseases with the greatest impact on animal health and productivity. Among the many possible strategies for controlling GIN infection, the enhancement of host resistance through the selection of resistant animals has been suggested by many authors. Because of the difficulty of routinely collecting phenotypic indicators of parasite resistance, information derived from molecular markers may be used to improve the efficiency of classical genetic breeding.

Methods

A total of 181 microsatellite markers evenly distributed along the 26 sheep autosomes were used in a genome scan analysis performed in a commercial population of Spanish Churra sheep to detect chromosomal regions associated with parasite resistance. Following a daughter design, we analysed 322 ewes distributed in eight half-sib families. The phenotypes studied included two faecal egg counts (LFEC0 and LFEC1), anti-Teladorsagia circumcincta LIV IgA levels (IgA) and serum pepsinogen levels (Peps).

Results

The regression analysis revealed one QTL at the 5% genome-wise significance level on chromosome 6 for LFEC1 within the marker interval BM4621-CSN3. This QTL was found to be segregating in three out of the eight families analysed. Four other QTL were identified at the 5% chromosome-wise level on chromosomes 1, 10 and 14. Three of these QTL influenced faecal egg count, and the other one had an effect on IgA levels.

Conclusion

This study has successfully identified segregating QTL for parasite resistance traits in a commercial population. For some of the QTL detected, we have identified interesting coincidences with QTL previously reported in sheep, although most of those studies have been focused on young animals. Some of these coincidences might indicate that some common underlying loci affect parasite resistance traits in different sheep breeds. The identification of new QTL may suggest the existence of complex host-parasite relationships that have unique features depending on the host-parasite combination, perhaps due to the different mechanisms underlying resistance in adult sheep (hypersensitivity reactions) and lambs (immunity). The most significant QTL identified on chromosome 6 for LFEC₁ may be the target for future fine-mapping research efforts.     

Regional Variation in Tissue Composition and Biomechanical Properties of Postmenopausal Ovine and Human Vagina

Objective

There are increasing numbers of reports describing human vaginal tissue composition in women with and without pelvic organ prolapse with conflicting results. The aim of this study was to compare ovine and human posterior vaginal tissue in terms of histological and biochemical tissue composition and to assess passive biomechanical properties of ovine vagina to further characterise this animal model for pelvic organ prolapse research.

Study Design

Vaginal tissue was collected from ovariectomised sheep (n = 6) and from postmenopausal women (n = 7) from the proximal, middle and distal thirds. Tissue histology was analyzed using Masson''s Trichrome staining; total collagen was quantified by hydroxyproline assays, collagen III/I+III ratios by delayed reduction SDS PAGE, glycosaminoglycans by dimethylmethylene blue assay, and elastic tissue associated proteins (ETAP) by amino acid analysis. Young''s modulus, maximum stress/strain, and permanent strain following cyclic loading were determined in ovine vagina.

Results

Both sheep and human vaginal tissue showed comparable tissue composition. Ovine vaginal tissue showed significantly higher total collagen and glycosaminoglycan values (p<0.05) nearest the cervix. No significant differences were found along the length of the human vagina for collagen, GAG or ETAP content. The proximal region was the stiffest (Young''s modulus, p<0.05), strongest (maximum stress, p<0.05) compared to distal region, and most elastic (permanent strain).

Conclusion

Sheep tissue composition and mechanical properties showed regional differences along the postmenopausal vaginal wall not apparent in human vagina, although the absolute content of proteins were similar. Knowledge of this baseline variation in the composition and mechanical properties of the vaginal wall will assist future studies using sheep as a model for vaginal surgery.     

Genome-wide association study combined with biological context can reveal more disease-related SNPs altering microRNA target seed sites

Background

Emerging studies demonstrate that single nucleotide polymorphisms (SNPs) resided in the microRNA recognition element seed sites (MRESSs) in 3′UTR of mRNAs are putative biomarkers for human diseases and cancers. However, exhaustively experimental validation for the causality of MRESS SNPs is impractical. Therefore bioinformatics have been introduced to predict causal MRESS SNPs. Genome-wide association study (GWAS) provides a way to detect susceptibility of millions of SNPs simultaneously by taking linkage disequilibrium (LD) into account, but the multiple-testing corrections implemented to suppress false positive rate always sacrificed the sensitivity. In our study, we proposed a method to identify candidate causal MRESS SNPs from 12 GWAS datasets without performing multiple-testing corrections. Alternatively, we used biological context to ensure credibility of the selected SNPs.

Results

In 11 out of the 12 GWAS datasets, MRESS SNPs were over-represented in SNPs with p-value ≤ 0.05 (odds ratio (OR) ranged from 1.1 to 2.4). Moreover, host genes of susceptible MRESS SNPs in each of the 11 GWAS dataset shared biological context with reported causal genes. There were 286 MRESS SNPs identified by our method, while only 13 SNPs were identified by multiple-testing corrections with a given threshold of 1 × 10⁻⁵, which is a common cutoff used in GWAS. 27 out of the 286 candidate SNPs have been reported to be deleterious while only 2 out of 13 multiple-testing corrected SNPs were documented in PubMed. MicroRNA-mRNA interactions affected by the 286 candidate SNPs were likely to present negatively correlated expression. These SNPs introduced greater alternation of binding free energy than other MRESS SNPs, especially when grouping by haplotypes (4210 vs. 4105 cal/mol by mean, 9781 vs. 8521 cal/mol by mean, respectively).

Conclusions

MRESS SNPs are promising disease biomarkers in multiple GWAS datasets. The method of integrating GWAS p-value and biological context is stable and effective for selecting candidate causal MRESS SNPs, it reduces the loss of sensitivity compared to multiple-testing corrections. The 286 candidate causal MRESS SNPs provide researchers a credible source to initialize their design of experimental validations in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-669) contains supplementary material, which is available to authorized users.     

Hyperfractionated high-dose rate brachytherapy in the treatment of oral tongue cancer

Background

Low-dose rate brachytherapy is a well established treatment modality of oral cancer. Data about high-dose rate (HDR) brachytherapy are still sparse with various fractionation schedules and heterogeneous results.

Aim

The aim of our retrospective study was to evaluate the results of HDR brachytherapy with doses of 3 Gy twice daily.

Patients and methods

Twenty patients with squamous cell tongue cancer were treated in the years 2001–2009 by exclusive HDR BT 18 × 3 Gy twice daily. The plastic tube technique was used. Median follow up was 47 months (7.8–118) since brachytherapy.

Results

The local and locoregional control was 85% and 68%, respectively. Bone necrosis developed in one case treated without mandibular shielding and soft tissue necrosis in 2 cases.

Conclusion

It can be concluded that HDR brachytherapy with 18 × 3 Gy twice daily is safe with promising local control. The risk of nodal recurrences is substantial.     

Nebulized perflubron and carbon dioxide rapidly dilate constricted airways in an ovine model of allergic asthma

Background

The low toxicity of perfluorocarbons (PFCs), their high affinity for respiratory gases and their compatibility with lung surfactant have made them useful candidates for treating respiratory diseases such as adult respiratory distress syndrome. We report results for treating acute allergic and non-allergic bronchoconstriction in sheep using S-1226 (a gas mixture containing carbon dioxide and small volumes of nebulized perflubron). The carbon dioxide, which is highly soluble in perflubron, was used to relax airway smooth muscle.

Methods

Sheep previously sensitized to house dust mite (HDM) were challenged with HDM aerosols to induce early asthmatic responses. At the maximal responses (characterised by an increase in lung resistance), the sheep were either not treated or treated with one of the following; nebulized S-1226 (perflubron + 12% CO₂), nebulized perflubron + medical air, 12% CO₂, salbutamol or medical air. Lung resistance was monitored for up to 20 minutes after cessation of treatment.In additional naïve sheep, a segmental bronchus was pre-contracted with methacholine (MCh) and treated with nebulized S-1226 administered via a bronchoscope catheter. Subsequent bronchodilatation was monitored by real time digital video recording.

Results

Treatment with S-1226 for 2 minutes following HDM challenge resulted in a more rapid, more profound and more prolonged decline in lung resistance compared with the other treatment interventions. Video bronchoscopy showed an immediate and complete (within 5 seconds) re-opening of MCh-constricted airways following treatment with S-1226.

Conclusions

S-1226 is a potent and rapid formulation for re-opening constricted airways. Its mechanism(s) of action are unknown. The formulation has potential as a rescue treatment for acute severe asthma.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0098-x) contains supplementary material, which is available to authorized users.     

A Common Variant in the SIAH2 Locus Is Associated with Estrogen Receptor-Positive Breast Cancer in the Chinese Han Population

Background

Genome-wide association studies (GWAS) have identified many loci associated with breast cancer risk. These studies have primarily been conducted in populations of European descent.

Objective

To determine whether previously reported susceptibility loci in other ethnic groups are also risk factors for breast cancer in a Chinese population.

Method

We genotyped 21 previously reported single nucleotide polymorphisms (SNPs) within a female Chinese cohort of 1203 breast cancer cases and 2525 healthy controls using the Sequenom iPlex platform. Fourteen SNPs passed the quality control test. These SNPs were subjected to statistical analysis for the entire cohort and were further analyzed for estrogen receptor (ER) status. The associations of the SNPs with disease susceptibility were assessed using logistic regression, adjusting for age. The Bonferroni correction was used to conservatively account for multiple testing, and the threshold for statistical significance was P<3.57×10⁻³ (0.05/14).

Result

Although none of the SNPs showed an overall association with breast cancer, an analysis of the ER status of the breast cancer patients revealed that the SIAH2 locus (rs6788895; P = 5.73×10⁻⁴, odds ratio [OR] = 0.81) is associated with ER-positive breast cancer.

Conclusion

A common variant in the SIAH2 locus is associated with ER-positive breast cancer in the Chinese Han population. The replication of published GWAS results in other ethnic groups provides important information regarding the genetic etiology of breast cancer.     

A genome-wide association study of limb bone length using a Large White?×?Minzhu intercross population

Background

In pig, limb bone length influences ham yield and body height to a great extent and has important economic implications for pig industry. In this study, an intercross population was constructed between the indigenous Chinese Minzhu pig breed and the western commercial Large White pig breed to examine the genetic basis for variation in limb bone length. The aim of this study was to detect potential genetic variants associated with porcine limb bone length.

Methods

A total of 571 F2 individuals from a Large White and Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for femur length (FL), humerus length (HL), hipbone length (HIPL), scapula length (SL), tibia length (TL), and ulna length (UL). A genome-wide association study was performed by applying the previously reported approach of genome-wide rapid association using mixed model and regression. Statistical significance of the associations was based on Bonferroni-corrected P-values.

Results

A total of 39 significant SNPs were mapped to a 11.93 Mb long region on pig chromosome 7 (SSC7). Linkage analysis of these significant SNPs revealed three haplotype blocks of 495 kb, 376 kb and 492 kb, respectively, in the 11.93 Mb region. Annotation based on the pig reference genome identified 15 genes that were located near or contained the significant SNPs in these linkage disequilibrium intervals. Conditioned analysis revealed that four SNPs, one on SSC2 and three on SSC4, showed significant associations with SL and HL, respectively.

Conclusions

Analysis of the 15 annotated genes that were identified in these three haplotype blocks indicated that HMGA1 and PPARD, which are expressed in limbs and influence chondrocyte cell growth and differentiation, could be considered as relevant biological candidates for limb bone length in pig, with potential applications in breeding programs. Our results may also be useful for the study of the mechanisms that underlie human limb length and body height.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0056-6) contains supplementary material, which is available to authorized users.