Genomic prediction for tuberculosis resistance in dairy cattle

Background: The increasing prevalence of bovine tuberculosis (bTB) in the UK and the limitations of the currently available diagnostic and control methods require the development of complementary approaches to assist in the sustainable control of the disease. One potential approach is the identification of animals that are genetically more resistant to bTB, to enable breeding of animals with enhanced resistance. This paper focuses on prediction of resistance to bTB. We explore estimation of direct genomic estimated breeding values (DGVs) for bTB resistance in UK dairy cattle, using dense SNP chip data, and test these genomic predictions for situations when disease phenotypes are not available on selection candidates. Methodology/Principal Findings: We estimated DGVs using genomic best linear unbiased prediction methodology, and assessed their predictive accuracies with a cross validation procedure and receiver operator characteristic (ROC) curves. Furthermore, these results were compared with theoretical expectations for prediction accuracy and area-under-the-ROC- curve (AUC). The dataset comprised 1151 Holstein-Friesian cows (bTB cases or controls). All individuals (592 cases and 559 controls) were genotyped for 727,252 loci (Illumina Bead Chip). The estimated observed heritability of bTB resistance was 0.23±0.06 (0.34 on the liability scale) and five-fold cross validation, replicated six times, provided a prediction accuracy of 0.33 (95% C.I.: 0.26, 0.40). ROC curves, and the resulting AUC, gave a probability of 0.58, averaged across six replicates, of correctly classifying cows as diseased or as healthy based on SNP chip genotype alone using these data. Conclusions/Significance: These results provide a first step in the investigation of the potential feasibility of genomic selection for bTB resistance using SNP data. Specifically, they demonstrate that genomic selection is possible, even in populations with no pedigree data and on animals lacking bTB phenotypes. However, a larger training population will be required to improve prediction accuracies. © 2014 Tsairidou et al.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated 1/42,000, 1/43,700 and 1/49,500 SNPs explained 1/421%, 1/424% and 1/429% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/I 2-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Digital Dermatitis in Dairy Cows: A Review of Risk Factors and Potential Sources of Between-Animal Variation in Susceptibility

Digital dermatitis (DD) is a bacterial disease that primarily affects the skin on the heels of cattle. It is a major cause of lameness in dairy cows and a significant problem for the dairy industry in many countries, causing reduced animal welfare and economic loss. A wide range of infection levels has been found on infected farms, prompting investigations into both farm level and animal level risk factors for DD occurrence. There also appears to be individual variation between animals in susceptibility to the disease. The identification of factors affecting individual variation in susceptibility to DD might allow changes in breeding policies or herd management which could be used to reduce DD prevalence. Factors mentioned in the literature as possibly influencing individual variation in susceptibility to DD include physical factors such as hoof conformation and properties of the skin, physiological factors such as the efficacy of the immune response, and behavioural factors such as standing half in cubicles. Further work is required to determine the influence of these factors, identify the genetic basis of variation, clarify the level of heritability of DD susceptibility and to determine how this is correlated with production and health traits currently used in breeding programmes.

Quantifying heritable variation in fitness-related traits of wild, farmed and hybrid Atlantic salmon families in a wild river environment

Farmed fish are typically genetically different from wild conspecifics. Escapees from fish farms may contribute one-way gene flow from farm to wild gene pools, which can depress population productivity, dilute local adaptations and disrupt coadapted gene complexes. Here, we reanalyse data from two experiments (McGinnity et al., 1997, 2003) where performance of Atlantic salmon (Salmo salar) progeny originating from experimental crosses between farm and wild parents (in three different cohorts) were measured in a natural stream under common garden conditions. Previous published analyses focussed on group-level differences but did not account for pedigree structure, as we do here using modern mixed-effect models. Offspring with one or two farm parents exhibited poorer survival in their first and second year of life compared with those with two wild parents and these group-level inferences were robust to excluding outlier families. Variation in performance among farm, hybrid and wild families was generally similar in magnitude. Farm offspring were generally larger at all life stages examined than wild offspring, but the differences were moderate (5–20%) and similar in magnitude in the wild versus hatchery environments. Quantitative genetic analyses conducted using a Bayesian framework revealed moderate heritability in juvenile fork length and mass and positive genetic correlations (>0.85) between these morphological traits. Our study confirms (using more rigorous statistical techniques) previous studies showing that offspring of wild fish invariably have higher fitness and contributes fresh insights into family-level variation in performance of farm, wild and hybrid Atlantic salmon families in the wild. It also adds to a small, but growing, number of studies that estimate key evolutionary parameters in wild salmonid populations. Such information is vital in modelling the impacts of introgression by escaped farm salmon.

Nuclear cataract shows significant familial aggregation in an older population after adjustment for possible shared environmental factors.

PURPOSE: To quantify the association between siblings in age-related nuclear cataract, after adjusting for known environmental and personal risk factors. METHODS: All participants (probands) in the Salisbury Eye Evaluation (SEE) project and their locally resident siblings underwent digital slit lamp photography and were administered a questionnaire to assess risk factors for cataract including: age, gender, lifetime sun exposure, smoking and diabetes history, and use of alcohol and medications such as estrogens and steroids. In addition, blood pressure, body mass index, and serum antioxidants were measured in all participants. Lens photographs were graded by trained observers masked to the subjects' identity, using the Wilmer Cataract Grading System. The odds ratio for siblings for affectedness with nuclear cataract and the sibling correlation of nuclear cataract grade, after adjusting for covariates, were estimated with generalized estimating equations. RESULTS: Among 307 probands (mean age, 77.6 +/- 4.5 years) and 434 full siblings (mean age, 72.4 +/- 7.4 years), the average sibship size was 2.7 per family. After adjustment for covariates, the probability of development of nuclear cataract was significantly increased (odds ratio [OR] = 2.07, 95% confidence interval [CI], 1.30-3.30) among individuals with a sibling with nuclear cataract (nuclear grade > or = 3.0). The final fitted model indicated a magnitude of heritability for nuclear cataract of 35.6% (95% CI: 21.0%-50.3%) after adjustment for the covariates. CONCLUSIONS: Findings in this study are consistent with a genetic effect for age-related nuclear cataract, a common and clinically significant form of lens opacity.

Cortical, but not posterior subcapsular, cataract shows significant familial aggregation in an older population after adjustment for possible shared environmental factors.

PURPOSE:

To quantify the risk for age-related cortical cataract and posterior subcapsular cataract (PSC) associated with having an affected sibling after adjusting for known environmental and personal risk factors.

DESIGN:

Sibling cohort study.

PARTICIPANTS:

Participants in the ongoing Salisbury Eye Evaluation (SEE) study (n = 321; mean age, 78.1+/-4.2 years) and their locally resident siblings (n = 453; mean age, 72.6+/-7.4 years) were recruited at the time of Rounds 3 and 4 of the SEE study. INTERVENTION/TESTING METHODS: Retroillumination photographs of the lens were graded for the presence of cortical cataract and PSC with the Wilmer grading system. The residual correlation between siblings' cataract grades was estimated after adjustment for a number of factors (age; gender; race; lifetime exposure to ultraviolet-B light; cigarette, alcohol, estrogen, and steroid use; serum antioxidants; history of diabetes; blood pressure; and body mass index) suspected to be associated with the presence of cataract.

RESULTS:

The average sibship size was 2.7 per family. Multivariate analysis revealed the magnitude of heritability (h(2)) for cortical cataract to be 24% (95% CI, 6%-42%), whereas that for PSC was not statistically significant (h(2) 4%; 95% CI, 0%-11%) after adjustment for the covariates. The model revealed that increasing age, female gender, a history of diabetes, and black race increased the odds of cortical cataract, whereas higher levels of provitamin A were protective. A history of diabetes and steroid use increased the odds for PSC.

CONCLUSIONS:

This study is consistent with a significant genetic effect for age-related cortical cataract but not PSC.

Contribution of genetic and environmental effects on lens thickness: the Guangzhou Twin Eye study.

PURPOSE:

This study investigated the heritability of lens thickness (LT) and relative lens thickness (LT/axial length, rLT) measured by Lenstar among Chinese children and adolescents in the Guangzhou Twin Eye study.

METHODS:

Twins aged 8 to 22 years were enrolled from the Guangzhou Twin Registry. A series of LT and axial length (AL) measurements using the Lenstar were taken for each twin. Zygosity was confirmed by genotyping in all same-sex twin pairs. Heritability was assessed by structural variance component genetic modeling, after adjustment for age and sex with the Mx program.

RESULTS:

Seven hundred sixty-eight twin pairs (482 monozygotic [MZ] and 286 dizygotic [DZ] twins) were available for data analysis. The mean (standard deviation) LT and rLT were 3.45 (0.18) mm and 0.142 (0.01), respectively. The intraclass correlation coefficients (ICCs) for LT were 0.90 for the MZ and 0.39 for the DZ twins; and those for rLT were 0.90 for the MZ and 0.40 for the DZ twins, respectively. The best-fitting model yielded 89.5% (95% CI: 87.8%-91.0%) of additive genetic effects and 10.5% (95% CI: 9.0%-12.2%) of unique environmental effects for LT, and 89.3% (95% CI: 89.2%-89.3%) of additive genetic effects and 10.7% (95% CI: 10.7%-11.4%) of unique environmental effects for rLT.

CONCLUSIONS:

This study confirms that the LT in young healthy subjects may be mainly affected by additive genetic factors. High heritability remains even when the data are corrected for the influence of AL with the use of rLT.

Genome-wide Association Study Identifies SESTD1 as a Novel Risk Gene for Lithium Responsive Bipolar Disorder

Lithium is the mainstay prophylactic treatment for bipolar disorder (BD), but treatment response varies considerably across individuals. Patients who respond well to lithium treatment might represent a relatively homogeneous subtype of this genetically and phenotypically diverse disorder. Here, we performed genome-wide association studies (GWAS) to identify (i) specific genetic variations influencing lithium response and (ii) genetic variants associated with risk for lithium-responsive BD. Patients with BD and controls were recruited from Sweden and the United Kingdom. GWAS were performed on 2698 patients with subjectively defined (self-reported) lithium response and 1176 patients with objectively defined (clinically documented) lithium response. We next conducted GWAS comparing lithium responders with healthy controls (1639 subjective responders and 8899 controls; 323 objective responders and 6684 controls). Meta-analyses of Swedish and UK results revealed no significant associations with lithium response within the bipolar subjects. However, when comparing lithium-responsive patients with controls, two imputed markers attained genome-wide significant associations, among which one was validated in confirmatory genotyping (rs116323614, P=2.74 × 10-8). It is an intronic single-nucleotide polymorphism (SNP) on chromosome 2q31.2 in the gene SEC14 and spectrin domains 1 (SESTD1), which encodes a protein involved in regulation of phospholipids. Phospholipids have been strongly implicated as lithium treatment targets. Furthermore, we estimated the proportion of variance for lithium-responsive BD explained by common variants ('SNP heritability') as 0.25 and 0.29 using two definitions of lithium response. Our results revealed a genetic variant in SESTD1 associated with risk for lithium-responsive BD, suggesting that the understanding of BD etiology could be furthered by focusing on this subtype of BD.

Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

MHC class I expression dependent on bacterial infection and parental factors in whitefish embryos (Salmonidae).

Ecological conditions can influence not only the expression of a phenotype, but also the heritability of a trait. As such, heritable variation for a trait needs to be studied across environments. We have investigated how pathogen challenge affects the expression of MHC genes in embryos of the lake whitefish Coregonus palaea. In order to experimentally separate paternal (i.e. genetic) from maternal and environmental effects, and determine whether and how stress affects the heritable variation for MHC expression, embryos were produced in full-factorial in vitro fertilizations, reared singly, and exposed at 208 degree days (late-eyed stage) to either one of two strains of Pseudomonas fluorescens that differ in their virulence characteristics (one increased mortality, while both delayed hatching time). Gene expression was assessed 48 h postinoculation, and virulence effects of the bacterial infection were monitored until hatching. We found no evidence of MHC class II expression at this stage of development. MHC class I expression was markedly down-regulated in reaction to both pseudomonads. While MHC expression could not be linked to embryo survival, the less the gene was expressed, the earlier the embryos hatched within each treatment group, possibly due to trade-offs between immune function and developmental rate or further factors that affect both hatching timing and MHC expression. We found significant additive genetic variance for MHC class I expression in some treatments. That is, changes in pathogen pressures could induce rapid evolution in MHC class I expression. However, we found no additive genetic variance in reaction norms in our study population.

Genome-wide profiling of blood pressure in adults and children.

Hypertension is an important determinant of cardiovascular morbidity and mortality and has a substantial heritability, which is likely of polygenic origin. The aim of this study was to assess to what extent multiple common genetic variants contribute to blood pressure regulation in both adults and children and to assess overlap in variants between different age groups, using genome-wide profiling. Single nucleotide polymorphism sets were defined based on a meta-analysis of genome-wide association studies on systolic blood pressure and diastolic blood pressure performed by the Cohort for Heart and Aging Research in Genome Epidemiology (n=29 136), using different P value thresholds for selecting single nucleotide polymorphisms. Subsequently, genetic risk scores for systolic blood pressure and diastolic blood pressure were calculated in an independent adult population (n=2072) and a child population (n=1034). The explained variance of the genetic risk scores was evaluated using linear regression models, including sex, age, and body mass index. Genetic risk scores, including also many nongenome-wide significant single nucleotide polymorphisms, explained more of the variance than scores based only on very significant single nucleotide polymorphisms in adults and children. Genetic risk scores significantly explained ≤1.2% (P=9.6*10(-8)) of the variance in adult systolic blood pressure and 0.8% (P=0.004) in children. For diastolic blood pressure, the variance explained was similar in adults and children (1.7% [P=8.9*10(-10)] and 1.4% [P=3.3*10(-5)], respectively). These findings suggest the presence of many genetic loci with small effects on blood pressure regulation both in adults and children, indicating also a (partly) common polygenic regulation of blood pressure throughout different periods of life.

Molecular autopsy in sudden cardiac death and its implication for families: discussion of the practical, legal and ethical aspects of the multidisciplinary collaboration.

Sudden cardiac death (SCD) is a major cause of premature death in young adults and children in developed countries. Standard forensic autopsy procedures are often unsuccessful in determining the cause of SCD. Post-mortem genetic testing, also called molecular autopsy, has revealed that a non-negligible number of these deaths are a result of inherited cardiac diseases, including arrhythmic disorders such as congenital long QT syndrome and Brugada syndrome. Due to the heritability of these diseases, the potential implications for living relatives must be taken into consideration. Advanced diagnostic analyses, genetic counselling, and interdisciplinary collaboration should be integral parts of clinical and forensic practice. In this article we present a multidisciplinary collaboration established in Lausanne, with the goal of properly informing families of these pathologies and their implications for surviving family members. In Switzerland, as in many other countries, legal guidelines for genetic testing do not address the use of molecular tools for post-mortem genetic analyses in forensic practice. In this article we present the standard practice guidelines established by our multidisciplinary team.

Systems Genetics of Metabolism: The Use of the BXD Murine Reference Panel for Multiscalar Integration of Traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.