THE ROLE OF AN EXPLICIT CAUSAL FRAMEWORK IN AFFECTED SIB PAIR DESIGNS WITH COVARIATES

The affected sib/relative pair (ASP/ARP) design is often used with covariates to find genes that can cause a disease in pathways other than through those covariates. However, such "covariates" can themselves have genetic determinants, and the validity of existing methods has so far only been argued under implicit assumptions. We propose an explicit causal formulation of the problem using potential outcomes and principal stratification. The general role of this formulation is to identify and separate the meaning of the different assumptions that can provide valid causal inference in linkage analysis. This separation helps to (a) develop better methods under explicit assumptions, and (b) show the different ways in which these assumptions can fail, which is necessary for developing further specific designs to test these assumptions and confirm or improve the inference. Using this formulation in the specific problem above, we show that, when the "covariate" (e.g., addiction to smoking) also has genetic determinants, then existing methods, including those previously thought as valid, can declare linkage between the disease and marker loci even when no such linkage exists. We also introduce design strategies to address the problem.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

A region on equine chromosome 13 is linked to recurrent airway obstruction in horses

REASONS FOR STUDY: Equine recurrent airway obstruction (RAO) is probably dependent on a complex interaction of genetic and environmental factors and shares many characteristic features with human asthma. Interleukin 4 receptor a chain (IL4RA) is a candidate gene because of its role in the development of human asthma, confirmation of this association is therefore required. METHODS: The equine BAC clone containing the IL4RA gene was localised to ECA13q13 by the FISH method. Microsatellite markers in this region were investigated for possible association and linkage with RAO in 2 large Warmblood halfsib families. Based on a history of clinical signs (coughing, nasal discharge, abnormal breathing and poor performance), horses were classified in a horse owner assessed respiratory signs index (HOARSI 1-4: from healthy, mild, moderate to severe signs). Four microsatellite markers (AHT133, LEX041, VHL47, ASB037) were analysed in the offspring of Sire 1 (48 unaffected HOARSI 1 vs. 59 affected HOARSI 2-4) and Sire 2 (35 HOARSI 1 vs. 50 HOARSI 2-4), age 07 years. RESULTS: For both sires haplotypes could be established in the order AHT133-LEXO47-VHL47-ASB37. The distances in this order were estimated to be 2.9, 0.9 and 2.3 centiMorgans, respectively. Haplotype association with mild to severe clinical signs of chronic lower airway disease (HOARSI 2-4) was significant in the offspring of Sire 1 (P = 0.026) but not significant for the offspring of Sire 2 (P = 0.32). Linkage analysis showed the ECA13q13 region containing IL4RA to be linked to equine chronic lower airway disease in one family (P<0.01), but not in the second family. CONCLUSIONS: This supports a genetic background for equine RAO and indicates that IL4RA is a candidate gene with possible locus heterogeneity for this disease. POTENTIAL RELEVANCE: Identification of major genes for RAO may provide a basis for breeding and individual prevention for this important disease.

Novel TULP1 mutation causing leber congenital amaurosis or early onset retinal degeneration

PURPOSE: To report a large, consanguineous Algerian family affected with Leber congenital amaurosis (LCA) or early-onset retinal degeneration (EORD). METHODS: All accessible family members underwent a complete ophthalmic examination, and blood was obtained for DNA extraction. Homozygosity mapping was performed with markers flanking 12 loci associated with LCA. The 15 exons of TULP1 were sequenced. RESULTS: Seven of 30 examined family members were affected, including five with EORD and two with LCA. All patients had nystagmus, hemeralopia, mild myopia, and low visual acuity without photophobia. Fundus features were variable among EORD patients: typical spicular retinitis pigmentosa or clumped pigmented retinopathy with age-dependent macular involvement. A salt-and-pepper retinopathy with midperipheral retinal pigment epithelium (RPE) atrophy was present in the older patients with LCA, whereas the retina appeared virtually normal in the younger ones. Both scotopic and photopic electroretinograms were nondetectable. Fundus imaging revealed a perifoveal ring of increased fundus autofluorescence (FAF) in the proband, and optical coherence tomography disclosed a thinned retina, mainly due to photoreceptor loss. Linkage analysis identified a region of homozygosity on chromosome 6, region p21.3, and mutation screening revealed a novel 6-base in-frame duplication, in the TULP1 gene. CONCLUSIONS: Mutation in the TULP1 gene is a rare cause of LCA/EORD, with only 14 mutations reported so far. The observed intrafamilial phenotypic variability could be attributed to disease progression or possibly modifier alleles. This study provides the first description of FAF and quantitative reflectivity profiles in TULP1-related retinopathy.

Study of a Swiss dopa-responsive dystonia family with a deletion in GCH1: redefining DYT14 as DYT5

OBJECTIVE: To report the study of a multigenerational Swiss family with dopa-responsive dystonia (DRD). METHODS: Clinical investigation was made of available family members, including historical and chart reviews. Subject examinations were video recorded. Genetic analysis included a genome-wide linkage study with microsatellite markers (STR), GTP cyclohydrolase I (GCH1) gene sequencing, and dosage analysis. RESULTS: We evaluated 32 individuals, of whom 6 were clinically diagnosed with DRD, with childhood-onset progressive foot dystonia, later generalizing, followed by parkinsonism in the two older patients. The response to levodopa was very good. Two additional patients had late onset dopa-responsive parkinsonism. Three other subjects had DRD symptoms on historical grounds. We found suggestive linkage to the previously reported DYT14 locus, which excluded GCH1. However, further study with more stringent criteria for disease status attribution showed linkage to a larger region, which included GCH1. No mutation was found in GCH1 by gene sequencing but dosage methods identified a novel heterozygous deletion of exons 3 to 6 of GCH1. The mutation was found in seven subjects. One of the patients with dystonia represented a phenocopy. CONCLUSIONS: This study rules out the previously reported DYT14 locus as a cause of disease, as a novel multiexonic deletion was identified in GCH1. This work highlights the necessity of an accurate clinical diagnosis in linkage studies as well as the need for appropriate allele frequencies, penetrance, and phenocopy estimates. Comprehensive sequencing and dosage analysis of known genes is recommended prior to genome-wide linkage analysis.

Extracellular matrix of the charophycean green algae

A comprehensive knowledge of cell wallstructure and function throughout the plant kingdom is essential to understanding cell wall evolution. The fundamental understanding of the charophycean green algal cell wall is broadening. The similarities and differences that exist between land plant and algal cell walls provide opportunities to understand plant evolution. A variety of polymers previously associated with higher plants were discovered in the charophycean green algae (CGA), including homogalacturonans, cross-linking glycans, arabinogalactan protein, β-glucans, and cellulose. The cellulose content of CGA cell walls ranged from 6% to 43%, with the higher valuescomparable to that found in the primary cell wall of land plants (20-30%). (1,3)β-glucans were found in the unicellular Chlorokybus atmophyticus, Penium margaritaceum, and Cosmarium turpini, the unbranched filamentous Klebsormidium flaccidum, and the multicellular Chara corallina. The discovery of homogalacturonan in Penium margaritaceum representsthe first confirmation of land plant-type pectinsin desmids and the second rigorous characterization of a pectin polymer from the charophycean algae. Homogalacturonan was also indicated from the basal species Chlorokybus atmophyticus and Klebsormidium flaccidum. There is evidence of branched pectins in Cosmarium turpini and linkage analysis suggests the presence of type I rhamnogalacturonan (RGI). Cross-linking β-glucans are associated with cellulose microfibrils during land plant cell growth, and were found in the cell wall of CGA. The evidence of mixed-linkage glucan (MLG) in the 11 charophytesis both suprising and significant given that MLG was once thought to be specific to some grasses. The organization and structure of Cosmarium turpini and Chara corallina MLG was found to be similar to that of Equisetumspp., whereas the basal species of the CGA, Chlorokybus atmophyticus and Klebsormidium flaccidum, have unique organization of alternating of 3- and 4-linkages. The significance of this result on the evolution of the MLG synthetic pathway has yet to be determined. The extracellular matrix (ECM) of Chlorokybus atmophyticus, Klebsormidium flaccidum, and Spirogyra spp. exhibits significant biochemical diversity, ranging from distinct “land plant” polymers to polysaccharides unique to these algae. The neutral sugar composition of Chlorokybus atmophyticus hot water extract and Spirogyra extracellular polymeric substance (EPS), combined with antibody labeling results, revealed the distinct possibility of an arabinogalactan protein in these organisms. Polysaccharide analysis of Zygnematales (desmid) EPS, indicated a probable range of different EPS backbones and substitution patterns upon the core portions of the molecules. Desmid EPS is predominately composed of a complex matrix of branched, uronic acid containing polysaccharides with ester sulfate substitutions and, as such, has an almost infinite capacity for various hydrogen bonding, hydrophobic interaction and ionic cross-bridging motifs, which characterize their unique function in biofilms. My observations support the hypothesis that members of the CGA represent the phylogenetic line that gave rise to vascular plants and that the primary cell wall of vascular plants many have evolved directly from structures typical of the cell wall of filamentous green algae found in the charophycean green algae.

Identification of novel mutations in X-linked retinitis pigmentosa families and implications for diagnostic testing

PURPOSE: The goal of this study was to identify mutations in X-chromosomal genes associated with retinitis pigmentosa (RP) in patients from Germany, The Netherlands, Denmark, and Switzerland. METHODS: In addition to all coding exons of RP2, exons 1 through 15, 9a, ORF15, 15a and 15b of RPGR were screened for mutations. PCR products were amplified from genomic DNA extracted from blood samples and analyzed by direct sequencing. In one family with apparently dominant inheritance of RP, linkage analysis identified an interval on the X chromosome containing RPGR, and mutation screening revealed a pathogenic variant in this gene. Patients of this family were examined clinically and by X-inactivation studies. RESULTS: This study included 141 RP families with possible X-chromosomal inheritance. In total, we identified 46 families with pathogenic sequence alterations in RPGR and RP2, of which 17 mutations have not been described previously. Two of the novel mutations represent the most 3'-terminal pathogenic sequence variants in RPGR and RP2 reported to date. In exon ORF15 of RPGR, we found eight novel and 14 known mutations. All lead to a disruption of open reading frame. Of the families with suggested X-chromosomal inheritance, 35% showed mutations in ORF15. In addition, we found five novel mutations in other exons of RPGR and four in RP2. Deletions in ORF15 of RPGR were identified in three families in which female carriers showed variable manifestation of the phenotype. Furthermore, an ORF15 mutation was found in an RP patient who additionally carries a 6.4 kbp deletion downstream of the coding region of exon ORF15. We did not identify mutations in 39 sporadic male cases from Switzerland. CONCLUSIONS: RPGR mutations were confirmed to be the most frequent cause of RP in families with an X-chromosomal inheritance pattern. We propose a screening strategy to provide molecular diagnostics in these families.

Arachnomelia in Brown Swiss cattle maps to chromosome 5

Arachnomelia in Brown Swiss cattle is a monogenic autosomal recessive inherited congenital disorder of the skeletal system giving affected calves a spidery look (OMIA ID 000059). Over a period of 20 years 15 cases were sampled in the Swiss and Italian Brown cattle population. Pedigree data revealed that all affected individuals trace back to a single acknowledged carrier founder sire. A genome scan using 240 microsatellites spanning the 29 bovine autosomes showed homozygosity at three adjacent microsatellite markers on bovine Chr 5 in all cases. Linkage analysis confirmed the localization of the arachnomelia mutation in the region of the marker ETH10. Fine-mapping and haplotype analysis using a total of 34 markers in this region refined the critical region of the arachnomelia locus to a 7.19-Mb interval on bovine Chr 5. The disease-associated IBD haplotype was shared by 36 proven carrier animals and allows marker-assisted selection. As the corresponding human and mouse chromosome segments do not contain any clear functional candidate genes for this disorder, the mutation causing arachnomelia in the Brown Swiss cattle might help to identify an unknown gene in bone development.

CHARACTERIZING A NOVEL GENETIC LOCUS ASSOCIATED WITH FAMILIAL CO-OCCURRENCE OF THORACIC AORTIC ANEURYSMS AND INTRACRANIAL ANEURYSMS

The Mendelian inheritance of genetic mutations can lead to adult-onset cardiovascular disease. Several genetic loci have been mapped for the familial form of Thoracic Aortic Aneurysms (TAA), and many causal mutations have been identified for this disease. Intracranial Aneurysms (ICA) also show linkage heterogeneity, but no mutations have been identified causing familial ICA alone. Here, we characterized a large family (TAA288) with an autosomal dominant pattern of inherited aneurysms. It is intriguing that female patients predominantly present with ICA and male patients predominantly with TAA in this family. To identify a causal mutation in this family, a genome-wide linkage analysis was previously performed on nine members of this family using the 50k GenChips Hind array from Affymetrix. This analysis eventually identified a single disease-segregating locus, on chromosome 5p15. We build upon this previous analysis in this study, hypothesizing that a genetic mutation inherited in this locus leads to the sex-specific phenotype of TAA and ICA in this family First we refined the boundaries of the 5p15 disease linked locus down to the genomic coordinates 5p15: 3,424,465- 6,312,925 (GRCh37/hg19 Assembly). This locus was named the TAA288 critical interval. Next, we sequenced candidate genes within the TAA288 critical interval. The selection of genes was simplified by the relatively small number of well-characterized genetic elements within the region. Seeking novel or rare disease-segregating variants, we initially observed a single point alteration in the metalloproteinase gene ADAMTS16 fulfilling this criteria. This variant was later classified as a low-frequency population polymorphism (rs72647757), but we continued to explore the potential role of the ADAMTS16 as the cause of disease in TAA288. We observed that fibroblasts cultured from TAA288 patients consistently upregulated the expression of this gene more strongly compared to matched control fibroblasts when treated with the cytokine TGF-β1, though there was some variation in the exact nature of this expression. We also observed evidence that this protein is expressed at elevated levels in aortic aneurysm tissue from patients with mutations in the gene TGFBR2 and Marfan syndrome, shown by immunohistochemical detection of this protein.

Genetic networks controlling retinal injury.

PURPOSE: The present study defines genomic loci underlying coordinate changes in gene expression following retinal injury. METHODS: A group of acute phase genes expressed in diverse nervous system tissues was defined by combining microarray results from injury studies from rat retina, brain, and spinal cord. Genomic loci regulating the brain expression of acute phase genes were identified using a panel of BXD recombinant inbred (RI) mouse strains. Candidate upstream regulators within a locus were defined using single nucleotide polymorphism databases and promoter motif databases. RESULTS: The acute phase response of rat retina, brain, and spinal cord was dominated by transcription factors. Three genomic loci control transcript expression of acute phase genes in brains of BXD RI mouse strains. One locus was identified on chromosome 12 and was highly correlated with the expression of classic acute phase genes. Within the locus we identified the inhibitor of DNA binding 2 (Id2) as a candidate upstream regulator. Id2 was upregulated as an acute phase transcript in injury models of rat retina, brain, and spinal cord. CONCLUSIONS: We defined a group of transcriptional changes associated with the retinal acute injury response. Using genetic linkage analysis of natural transcript variation, we identified regulatory loci and candidate regulators that control transcript levels of acute phase genes.

A study of reduced penetrance retinoblastoma

Retinoblastoma is a pediatric tumor which is associated with somatic and inherited mutations at the retinoblastoma susceptibility locus, RB1. Although most cases of retinoblastoma fit the previously described 'two hit' model of oncogenesis, the molecular mechanisms underlying rare instances of familial retinoblastoma with reduced penetrance are not well understood. To better understand this phenomenon, a study was undertaken to uncover the molecular cause of low penetrance retinoblastoma in a limited number of families. In one case, a unique cryptic splicing alteration was discovered in the RB1 gene and demonstrated to reduce the level of normal RB1 mRNA produced. Penetrance in the large family known to carry this mutation is less than 50%. Data about the mutation supports a theory that reduced penetrance retinoblastoma is caused by partially functional mutations in RB1. In another family, three independent causes of retinoblastoma or the related phenotype of retinoma were indicated by linkage analysis, a finding unique in retinoblastoma research. A novel polymorphism restricted to Asian populations was also described during the course of this study. ^

Defining the regions of homology between human chromosome 16p12--13 and mouse chromosomes

In this study, the evolutionary relationship between human chromosome 16p12-p13 and mouse chromosomes was investigated by determining the order of marker loci in the region and then identifying the chromosomal locations of the homologous loci in mice. Eighteen genes from human 16 were mapped to fifteen subchromosomal regions by a variety of mapping approaches.^ Thirteen of the genes were mapped in the mouse. Linkage analysis with backcross mice and segregation analysis in a mouse - Chinese Hamster Ovary (CHO) somatic cell hybrid panel informative for different regions of mouse genome were used. The results assigned the thirteen genes to three different mouse chromosomes.^ A group of six genes on mouse 16 was found to be closely linked to Scid. The order of Myh11 and Mrp remains ambiguous since no recombination was detected in backcross analysis. Their relative position in human is also uncertain since they were shown to be very close to each other. For the other mouse loci, an unambiguous gene order could be determined and was found to be identical to that in human. Therefore, they comprise a new conserved linkage group between the two species. The orientation of the group was inverted relative to the centromeres, i.e. the proximal loci in one species become distal in another. The size of the group was estimated to be from 4.4 to 8 Mb and 10 to 32 cM in human. In mouse, it was about 21 cM in the backcross analysis. The two boundaries of the conserved linkage were defined within a 1 Mb range. It is now possible to predict the locations of mouse homologs for some human disease genes based on their locations on human 16p.^ The six human 16p genes that map to MMU7 showed a different gene order in mouse than in human. No recombination was found between Crym and Umod while Crym was distal to D16S79A and proximal to D16S92. The location of Stp and Cdr2 with respect to the above four loci was not determined since they were not mapped in the same set of backcross mice. These genes greatly expanded an existing conserved synteny group between the human 16p12-p13 region and the MMU7. It now consists of eleven loci that span a region of probably more than 10 Mb in human. The gene order derived from this study provided further evidence for chromosomal rearrangements within the conserved synteny. (Abstract shortened by UMI.) ^

Contribution of the genes of the renin-angiotensin system to interindividual differences in blood pressure levels in Caucasians from Rochester, Minnesota

The Renin-Angiotensin system (RAS) regulates blood pressure through its effects on vascular tone, renal hemodynamics, and renal sodium and fluid balance. The genes encoding the four major components of the RAS, angiotensinogen, renin, angiotensin I-converting enzyme (ACE), and angiotensin II receptor type 1 (AT1), have been investigated as candidate genes in the pathogenesis of essential hypertension. However, studies have primarily focused on small samples of diseased individuals, and, therefore, have provided little information about the determinants of interindividual variation in blood pressure (BP) in the general population.^ Using data from a large population-based sample from Rochester, MN, I have evaluated the contribution of variation in the region of the RAS genes to interindividual variation in systolic, diastolic, and mean arterial pressure in the population-at-large. Marker genotype data from four polymorphisms located within or very near these genes were first collected on 3,974 individuals from 583 randomly ascertained three-generation pedigrees. Haseman-Elston regression and variance component methods of linkage analysis were then carried out to estimate the proportion of interindividual variance in BP attributable to the effects of variation at these four measured loci.^ A significant effect of the ACE locus on interindividual variation in mean arterial pressure (MAP) was detected in a sample of siblings belonging to the youngest generation. After allowing for measured covariates, this effect accounted for 15-25% of the interindividual variance in MAP, and was even greater in a subset with a positive family history of hypertension. When gender-specific analyses were carried out, this effect was significant in males but not in females. Extended pedigree analyses also provided evidence for an effect of the ACE locus on interindividual variation in MAP, but no difference between males and females was observed. Circumstantial evidence suggests that the ACE gene itself may be responsible for the observed effects on BP, although the possibility that other genes in the region may be at play cannot be excluded.^ No definitive evidence for an effect of the renin, angiotensinogen, or AT1 loci on interindividual variation in BP was obtained in this study, suggesting that the impact of these genes on BP may not be great in the Caucasian population-at-large. However, this does not preclude a larger effect of these genes in some subsets of individuals, especially among those with clinically manifest hypertension or coronary heart disease, or in other populations. ^

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

Identification of a missense mutation in the bovine ATP2A1 gene in congenital pseudomyotonia of Chianina cattle: an animal model of human Brody disease

Congenital pseudomyotonia in Chianina cattle is a muscle function disorder very similar to that of Brody disease in humans. Mutations in the human ATP2A1 gene, encoding SERCA1, cause Brody myopathy. The analysis of the collected Chianina pedigree data suggested monogenic autosomal recessive inheritance and revealed that all 17 affected individuals traced back to a single founder. A deficiency of SERCA1 function in skeletal muscle of pseudomyotonia affected Chianina cattle was observed as SERCA1 activity in affected animals was decreased by about 70%. Linkage analysis showed that the mutation was located in the ATP2A1 gene region on BTA25 and subsequent mutation analysis of the ATP2A1 exons revealed a perfectly associated missense mutation in exon 6 (c.491G>A) leading to a p.Arg164His substitution. Arg164 represents a functionally important and strongly conserved residue of SERCA1. This study provides a suitable large animal model for human Brody disease.