Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

LOH at 6p21.3 region and HLA class altered phenotypes in bladder carcinomas

Alterations in HLA class I antigen expression have been frequently described in different epithelial tumors and are thought to favor tumor immune escape from T lymphocyte recognition. Multiple molecular mechanisms are responsible for these altered HLA class I tumor phenotypes. Some are structural defects that produce unresponsiveness to treatment with interferons. Others include alterations in regulatory mechanisms that can be switched on by treatment of tumor cells with different cytokines. One important mechanism belonging to the first group is loss of heterozygosity (LOH) at chromosome region 6p21.3, which can lead to HLA haplotype loss. In this investigation, the frequency of LOH at 6p21 chromosome region was studied in 69 bladder carcinomas. Short tandem repeat analysis showed that 35% of cases had LOH in this chromosome region. By considering these results together with immunohistological findings previously published by our group, we identified a distribution pattern of HLA class I altered phenotypes in bladder cancer. The most frequently altered phenotype in bladder carcinomas was total loss of HLA class I expression (17 cases, 25%), followed by phenotype II associated with HLA haplotype loss (12 cases, 17.5%), and HLA allelic loss (ten cases, 14.5%). Nine cases (13%) were classified as having a compound phenotype, five cases (7%) as having HLA locus loss, and in 16 cases (23%) no alteration in HLA expression was detected. An important conclusion of this report is that a combination of different molecular and immunohistological techniques is required to precisely define which HLA alleles are lost during tumor progression and to characterize the underlying mechanisms of these losses. These studies should be performed when a cancer patient is to be included in an immunotherapy protocol that aims to stimulate different immune effector mechanisms.     

A Genome Wide Study of Copy Number Variation Associated with Nasopharyngeal Carcinoma in Malaysian Chinese Identifies CNVs at 11q14.3 and 6p21.3 as Candidate Loci

Background

Nasopharyngeal carcinoma (NPC) is a neoplasm of the epithelial lining of the nasopharynx. Despite various reports linking genomic variants to NPC predisposition, very few reports were done on copy number variations (CNV). CNV is an inherent structural variation that has been found to be involved in cancer predisposition.

Methods

A discovery cohort of Malaysian Chinese descent (NPC patients, n = 140; Healthy controls, n = 256) were genotyped using Illumina^® HumanOmniExpress BeadChip. PennCNV and cnvPartition calling algorithms were applied for CNV calling. Taqman CNV assays and digital PCR were used to validate CNV calls and replicate candidate copy number variant region (CNVR) associations in a follow-up Malaysian Chinese (NPC cases, n = 465; and Healthy controls, n = 677) and Malay cohort (NPC cases, n = 114; Healthy controls, n = 124).

Results

Six putative CNVRs overlapping GRM5, MICA/HCP5/HCG26, LILRB3/LILRA6, DPY19L2, RNase3/RNase2 and GOLPH3 genes were jointly identified by PennCNV and cnvPartition. CNVs overlapping GRM5 and MICA/HCP5/HCG26 were subjected to further validation by Taqman CNV assays and digital PCR. Combined analysis in Malaysian Chinese cohort revealed a strong association at CNVR on chromosome 11q14.3 (P_combined = 1.54x10^-5; odds ratio (OR) = 7.27; 95% CI = 2.96–17.88) overlapping GRM5 and a suggestive association at CNVR on chromosome 6p21.3 (P_combined = 1.29x10^-3; OR = 4.21; 95% CI = 1.75–10.11) overlapping MICA/HCP5/HCG26 genes.

Conclusion

Our results demonstrated the association of CNVs towards NPC susceptibility, implicating a possible role of CNVs in NPC development.     

A Genome-wide Association Study of Lung Cancer Identifies a Region of Chromosome 5p15 Associated with Risk for Adenocarcinoma

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR] = 1.23, 95% confidence interval [CI] = 1.13–1.33, p = 3.02 × 10⁻⁷), but not with other histologic types (OR = 1.01, p = 0.84 and OR = 1.00, p = 0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR = 1.24, 95% CI = 1.17–1.31, p = 3.74 × 10⁻¹⁴ for AD; OR = 0.99, p = 0.69 and OR = 0.97, p = 0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.     

A cosmid library specific for human Chromosome regions 6p21.3 and 6q27

We have explored the potential of irradiation-fusion gene transfer (IFGT) hybrids as a source of well-defined human chromosome fragments from which probes can be derived. Extensive characterization of the IFGT hybrid 4J4 with a full panel of markers from Chromosome (Chr) 6 showed that the human DNA content derives largely from 6p21.3 and 6q27. A cosmid library has been constructed from 4J4 DNA, and 370 recombinants containing human DNA have been isolated and overlapping clones ordered into 20 contigs. Regional localization of representative clones from each contig, determined by fluorescent in situ hybridization (FISH), places 13 contigs in 6q27 and 6 in 6p21.3. Preliminary screening of cDNA libraries with selected cosmids has identified two expressed sequences. Since there are a number of medically important genes in both these regions of human Chr 6 with several disease loci linked to the HLA-A region in 6p21.3 and various tumor suppressor genes to 6q27, this library will provide a valuable resource to aid the isolation of candidate genes for these diseases. In addition, unique markers for detailed physical and genetic mapping of these regions of human Chr 6 can be easily obtained.     

Genome-Wide Association Study Reveals Novel Quantitative Trait Loci Associated with Resistance to Multiple Leaf Spot Diseases of Spring Wheat

Accelerated wheat development and deployment of high-yielding, climate resilient, and disease resistant cultivars can contribute to enhanced food security and sustainable intensification. To facilitate gene discovery, we assembled an association mapping panel of 528 spring wheat landraces of diverse geographic origin for a genome-wide association study (GWAS). All accessions were genotyped using an Illumina Infinium 9K wheat single nucleotide polymorphism (SNP) chip and 4781 polymorphic SNPs were used for analysis. To identify loci underlying resistance to the major leaf spot diseases and to better understand the genomic patterns, we quantified population structure, allelic diversity, and linkage disequilibrium. Our results showed 32 loci were significantly associated with resistance to the major leaf spot diseases. Further analysis identified QTL effective against major leaf spot diseases of wheat which appeared to be novel and others that were previously identified by association analysis using Diversity Arrays Technology (DArT) and bi-parental mapping. In addition, several identified SNPs co-localized with genes that have been implicated in plant disease resistance. Future work could aim to select the putative novel loci and pyramid them in locally adapted wheat cultivars to develop broad-spectrum resistance to multiple leaf spot diseases of wheat via marker-assisted selection (MAS).     

Genome-wide Association Study in a High-Risk Isolate for Multiple Sclerosis Reveals Associated Variants in STAT3 Gene

Genetic risk for multiple sclerosis (MS) is thought to involve both common and rare risk alleles. Recent GWAS and subsequent meta-analysis have established the critical role of the HLA locus and identified new common variants associated to MS. These variants have small odds ratios (ORs) and explain only a fraction of the genetic risk. To expose potentially rare, high-impact alleles, we conducted a GWAS of 68 distantly related cases and 136 controls from a high-risk internal isolate of Finland with increased prevalence and familial occurrence of MS. The top 27 loci with p < 10⁻⁴ were tested in 711 cases and 1029 controls from Finland, and the top two findings were validated in 3859 cases and 9110 controls from more heterogeneous populations. SNP (rs744166) within the STAT3 gene was associated to MS (p = 2.75 × 10⁻¹⁰, OR 0.87, confidence interval 0.83–0.91). The protective haplotype for MS in STAT3 is a risk allele for Crohn disease, implying that STAT3 represents a shared risk locus for at least two autoimmune diseases. This study also demonstrates the potential of special isolated populations in search for variants contributing to complex traits.     

Deletion Mapping of Four Loci Defined by N-Ethyl-N-Nitrosourea-Induced Postimplantation-Lethal Mutations within the Pid-Hbb Region of Mouse Chromosome 7

As part of a long-term effort to refine the physical and functional maps of the Fes-Hbb region of mouse chromosome 7, four loci [l(7)1Rn, l(7)2Rn, l(7)3Rn, l(7)4Rn] defined by N-ethyl-N-nitrosourea (ENU)-induced, prenatally lethal mutations were mapped by means of trans complementation crosses to mice carrying lethal deletions of the mouse chromosome-7 albino (c) locus. Each locus was assigned to a defined subregion of the deletion map at the distal end of the Fes-Hbb interval. Of particular use for this mapping were preimplantation-lethal deletions having distal breakpoints localized between pid and Omp. Hemizygosity or homozygosity for each of the ENU-induced lethals was found to arrest development after uterine implantation; the specific time of postimplantation death varied, and depended on both the mutation itself and on whether it was hemizygous or homozygous. Based on their map positions outside of and distal to deletions that cause death at preimplantation stages, these ENU-induced mutations identify loci, necessary for postimplantation development, that could not have been discovered by phenotypic analyses of mice homozygous for any albino deletion. The mapping of these loci to specific genetic intervals defined by deletion breakpoints suggests a number of positional-cloning strategies for the molecular isolation of these genes. Phenotypic and genetic analyses of these mutations should provide useful information on the functional composition of the corresponding segment of the human genome (perhaps human 11q13.5).     

A Locus on Chromosome 1p36 Is Associated with Thyrotropin and Thyroid Function as Identified by Genome-wide Association Study

Thyroid hormones are key regulators of cellular growth, development, and metabolism, and thyroid disorders are a common cause of ill health in the community. Circulating concentrations of thyrotropin (TSH), thyroxine (T4) and triiodothyronine (T3) have a strong heritable component and are thought to be under polygenic control, but the genes responsible are mostly unknown. In order to identify genetic loci associated with these metabolic phenotypes, we performed a genome-wide association study of 2,120,505 SNPs in 2014 female twins from the TwinsUK study and found a significant association between rs10917469 on chromosome 1p36.13 and serum TSH (p = 3.2 × 10⁻⁸). The association of rs10917469 with serum TSH was replicated (p = 2.0 × 10⁻⁴) in an independent community-based sample of 1154 participants in the Busselton Health Study. This SNP is located near CAPZB, which might be a regulator of TSH secretion and thus of pituitary-thyroid axis function. Twenty-nine percent of white individuals carry the variant, and the difference in mean TSH concentrations between wild-type individuals and those homozygous for the minor G allele was 0.5 mU/l, which is likely to be clinically relevant. We also provide evidence of suggestive association (p < 5.0 × 10⁻⁶) of other SNPs with serum TSH, free T4, and free T3 concentrations, and these SNPs might be good targets for further studies. These results advance understanding of the genetic basis of pituitary-thyroid axis function and metabolic regulation.     

Detection of a polymorphism within the pepsinogen C gene with PCR: construction of a linkage map around PGC from 6p11-6p21.3.

An insertion/deletion polymorphism between exons 7 and 8 of the pepsinogen C gene (PGC), previously detectable with Southern analysis, was formatted for detection with PCR. Alleles were rapidly typed by UV irradiation of ethidium bromide-stained agarose gels. Whereas Southern analysis revealed two alleles, the smaller fragments generated with PCR allowed the resolution of three alleles that were previously scored as a single allele and increased the heterozygosity of the system from 0.20 to 0.53. After a set of reference families was genotyped with the PCR-based polymorphism, a linkage map around the PGC gene on chromosome 6 was constructed. This included the HLA cluster and the highly informative D6S223 locus. PGC lies 22 cM proximal to HLA-DPB and between D6S5 and D6S4 at distances of 4.5 and 13.1 cM, respectively.     

A Large Cohort Study Reveals the Association of Elevated Peripheral Blood Lymphocyte-to-Monocyte Ratio with Favorable Prognosis in Nasopharyngeal Carcinoma

Background

Nasopharyngeal carcinoma (NPC) is an endemic neoplasm in southern China. Although NPC sufferers are sensitive to radiotherapy, 20–30% of patients finally progress with recurrence and metastases. Elevated lymphocyte-to-monocyte ratio (LMR) has been reported to be associated with favorable prognosis in some hematology malignancies, but has not been studied in NPC. The aim of this study was to evaluate whether LMR could predict the prognosis of NPC patients.

Methods

A retrospective cohort of 1,547 non-metastatic NPC patients was recruited between January 2005 and June 2008. The counts for peripheral lymphocyte and monocyte were retrieved, and the LMR was calculated. Receiver operating characteristic curve analysis, univariate and multivariate COX proportional hazards analyses were applied to evaluate the associations of LMR with overall survival (OS), disease-free survival (DFS), distant metastasis-free survival (DMFS) and loco-regional recurrence-free survival (LRRFS), respectively.

Results

Univariate analysis revealed that higher LMR level (≥5.220) was significantly associated with superior OS, DFS and DMFS (P values <0.001). The higher lymphocyte count (≥2.145×10⁹/L) was significantly associated with better OS (P = 0.002) and DMFS (P = 0.031), respectively, while the lower monocyte count (<0.475×10⁹/L) was associated with better OS (P = 0.012), DFS (P = 0.011) and DMFS (P = 0.003), respectively. Multivariate Cox proportional hazard analysis showed that higher LMR level was a significantly independent predictor for superior OS (hazard ratio or HR  = 0.558, 95% confidence interval or 95% CI  = 0.417–0.748; P<0.001), DFS (HR  = 0.669, 95% CI  = 0.535–0.838; P<0.001) and DMFS (HR = 0.543, 95% CI  = 0.403–0.732; P<0.001), respectively. The advanced T and N stages were also independent indicators for worse OS, DFS, and DMFS, except that T stage showed borderline statistical significance for DFS (P = 0.053) and DMFS (P = 0.080).

Conclusions

The elevated pretreatment peripheral LMR level was a significant favorable factor for NPC prognosis and this easily accessed variable may serve as a potent marker to predict the outcomes of NPC patients.     

A Genome Wide Association Study of Mathematical Ability Reveals an Association at Chromosome 3q29, a Locus Associated with Autism and Learning Difficulties: A Preliminary Study

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10⁻⁵, 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10⁻⁶). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10⁻⁴). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.     

A gene for hypotrichosis simplex of the scalp maps to chromosome 6p21.3

Hypotrichosis simplex of the scalp (HSS) is an autosomal dominant form of isolated alopecia causing almost complete loss of scalp hair, with onset in childhood. After exclusion of candidate regions previously associated with hair-loss disorders, we performed a genomewide linkage analysis in two Danish families and localized the gene to chromosome 6p21.3. This was confirmed in a Spanish family, with a total LOD score of 11.97 for marker D6S1701 in all families. The combined haplotype data identify a critical interval of 14.9 cM between markers D6S276 and D6S1607. Localization of the locus for HSS to 6p21.3 is a first step toward identification of the gene. The gene will give important insights into the molecular and cellular basis of hair growth on the scalp.