The NF-{kappa}B negative regulator TNFAIP3 (A20) is inactivated by somatic mutations and genomic deletions in marginal zone lymphomas

Unique and shared cytogenetic abnormalities have been documented for marginal zone lymphomas (MZLs) arising at different sites. Recently, homozygous deletions of the chromosomal band 6q23, involving the tumor necrosis factor alpha-induced protein 3 (TNFAIP3, A20) gene, a negative regulator of NF-kappaB, were described in ocular adnexal MZL, suggesting a role for A20 as a tumor suppressor in this disease. Here, we investigated inactivation of A20 by DNA mutations or deletions in a panel of extranodal MZL (EMZL), nodal MZL (NMZL), and splenic MZL (SMZL). Inactivating mutations encoding truncated A20 proteins were identified in 6 (19%) of 32 MZLs, including 2 (18%) of 11 EMZLs, 3 (33%) of 9 NMZLs, and 1 (8%) of 12 SMZLs. Two additional unmutated nonsplenic MZLs also showed monoallelic or biallelic A20 deletions by fluorescent in situ hybridization (FISH) and/or SNP-arrays. Thus, A20 inactivation by either somatic mutation and/or deletion represents a common genetic aberration across all MZL subtypes, which may contribute to lymphomagenesis by inducing constitutive NF-kappaB activation.

Genotyping of Francisella tularensis strains by pulsed-field gel electrophoresis, amplified fragment length polymorphism fingerprinting, and 16S rRNA gene sequencing

We evaluated three molecular methods for identification of Francisella strains: pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphism (AFLP) analysis, and 16S rRNA gene sequencing. The analysis was performed with 54 Francisella tularensis subsp. holarctica, 5 F. tularensis subsp. tularensis, 2 F. tularensis subsp. novicida, and 1 F. philomiragia strains. On the basis of the combination of results obtained by PFGE with the restriction enzymes XhoI and BamHI, PFGE revealed seven pulsotypes, which allowed us to discriminate the strains to the subspecies level and which even allowed us to discriminate among some isolates of F. tularensis subsp. holarctica. The AFLP analysis technique produced some degree of discrimination among F. tularensis subsp. holarctica strains (one primary cluster with three major subclusters and minor variations within subclusters) when EcoRI-C and MseI-A, EcoRI-T and MseI-T, EcoRI-A and MseI-C, and EcoRI-0 and MseI-CA were used as primers. The degree of similarity among the strains was about 94%. The percent similarities of the AFLP profiles of this subspecies compared to those of F. tularensis subsp. tularensis, F. tularensis subsp. novicida, and F. philomiragia were less than 90%, about 72%, and less than 24%, respectively, thus permitting easy differentiation of this subspecies. 16S rRNA gene sequencing revealed 100% similarity for all F. tularensis subsp. holarctica isolates compared in this study. These results suggest that although limited genetic heterogeneity among F. tularensis subsp. holarctica isolates was observed, PFGE and AFLP analysis appear to be promising tools for the diagnosis of infections caused by different subspecies of F. tularensis and suitable techniques for the differentiation of individual strains.

Equine cytochrome P450 2B6--genomic identification, expression and functional characterization with ketamine

Ketamine is an anesthetic and analgesic regularly used in veterinary patients. As ketamine is almost always administered in combination with other drugs, interactions between ketamine and other drugs bear the risk of either adverse effects or diminished efficacy. Since cytochrome P450 enzymes (CYPs) play a pivotal role in the phase I metabolism of the majority of all marketed drugs, drug-drug interactions often occur at the active site of these enzymes. CYPs have been thoroughly examined in humans and laboratory animals, but little is known about equine CYPs. The characterization of equine CYPs is essential for a better understanding of drug metabolism in horses. We report annotation, cloning and heterologous expression of the equine CYP2B6 in V79 Chinese hamster fibroblasts. After computational annotation of all CYP2B genes, the coding sequence (CDS) of equine CYP2B6 was amplified by RT-PCR from horse liver total RNA and revealed an amino acid sequence identity of 77% and a similarity of 93.7% to its human ortholog. A non-synonymous variant c.226G>A in exon 2 of the equine CYP2B6 was detected in 97 horses. The mutant A-allele showed an allele frequency of 82%. Two further variants in exon 3 were detected in one and two horses of this group, respectively. Transfected V79 cells were incubated with racemic ketamine and norketamine as probe substrates to determine metabolic activity. The recombinant equine CYP2B6 N-demethylated ketamine to norketamine and produced metabolites of norketamine, such as hydroxylated norketamines and 5,6-dehydronorketamine. V(max) for S-/and R-norketamine formation was 0.49 and 0.45nmol/h/mg cellular protein and K(m) was 3.41 and 2.66μM, respectively. The N-demethylation of S-/R-ketamine was inhibited concentration-dependently with clopidogrel showing an IC(50) of 5.63 and 6.26μM, respectively. The functional importance of the recorded genetic variants remains to be explored. Equine CYP2B6 was determined to be a CYP enzyme involved in ketamine and norketamine metabolism, thus confirming results from inhibition studies with horse liver microsomes. Clopidogrel seems to be a feasible inhibitor for equine CYP2B6. The specificity still needs to be established with other single equine CYPs. Heterologous expression of single equine CYP enzymes opens new possibilities to substantially improve the understanding of drug metabolism and drug interactions in horses.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Sarcoid-derived fibroblasts: links between genomic instability, energy metabolism and senescence.

Bovine papillomavirus 1 (BPV-1) is a well recognized etiopathogenetic factor in a cancer-like state in horses, namely equine sarcoid disease. Nevertheless, little is known about BPV-1-mediated cell transforming effects. It was shown that BPV-1 triggers genomic instability through DNA hypomethylation and oxidative stress. In the present study, we further characterized BPV-1-positive fibroblasts derived from sarcoid tumors. The focus was on cancer-like features of sarcoid-derived fibroblasts, including cell cycle perturbation, comprehensive DNA damage analysis, end-replication problem, energy metabolism and oncogene-induced premature senescence. The S phase of the cell cycle, polyploidy events, DNA double strand breaks (DSBs) and DNA single strand breaks (SSBs) were increased in BPV-1-positive cells compared to control fibroblasts. BPV-1-mediated oxidative stress may contribute to telomere dysfunction in sarcoid-derived fibroblasts. Loss of mitochondrial membrane potential and concurrent elevation in intracellular ATP production may be a consequence of changes in energy-supplying pathways in BPV-1-positive cells which is also typical for cancer cells. Shifts in energy metabolism may support rapid proliferation in cells infected by BPV-1. Nevertheless, sarcoid-derived fibroblasts representing a heterogeneous cell fraction vary in some aspects of metabolic phenotype due to a dual role of BPV-1 in cell transformation and oncogene-induced premature senescence. This was shown with increased senescence-associated β-galactosidase (SA-β-gal) activity. Taken together, metabolic phenotypes in sarcoid-derived fibroblasts are plastic, which are similar to greater plasticity of cancer tissues than normal tissues.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Genomic screening identifies novel linkages and provides further evidence for a role of MYH9 in nonsyndromic cleft lip and palate.

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth anomaly that requires prolonged multidisciplinary rehabilitation. Although variation in several genes has been identified as contributing to NSCLP, most of the genetic susceptibility loci have yet to be defined. To identify additional contributory genes, a high-throughput genomic scan was performed using the Illumina Linkage IVb Panel platform. We genotyped 6008 SNPs in nine non-Hispanic white NSCLP multiplex families and a single large African-American NSCLP multiplex family. Fourteen chromosomal regions were identified with LOD>1.5, including six regions not previously reported. Analysis of the data from the African-American and non-Hispanic white families revealed two likely chromosomal regions: 8q21.3-24.12 and 22q12.2-12.3 with LOD scores of 2.98 and 2.66, respectively. On the basis of biological function, syndecan 2 (SDC2) and growth differentiation factor 6 (GDF6) in 8q21.3-24.12 and myosin heavy-chain 9, non-muscle (MYH9) in 22q12.2-12.3 were selected as candidate genes. Association analyses from these genes yielded marginally significant P-values for SNPs in SDC2 and GDF6 (0.01

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

A conserved motif at the 3$\sp\prime$ end of genomic RNA of mouse hepatitis virus required for host protein binding and viral RNA replication

The initial step in coronavirus-mouse hepatitis virus (MHV) replication is the synthesis of negative strand RNA from a positive strand genomic RNA template. Our approach to studying MHV RNA replication is to identify the cis-acting signals for RNA synthesis and the protein(s) which recognizes these signals at the 3$\sp\prime$ end of genomic RNA of MHV. To determine whether host cellular and/or virus-specific proteins interact with the 3$\sp\prime$ end of the coronavirus genome, an RNase T$\sb1$ protection/gel mobility shift electrophoresis assay was used to examine cytoplasmic extracts from either mock- or MHV-JHM-infected 17Cl-1 murine cells for the ability to form complexes with defined regions of the genomic RNA. A conserved 11 nucleotide sequence UGAAUGAAGUU at nucleotide positions 36 to 26 from the 3$\sp\prime$ end of genomic RNA was identified to be responsible for the specific binding of host proteins, by using a series of RNA probes with deletions and mutations in this region. The RNA probe containing the 11 nucleotide sequence bound approximately four host cellular proteins with a highly labeled 120 kDa and three minor species with sizes of 103, 81 and 55 kDa, assayed by UV-induced covalent cross-linking. Mutation of the 11 nucleotide motif strongly inhibited cellular protein binding, and decreased the amount of the 103 and 81 kDa proteins in the complex to undetectable levels and strongly reduced the binding of the 120 kDa protein. Less extensive mutations within this 11 nucleotide motif resulted in variable decreases in RNA-protein complex formation depending on each probe tested. The RNA-protein complexes observed with cytoplasmic extracts from MHV-JHM-infected cells in both RNase protection/gel mobility shift and UV cross-linking assays were indistinguishable to those observed with extracts from uninfected cells.^ To investigate the possible role of this 3$\sp\prime$ protein binding element in viral RNA replication in vivo, defective interfering RNA molecules with complete or partial mutations of the 11 nucleotide conserved sequence were transcribed in vitro, transfected to host 17Cl-1 cells in the presence of helper virus MHV-JHM and analyzed by agarose gel electrophoresis, competitive RT-PCR and direct sequencing of the RT-PCR products. Both negative strand synthesis and positive strand replication of DI RNA were affected by mutation that disrupts RNA-protein complex formation, even though the 11 mutated nucleotides were converted to wild type sequence, presumably by recombination with helper virus. Kinetic analysis indicated that recombination between DI RNA and helper virus occurred 5.5 to 7.5 hours post infection when replication of positive strand DI RNA was barely observed. Replication of positive strand DI RNAs carrying partial mutations within the 11 nucleotide motif was dependent upon recombination events after transfection. Replication was strongly inhibited when reversion to wild type sequence did not occur, and after recombination, reached similar levels as wild type DI RNA. A DI RNA with mutation upstream of the protein binding motif replicated as efficiently as wild type without undergoing recombination. Thus the conserved 11 nucleotide host protein binding motif appears to play an important role in viral RNA replication. ^

TUMOR PROGRESSION, METASTATIC POTENTIAL AND GENOMIC INSTABILITY

To investigate the hypothesis that increased malignant potential correlates with increased levels of genetic instability, the following parameters of instability were measured: (1) spontaneous mutation rates for ouabain resistance in murine cell lines of different malignant potentials, (2) the background prevalence of 6-thioguanine (6-TG) resistance in clone 4 (highly metastatic) and clone 19 (poorly metastatic) of the K1735 murine melanoma, (3) the prevalence of ouabain resistant variants in three murine cell lines and their variants after exposure to the mutagen MNNG, (4) the rate of generation of major karyotypic abnormalities in B16 F1 (poorly metastatic) and B16 F10 (highly metastatic) murine melanoma, and (5) analysis of the G-banded karyotypes of cloned B16 F1 and B16 F10 melanoma.^ No correlation of increased spontaneous mutation rates with increased malignant potential was found in repeated experiments with three murine cell lines and their variants of different malignant potential. The background prevalence of g-TG resistance was not significantly different for the poorly and highly metastatic clones of K1735 melanoma. The studies with MNNG-induced mutation showed no increased sensitivity of the highly metastatic variants of the three murine cell lines to mutagenesis. Neither did the rate of generation of major karyotypic abnormalities correlate with malignant potential. However, certain karyotypic differences were demonstrated after G-banding of the B16 F1 and F10 melanomas.^ One hypothesis which is consistent with these results is that the rate of generation of genetic abnormalities need not be strongly related to the degree of malignant potential. An increased prevalence of genetic changes may merely reflect the accumulation of abnormalities while their rate of production remains constant. The presence of specific nonrandom changes likely is the main determinant of malignant potential rather than the rate of production of random changes. ^