Two Loci on Chromosome 5 Are Associated with Serum IgE Levels in Labrador Retrievers

Crosslinking of immunoglobulin E antibodies (IgE) bound at the surface of mast cells and subsequent mediator release is considered the most important trigger for allergic reactions. Therefore, the genetic control of IgE levels is studied in the context of allergic diseases, such as asthma, atopic rhinitis, or atopic dermatitis (AD). We performed genome-wide association studies in 161 Labrador Retrievers with regard to total and allergen-specific immunoglobulin E (IgE) levels. We identified a genome-wide significant association on CFA 5 with the antigen-specific IgE responsiveness to Acarus siro. We detected a second genome-wide significant association with respect to the antigen-specific IgE responsiveness to Tyrophagus putrescentiae at a different locus on chromosome 5. A. siro and T. putrescentiae both belong to the family Acaridae and represent so-called storage or forage mites. These forage mites are discussed as major allergen sources in canine AD. No obvious candidate gene for the regulation of IgE levels is located under the two association signals. Therefore our studies offer a chance of identifying a novel mechanism controlling the host's IgE response.

Congenital Hepatic Fibrosis in the Franches-Montagnes Horse Is Associated with the Polycystic Kidney and Hepatic Disease 1 (PKHD1) Gene

Congenital hepatic fibrosis has been described as a lethal disease with monogenic autosomal recessive inheritance in the Swiss Franches-Montagnes horse breed. We performed a genome-wide association study with 5 cases and 12 controls and detected an association on chromosome 20. Subsequent homozygosity mapping defined a critical interval of 952 kb harboring 10 annotated genes and loci including the polycystic kidney and hepatic disease 1 (autosomal recessive) gene (PKHD1). PKHD1 represents an excellent functional candidate as variants in this gene were identified in human patients with autosomal recessive polycystic kidney and hepatic disease (ARPKD) as well as several mouse and rat mutants. Whereas most pathogenic PKHD1 variants lead to polycystic defects in kidney and liver, a small subset of the human ARPKD patients have only liver symptoms, similar to our horses with congenital hepatic fibrosis. The PKHD1 gene is one of the largest genes in the genome with multiple alternative transcripts that have not yet been fully characterized. We sequenced the genomes of an affected foal and 46 control horses to establish a comprehensive list of variants in the critical interval. We identified two missense variants in the PKHD1 gene which were strongly, but not perfectly associated with congenital hepatic fibrosis. We speculate that reduced penetrance and/or potential epistatic interactions with hypothetical modifier genes may explain the imperfect association of the detected PKHD1 variants. Our data thus indicate that horses with congenital hepatic fibrosis represent an interesting large animal model for the liver-restricted subtype of human ARPKD.

A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH)

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p(raw) = 1.0×10(-13)) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p(raw) = 6.9×10(-10)). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.

The brown coat colour of Coppernecked goats is associated with a non-synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome-wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome-wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non-synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss-of-function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the (496) Asp allele might possibly act in a dominant manner. The (496) Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.

Deletion in the EVC2 gene causes chondrodysplastic dwarfism in Tyrolean Grey cattle.

During the summer of 2013 seven Italian Tyrolean Grey calves were born with abnormally short limbs. Detailed clinical and pathological examination revealed similarities to chondrodysplastic dwarfism. Pedigree analysis showed a common founder, assuming autosomal monogenic recessive transmission of the defective allele. A positional cloning approach combining genome wide association and homozygosity mapping identified a single 1.6 Mb genomic region on BTA 6 that was associated with the disease. Whole genome re-sequencing of an affected calf revealed a single candidate causal mutation in the Ellis van Creveld syndrome 2 (EVC2) gene. This gene is known to be associated with chondrodysplastic dwarfism in Japanese Brown cattle, and dwarfism, abnormal nails and teeth, and dysostosis in humans with Ellis-van Creveld syndrome. Sanger sequencing confirmed the presence of a 2 bp deletion in exon 19 (c.2993_2994ACdel) that led to a premature stop codon in the coding sequence of bovine EVC2, and was concordant with the recessive pattern of inheritance in affected and carrier animals. This loss of function mutation confirms the important role of EVC2 in bone development. Genetic testing can now be used to eliminate this form of chondrodysplastic dwarfism from Tyrolean Grey cattle.

Looking the cow in the eye: deletion in the NID1 gene is associated with recessive inherited cataract in Romagnola cattle

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.

A Deletion in the VLDLR Gene in Eurasier Dogs with Cerebellar Hypoplasia Resembling a Dandy-Walker-Like Malformation (DWLM).

Dandy-Walker-like malformation (DWLM) is the result of aberrant brain development and mainly characterized by cerebellar hypoplasia. DWLM affected dogs display a non-progressive cerebellar ataxia. Several DWLM cases were recently observed in the Eurasier dog breed, which strongly suggested a monogenic autosomal recessive inheritance in this breed. We performed a genome-wide association study (GWAS) with 9 cases and 11 controls and found the best association of DWLM with markers on chromosome 1. Subsequent homozygosity mapping confirmed that all 9 cases were homozygous for a shared haplotype in this region, which delineated a critical interval of 3.35 Mb. We sequenced the genome of an affected Eurasier and compared it with the Boxer reference genome and 47 control genomes of dogs from other breeds. This analysis revealed 4 private non-synonymous variants in the critical interval of the affected Eurasier. We genotyped these variants in additional dogs and found perfect association for only one of these variants, a single base deletion in the VLDLR gene encoding the very low density lipoprotein receptor. This variant, VLDLR:c.1713delC is predicted to cause a frameshift and premature stop codon (p.W572Gfs*10). Variants in the VLDLR gene have been shown to cause congenital cerebellar ataxia and mental retardation in human patients and Vldlr knockout mice also display an ataxia phenotype. Our combined genetic data together with the functional knowledge on the VLDLR gene from other species thus strongly suggest that VLDLR:c.1713delC is indeed causing DWLM in Eurasier dogs.

The interplay between host genetic variation, viral replication and microbial translocation in untreated HIV-infected individuals.

Systemic immune activation, a major determinant of HIV disease progression, is the result of a complex interplay between viral replication, dysregulation of the immune system, and microbial translocation due to gut mucosal damage. While human genetic variants influencing HIV viral load have been identified, it is unknown to what extent the host genetic background contributes to inter-individual differences in other determinants of HIV pathogenesis like gut damage and microbial translocation. Using samples and data from 717 untreated participants in the Swiss HIV Cohort Study and a genome-wide association study design, we searched for human genetic determinants of plasma levels of intestinal fatty-acid binding protein (I-FABP/FABP2), a marker of gut damage, and of soluble sCD14 (sCD14), a marker of LPS bioactivity and microbial translocation. We also assessed the correlations between HIV viral load, sCD14 and I-FABP. While we found no genome-wide significant determinant of the tested plasma markers, we observed strong associations between sCD14 and both HIV viral load and I-FABP, shedding new light on the relationships between processes that drive progression of untreated HIV infection.

Disease-modifying genes and monogenic disorders: experience in cystic fibrosis

The mechanisms responsible for the determination of phenotypes are still not well understood; however, it has become apparent that modifier genes must play a considerable role in the phenotypic heterogeneity of Mendelian disorders. Significant advances in genetic technologies and molecular medicine allow huge amounts of information to be generated from individual samples within a reasonable time frame. This review focuses on the role of modifier genes using the example of cystic fibrosis, the most common lethal autosomal recessive disorder in the white population, and discusses the advantages and limitations of candidate gene approaches versus genome-wide association studies. Moreover, the implications of modifier gene research for other monogenic disorders, as well as its significance for diagnostic, prognostic, and therapeutic approaches are summarized. Increasing insight into modifying mechanisms opens up new perspectives, dispelling the idea of genetic disorders being caused by one single gene.

Wattles in goats are associated with the FMN1/GREM1 region on chromosome 10.

The presence of congenital appendages (wattles) on the throat of goats is supposed to be under genetic control with a dominant mode of inheritance. Wattles contain a cartilaginous core covered with normal skin resembling early stages of extremities. To map the dominant caprine wattles (W) locus, we collected samples of 174 goats with wattles and 167 goats without wattles from nine different Swiss goat breeds. The samples were genotyped with the 53k goat SNP chip for a subsequent genome-wide association study. We obtained a single strong association signal on chromosome 10 in a region containing functional candidate genes for limb development and outgrowth. We sequenced the whole genomes of an informative family trio containing an offspring without wattles and its heterozygous parents with wattles. In the associated goat chromosome 10 region, a total of 1055 SNPs and short indels perfectly co-segregate with the W allele. None of the variants were perfectly associated with the phenotype after analyzing the genome sequences of eight additional goats. We speculate that the causative mutation is located in one of the numerous gaps in the current version of the goat reference sequence and/or represents a larger structural variant which influences the expression of the FMN1 and/or GREM1 genes. Also, we cannot rule out possible genetic or allelic heterogeneity. Our genetic findings support earlier assumptions that wattles are rudimentary developed extremities.

Impact of common risk factors of fibrosis progression in chronic hepatitis C.

OBJECTIVE The natural course of chronic hepatitis C varies widely. To improve the profiling of patients at risk of developing advanced liver disease, we assessed the relative contribution of factors for liver fibrosis progression in hepatitis C. DESIGN We analysed 1461 patients with chronic hepatitis C with an estimated date of infection and at least one liver biopsy. Risk factors for accelerated fibrosis progression rate (FPR), defined as ≥0.13 Metavir fibrosis units per year, were identified by logistic regression. Examined factors included age at infection, sex, route of infection, HCV genotype, body mass index (BMI), significant alcohol drinking (≥20 g/day for ≥5 years), HIV coinfection and diabetes. In a subgroup of 575 patients, we assessed the impact of single nucleotide polymorphisms previously associated with fibrosis progression in genome-wide association studies. Results were expressed as attributable fraction (AF) of risk for accelerated FPR. RESULTS Age at infection (AF 28.7%), sex (AF 8.2%), route of infection (AF 16.5%) and HCV genotype (AF 7.9%) contributed to accelerated FPR in the Swiss Hepatitis C Cohort Study, whereas significant alcohol drinking, anti-HIV, diabetes and BMI did not. In genotyped patients, variants at rs9380516 (TULP1), rs738409 (PNPLA3), rs4374383 (MERTK) (AF 19.2%) and rs910049 (major histocompatibility complex region) significantly added to the risk of accelerated FPR. Results were replicated in three additional independent cohorts, and a meta-analysis confirmed the role of age at infection, sex, route of infection, HCV genotype, rs738409, rs4374383 and rs910049 in accelerating FPR. CONCLUSIONS Most factors accelerating liver fibrosis progression in chronic hepatitis C are unmodifiable.

A conserved role for p48 homologs in protecting dopaminergic neurons from oxidative stress.

Parkinson's disease (PD) is the most common neurodegenerative movement disorder characterized by the progressive loss of dopaminergic (DA) neurons. Both environmental and genetic factors are thought to contribute to the pathogenesis of PD. Although several genes linked to rare familial PD have been identified, endogenous risk factors for sporadic PD, which account for the majority of PD cases, remain largely unknown. Genome-wide association studies have identified many single nucleotide polymorphisms associated with sporadic PD in neurodevelopmental genes including the transcription factor p48/ptf1a. Here we investigate whether p48 plays a role in the survival of DA neurons in Drosophila melanogaster and Caenorhabditis elegans. We show that a Drosophila p48 homolog, 48-related-2 (Fer2), is expressed in and required for the development and survival of DA neurons in the protocerebral anterior medial (PAM) cluster. Loss of Fer2 expression in adulthood causes progressive PAM neuron degeneration in aging flies along with mitochondrial dysfunction and elevated reactive oxygen species (ROS) production, leading to the progressive locomotor deficits. The oxidative stress challenge upregulates Fer2 expression and exacerbates the PAM neuron degeneration in Fer2 loss-of-function mutants. hlh-13, the worm homolog of p48, is also expressed in DA neurons. Unlike the fly counterpart, hlh-13 loss-of-function does not impair development or survival of DA neurons under normal growth conditions. Yet, similar to Fer2, hlh-13 expression is upregulated upon an acute oxidative challenge and is required for the survival of DA neurons under oxidative stress in adult worms. Taken together, our results indicate that p48 homologs share a role in protecting DA neurons from oxidative stress and degeneration, and suggest that loss-of-function of p48 homologs in flies and worms provides novel tools to study gene-environmental interactions affecting DA neuron survival.

A coding variant in RARG confers susceptibility to anthracycline-induced cardiotoxicity in childhood cancer

Anthracyclines are used in over 50% of childhood cancer treatment protocols, but their clinical usefulness is limited by anthracycline-induced cardiotoxicity (ACT) manifesting as asymptomatic cardiac dysfunction and congestive heart failure in up to 57% and 16% of patients, respectively. Candidate gene studies have reported genetic associations with ACT, but these studies have in general lacked robust patient numbers, independent replication or functional validation. Thus, the individual variability in ACT susceptibility remains largely unexplained. We performed a genome-wide association study in 280 patients of European ancestry treated for childhood cancer, with independent replication in similarly treated cohorts of 96 European and 80 non-European patients. We identified a nonsynonymous variant (rs2229774, p.Ser427Leu) in RARG highly associated with ACT (P = 5.9 × 10(-8), odds ratio (95% confidence interval) = 4.7 (2.7-8.3)). This variant alters RARG function, leading to derepression of the key ACT genetic determinant Top2b, and provides new insight into the pathophysiology of this severe adverse drug reaction.

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.