Comprehensive evaluation of imputation performance in African Americans.

Imputation of genome-wide single-nucleotide polymorphism (SNP) arrays to a larger known reference panel of SNPs has become a standard and an essential part of genome-wide association studies. However, little is known about the behavior of imputation in African Americans with respect to the different imputation algorithms, the reference population(s) and the reference SNP panels used. Genome-wide SNP data (Affymetrix 6.0) from 3207 African American samples in the Atherosclerosis Risk in Communities Study (ARIC) was used to systematically evaluate imputation quality and yield. Imputation was performed with the imputation algorithms MACH, IMPUTE and BEAGLE using several combinations of three reference panels of HapMap III (ASW, YRI and CEU) and 1000 Genomes Project (pilot 1 YRI June 2010 release, EUR and AFR August 2010 and June 2011 releases) panels with SNP data on chromosomes 18, 20 and 22. About 10% of the directly genotyped SNPs from each chromosome were masked, and SNPs common between the reference panels were used for evaluating the imputation quality using two statistical metrics-concordance accuracy and Cohen's kappa (κ) coefficient. The dependencies of these metrics on the minor allele frequencies (MAF) and specific genotype categories (minor allele homozygotes, heterozygotes and major allele homozygotes) were thoroughly investigated to determine the best panel and method for imputation in African Americans. In addition, the power to detect imputed SNPs associated with simulated phenotypes was studied using the mean genotype of each masked SNP in the imputed data. Our results indicate that the genotype concordances after stratification into each genotype category and Cohen's κ coefficient are considerably better equipped to differentiate imputation performance compared with the traditionally used total concordance statistic, and both statistics improved with increasing MAF irrespective of the imputation method. We also find that both MACH and IMPUTE performed equally well and consistently better than BEAGLE irrespective of the reference panel used. Of the various combinations of reference panels, for both HapMap III and 1000 Genomes Project reference panels, the multi-ethnic panels had better imputation accuracy than those containing only single ethnic samples. The most recent 1000 Genomes Project release June 2011 had substantially higher number of imputed SNPs than HapMap III and performed as well or better than the best combined HapMap III reference panels and previous releases of the 1000 Genomes Project.

Vif of feline immunodeficiency virus from domestic cats protects against APOBEC3 restriction factors from many felids.

To get more insight into the role of APOBEC3 (A3) cytidine deaminases in the species-specific restriction of feline immunodeficiency virus (FIV) of the domestic cat, we tested the A3 proteins present in big cats (puma, lion, tiger, and lynx). These A3 proteins were analyzed for expression and sensitivity to the Vif protein of FIV. While A3Z3s and A3Z2-Z3s inhibited Deltavif FIV, felid A3Z2s did not show any antiviral activity against Deltavif FIV or wild-type (wt) FIV. All felid A3Z3s and A3Z2-Z3s were sensitive to Vif of the domestic cat FIV. Vif also induced depletion of felid A3Z2s. Tiger A3s showed a moderate degree of resistance against the Vif-mediated counter defense. These findings may imply that the A3 restriction system does not play a major role to prevent domestic cat FIV transmission to other Felidae. In contrast to the sensitive felid A3s, many nonfelid A3s actively restricted wt FIV replication. To test whether Vif(FIV) can protect also the distantly related human immunodeficiency virus type 1 (HIV-1), a chimeric HIV-1.Vif(FIV) was constructed. This HIV-1.Vif(FIV) was replication competent in nonpermissive feline cells expressing human CD4/CCR5 that did not support the replication of wt HIV-1. We conclude that the replication of HIV-1 in some feline cells is inhibited only by feline A3 restriction factors and the absence of the appropriate receptor or coreceptor.

Molecular genetics and evolution of melanism in the cat family.

Melanistic coat coloration occurs as a common polymorphism in 11 of 37 felid species and reaches high population frequency in some cases but never achieves complete fixation. To investigate the genetic basis, adaptive significance, and evolutionary history of melanistic variants in the Felidae, we mapped, cloned, and sequenced the cat homologs of two putative candidate genes for melanism (ASIP [agouti] and MC1R) and identified three independent deletions associated with dark coloration in three different felid species. Association and transmission analyses revealed that a 2 bp deletion in the ASIP gene specifies black coloration in domestic cats, and two different "in-frame" deletions in the MC1R gene are implicated in melanism in jaguars and jaguarundis. Melanistic individuals from five other felid species did not carry any of these mutations, implying that there are at least four independent genetic origins for melanism in the cat family. The inferred multiple origins and independent historical elevation in population frequency of felid melanistic mutations suggest the occurrence of adaptive evolution of this visible phenotype in a group of related free-ranging species.

Phylogeography and genetic ancestry of tigers (Panthera tigris).

Eight traditional subspecies of tiger (Panthera tigris),of which three recently became extinct, are commonly recognized on the basis of geographic isolation and morphological characteristics. To investigate the species' evolutionary history and to establish objective methods for subspecies recognition, voucher specimens of blood, skin, hair, and/or skin biopsies from 134 tigers with verified geographic origins or heritage across the whole distribution range were examined for three molecular markers: (1) 4.0 kb of mitochondrial DNA (mtDNA) sequence; (2) allele variation in the nuclear major histocompatibility complex class II DRB gene; and (3) composite nuclear microsatellite genotypes based on 30 loci. Relatively low genetic variation with mtDNA,DRB,and microsatellite loci was found, but significant population subdivision was nonetheless apparent among five living subspecies. In addition, a distinct partition of the Indochinese subspecies P. t. corbetti in to northern Indochinese and Malayan Peninsula populations was discovered. Population genetic structure would suggest recognition of six taxonomic units or subspecies: (1) Amur tiger P. t. altaica; (2) northern Indochinese tiger P. t. corbetti; (3) South China tiger P. t. amoyensis; (4) Malayan tiger P. t. jacksoni, named for the tiger conservationist Peter Jackson; (5) Sumatran tiger P. t. sumatrae; and (6) Bengal tiger P. t. tigris. The proposed South China tiger lineage is tentative due to limited sampling. The age of the most recent common ancestor for tiger mtDNA was estimated to be 72,000-108,000 y, relatively younger than some other Panthera species. A combination of population expansions, reduced gene flow, and genetic drift following the last genetic diminution, and the recent anthropogenic range contraction, have led to the distinct genetic partitions. These results provide an explicit basis for subspecies recognition and will lead to the improved management and conservation of these recently isolated but distinct geographic populations of tigers.

Differential gene and microRNA expression between etoposide resistant and etoposide sensitive MCF7 breast cancer cell lines.

In order to develop targeted strategies for combating drug resistance it is essential to understand it's basic molecular mechanisms. In an exploratory study we have found several possible indicators of etoposide resistance operating in MCF7VP cells, including up-regulation of ABC transporter genes, modulation of miRNA, and alteration in copy numbers of genes.

Sequences, annotation and single nucleotide polymorphism of the major histocompatibility complex in the domestic cat.

Two sequences of major histocompatibility complex (MHC) regions in the domestic cat, 2.976 and 0.362 Mbps, which were separated by an ancient chromosome break (55-80 MYA) and followed by a chromosomal inversion were annotated in detail. Gene annotation of this MHC was completed and identified 183 possible coding regions, 147 human homologues, possible functional genes and 36 pseudo/unidentified genes) by GENSCAN and BLASTN, BLASTP RepeatMasker programs. The first region spans 2.976 Mbp sequence, which encodes six classical class II antigens (three DRA and three DRB antigens) lacking the functional DP, DQ regions, nine antigen processing molecules (DOA/DOB, DMA/DMB, TAPASIN, and LMP2/LMP7,TAP1/TAP2), 52 class III genes, nineteen class I genes/gene fragments (FLAI-A to FLAI-S). Three class I genes (FLAI-H, I-K, I-E) may encode functional classical class I antigens based on deduced amino acid sequence and promoter structure. The second region spans 0.362 Mbp sequence encoding no class I genes and 18 cross-species conserved genes, excluding class I, II and their functionally related/associated genes, namely framework genes, including three olfactory receptor genes. One previously identified feline endogenous retrovirus, a baboon retrovirus derived sequence (ECE1) and two new endogenous retrovirus sequences, similar to brown bat endogenous retrovirus (FERVmlu1, FERVmlu2) were found within a 140 Kbp interval in the middle of class I region. MHC SNPs were examined based on comparisons of this BAC sequence and MHC homozygous 1.9x WGS sequences and found that 11,654 SNPs in 2.84 Mbp (0.00411 SNP per bp), which is 2.4 times higher rate than average heterozygous region in the WGS (0.0017 SNP per bp genome), and slightly higher than the SNP rate observed in human MHC (0.00337 SNP per bp).

Functions, structure, and read-through alternative splicing of feline APOBEC3 genes.

BACKGROUND: Over the past years a variety of host restriction genes have been identified in human and mammals that modulate retrovirus infectivity, replication, assembly, and/or cross-species transmission. Among these host-encoded restriction factors, the APOBEC3 (A3; apolipoprotein B mRNA-editing catalytic polypeptide 3) proteins are potent inhibitors of retroviruses and retrotransposons. While primates encode seven of these genes (A3A to A3H), rodents carry only a single A3 gene. RESULTS: Here we identified and characterized several A3 genes in the genome of domestic cat (Felis catus) by analyzing the genomic A3 locus. The cat genome presents one A3H gene and three very similar A3C genes (a-c), probably generated after two consecutive gene duplications. In addition to these four one-domain A3 proteins, a fifth A3, designated A3CH, is expressed by read-through alternative splicing. Specific feline A3 proteins selectively inactivated only defined genera of feline retroviruses: Bet-deficient feline foamy virus was mainly inactivated by feA3Ca, feA3Cb, and feA3Cc, while feA3H and feA3CH were only weakly active. The infectivity of Vif-deficient feline immunodeficiency virus and feline leukemia virus was reduced only by feA3H and feA3CH, but not by any of the feA3Cs. Within Felidae, A3C sequences show significant adaptive selection, but unexpectedly, the A3H sequences present more sites that are under purifying selection. CONCLUSION: Our data support a complex evolutionary history of expansion, divergence, selection and individual extinction of antiviral A3 genes that parallels the early evolution of Placentalia, becoming more intricate in taxa in which the arms race between host and retroviruses is harsher.

Comparative genomic structure of human, dog, and cat MHC: HLA, DLA, and FLA.

Comparisons of the genomic structure of 3 mammalian major histocompatibility complexes (MHCs), human HLA, canine DLA, and feline FLA revealed remarkable structural differences between HLA and the other 2 MHCs. The 4.6-Mb HLA sequence was compared with the 3.9-Mb DLA sequence from 2 supercontigs generated by 7x whole-genome shotgun assembly and 3.3-Mb FLA draft sequence. For FLA, we confirm that 1) feline FLA was split into 2 pieces within the TRIM (member of the tripartite motif) gene family found in human HLA, 2) class II, III, and I regions were placed in the pericentromeric region of the long arm of chromosome B2, and 3) the remaining FLA was located in subtelomeric region of the short arm of chromosome B2. The exact same chromosome break was found in canine DLA structure, where class II, III, and I regions were placed in a pericentromeric region of chromosome 12 whereas the remaining region was located in a subtelomeric region of chromosome 35, suggesting that this chromosome break occurred once before the split of felid and canid more than 55 million years ago. However, significant differences were found in the content of genes in both pericentromeric and subtelomeric regions in DLA and FLA, the gene number, and amplicon structure of class I genes plus 2 other class I genes found on 2 additional chromosomes; canine chromosomes 7 and 18 suggest the dynamic nature in the evolution of MHC class I genes.

Gene conversion between mammalian CCR2 and CCR5 chemokine receptor genes: a potential mechanism for receptor dimerization.

The chemokine receptor genes of the CCR cluster on human chromosome 3p21 play important roles in humoral and cellular immune responses. Several of these receptors have been shown to influence human immunodeficiency virus infection and progression to AIDS, and their homologues may play a role in feline immunodeficiency virus infection. We report the isolation and sequencing of a 150-kb domestic cat BAC clone containing the feline CCR genes CCR1, CCR2, CCR3, and CCR5 to further analyze these four receptor genes within the family Felidae. Comparative and phylogenetic analyses reveal evidence for historic gene conversion between the adjacent CCR2 and CCR5 genes in the Felidae and in three independent mammalian orders (Primates, Cetartiodactyla, and Rodentia), resulting in higher than expected levels of sequence similarity between the two paralogous genes within each order. The gene conversion was restricted to the structural (transmembrane) domains of the CCR2 and CCR5 genes. We also discovered a recent gene conversion event between the third extracellular loop of CCR2 and CCR5 genes that was fixed in Asian lions and found at low frequency in African lions (Panthera leo), suggesting that this domain may have an important functional role. Our results suggest that ongoing parallel gene conversion between CCR2 and CCR5 promotes receptor heterodimerization in independent evolutionary lineages and offers an effective adaptive strategy for gene editing and coevolution among interactive immune response genes in mammals.

A complex rearrangement in GBE1 causes both perinatal hypoglycemic collapse and late-juvenile-onset neuromuscular degeneration in glycogen storage disease type IV of Norwegian forest cats.

Deficiency of glycogen branching enzyme (GBE) activity causes glycogen storage disease type IV (GSD IV), an autosomal recessive error of metabolism. Abnormal glycogen accumulates in myocytes, hepatocytes, and neurons, causing variably progressive, benign to lethal organ dysfunctions. A naturally occurring orthologue of human GSD IV was described previously in Norwegian forest cats (NFC). Here, we report that while most affected kittens die at or soon after birth, presumably due to hypoglycemia, survivors of the perinatal period appear clinically normal until onset of progressive neuromuscular degeneration at 5 months of age. Molecular investigation of affected cats revealed abnormally spliced GBE1 mRNA products and lack of GBE cross-reactive material in liver and muscle. Affected cats are homozygous for a complex rearrangement of genomic DNA in GBE1, constituted by a 334 bp insertion at the site of a 6.2 kb deletion that extends from intron 11 to intron 12 (g. IVS11+1552_IVS12-1339 del6.2kb ins334 bp), removing exon 12. An allele-specific, PCR-based test demonstrates that the rearrangement segregates with the disease in the GSD IV kindred and is not found in unrelated normal cats. Screening of 402 privately owned NFC revealed 58 carriers and 4 affected cats. The molecular characterization of feline GSD IV will enhance further studies of GSD IV pathophysiology and development of novel therapies in this unique animal model.

Development of Y chromosome intraspecific polymorphic markers in the Felidae.

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.

Initial sequence and comparative analysis of the cat genome.

The genome sequence (1.9-fold coverage) of an inbred Abyssinian domestic cat was assembled, mapped, and annotated with a comparative approach that involved cross-reference to annotated genome assemblies of six mammals (human, chimpanzee, mouse, rat, dog, and cow). The results resolved chromosomal positions for 663,480 contigs, 20,285 putative feline gene orthologs, and 133,499 conserved sequence blocks (CSBs). Additional annotated features include repetitive elements, endogenous retroviral sequences, nuclear mitochondrial (numt) sequences, micro-RNAs, and evolutionary breakpoints that suggest historic balancing of translocation and inversion incidences in distinct mammalian lineages. Large numbers of single nucleotide polymorphisms (SNPs), deletion insertion polymorphisms (DIPs), and short tandem repeats (STRs), suitable for linkage or association studies were characterized in the context of long stretches of chromosome homozygosity. In spite of the light coverage capturing approximately 65% of euchromatin sequence from the cat genome, these comparative insights shed new light on the tempo and mode of gene/genome evolution in mammals, promise several research applications for the cat, and also illustrate that a comparative approach using more deeply covered mammals provides an informative, preliminary annotation of a light (1.9-fold) coverage mammal genome sequence.

The feline major histocompatibility complex is rearranged by an inversion with a breakpoint in the distal class I region.

In order to determine the genomic organization of the major histocompatibility complex (MHC) of the domestic cat (Felis catus), DNA probes for 61 markers were designed from human MHC reference sequences and used to construct feline MHC BAC contig map spanning ARE1 in the class II region to the olfactory receptor complex in the extended class I region. Selected BAC clones were then used to identify feline-specific probes for the three regions of the mammalian MHC (class II-class III-class I) for radiation hybrid mapping and fluorescent in situ hybridization to refine the organization of the domestic cat MHC. The results not only confirmed that the p-arm of domestic cat B2 is inverted relative to human Chromosome 6, but also demonstrated that one inversion breakpoint localized to the distal segment of the MHC class I between TRIM39 and TRIM26. The inversion thus disjoined the approximately 2.85 Mb of MHC containing class II-class III-class I (proximal region) from the approximately 0.50 Mb of MHC class I/extended class I region, such that TRIM39 is adjacent to the Chromosome B2 centromere and TRIM26 is adjacent to the B2 telomere in the domestic cat.

The Feline Genome Project.

The compilation of a dense gene map and eventually a whole genome sequence (WGS) of the domestic cat holds considerable value for human genome annotation, for veterinary medicine, and for insight into the evolution of genome organization among mammals. Human association and veterinary studies of the cat, its domestic breeds, and its charismatic wild relatives of the family Felidae have rendered the species a powerful model for human hereditary diseases, for infectious disease agents, for adaptive evolutionary divergence, for conservation genetics, and for forensic applications. Here we review the advantages, rationale, and present strategy of a feline genome project, and we describe the disease models, comparative genomics, and biological applications posed by the full resolution of the cat's genome.

Quantitative polymerase chain reaction-based assay for estimating DNA yield extracted from domestic cat specimens.

A quantitative polymerase chain reaction (PCR) assay has been developed for the quantification of genomic DNA extracted from domestic cat samples. The assay, which targets highly repetitive genomic short interspersed nuclear elements (SINE), can be performed rapidly and is highly sensitive, detecting as little as 10 fg of feline genomic DNA. The assay was linear over a 10(6) dilution range. We have recently developed a short tandem repeat (STR) multiplex panel for forensic analysis of feline specimens. The SINE assay is an integral part of the forensic typing system. The sensitivity of the assay will enable forensic examiners to determine the likelihood of success of genotyping sample extracts with the STR panel without sacrificing valuable DNA necessary to perform genotyping of samples.

Comparative genome organization of human, murine, and feline MHC class II region.

To study comparative molecular dynamics in the genesis of the major histocompatibility complex (MHC), we determined a complete nucleotide sequence spanning 758,291 bp of the domestic cat (Felis catus) extended and classical class II region. The feline class II MHC includes 44 genes (31 predicted to be expressed) which display DNA sequence homology and ordered gene synteny with human HLA and mouse H2, in extended class II and centromere proximal regions (DM to DO) of the classical class II region. However, remarkable genomic alterations including gene gain and loss plus size differentials of 250 kb are evident in comparisons of the cat class II with those of human and mouse. The cat MHC lacks the entire DQ region and retains only relict pseudogene homologs of DP genes, compensated by expansion and reorganization of seven modern DR genes. Repetitive gene families within the feline MHC comprise 35% of the feline MHC with very different density and abundance of GC levels, SINES, LINES, STRs, and retro-elements from the same repeats in human and mouse MHC. Comparison of the feline MHC with the murine and human MHC offers a detailed view of the consequences of genome organization in three mammalian lineages.