Evaluating Prophylactic Effect of Bovine Colostrum on Intestinal Barrier Function in Zonulin Transgenic Mice: A Transcriptomic Study.

The intestinal barrier comprises a single layer of epithelial cells tightly joined to form a physical barrier. Disruption or compromise of the intestinal barrier can lead to the inadvertent activation of immune cells, potentially causing an increased risk of chronic inflammation in various tissues. Recent research has suggested that specific dietary components may influence the function of the intestinal barrier, potentially offering a means to prevent or mitigate inflammatory disorders. However, the precise mechanism underlying these effects remains unclear. Bovine colostrum (BC), the first milk from cows after calving, is a natural source of nutrients with immunomodulatory, anti-inflammatory, and gut-barrier fortifying properties. This novel study sought to investigate the transcriptome in BC-treated Zonulin transgenic mice (Ztm), characterized by dysbiotic microbiota, intestinal hyperpermeability, and mild hyperactivity, applying RNA sequencing. Seventy-five tissue samples from the duodenum, colon, and brain of Ztm and wild-type (WT) mice were dissected, processed, and RNA sequenced. The expression profiles were analyzed and integrated to identify differentially expressed genes (DEGs) and differentially expressed transcripts (DETs). These were then further examined using bioinformatics tools. RNA-seq analysis identified 1298 DEGs and 20,952 DETs in the paired (Ztm treatment vs. Ztm control) and reference (WT controls) groups. Of these, 733 DEGs and 10,476 DETs were upregulated, while 565 DEGs and 6097 DETs were downregulated. BC-treated Ztm female mice showed significant upregulation of cingulin (Cgn) and claudin 12 (Cldn12) duodenum and protein interactions, as well as molecular pathways and interactions pertaining to tight junctions, while BC-treated Ztm males displayed an upregulation of transcripts like occludin (Ocln) and Rho/Rac guanine nucleotide exchange factor 2 (Arhgf2) and cellular structures and interfaces, protein-protein interactions, and organization and response mechanisms. This comprehensive analysis reveals the influence of BC treatment on tight junctions (TJs) and Nuclear factor kappa-light-chain-enhancer of activated B cells (NF-kB) signaling pathway gene expressions. The present study is the first to analyze intestinal and brain samples from BC-treated Ztm mice applying high-throughput RNA sequencing. This study revealed molecular interaction in intestinal barrier function and identified hub genes and their functional pathways and biological processes in response to BC treatment in Ztm mice. Further research is needed to validate these findings and explore their implications for dietary interventions aimed at improving intestinal barrier integrity and function. The MGH Institutional Animal Care and Use Committee authorized the animal study (2013N000013).

Aroylated phenylenediamine HO53 modulates innate immunity, histone acetylation and metabolism.

In the current context of antibiotic resistance, the need to find alternative treatment strategies is urgent. Our research aimed to use synthetized aroylated phenylenediamines (APDs) to induce the expression of cathelicidin antimicrobial peptide gene (CAMP) to minimize the necessity of antibiotic use during infection. One of these compounds, HO53, showed promising results in inducing CAMP expression in bronchial epithelium cells (BCi-NS1.1 hereafter BCi). Thus, to decipher the cellular effects of HO53 on BCi cells, we performed RNA sequencing (RNAseq) analysis after 4, 8 and 24 h treatment of HO53. The number of differentially expressed transcripts pointed out an epigenetic modulation. Yet, the chemical structure and in silico modeling indicated HO53 as a histone deacetylase (HDAC) inhibitor. When exposed to a histone acetyl transferase (HAT) inhibitor, BCi cells showed a decreased expression of CAMP. Inversely, when treated with a specific HDAC3 inhibitor (RGFP996), BCi cells showed an increased expression of CAMP, indicating acetylation status in cells as determinant for the induction of the expression of the gene CAMP expression. Interestingly, a combination treatment with both HO53 and HDAC3 inhibitor RGFP966 leads to a further increase of CAMP expression. Moreover, HDAC3 inhibition by RGFP966 leads to increased expression of STAT3 and HIF1A, both previously demonstrated to be involved in pathways regulating CAMP expression. Importantly, HIF1α is considered as a master regulator in metabolism. A significant number of genes of metabolic enzymes were detected in our RNAseq data with enhanced expression conveying a shift toward enhanced glycolysis. Overall, we are demonstrating that HO53 might have a translational value against infections in the future through a mechanism leading to innate immunity strengthening involving HDAC inhibition and shifting the cells towards an immunometabolism, which further favors innate immunity activation.

Ventilator-induced lung injury results in oxidative stress response and mitochondrial swelling in a mouse model.

BACKGROUND: Mechanical ventilation is a life-saving therapy for critically ill patients, providing rest to the respiratory muscles and facilitating gas exchange in the lungs. Ventilator-induced lung injury (VILI) is an unfortunate side effect of mechanical ventilation that may lead to serious consequences for the patient and increase mortality. The four main injury mechanisms associated with VILI are: baro/volutrauma caused by overstretching the lung tissues; atelectrauma, caused by repeated opening and closing of the alveoli resulting in shear stress; oxygen toxicity due to use of high ratio of oxygen in inspired air, causing formation of free radicals; and biotrauma, the resulting biological response to tissue injury, that leads to a cascade of events due to excessive inflammatory reactions and may cause multi-organ failure. An often-overlooked part of the inflammatory reaction is oxidative stress. In this research, a mouse model of VILI was set up with three tidal volume settings (10, 20 and 30 mL/kg) at atmospheric oxygen level. Airway pressures and heart rate were monitored and bronchoalveolar lavage fluid (BALF) and lung tissue samples were taken. RESULTS: We show a correlation between increased inflammation and barrier failure, and higher tidal volumes, evidenced by increased IL-6 expression, high concentration of proteins in BALF along with changes in expression of adhesion molecules. Furthermore, swelling of mitochondria in alveolar type II cells was seen indicating their dysfunction and senescence-like state. RNA sequencing data present clear increases in inflammation, mitochondrial biogenesis and oxidative stress as tidal volume is increased, supported by degradation of Keap1, a redox-regulated substrate adaptor protein. CONCLUSIONS: Oxidative stress seems to be a more prominent mechanism of VILI than previously considered, indicating that possible treatment methods against VILI might be identified by impeding oxidative pathways.

The Novel Inducer of Innate Immunity HO53 Stimulates Autophagy in Human Airway Epithelial Cells.

Aroylated phenylenediamines (APDs) are novel modulators of innate immunity with respect to enhancing the expression of antimicrobial peptides and maintaining epithelial barrier integrity. Here, we present a new study on induction of autophagy in human lung epithelial cells by the APD HO53. Interestingly, HO53 affected autophagy in a dose-dependent manner, demonstrated by increased microtubule-associated proteins 1A/1B light-chain 3B (LC3B) processing in mature polarized bronchial epithelial cells. The quantification of LC3B puncta showed increased autophagy flux and formation of autophagosomes visualized by transmission electron microscopy. The phenotypic changes indicated that autophagy induction was associated with activation of 5' adenosine monophosphate-activated protein kinase (AMPK), nuclear translocation of transcription factor EB (TFEB), and changes in expression of autophagy-related genes. The kinetics of the explored signaling pathways indicated on activation of AMPK followed by the nuclear translocation of TFEB. Moreover, our data suggest that HO53 modulates epigenetic changes related to induction of autophagy manifested by transcriptional regulation of histone-modifying enzymes. These changes were reflected by decreased ubiquitination of histone 2B at the lysine 120 residue that is associated with autophagy induction. Taken together, HO53 modulates autophagy, a part of the host defense system, through a complex mechanism involving several pathways and epigenetic events.

Metabolic and Transcriptional Changes across Osteogenic Differentiation of Mesenchymal Stromal Cells.

Mesenchymal stromal cells (MSCs) are multipotent post-natal stem cells with applications in tissue engineering and regenerative medicine. MSCs can differentiate into osteoblasts, chondrocytes, or adipocytes, with functional differences in cells during osteogenesis accompanied by metabolic changes. The temporal dynamics of these metabolic shifts have not yet been fully characterized and are suspected to be important for therapeutic applications such as osteogenesis optimization. Here, our goal was to characterize the metabolic shifts that occur during osteogenesis. We profiled five key extracellular metabolites longitudinally (glucose, lactate, glutamine, glutamate, and ammonia) from MSCs from four donors to classify osteogenic differentiation into three metabolic stages, defined by changes in the uptake and secretion rates of the metabolites in cell culture media. We used a combination of untargeted metabolomic analysis, targeted analysis of 13C-glucose labelled intracellular data, and RNA-sequencing data to reconstruct a gene regulatory network and further characterize cellular metabolism. The metabolic stages identified in this proof-of-concept study provide a framework for more detailed investigations aimed at identifying biomarkers of osteogenic differentiation and small molecule interventions to optimize MSC differentiation for clinical applications.

Azithromycin ameliorates sulfur dioxide-induced airway epithelial damage and inflammatory responses.

BACKGROUND: The airway epithelium (AE) forms the first line of defence against harmful particles and pathogens. Barrier failure of the airway epithelium contributes to exacerbations of a range of lung diseases that are commonly treated with Azithromycin (AZM). In addition to its anti-bacterial function, AZM has immunomodulatory effects which are proposed to contribute to its clinical effectiveness. In vitro studies have shown the AE barrier-enhancing effects of AZM. The aim of this study was to analyze whether AE damage caused by inhalation of sulfur dioxide (SO2) in a murine model could be reduced by pre-treatment with AZM. METHODS: The leakiness of the AE barrier was evaluated after SO2 exposure by measuring levels of human serum albumin (HSA) in bronchoalveolar lavage fluid (BALF). Protein composition in BALF was also assessed and lung tissues were evaluated across treatments using histology and gene expression analysis. RESULTS: AZM pre-treatment (2 mg/kg p.o. 5 times/week for 2 weeks) resulted in reduced glutathione-S-transferases in BALF of SO2 injured mice compared to control (without AZM treatment). AZM treated mice had increased intracellular vacuolization including lamellar bodies and a reduction in epithelial shedding after injury in addition to a dampened SO2-induced inflammatory response. CONCLUSIONS: Using a mouse model of AE barrier dysfunction we provide evidence for the protective effects of AZM in vivo, possibly through stabilizing the intracellular microenvironment and reducing inflammatory responses. Our data provide insight into the mechanisms contributing to the efficacy of AZM in the treatment of airway diseases.

Expression of ncRNAs on the DLK1-DIO3 Locus Is Associated With Basal and Mesenchymal Phenotype in Breast Epithelial Progenitor Cells.

Epithelial-to-mesenchymal transition (EMT) and its reversed process mesenchymal-to-epithelial transition (MET) play a critical role in epithelial plasticity during development and cancer progression. Among important regulators of these cellular processes are non-coding RNAs (ncRNAs). The imprinted DLK1-DIO3 locus, containing numerous maternally expressed ncRNAs including the lncRNA maternally expressed gene 3 (MEG3) and a cluster of over 50 miRNAs, has been shown to be a modulator of stemness in embryonic stem cells and in cancer progression, potentially through the tumor suppressor role of MEG3. In this study we analyzed the expression pattern and functional role of ncRNAs from the DLK1-DIO3 locus in epithelial plasticity of the breast. We studied their expression in various cell types of breast tissue and revisit the role of the locus in EMT/MET using a breast epithelial progenitor cell line (D492) and its isogenic mesenchymal derivative (D492M). Marked upregulation of ncRNAs from the DLK1-DIO3 locus was seen after EMT induction in two cell line models of EMT. In addition, the expression of MEG3 and the maternally expressed ncRNAs was higher in stromal cells compared to epithelial cell types in primary breast tissue. We also show that expression of MEG3 is concomitant with the expression of the ncRNAs from the DLK1-DIO3 locus and its expression is therefore likely indicative of activation of all ncRNAs at the locus. MEG3 expression is correlated with stromal markers in normal tissue and breast cancer tissue and negatively correlated with the survival of breast cancer patients in two different cohorts. Overexpression of MEG3 using CRISPR activation in a breast epithelial cell line induced partial EMT and enriched for a basal-like phenotype. Conversely, knock down of MEG3 using CRISPR inhibition in a mesenchymal cell line reduced the mesenchymal and basal-like phenotype of the cell line. In summary our study shows that maternally expressed ncRNAs are markers of EMT and suggests that MEG3 is a novel regulator of EMT/MET in breast tissue. Nevertheless, further studies are needed to fully dissect the molecular pathways influenced by non-coding RNAs at the DLK1-DIO3 locus in breast tissue.

Azithromycin has lung barrier protective effects in a cell model mimicking ventilator-induced lung injury.

Azithromycin (AZM) is a broad-spectrum antibiotic widely used to treat infections. AZM also has been shown to have anti-inflammatory and immunomodulatory functions unrelated to its antibacterial activity that contribute to the effectiveness of this drug in chronic respiratory diseases. The mechanisms behind these beneficial effects are not yet fully elucidated. We have previously shown that AZM enhances barrier integrity of bronchial epithelial cells and directs them towards epidermal differentiation. In this study, we analyzed the effect of AZM pre-treatment of human bronchial and alveolar derived cell lines on mechanical stress in a cyclical pressure air-liquid interface device (CPAD) that models the disruption of the epithelial barrier with increased inflammatory response in lung tissue, which is associated with ventilator-induced lung injury (VILI). Immunostaining and electron microscopy showed that barrier integrity of the epithelium was compromised by cyclically stressing the cells but maintained when cells had been pre-treated with AZM. Lamellar body formation was revealed in AZM pre-treated cells, possibly further supporting the barrier-enhancing effects. RNA sequencing showed that the inflammatory response was attenuated by AZM treatment before cyclical stress. YKL-40, an emerging inflammatory marker, increased both due to cyclical stress and upon AZM treatment. These data confirm the usefulness of the CPAD to model ventilator-induced lung injury and suggest that AZM has barrier protective and immunomodulatory effects, attenuating the inflammatory response during mechanical stress, and might therefore be lung protective during mechanical ventilation. The model could be used to assess further drug candidates that influence barrier integrity and modulate inflammatory response.

Common and Rare Sequence Variants Influencing Tumor Biomarkers in Blood.

BACKGROUND: Alpha-fetoprotein (AFP), cancer antigens 15.3, 19.9, and 125, carcinoembryonic antigen, and alkaline phosphatase (ALP) are widely measured in attempts to detect cancer and to monitor treatment response. However, due to lack of sensitivity and specificity, their utility is debated. The serum levels of these markers are affected by a number of nonmalignant factors, including genotype. Thus, it may be possible to improve both sensitivity and specificity by adjusting test results for genetic effects. METHODS: We performed genome-wide association studies of serum levels of AFP (N = 22,686), carcinoembryonic antigen (N = 22,309), cancer antigens 15.3 (N = 7,107), 19.9 (N = 9,945), and 125 (N = 9,824), and ALP (N = 162,774). We also examined the correlations between levels of these biomarkers and the presence of cancer, using data from a nationwide cancer registry. RESULTS: We report a total of 84 associations of 79 sequence variants with levels of the six biomarkers, explaining between 2.3% and 42.3% of the phenotypic variance. Among the 79 variants, 22 are cis (in- or near the gene encoding the biomarker), 18 have minor allele frequency less than 1%, 31 are coding variants, and 7 are associated with gene expression in whole blood. We also find multiple conditions associated with higher biomarker levels. CONCLUSIONS: Our results provide insights into the genetic contribution to diversity in concentration of tumor biomarkers in blood. IMPACT: Genetic correction of biomarker values could improve prediction algorithms and decision-making based on these biomarkers.

Azithromycin induces epidermal differentiation and multivesicular bodies in airway epithelia.

BACKGROUND: Azithromycin (Azm) is a macrolide recognized for its disease-modifying effects and reduction in exacerbation of chronic airway diseases. It is not clear whether the beneficial effects of Azm are due to its anti-microbial activity or other pharmacological actions. We have shown that Azm affects the integrity of the bronchial epithelial barrier measured by increased transepithelial electrical resistance. To better understand these effects of Azm on bronchial epithelia we have investigated global changes in gene expression. METHODS: VA10 bronchial epithelial cells were treated with Azm and cultivated in air-liquid interface conditions for up to 22 days. RNA was isolated at days 4, 10 and 22 and analyzed using high-throughput RNA sequencing. qPCR and immunostaining were used to confirm key findings from bioinformatic analyses. Detailed assessment of cellular changes was done using microscopy, followed by characterization of the lipidomic profiles of the multivesicular bodies present. RESULTS: Bioinformatic analysis revealed that after 10 days of treatment genes encoding effectors of sterol and cholesterol metabolism were prominent. Interestingly, expression of genes associated with epidermal barrier differentiation, KRT1, CRNN, SPINK5 and DSG1, increased significantly at day 22. Together with immunostaining, these results suggest an epidermal differentiation pattern. We also found that Azm induced the formation of multivesicular and lamellar bodies in two different airway epithelial cell lines. Lipidomic analysis revealed that Azm was entrapped in multivesicular bodies linked to different types of lipids, most notably palmitate and stearate. Furthermore, targeted analysis of lipid species showed accumulation of phosphatidylcholines, as well as ceramide derivatives. CONCLUSIONS: Taken together, we demonstrate how Azm might confer its barrier enhancing effects, via activation of epidermal characteristics and changes to intracellular lipid dynamics. These effects of Azm could explain the unexpected clinical benefit observed during Azm-treatment of patients with various lung diseases affecting barrier function.

A homozygous loss-of-function mutation leading to CYBC1 deficiency causes chronic granulomatous disease.

Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase's main subunit (gp91phox) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chip-genotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10-8; OR = 67.6), as well as reduced height (P = 3.3 × 10-4; -8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction.

Fourteen sequence variants that associate with multiple sclerosis discovered by meta-analysis informed by genetic correlations.

A meta-analysis of publicly available summary statistics on multiple sclerosis combined with three Nordic multiple sclerosis cohorts (21,079 cases, 371,198 controls) revealed seven sequence variants associating with multiple sclerosis, not reported previously. Using polygenic risk scores based on public summary statistics of variants outside the major histocompatibility complex region we quantified genetic overlap between common autoimmune diseases in Icelanders and identified disease clusters characterized by autoantibody presence/absence. As multiple sclerosis-polygenic risk scores captures the risk of primary biliary cirrhosis and vice versa (P = 1.6 × 10-7, 4.3 × 10-9) we used primary biliary cirrhosis as a proxy-phenotype for multiple sclerosis, the idea being that variants conferring risk of primary biliary cirrhosis have a prior probability of conferring risk of multiple sclerosis. We tested 255 variants forming the primary biliary cirrhosis-polygenic risk score and found seven multiple sclerosis-associating variants not correlated with any previously established multiple sclerosis variants. Most of the variants discovered are close to or within immune-related genes. One is a low-frequency missense variant in TYK2, another is a missense variant in MTHFR that reduces the function of the encoded enzyme affecting methionine metabolism, reported to be dysregulated in multiple sclerosis brain.

Effect of sequence variants on variance in glucose levels predicts type 2 diabetes risk and accounts for heritability.

Sequence variants that affect mean fasting glucose levels do not necessarily affect risk for type 2 diabetes (T2D). We assessed the effects of 36 reported glucose-associated sequence variants on between- and within-subject variance in fasting glucose levels in 69,142 Icelanders. The variant in TCF7L2 that increases fasting glucose levels increases between-subject variance (5.7% per allele, P = 4.2 × 10-10), whereas variants in GCK and G6PC2 that increase fasting glucose levels decrease between-subject variance (7.5% per allele, P = 4.9 × 10-11 and 7.3% per allele, P = 7.5 × 10-18, respectively). Variants that increase mean and between-subject variance in fasting glucose levels tend to increase T2D risk, whereas those that increase the mean but reduce variance do not (r2 = 0.61). The variants that increase between-subject variance increase fasting glucose heritability estimates. Intuitively, our results show that increasing the mean and variance of glucose levels is more likely to cause pathologically high glucose levels than increase in the mean offset by a decrease in variance.

Sequence variant at 4q25 near PITX2 associates with appendicitis.

Appendicitis is one of the most common conditions requiring acute surgery and can pose a threat to the lives of affected individuals. We performed a genome-wide association study of appendicitis in 7,276 Icelandic and 1,139 Dutch cases and large groups of controls. In a combined analysis of the Icelandic and Dutch data, we detected a single signal represented by an intergenic variant rs2129979 [G] close to the gene PITX2 associating with increased risk of appendicitis (OR = 1.15, P = 1.8 × 10-11). We only observe the association in patients diagnosed in adulthood. The marker is close to, but distinct from, a set of markers reported to associate with atrial fibrillation, which have been linked to PITX2. PITX2 has been implicated in determination of right-left symmetry during development. Anomalies in organ arrangement have been linked to increased prevalence of gastrointestinal and intra-abdominal complications, which may explain the effect of rs2129979 on appendicitis risk.

Sequence variants in ARHGAP15, COLQ and FAM155A associate with diverticular disease and diverticulitis.

Diverticular disease is characterized by pouches (that is, diverticulae) due to weakness in the bowel wall, which can become infected and inflamed causing diverticulitis, with potentially severe complications. Here, we test 32.4 million sequence variants identified through whole-genome sequencing (WGS) of 15,220 Icelanders for association with diverticular disease (5,426 cases) and its more severe form diverticulitis (2,764 cases). Subsequently, 16 sequence variants are followed up in a diverticular disease sample from Denmark (5,970 cases, 3,020 controls). In the combined Icelandic and Danish data sets we observe significant association of intronic variants in ARHGAP15 (Rho GTPase-activating protein 15; rs4662344-T: P=1.9 × 10-18, odds ratio (OR)=1.23) and COLQ (collagen-like tail subunit of asymmetric acetylcholinesterase; rs7609897-T: P=1.5 × 10-10, OR=0.87) with diverticular disease and in FAM155A (family with sequence similarity 155A; rs67153654-A: P=3.0 × 10-11, OR=0.82) with diverticulitis. These are the first loci shown to associate with diverticular disease in a genome-wide study.

Rare mutations associating with serum creatinine and chronic kidney disease.

Chronic kidney disease (CKD) is a complex disorder with a strong genetic component. A number of common sequence variants have been found to associate with serum creatinine (SCr), estimated glomerular filtration rate (eGFR) and/or CKD. We imputed 24 million single-nucleotide polymorphisms and insertions/deletions identified by whole-genome sequencing of 2230 Icelanders into 81 656 chip-typed individuals and 112 630 relatives of genotyped individuals over the age of 18 with SCr measurements. The large set of sequenced individuals allowed accurate imputation of variants to a minor allele frequency (MAF) of 0.1%. We tested the imputed variants for association with SCr. In addition to replicating established loci, we discovered missense and loss-of-function variants associating with SCr in three solute carriers (SLC6A19, SLC25A45 and SLC47A1) and two E3 ubiquitin ligases (RNF186 and RNF128). All the variants are within coding sequences and all but one are rare (MAF <2%) with SCr effects between 0.085 and 0.129 standard deviations. These rare variants have a larger effect on SCr than previously reported common variants, explaining 0.5% of the variability of SCr in Icelanders in addition to the 1% already accounted for. We tested the five variants associating with SCr for association with CKD in an Icelandic sample of 15 594 cases and 291 428 controls. Three of the variants also associated with CKD. These variants may either affect kidney function or creatinine synthesis and excretion. Of note were four mutations in SLC6A19 that associate with reduced SCr, three of which have been shown to cause Hartnup disease.

Candidate genes and functional noncoding variants identified in a canine model of obsessive-compulsive disorder.

BACKGROUND: Obsessive-compulsive disorder (OCD), a severe mental disease manifested in time-consuming repetition of behaviors, affects 1 to 3% of the human population. While highly heritable, complex genetics has hampered attempts to elucidate OCD etiology. Dogs suffer from naturally occurring compulsive disorders that closely model human OCD, manifested as an excessive repetition of normal canine behaviors that only partially responds to drug therapy. The limited diversity within dog breeds makes identifying underlying genetic factors easier. RESULTS: We use genome-wide association of 87 Doberman Pinscher cases and 63 controls to identify genomic loci associated with OCD and sequence these regions in 8 affected dogs from high-risk breeds and 8 breed-matched controls. We find 119 variants in evolutionarily conserved sites that are specific to dogs with OCD. These case-only variants are significantly more common in high OCD risk breeds compared to breeds with no known psychiatric problems. Four genes, all with synaptic function, have the most case-only variation: neuronal cadherin (CDH2), catenin alpha2 (CTNNA2), ataxin-1 (ATXN1), and plasma glutamate carboxypeptidase (PGCP). In the 2 Mb gene desert between the cadherin genes CDH2 and DSC3, we find two different variants found only in dogs with OCD that disrupt the same highly conserved regulatory element. These variants cause significant changes in gene expression in a human neuroblastoma cell line, likely due to disrupted transcription factor binding. CONCLUSIONS: The limited genetic diversity of dog breeds facilitates identification of genes, functional variants and regulatory pathways underlying complex psychiatric disorders that are mechanistically similar in dogs and humans.

Canine hereditary ataxia in old english sheepdogs and gordon setters is associated with a defect in the autophagy gene encoding RAB24.

Old English Sheepdogs and Gordon Setters suffer from a juvenile onset, autosomal recessive form of canine hereditary ataxia primarily affecting the Purkinje neuron of the cerebellar cortex. The clinical and histological characteristics are analogous to hereditary ataxias in humans. Linkage and genome-wide association studies on a cohort of related Old English Sheepdogs identified a region on CFA4 strongly associated with the disease phenotype. Targeted sequence capture and next generation sequencing of the region identified an A to C single nucleotide polymorphism (SNP) located at position 113 in exon 1 of an autophagy gene, RAB24, that segregated with the phenotype. Genotyping of six additional breeds of dogs affected with hereditary ataxia identified the same polymorphism in affected Gordon Setters that segregated perfectly with phenotype. The other breeds tested did not have the polymorphism. Genome-wide SNP genotyping of Gordon Setters identified a 1.9 MB region with an identical haplotype to affected Old English Sheepdogs. Histopathology, immunohistochemistry and ultrastructural evaluation of the brains of affected dogs from both breeds identified dramatic Purkinje neuron loss with axonal spheroids, accumulation of autophagosomes, ubiquitin positive inclusions and a diffuse increase in cytoplasmic neuronal ubiquitin staining. These findings recapitulate the changes reported in mice with induced neuron-specific autophagy defects. Taken together, our results suggest that a defect in RAB24, a gene associated with autophagy, is highly associated with and may contribute to canine hereditary ataxia in Old English Sheepdogs and Gordon Setters. This finding suggests that detailed investigation of autophagy pathways should be undertaken in human hereditary ataxia.

Genome-wide analyses implicate 33 loci in heritable dog osteosarcoma, including regulatory variants near CDKN2A/B.

BACKGROUND: Canine osteosarcoma is clinically nearly identical to the human disease, but is common and highly heritable, making genetic dissection feasible. RESULTS: Through genome-wide association analyses in three breeds (greyhounds, Rottweilers, and Irish wolfhounds), we identify 33 inherited risk loci explaining 55% to 85% of phenotype variance in each breed. The greyhound locus exhibiting the strongest association, located 150 kilobases upstream of the genes CDKN2A/B, is also the most rearranged locus in canine osteosarcoma tumors. The top germline candidate variant is found at a >90% frequency in Rottweilers and Irish wolfhounds, and alters an evolutionarily constrained element that we show has strong enhancer activity in human osteosarcoma cells. In all three breeds, osteosarcoma-associated loci and regions of reduced heterozygosity are enriched for genes in pathways connected to bone differentiation and growth. Several pathways, including one of genes regulated by miR124, are also enriched for somatic copy-number changes in tumors. CONCLUSIONS: Mapping a complex cancer in multiple dog breeds reveals a polygenic spectrum of germline risk factors pointing to specific pathways as drivers of disease.

Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing.

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (∼1.5-5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing.