Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency

A retrospective multi-center study of treatment, outcome, and prognostic factors in 34 dogs with disseminated aspergillosis in Australia

BACKGROUND

Disseminated aspergillosis (DA) in dogs has a guarded prognosis and there is a lack of a gold standard treatment protocol.

OBJECTIVE

To retrospectively assess survival times and factors influencing survival times.

ANIMALS

Dogs diagnosed with DA from January 2007 to June 2017.

METHODS

Disseminated aspergillosis case data were retrieved from 13 Australian veterinary referral centers, with a diagnosis confirmed with culture or PCR. Factors influencing survival time after diagnosis were quantified using a Cox proportional hazards regression model.

RESULTS

Thirty-four dogs met the study inclusion criteria. Twenty-two dogs were treated with antifungal treatment and 12 dogs received no antifungal treatment. Accounting for censoring of dogs that were either still alive on the date of data collection or were loss to follow-up, dogs treated with itraconazole alone (n = 8) had a median survival time (MST) of 63 (95% CI: 20-272) days compared to 830 (95% CI: 267-1259) days for the n = 14 dogs that received multimodal antifungal therapy ( χ 2 test statistic 8.6; df = 1; P < .01). The daily hazard of death (DHOD) for dogs with abnormally high serum creatinine concentration at the time of diagnosis was 7.4 (95% CI: 1.9-29) times that of dogs with serum creatinine within the reference interval.

CONCLUSION AND CLINICAL IMPORTANCE

Serum creatinine concentration at the time of diagnosis is a useful prognostic indicator for survival after a diagnosis of DA. The MST for dogs treated with multimodal antifungal therapy is longer than itraconazole alone and warrant further investigation (P < .01).

Group Feature Screening via the F Statistic

Feature screening is crucial in the analysis of ultrahigh dimensional data, where the number of variables (features) is in an exponential order of the number of observations. In various ultrahigh dimensional data, variables are naturally grouped, giving us a good rationale to develop a screening method using joint effect of multiple variables. In this article, we propose a group screening procedure via the F-test statistic. The proposed method is a direct extension of the original sure independence screening procedure, when the group information is known, for example, from prior knowledge. Under certain regularity conditions, we prove that the proposed group screening procedure possesses the sure screening property that selects all effective groups with a probability approaching one at an exponential rate. We use simulations to demonstrate the advantages of the proposed method and show its application in a genome-wide association study. We conclude that the grouping method is very useful in the analysis of ultrahigh dimensional data, as the optimal F-test can detect true signals with desired properties.

Genomic association and further characterisation of faecal immunoglobulin A deficiency in German Shepherd dogs

BACKGROUND

Immunoglobulin A (IgA) deficiency, chronic enteropathies and exocrine pancreatic insufficiency (EPI) have a high prevalence in German Shepherd dogs (GSD). This prospective study determined the prevalence of faecal IgA deficiency (IgAD) in GSD and investigated several candidate genes and the canine genome for a region or locus co-segregating with IgAD in GSD. Faecal IgA concentrations were quantified and genomic DNA was extracted from 8 GSD with an undetectable faecal IgA (classified as IgAD) and 80 non-IgAD GSD. The canine minimal screening set II microsatellite markers were genotyped, with evidence of an association at p < 1.0 × 10-3 . Faecal IgA concentrations were also tested for an association with patient clinical and biochemical variables.

RESULTS

Allele frequencies observed using the candidate gene approach were not associated with faecal IgAD in GSD. In the genome-wide association study (GWAS), the microsatellite marker FH2361 on canine chromosome 33 approached statistical significance for a link with IgAD in GSD (p = 1.2 × 10-3 ). A subsequent GWAS in 11 GSD with EPI and 80 control GSD revealed a significant association between EPI and FH2361 (p = 8.2 × 10-4 ).

CONCLUSIONS

The lack of an association with the phenotype of faecal IgAD in GSD using the candidate gene approach and GWAS might suggests that faecal IgAD in GSD is a relative or transient state of deficiency. However, the prevalence of faecal IgAD in GSD appears to be low (<3%). The relationship between faecal IgAD, EPI and loci close to FH2361 on canine chromosome 33 in GSD warrants further investigation.

Innate Mechanisms in Selective IgA Deficiency

Selective IgA deficiency (SIgAD), characterized by a serum IgA level below 0.07 mg/ml, while displaying normal serum levels of IgM and IgG antibodies, is the most frequently occurring primary immunodeficiency that reveals itself after the first four years after birth. These individuals with SIgAD are for the majority healthy and even when they are identified they are usually not investigated further or followed up. However, recent studies show that newborns and young infants already display clinical manifestations of this condition due to aberrancies in their immune defense. Interestingly, there is a huge heterogeneity in the clinical symptoms of the affected individuals. More than 50% of the affected individuals do not have clinical symptoms, while the individuals that do show clinical symptoms can suffer from mild to severe infections, allergies and autoimmune diseases. However, the reason for this heterogeneity in the manifestation of clinical symptoms of the individuals with SIgAD is unknown. Therefore, this review focusses on the characteristics of innate immune system driving T-cell independent IgA production and providing a mechanism underlying the development of SIgAD. Thereby, we focus on some important genes, including TNFRSF13B (encoding TACI), associated with SIgAD and the involvement of epigenetics, which will cover the methylation degree of TNFRSF13B, and environmental factors, including the gut microbiota, in the development of SIgAD. Currently, no specific treatment for SIgAD exists and novel therapeutic strategies could be developed based on the discussed information.

SERPINA9 and SERPINB2: Novel Cartilage Lineage Differentiation Markers of Human Mesenchymal Stem Cells with Kartogenin

OBJECTIVE

Human mesenchymal stem cells (hMSCs) are a promising source for regenerative medicine, especially mesodermal lineages. Clinical applications require an understanding of the mechanisms for transcriptional control to maintain the desired cell type. The aim of this study was to identify novel markers for differentiation of hMSCs into bone or cartilage with the use of Kartogenin, by RNA analysis using microarray technology, and explore the role of RhoA-Rho associated protein kinase (ROCK) inhibition in these.

METHODS

Commercial human bone marrow derived primary mesenchymal stem cells were purchased from ATCC. Cells were differentiated in vitro in 2-dimensional cultures using Kartogenin as the main cartilage inducer and bone morphogenetic protein 2 for bone differentiation; cells were cultured with and without ROCK inhibitor Y-27632. After 21 days of culture, whole RNA was extracted and analyzed via Affimetrix microarrays. The most significant hits were validated by quantitative polymerase chain reaction.

RESULTS

We found a total of 1,757 genes that were either up- or downregulated on differentiation, when compared to P1 hMSC (control) at day 0 of differentiation. Two members of the Serpin superfamily, SERPINA9 and SERPINB2, were significantly upregulated in the cartilage groups, whereas they were unchanged in the bone groups with and without ROCK inhibition.

CONCLUSIONS

SERPINA9 and SERPINB2 are novel differentiation markers, and molecular regulator candidates for hMSC lineage commitment toward bone and cartilage, providing a new tool for regenerative medicine. Our study highlights the roles of these 2 genes, with significant upregulation of both in cell cultures stimulated with Kartogenin.

Chromosomal segments may explain the antibody response cooperation for canine leishmaniasis pathogenesis

Visceral leishmaniasis (VL) is marked by hyperactivation of a humoral response secreting high quantity of immunoglobulins (Igs) that are inaccessible to intracellular parasites. Here we investigated the contributions of the antibody response to the canine leishmaniasis pathogenesis. Using correlation and genome-wide association analysis, we investigated the relationship of anti-Leishmania infantum immunoglobulin classes levels with parasite burden, clinical response, renal/hepatic biochemical, and oxidative stress markers in dogs from endemic areas of VL. Immunoglobulin G (IgG) and IgA were positively correlated with parasite burden on lymph node and blood. Increased IgG, IgA and IgE levels were associated with severe canine leishmaniasis (CanL) whereas IgM was elevated in uninfected exposed dogs. Correlations of IgM, IgG and IgA with creatinine, urea, AST and ALT levels in the serum were suggested an involvement of those Igs with renal and hepatic changes. The correlogram of oxidative radicals and antioxidants revealed a likely relationship of IgM, IgG and IgA with oxidative stress and lipid peroxidation in the blood, suggested as mechanisms mediating tissue damage and CanL worsening. The gene mapping on chromosomal segments associated with the quantitative variation of immunoglobulin classes identified genetic signatures involved with reactive oxygen species generation, phagolysosome maturation and rupture, free iron availability, Th1/Th2 differenciation and, immunoglobulin clearance. The findings demonstrated the roles of the antibody response as resistance or susceptibility markers and mediators of CanL pathogenesis. In addition we pinpointed candidate genes as potential targets for the therapy against the damage caused by exacerbated antibody response and parasitism in VL.

Quality of DNA extracted from formalin-fixed, paraffin-embedded canine tissues

Veterinary pathology tissue banks are valuable resources for genetic studies. However, limited data exist as to whether quality DNA can be extracted from these tissues for use in canine genotyping studies. We extracted DNA from 44 formalin-fixed, paraffin-embedded (FFPE) tissue blocks from dogs; 9 of these dogs had DNA available from whole blood samples that had been banked. We genotyped DNA from 30 of 44 tissue blocks and 9 whole blood samples on the Illumina CanineHD BeadChip; DNA quality was insufficient in 14 of 44 samples from tissue blocks. There was significant correlation between the 260/280 ratio and single-nucleotide variation (SNV) call rate (p = 0.0276; r² = 0.162); 23 of 30 samples from FFPE were genotyped with > 65% call rates. Median pairwise identical-by-state (IBS) analysis was 0.99 in 8 pairs of dogs with call rates > 65%. Neither age of tissue block nor specific tissue types were associated with significant differences in DNA concentration, 260/280 ratio, or SNV call rate. DNA extracted from tissue blocks can have variable quality, although comparable levels of homozygosity suggest that extracts from FFPE with call rates > 65% might provide similar results to samples from whole blood when analyzed on the Illumina CanineHD BeadChip.

Whole-genome genotyping and resequencing reveal the association of a deletion in the complex interferon alpha gene cluster with hypothyroidism in dogs

Background

Hypothyroidism is a common complex endocrinopathy that typically has an autoimmune etiology, and it affects both humans and dogs. Genetic and environmental factors are both known to play important roles in the disease development. In this study, we sought to identify the genetic risk factors potentially involved in the susceptibility to the disease in the high-risk Giant Schnauzer dog breed.

Results

By employing genome-wide association followed by fine-mapping (top variant p-value = 5.7 × 10^− 6), integrated with whole-genome resequencing and copy number variation analysis, we detected a ~ 8.9 kbp deletion strongly associated (p-value = 0.0001) with protection against development of hypothyroidism. The deletion is located between two predicted Interferon alpha (IFNA) genes and it may eliminate functional elements potentially involved in the transcriptional regulation of these genes. Remarkably, type I IFNs have been extensively associated to human autoimmune hypothyroidism and general autoimmunity. Nonetheless, the extreme genomic complexity of the associated region on CFA11 warrants further long-read sequencing and annotation efforts in order to ascribe functions to the identified deletion and to characterize the canine IFNA gene cluster in more detail.

Conclusions

Our results expand the current knowledge on genetic determinants of canine hypothyroidism by revealing a significant link with the human counterpart disease, potentially translating into better diagnostic tools across species, and may contribute to improved canine breeding strategies.

Canine IgA and IgA deficiency: Implications for immunization against respiratory pathogens

Immunoglobulin A (IgA) is widely recognized as the important antibody isotype involved in protective responses on mucosal surfaces, where it acts primarily by effectuating immune exclusion of foreign material. Selective IgA deficiency (SIgAD) is the most common immunodeficiency disease in dogs and humans and has consequences for mucosal immunity. This review is a comparative look at the biology of IgA and SIgAD with a focus on how this branch of immunology relates to vaccine selection and efficacy for canine infectious respiratory disease.

Humoral immune responses against gut bacteria in dogs with inflammatory bowel disease

Inflammatory bowel disease (IBD) in dogs is associated with clinical signs of intestinal dysfunction, as well as abnormal lymphocytic and myeloid cell infiltrates in the small and/or large intestine. Thus, in many respects IBD in dogs resembles IBD in humans. However, the factors that trigger intestinal inflammation in dogs with IBD are not well understood and have been variously attributed to immune responses against dietary antigens or intestinal antigens. Previous studies in humans with IBD have documented increased production of IgG and IgA antibodies specific to intestinal bacteria, and this abnormal immune response has been linked to disease pathogenesis. Therefore, we investigated the humoral immune response against gut bacteria in dogs with IBD, using flow cytometry to quantitate IgG and IgA binding. Studies were also done to investigate the source of these antibodies (locally produced versus systemic production) and whether greater antibody binding to bacteria is associated with increased inflammatory responses. We found that dogs with IBD had significantly higher percentages and overall amounts of IgG bound to their intestinal bacteria compared to healthy dogs. Similarly, significantly higher percentages of bacteria were IgA+ bacteria were also found in dogs with IBD. Serum antibody recognition of gut bacteria was not different between healthy dogs and dogs with IBD, suggesting that anti-bacterial antibodies were primarily produced locally in the gut rather than systemically. Importantly, bacteria in the Actinobacteria phylum and in particular the genus Collinsella had significantly greater levels of antibody binding in dogs with IBD. Based on these findings, we concluded that antibody binding to commensal gut bacteria was significantly increased in dogs with IBD, that particular phyla were preferential targets for gut antibodies, and that anti-bacterial antibody responses may play an important role in regulating gut inflammation.

The long and winding road to IgA deficiency: causes and consequences

INTRODUCTION

The most common humoral immunodeficiency is IgA deficiency. One of the first papers addressing the cellular and molecular mechanisms underlying IgA deficiency indicated that immature IgA-positive B-lymphocytes are present in these patients. This suggests that the genetic background for IgA is still intact and that class switching can take place. At this moment, it cannot be ruled out that genetic as well as environmental factors are involved. Areas covered: A clinical presentation, the biological functions of IgA, and the management of IgA deficiency are reviewed. In some IgA deficient patients, a relationship with a loss-of-function mutation in the TACI (transmembrane activator and calcium-modulating cyclophilin ligand interaction) gene has been found. Many other genes also have been associated. Gut microbiota are an important environmental trigger for IgA synthesis. Expert commentary: Expression of IgA deficiency is due to both genetic and environmental factors and a role for gut microbiota cannot be excluded.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs

Background

Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2.

Results

Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10⁻⁷) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10⁻⁵), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (p_difference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region.

Conclusions

Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0404-3) contains supplementary material, which is available to authorized users.