Toward the most ideal case–control design with related and unrelated dogs in whole‐exome sequencing studies

With the recent development of whole‐exome sequencing enrichment designs for the dog, a novel tool for disease‐association studies became available. The aim of disease‐association studies is to identify one or a very limited number of putative causal variants or genes from the large pool of genetic variation. To maximize the efficiency of these studies and to provide some directions of what to expect, we evaluated the effect on variant reduction for various combinations of cases and controls for both dominant and recessive types of inheritance assuming variable degrees of penetrance and detectance. In this study, variant data of 14 dogs (13 Labrador Retrievers and one Dogue de Bordeaux), obtained by whole‐exome sequencing, were analyzed. In the filtering process, we found that unrelated dogs from the same breed share up to 70% of their variants, which is likely a consequence of the breeding history of the dog. For the designs tested with unrelated dogs, combining two cases and two controls gave the best result. These results were improved further by adding closely related dogs. Reduced penetrance and/or detectance has a drastic effect on the efficiency and is likely to have a profound effect on the sample size needed to elucidate the causal variant. Overall, we demonstrated that sequencing a small number of dogs results in a marked reduction of variants that are likely sufficient to pinpoint causal variants or genes.     

Antisense PMO Found in Dystrophic Dog Model Was Effective in Cells from Exon 7-Deleted DMD Patient

Background

Antisense oligonucleotide-induced exon skipping is a promising approach for treatment of Duchenne muscular dystrophy (DMD). We have systemically administered an antisense phosphorodiamidate morpholino oligomer (PMO) targeting dystrophin exons 6 and 8 to a dog with canine X-linked muscular dystrophy in Japan (CXMD_J) lacking exon 7 and achieved recovery of dystrophin in skeletal muscle. To date, however, antisense chemical compounds used in DMD animal models have not been directly applied to a DMD patient having the same type of exon deletion. We recently identified a DMD patient with an exon 7 deletion and tried direct translation of the antisense PMO used in dog models to the DMD patient''s cells.

Methodology/Principal Findings

We converted fibroblasts of CXMD_J and the DMD patient to myotubes by FACS-aided MyoD transduction. Antisense PMOs targeting identical regions of dog and human dystrophin exons 6 and 8 were designed. These antisense PMOs were mixed and administered as a cocktail to either dog or human cells in vitro. In the CXMD_J and human DMD cells, we observed a similar efficacy of skipping of exons 6 and 8 and a similar extent of dystrophin protein recovery. The accompanying skipping of exon 9, which did not alter the reading frame, was different between cells of these two species.

Conclusion/Significance

Antisense PMOs, the effectiveness of which has been demonstrated in a dog model, achieved multi-exon skipping of dystrophin gene on the FACS-aided MyoD-transduced fibroblasts from an exon 7-deleted DMD patient, suggesting the feasibility of systemic multi-exon skipping in humans.     

Comparative sequence analysis and radiation hybrid mapping of the canine keratin 10 gene.

The type I keratin, K10, is expressed in epidermal keratinocytes undergoing terminal differentiation to form the stratum corneum, a barrier essential for life. In order to facilitate the study of keratinization disorders in the dog, the sequence and mapping of KRT10 is reported. The coding region of KRT10 is 1707 bp and is comprised of eight exons. Although the length of KRT10 has been reported to be polymorphic in humans, this was not observed in the eight domestic dog breeds studied, although one wild canid displayed a size difference. The structure and sequence of this gene is highly conserved across mammalian species. Canine K10 had an 86% amino acid identity with the human gene. KRT10 was localized to the on-going canine radiation hybrid map to chromosome 9 in the type I keratin gene cluster.     

Structure and sequence of the human homeobox gene HOX7.

A cosmid containing the human sequence HOX7, homologous to the murine Hox-7 gene, was isolated from a genomic library, and the positions of the coding sequences were determined by hybridization. DNA sequence analysis demonstrated two exons that code for a homeodomain-containing protein of 297 amino acids. The open reading frame is interrupted by a single intron of approximately 1.6 kb, the splice donor and acceptor sites of which conform to known consensus sequences. The human HOX7 coding sequence has a very high degree of identity with the murine Hox-7 cDNA. Within the homeobox, the two sequences share 94% identity at the DNA level, all substitutions being silent. This high level of sequence similarity is not confined to the homeodomain; overall the human and murine HOX7 gene products show 80% identity at the amino acid level. Both the 5' and 3' untranslated regions also show significant similarity to the murine gene, with 79 and 70% sequence identity, respectively. The sequence upstream of the coding sequence of exon 1 contains a GC-rich putative promoter region. There is no TATA box, but a CCAAT and numerous GC boxes are present. The region encompassing the promoter region, exon 1, and the 5' region of exon 2 have a higher than expected frequency of CpG dinucleotides; numerous sites for rare-cutter restriction enzymes are present, a characteristic of HTF islands.     

A comparative analysis of MC4R gene sequence, polymorphism, and chromosomal localization in Chinese raccoon dog and Arctic fox

Numerous mutations of the human melanocortin receptor type 4 (MC4R) gene are responsible for monogenic obesity, and some of them appear to be associated with predisposition or resistance to polygenic obesity. Thus, this gene is considered a functional candidate for fat tissue accumulation and body weight in domestic mammals. The aim of the study was comparative analysis of chromosome localization, nucleotide sequence, and polymorphism of the MC4R gene in two farmed species of the Canidae family, namely the Chinese raccoon dog (Nycterutes procyonoides procyonoides) and the arctic fox (Alopex lagopus). The whole coding sequence, including fragments of 3'UTR and 5'UTR, shows 89% similarity between the arctic fox (1276 bp) and Chinese raccoon dog (1213 bp). Altogether, 30 farmed Chinese raccoon dogs and 30 farmed arctic foxes were searched for polymorphisms. In the Chinese raccoon dog, only one silent substitution in the coding sequence was identified; whereas in the arctic fox, four InDels and two single-nucleotide polymorphisms (SNPs) in the 5'UTR and six silent SNPs in the exon were found. The studied gene was mapped by FISH to the Chinese raccoon dog chromosome 9 (NPP9q1.2) and arctic fox chromosome 24 (ALA24q1.2-1.3). The obtained results are discussed in terms of genome evolution of species belonging to the family Canidae and their potential use in animal breeding.     

大鼠FMR1同源基因cDNA的克隆、测序与选择剪接的初步分析

应用ＲＴ－ＰＣＲ方法,从新生大鼠脑组织总ＲＮＡ扩增大鼠ＦＭＲ１同源基因的ｃＤＮＡ片段,克降至ｐＵＣ１８质粒中进行序列分析．获得从终止密码子起共１６８１ｂｐ的编码序列,尚缺少约２００ｂｐ的５′序列．所克隆的这部分大鼠ＦＭＲ１ｃＤＮＡ,不含有对应于人ＦＭＲ１基因的外显子１２及外显子１７第一和第三剪接受点之间的序列,提示大鼠ＦＭＲ１基因也有选择剪接表达．同源性分析显示,大鼠ＦＭＲ１与小鼠ＦＭＲ１基因的同源性为９７．７％,与人ＦＭＲ１基因的同源性为９４．９％;与小鼠ＦＭＲＰ（ＦＭＲ１蛋白）的氨基酸序列同源性为９８．４％,与人ＦＭＲＰ的氨基酸序列同源性为９７．９％．以大鼠ＦＭＲ１ｃＤＮＡ片段为探针检测到大鼠不同组织中ＦＭＲ１基因的选择剪接表达．上述结果为以大鼠为动物模型深入研究ＦＭＲ１基因功能奠定了基础．     

Amplification of CFTR Exon 9 Sequences to Multiple Locations in the Human Genome

Cloning and characterization of the cystic fibrosis transmembrane conductance regulator (CFTR) gene led to the identification and isolation of cDNA and genomic sequences that cross-hybridized to the first nucleotide binding fold of CFTR. DNA sequence analysis of these clones showed that the cross-hybridizing sequences corresponded to CFTR exon 9 and its flanking introns, juxtapositioned with two segments of LINE1 sequences. The CFTR sequence appeared to have been transcribed from the opposite direction of the gene, reversely transcribed, and co-integrated with the L1 sequences into a chromosome location distinct from that of the CFTR locus. Based on hybridization intensity and complexity of the restriction fragments, it was estimated that there were at least 10 copies of the “amplified” CFTR exon 9 sequences in the human genome. Furthermore, when DNA segments adjacent to the insertion site were used in genomic DNA blot hybridization analysis, multiple copies were also detected. The overall similarity between these CFTR exon 9-related sequences suggested that they were derived from a single retrotransposition event and subsequent sequence amplification. The amplification unit appeared to be greater than 30 kb. Physical mapping studies includingin situhybridization to human metaphase chromosomes showed that multiple copies of these amplified sequences (with and without the CFTR exon 9 insertion) were dispersed throughout the genome. These findings provide insight into the structure and evolution of the human genome.     

Keratin 9 gene mutational heterogeneity in patients with epidermolytic palmoplantar keratoderma

Mutations in the human keratin 9 gene have recently been shown to be involved in the etiology of palmoplantar keratoderma (PPK). We have investigated eleven unrelated German kindreds with the epidermolytic variant of PPK (EPPK) for mutations in the keratin 9 gene. We have identified two novel mutations, M156V and Q171P, both in the coil 1A segment of keratin 9. Mutation M156V was detected in two unrelated patients with EPPK, and mutation Q171P was shown to cosegregate with the disease in a large four-generation family. These findings confirm the functional importance of coil 1A integrity for heterodimerisation in keratins and for intermediate filament assembly. Our results provide further evidence for mutational heterogeneity in EPPK, and for the involvement of keratins in diseases of hyperkeratinisation and epidermolysis.     

Relationships between the gene and protein structure in human complement component C9

Human complement component C9 is a multidomain protein for which a large number of surface topographical features have been determined. We have analyzed the exon-intron boundaries of the human C9 gene and find a good correlation between splice sites and surface features of the protein but little correlation with the putative protein domain structure, even in the cysteine-rich sequence homology with the low-density lipoprotein (LDL) receptor which is likely to be an independently folded structural motif. This is surprising because in the LDL receptor the same sequence is precisely bounded by introns, and it has been assumed that this sequence is present in both proteins as a result of exon shuffling. We deduce that substantial rearrangement of the exon-intron structure of the C9 gene must have occurred before the exchange of cysteine-rich domains, possibly linked to the process of exon duplication which was required to generate the repeats in the LDL receptor.     

Splicing factors induce cystic fibrosis transmembrane regulator exon 9 skipping through a nonevolutionary conserved intronic element

In monosymptomatic forms of cystic fibrosis such as congenital bilateral absence of vas deferens, variations in the TG(m) and T(n) polymorphic repeats at the 3' end of intron 8 of the cystic fibrosis transmembrane regulator (CFTR) gene are associated with the alternative splicing of exon 9, which results in a nonfunctional CFTR protein. Using a minigene model system, we have previously shown a direct relationship between the TG(m)T(n) polymorphism and exon 9 splicing. We have now evaluated the role of splicing factors in the regulation of the alternative splicing of this exon. Serine-arginine-rich proteins and the heterogeneous nuclear ribonucleoprotein A1 induced exon skipping in the human gene but not in its mouse counterpart. The effect of these proteins on exon 9 exclusion was strictly dependent on the composition of the TG(m) and T(n) polymorphic repeats. The comparative and functional analysis of the human and mouse CFTR genes showed that a region of about 150 nucleotides, present only in the human intron 9, mediates the exon 9 splicing inhibition in association with exonic regulatory elements. This region, defined as the CFTR exon 9 intronic splicing silencer, is a target for serine-arginine-rich protein interactions. Thus, the nonevolutionary conserved CFTR exon 9 alternative splicing is modulated by the TG(m) and T(n) polymorphism at the 3' splice region, enhancer and silencer exonic elements, and the intronic splicing silencer in the proximal 5' intronic region. Tissue levels and individual variability of splicing factors would determine the penetrance of the TG(m)T(n) locus in monosymptomatic forms of cystic fibrosis.     

Inactivation of human keratin genes: the spectrum of mutations in the sequence of an acidic keratin pseudogene

Keratins are cytoskeletal proteins encoded by a multigene family. We have identified the first human keratin pseudogene and determined its complete nucleotide sequence. Sequence comparisons indicate that the pseudogene arose from a very recent duplication of the 50-kd keratin (K14) gene. The coding and the intron sequences of the two genes are 95% and 93% identical, respectively. Although the sequence of the regulatory region in the pseudogene is virtually identical to that in the 50-kd functional gene, several deleterious mutations have been identified in the pseudogene. There are three frameshifts in the coding regions, one of which is a perfect 8-bp duplication. A single-base-pair deletion in the first exon and a single-base-pair insertion in the penultimate exon also result in frameshifts. The three remaining deleterious mutations interfere with the mRNA processing signals: two alter the intron/exon boundaries, and the third disrupts the polyadenylation signal. These mutations clearly identify the sequence as a human keratin pseudogene.     

Genomic structure and expression of human guanosine monophosphate reductase

Summary In vitro translation in the rabbit reticulocyte system and transient expression in Cos7 cells were performed to characterize the protein encoded by a chromosome 6-linked human cDNA clone, whose nucleotide sequence is homologous to that of Escherichia coli guanosine monophosphate reductase (GMP reductase) cDNA. The molecular weight of the peptide produced by the cDNA was about 37,000 Dalton, and the protein produced in the Cos7 cells exhibited GMP reductase activity, substantiating that the cDNA is for human GMP reductase. The corresponding genomic clones were obtained from two human genomic libraries. The gene spans about 50 Kb and is composed of 9 exons, which encode 345 amino acid residues. Organization of exons and introns was established by DNA sequencing of each exon and splicing junctions. The gene contains two potential SpI binding sites within exon 1, and a functional atypical polyadenylation signal in exon 9.