Integrating understanding of epidemiology and genomics in B-cell non-Hodgkin lymphoma as a pathway to novel management strategies.

Non-Hodgkin lymphomas include a biologically and clinically heterogeneous group of cancers distinguished by genetics, histology, and treatment outcomes. New discoveries regarding the genomic alterations and epidemiological exposures associated with these lymphomas have enhanced our understanding of factors that contribute to lymphomagenesis for specific subtypes. We explore the impact of normal B-cell biology engineered for recognizing a wide variety of antigens on the development of specific lymphoma subtypes, review lymphoma genetics, and examine the epidemiology of B-cell NHLs including recent investigations of risk factors for particular lymphoma subtypes based on large pooled analyses. Burkitt lymphoma, an aggressive form of B-cell NHL involving translocation of the MYC gene and an immunoglobulin gene has been associated with a history of eczema, hepatitis C, and occupation as a cleaner. Increased risk of diffuse large B-cell lymphoma has been associated with increased young adult body mass index, history of B-cell-activating autoimmune diseases, hepatitis C, and several single nucleotide variants involving the human leukocyte antigen (HLA) region of chromosome 6 and non-HLA loci near EXOC2, PVT1, MYC, and NCOA1. Tumor sequencing studies suggest that multiple pathways are involved in the development of DLBCL. Additional studies of epidemiological exposures, genome wide associations, and tumor sequencing in follicular, lymphoplasmacytic, marginal zone, and mantle cell lymphoma demonstrate overlapping areas of increased risk factors and unique factors for specific subtypes. Integrating these findings is important for constructing comprehensive models of NHL pathogenesis, which could yield novel targets for therapy and strategies for lymphoma prevention in certain populations.

A male-specific QTL for social interaction behavior in mice mapped with automated pattern detection by a hidden Markov model incorporated into newly developed freeware.

BACKGROUND: Owing to their complex nature, social interaction tests normally require the observation of video data by a human researcher, and thus are difficult to use in large-scale studies. We previously established a statistical method, a hidden Markov model (HMM), which enables the differentiation of two social states ("interaction" and "indifference"), and three social states ("sniffing", "following", and "indifference"), automatically in silico. NEW METHOD: Here, we developed freeware called DuoMouse for the rapid evaluation of social interaction behavior. This software incorporates five steps: (1) settings, (2) video recording, (3) tracking from the video data, (4) HMM analysis, and (5) visualization of the results. RESULTS: Using DuoMouse, we mapped a genetic locus related to social interaction. We previously reported that a consomic strain, B6-Chr6C(MSM), with its chromosome 6 substituted for one from MSM/Ms, showed more social interaction than C57BL/6 (B6). We made four subconsomic strains, C3, C5, C6, and C7, each of which has a shorter segment of chromosome 6 derived from B6-Chr6C, and conducted social interaction tests on these strains. DuoMouse indicated that C6, but not C3, C5, and C7, showed higher interaction, sniffing, and following than B6, specifically in males. COMPARISON WITH EXISTING METHOD: The data obtained by human observation showed high concordance to those from DuoMouse. The results indicated that the MSM-derived chromosomal region present in C6-but not in C3, C5, and C7-associated with increased social behavior. CONCLUSIONS: This method to analyze social interaction will aid primary screening for difference in social behavior in mice.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

A hidden Markov random field model for genome-wide association studies.

Genome-wide association studies (GWAS) are increasingly utilized for identifying novel susceptible genetic variants for complex traits, but there is little consensus on analysis methods for such data. Most commonly used methods include single single nucleotide polymorphism (SNP) analysis or haplotype analysis with Bonferroni correction for multiple comparisons. Since the SNPs in typical GWAS are often in linkage disequilibrium (LD), at least locally, Bonferroni correction of multiple comparisons often leads to conservative error control and therefore lower statistical power. In this paper, we propose a hidden Markov random field model (HMRF) for GWAS analysis based on a weighted LD graph built from the prior LD information among the SNPs and an efficient iterative conditional mode algorithm for estimating the model parameters. This model effectively utilizes the LD information in calculating the posterior probability that an SNP is associated with the disease. These posterior probabilities can then be used to define a false discovery controlling procedure in order to select the disease-associated SNPs. Simulation studies demonstrated the potential gain in power over single SNP analysis. The proposed method is especially effective in identifying SNPs with borderline significance at the single-marker level that nonetheless are in high LD with significant SNPs. In addition, by simultaneously considering the SNPs in LD, the proposed method can also help to reduce the number of false identifications of disease-associated SNPs. We demonstrate the application of the proposed HMRF model using data from a case-control GWAS of neuroblastoma and identify 1 new SNP that is potentially associated with neuroblastoma.

Assessment of transmission distortion on chromosome 6p in healthy individuals using tagSNPs.

The best-documented example for transmission distortion (TD) to normal offspring are the t haplotypes on mouse chromosome 17. In healthy humans, TD has been described for whole chromosomes and for particular loci, but multiple comparisons have presented a statistical obstacle in wide-ranging analyses. Here we provide six high-resolution TD maps of the short arm of human chromosome 6 (Hsa6p), based on single-nucleotide polymorphism (SNP) data from 60 trio families belonging to two ethnicities that are available through the International HapMap Project. We tested all approximately 70,000 previously genotyped SNPs within Hsa6p by the transmission disequilibrium test. TagSNP selection followed by permutation testing was performed to adjust for multiple testing. A statistically significant evidence for TD was observed among male parents of European ancestry, due to strong and wide-ranging skewed segregation in a 730 kb long region containing the transcription factor-encoding genes SUPT3H and RUNX2, as well as the microRNA locus MIRN586. We also observed that this chromosomal segment coincides with pronounced linkage disequilibrium (LD), suggesting a relationship between TD and LD. The fact that TD may be taking place in samples not selected for a genetic disease implies that linkage studies must be assessed with particular caution in chromosomal segments with evidence of TD.

Stochastic search gene suggestion: a Bayesian hierarchical model for gene mapping.

Mapping the genes for a complex disease, such as diabetes or rheumatoid arthritis (RA), involves finding multiple genetic loci that may contribute to the onset of the disease. Pairwise testing of the loci leads to the problem of multiple testing. Looking at haplotypes, or linear sets of loci, avoids multiple tests but results in a contingency table with sparse counts, especially when using marker loci with multiple alleles. We propose a hierarchical Bayesian model for case-parent triad data that uses a conditional logistic regression likelihood to model the probability of transmission to a diseased child. We define hierarchical prior distributions on the allele main effects to model the genetic dependencies present in the human leukocyte antigen (HLA) region of chromosome 6. First, we add a hierarchical level for model selection that accounts for both locus and allele selection. This allows us to cast the problem of identifying genetic loci relevant to the disease into a problem of Bayesian variable selection. Second, we attempt to include linkage disequilibrium as a covariance structure in the prior for model coefficients. We evaluate the performance of the procedure with some simulated examples and then apply our procedure to identifying genetic markers in the HLA region that influence risk for RA. Our software is available on the website http://www.epigenetic.org/Linkage/ssgs-public/.

Assessment of the role of common genetic variation in the transient neonatal diabetes mellitus (TNDM) region in type 2 diabetes: a comparative genomic and tagging single nucleotide polymorphism approach.

Recent evidence supports the strong overlap between genes implicated in monogenic diabetes and susceptibility to type 2 diabetes. Transient neonatal diabetes mellitus (TNDM) is a rare disorder associated with overexpression of genes at a paternally expressed imprinted locus on chromosome 6q24. There are two overlapping genes in this region: the transcription factor zinc finger protein associated with cell cycle control and apoptosis (ZAC also known as PLAGL1) and HYMA1, which encodes an untranslated mRNA. Several type 2 diabetes linkage studies have reported linkage to chromosome 6q22-25. We hypothesized that common genetic variation at this TNDM region influences type 2 diabetes susceptibility. In addition to the coding regions, we used comparative genomic analysis to identify conserved noncoding regions, which were resequenced for single nucleotide polymorphism (SNP) discovery in 47 individuals. Twenty-six SNPs were identified. Fifteen tag SNPs (tSNPs) were successfully genotyped in a large case-control (n = 3,594) and family-based (n = 1,654) study. We did not find any evidence of association or overtransmission of any tSNP to affected offspring or of a parent-of-origin effect. Using a study sufficiently powered to detect odds ratios of <1.2, we conclude that common variation in the TNDM region does not play an important role in the genetic susceptibility to type 2 diabetes.

Phenotypic and molecular assessment of seven patients with 6p25 deletion syndrome: relevance to ocular dysgenesis and hearing impairment.

BACKGROUND: Thirty-nine patients have been described with deletions involving chromosome 6p25. However, relatively few of these deletions have had molecular characterization. Common phenotypes of 6p25 deletion syndrome patients include hydrocephalus, hearing loss, and ocular, craniofacial, skeletal, cardiac, and renal malformations. Molecular characterization of deletions can identify genes that are responsible for these phenotypes. METHODS: We report the clinical phenotype of seven patients with terminal deletions of chromosome 6p25 and compare them to previously reported patients. Molecular characterization of the deletions was performed using polymorphic marker analysis to determine the extents of the deletions in these seven 6p25 deletion syndrome patients. RESULTS: Our results, and previous data, show that ocular dysgenesis and hearing impairment are the two most highly penetrant phenotypes of the 6p25 deletion syndrome. While deletion of the forkhead box C1 gene (FOXC1) probably underlies the ocular dysgenesis, no gene in this region is known to be involved in hearing impairment. CONCLUSIONS: Ocular dysgenesis and hearing impairment are the two most common phenotypes of 6p25 deletion syndrome. We conclude that a locus for dominant hearing loss is present at 6p25 and that this locus is restricted to a region distal to D6S1617. Molecular characterization of more 6p25 deletion patients will aid in refinement of this locus and the identification of a gene involved in dominant hearing loss.

Assessment of killer cell immunoglobulinlike receptor expression and corresponding HLA class I phenotypes demonstrates heterogenous KIR expression independent of anticipated HLA class I ligands.

Natural killer (NK) cell-mediated cytolysis is stimulated and downregulated through the interaction of distinct human leukocyte antigen (HLA) class I molecules on target cells with specific killer cell immunoglobulinlike receptors (KIRs) on NK cells. Killer cell immunoglobulinlike receptors are highly polymorphic and are clonally distributed on NK cell populations within individuals. However, the regulation of KIR expression by individual HLA class I phenotypes is not well understood. To examine a potential influence of the HLA class I phenotype on KIR expression patterns we studied the KIR expression in individuals that were subgrouped according to the major HLA-C encoded KIR-epitopes (group C1 versus C2). In these individuals, NK cells were analyzed for KIR expression using flow cytometry and RNA-based expression analysis. Our results demonstrate that KIR genes are transmitted very heterogeneously with two main patterns of KIR genotypes as previously described; group A and group B (with 21 different genotypes). There are distinct populations exhibiting different densities of CD158a and/or CD158b positive NK cells that coexist in all individuals. A clear correlation between KIR expression and the currently known HLA class I ligands was not observed. In conclusion, the surface expression of KIRs in individuals with different HLA class I genotypes indicates that other non-HLA class I encoded factors contribute to the shaping of the KIR repertoire.

Nonsyndromic cleft lip with or without cleft palate in China: assessment of candidate regions.

OBJECTIVE: Although Asians have the highest birth prevalence of oral-facial clefts, the majority of gene mapping studies of cleft lip with or without cleft palate (CL/P) have been in European or American Caucasians. Therefore, the objective of this study of Chinese families was to evaluate linkage and association between CL/P and 10 genetic markers in five chromosomal regions that have shown positive results in Caucasians. SETTING: Families were ascertained through nonsyndromic CL/P surgical probands from hospitals throughout Shanghai, China. PARTICIPANTS: Study participants included 671 individuals from 60 families with two or more members affected with oral-facial clefts. Of the 671 total individuals, 145 were affected. RESULTS: Ten markers from chromosomes 2, 4, 6, 17, and 19 were assessed (TGFA, MSX1, D4S194, D4S175, F13A1, GATA185H, D17S250, D17S579, D19S49, APOC2). LOD scores were calculated between each of the 10 markers and CL/P as well as model-free statistics of linkage (SimIBD) and association (TDT). None of the markers showed significantly positive LOD scores with CL/P. A significantly positive result (p =.01) was seen using SimIBD for APOC2 on chromosome 19, and a positive TDT result (p =.004) was obtained for D19S49, near APOC2. CONCLUSIONS: This is the first gene mapping study of CL/P in China. These results indicate that most of the genetic regions with positive results in Caucasian families may not be involved in CL/P found in China, although there is some positive evidence for the candidate region on chromosome 19.