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

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

A Nonsense Variant in COL6A1 in Landseer Dogs with Muscular Dystrophy.

A novel canine muscular dystrophy in Landseer dogs was observed. We had access to five affected dogs from two litters. The clinical signs started at a few weeks of age and the severe progressive muscle weakness led to euthanasia between 5 and 15 months of age. The pedigrees of the affected dogs suggested a monogenic autosomal recessive inheritance of the trait. Linkage and homozygosity mapping indicated two potential genome segments for the causative variant on chromosomes 10 and 31 harboring a total of 4.8 Mb of DNA or 0.2% of the canine genome. Using the illumina sequencing technology we obtained a whole genome sequence from one affected Landseer. Variants were called with respect to the dog reference genome and compared to the genetic variants of 170 control dogs from other breeds. The affected Landseer dog was homozygous for a single private non-synonymous variant in the critical intervals, a nonsense variant in the COL6A1 gene (Chr31:39,303,964G>T; COL6A1:c.289G>T; p.E97*). Genotypes at this variant showed perfect concordance with the muscular dystrophy phenotype in all five cases and more than one thousand control dogs. Variants in the human COL6A1 gene cause Bethlem myopathy or Ullrich congenital muscular dystrophy. We therefore conclude that the identified canine COL6A1 variant is most likely causative for the observed muscular dystrophy in Landseer dogs. Based on the nature of the genetic variant in Landseer dogs and their severe clinical phenotype these dogs represent a model for human Ullrich congenital muscular dystrophy.

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

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

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

AsaGEI2b: a new variant of a genomic island identified in the Aeromonas salmonicida subsp. salmonicida JF3224 strain isolated from a wild fish in Switzerland.

Aeromonas salmonicida subsp. salmonicida is the causal agent of furunculosis in salmonids. We recently identified a group of genomic islands (AsaGEI) in this bacterium. AsaGEI2a, one of these genomic islands, has almost exclusively been identified in isolates from North America. To date, Aeromonas salmonicida subsp. salmonicida JF3224, a strain isolated from a wild brown trout (Salmo trutta) caught in Switzerland, was the only European isolate that appeared to bear AsaGEI2a. We analyzed the genome of JF3224 and showed that the genomic island in JF3224 is a new variant of AsaGEI, which we have called AsaGEI2b. While AsaGEI2b shares the same integrase gene and insertion site as AsaGEI2a, it is very different in terms of many other features. Additional genomic investigations combined with PCR genotyping revealed that JF3224 is sensitive to growth at 25°C, leading to insertion sequence-dependent rearrangement of the locus on the pAsa5 plasmid that encodes a type three secretion system, which is essential for the virulence of the bacterium. The analysis of the JF3224 genome confirmed that AsaGEIs are accurate indicators of the geographic origins of A. salmonicida subsp. salmonicida isolates and is another example of the susceptibility of the pAsa5 plasmid to DNA rearrangements.

A variant of DNA polymerase β acts as a dominant negative mutant

In eukaryotic cells, DNA polymerase β (polβ) carries out base-excision repair (BER) required for DNA maintenance, replication, recombination, and drug resistance. A specific deletion in one allele in the coding sequence of the polβ gene occurs in colorectal and breast carcinomas. The 87-bp deleted region encodes amino acid residues 208–236 in the catalytic domain of the enzyme. Here, we report evidence for expression of the wild-type (WT) and the truncated polβ proteins in colorectal tumors. To elucidate the potential functional consequences of polβ truncation, stable HeLa cell lines were established from cloned WT and variant polβΔ208–236. Cells expressing the variant protein exhibited substantially decreased BER activity. To test our hypothesis that truncated polβ may disrupt the function of the WT enzyme, we stably transfected mouse embryonic fibroblast 16.3 cells with polβΔ208–236 cDNA. Reverse transcription–PCR and Western blot analyses showed that the new cell line, 16.3ΔP, expresses the WT and the truncated polβ mRNA and proteins. BER and binding activities were undetectable in these cells. Furthermore, in vivo the 16.3ΔP cells were more sensitive to N-methyl-N′-nitro-N-nitrosoguanidine than the 16.3 cells. On adding increasing amounts of 16.3ΔP protein extracts, the BER and DNA binding activities of extracts of the parent 16.3 cell line progressively declined. These results strongly suggest that truncated polβ acts as a dominant negative mutant. The defective polβ may facilitate accumulation of mutations, leading to the expression of a mutator phenotype in tumor cells.

Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Genomic and cDNA sequence analysis of the cell matrix adhesion regulator gene

The cell matrix adhesion regulator (CMAR) gene has been suggested to be a signal transduction molecule influencing cell adhesion to collagen and, through this, possibly involved in tumor suppression. The originally reported CMAR cDNA was 464 bp long with a tyrosine phosphorylation site at the extreme 3′ end, which mutagenesis studies had shown to be central to the function of this gene. Since the discovery of a 4-bp insertion polymorphism within the originally reported coding region, further sequence information has been obtained. The cDNA has been extended 5′ by ≈2 kb revealing a 559-bp region showing strong homology to the proposed 5′ untranslated sequence of a murine protein kinase receptor family member, variant in kinase (vik). CMAR genomic sequencing has shown the presence of an intron, the intron/exon boundary lying within this region of homology. An RNA transcript for CMAR of ≈2.5 kb has also been identified. The data suggest complex mechanisms for control of expression of two closely associated genes, CMAR and the vik- associated sequence.

Variant antigens and endothelial receptor adhesion in Plasmodium falciparum.

Parasite-derived proteins expressed on the surface of erythrocytes infected with Plasmodium falciparum are important virulence factors, since they mediate binding of infected cells to diverse receptors on vascular endothelium and are targets of a protective immune response. They are difficult to study because they undergo rapid clonal antigenic variation in vitro, which precludes the derivation of phenotypically homogeneous cultures. Here we have utilized sequence-specific proteases to dissect the role of defined antigenic variants in binding to particular receptors. By selection of protease-resistant subpopulations of parasites on defined receptors we (i) confirm the high rate of antigenic variation in vitro; (ii) demonstrate that a single infected erythrocyte can bind to intercellular adhesion molecule 1, CD36, and thrombospondin; (iii) show that binding to intercellular adhesion molecule 1 and CD36 are functions of the variant antigen; and (iv) suggest that binding to thrombospondin may be mediated by other components of the infected erythrocyte surface.

Vascular variant of prion protein cerebral amyloidosis with tau-positive neurofibrillary tangles: the phenotype of the stop codon 145 mutation in PRNP.

Deposition of PrP amyloid in cerebral vessels in conjunction with neurofibrillary lesions is the neuropathologic hallmark of the dementia associated with a stop mutation at codon 145 of PRNP, the gene encoding the prion protein (PrP). In this disorder, the vascular amyloid in tissue sections and the approximately 7.5-kDa fragment extracted from amyloid are labeled by antibodies to epitopes located in the PrP sequence including amino acids 90-147. Amyloid-laden vessels are also labeled by antibodies against the C terminus, suggesting that PrP from the normal allele is involved in the pathologic process. Abundant neurofibrillary lesions are present in the cerebral gray matter. They are composed of paired helical filaments, are labeled with antibodies that recognize multiple phosphorylation sites in tau protein, and are similar to those observed in Alzheimer disease. A PrP cerebral amyloid angiopathy has not been reported in diseases caused by PRNP mutations or in human transmissible spongiform encephalopathies; we propose to name this phenotype PrP cerebral amyloid angiopathy (PrP-CAA).

Cell cycle regulation of the cyclin A gene promoter is mediated by a variant E2F site.

Cyclin A is involved in the control of S phase and mitosis in mammalian cells. Expression of the cyclin A gene in nontransformed cells is characterized by repression of its promoter during the G1 phase of the cell cycle and its induction at S-phase entry. We show that this mode of regulation is mediated by the transcription factor E2F, which binds to a specific site in the cyclin A promoter. It differs from the prototype E2F site in nucleotide sequence and protein binding; it is bound by E2F complexes containing cyclin E and p107 but not pRB. Ectopic expression of cyclin D1 triggers premature activation of the cyclin A promoter by E2F, and this effect is blocked by the tumor suppressor protein p16INK4.

Use of a novel DNA melting profile assay for the identification of PCR-amplified genomic sequences encoding for variant-specific surface proteins from the clonal GS/M-83-H7 line of Giardia lamblia.

During infections, Giardia lamblia undergoes a continuous change of its major surface antigens, the variant-specific surface proteins (VSPs). Many studies on antigenic variation have been performed using G. lamblia clone GS/M-83-H7, which expresses surface antigen VSP H7. The present study was focused on the identification and characterization of vsp gene sequences within the genome of the clonal G. lamblia GS/M-83-H7 line. For this purpose, we applied a PCR which specifically amplified truncated sequences from the 3'-terminal region of the vsp genes. Upon cloning, most of the vsp gene amplification products were shown to be approximately identical in size and thus could not be distinguished from each other by conventional gel electrophoresis. In order to pre-estimate the sequence complexity within the large panel of vsp clones isolated, we elaborated a novel concept which facilitated our large-scale genetic screening approach: PCR products from cloned DNA molecules were generated and then subjected to a DNA melting profile assay based on the use of the LightCycler Instrument. This high-throughput assay system proved to be well suited to monitor sequence differences between the amplification products from closely related vsp genes and thus could be used for the primary, sequence-related discrimination of the corresponding clones. After testing 50 candidates, vsp clones could be divided into five groups, each characterized by an individual DNA melting profile of the corresponding amplification products. Sequence analysis of some of these 50 candidates confirmed data from the aforementioned assay in that clones were demonstrated to be identical within, but different between, the distinct groups. The nucleotide and deduced amino acid sequences of five representative vsp clones showed high similarities both among each other and also with the corresponding gene segment of the variant-specific surface antigen (VSP H7) expressed by the original GS/M-83-H7 variant type. Furthermore, three of the genomic vsp sequences turned out to be identical to vsp sequences that represented previously characterized transcription products from in vivo- or in vitro-switched GS/M-83-H7 trophozoites. In conclusion, the DNA melting profile assay seems to be a versatile tool for the PCR-based genotyping of moderately or highly diversified sequence orthologues.

Activation of the cAMP pathway by variant human MC1R alleles expressed in HEK and in melanoma cells

alpha-Melanocyte-stimulating hormone (alpha-MSH) activates the melanocortin-1 receptor (MC1R) on melanocytes to promote a switch from red/yellow pheomelanin synthesis to darker eumelanins via positive coupling to adenylate cyclase. The human MC1R locus is highly polymorphic with the specific variants associated with red hair and fair skin (RHC phenotype) postulated to be loss-of-function receptors. We have examined the ability of MC1R variants to activate the cAMP pathway in stably transfected REK293 cells. The RHC associated variants, Arg151Cys, Arg160Trp and Asp294His, demonstrated agonist-mediated increases in cAMP and phosphorylation of cAMP-responsive element-binding protein (CREB). Whereas the Asp294His variant showed severely impaired functional responses, the Arg151Cys and Arg160Trp variants retained considerable signaling capacity. Melanoma cells homozygous for either the Arg151Cys variant or consensus sequence both elicited CREB phosphorylation in response to alpha-MSH in the presence of IBMX. The common RHC alleles, Arg151Cys, Arg160Trp and Asp294His, are neither complete loss-of-function receptors nor are they functionally equivalent. (c) 2005 Elsevier Inc. All rights reserved.

FRAP analysis of nucleocytoplasmic dynamics of the vitamin D receptor splice variant VDRB1: preferential targeting to nuclear speckles

Although the key components of the cellular nuclear transport machinery have largely been characterized through extensive efforts in recent years, in vivo measurements of the kinetics of nuclear protein import/export are patently few. The present study applies the approach of FRAP (fluorescence recovery after photobleaching) to examine the nucleocytoplasmic flux of a novel human VDRB1 (vitamin D receptor B I) isoform in living cells. Through an N-terminal extension containing a consensus nuclear targeting sequence, VDRB1 is capable of localizing in nuclear speckles adjacent to SC-35 (35 kDa splicing component)containing speckles as well as in the nucleoplasm, dependent on ligand. Investigation of VDRB1 nucleocytoplasmic transport using FRAP indicates for the first time that the VDRB1 has a serum-modulated, active nuclear-import mechanism. There is no evidence of an efficient, active export mechanism for VDRB1, probably as a result of nuclear retention. VDRB1 nuclear import in the absence of serum occurred more rapidly and to a greater extent to nuclear speckles compared with import to other nuclear sites. This preferential transport from the cytoplasm to and accumulation within nuclear speckles is consistent with the idea that the latter represent dynamic centres of VDRB1 interaction with other nuclear proteins. The results are consistent with the existence of specialized pathways to target proteins to nuclear subdomains.

The variant rpoS allele of S. enteritidis strain 27655R does not affect virulence in a chick model nor constitutive curliation but does generate a cold-sensitive phenotype

Adherence of pathogenic Escherichia coli and Salmonella spp. to host cells is in part mediated by curli fimbriae which, along with other virulence determinants, are positively regulated by RpoS. Interested in the role and regulation of curli (SEF17) fimbriae of Salmonella enteritidis in poultry infection, we tested the virulence of naturally occurring S. enteritidis PT4 strains 27655R and 27655S which displayed constitutive and null expression of curli (SEF17) fimbriae, respectively, in a chick invasion assay and analysed their rpoS alleles. Both strains were shown to be equally invasive and as invasive as a wild-type phage type 4 strain and an isogenic derivative defective for the elaboration of curli. We showed that the rpoS allele of 27655S was intact even though this strain was non-curliated and we confirmed that a S. enteritidis rpoS::strr null mutant was unable to express curli, as anticipated. Strain 27655R, constitutively curliated, possessed a frameshift mutation at position 697 of the rpoS coding sequence which resulted in a truncated product and remained curliated even when transduced to rpoS::strr. Additionally, rpoS mutants are known to be cold-sensitive, a phenotype confirmed for strain 27655R. Collectively, these data indicated that curliation was not a significant factor for pathogenesis of S. enteritidis in this model and that curliation of strains 27655R and 27655S was independent of RpoS. Significantly, strain 27655R possessed a defective rpoS allele and remained virulent. Here was evidence that supported the concept that different naturally occurring rpoS alleles may generate varying virulence phenotypic traits.