Genome-wide Study Identifies Association between HLA-B∗55:01 and Self-Reported Penicillin Allergy.

Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B∗55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B∗55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.

Maternal uniparental isodisomy of chromosome 6 unmasks a novel variant in TULP1 in a patient with early onset retinal dystrophy.

Purpose: Inherited retinal dystrophies are a clinically and genetically heterogeneous group of disorders. Molecular diagnosis has proven utility for affected individuals. In this study, we report an individual enrolled in the Australian Inherited Retinal Disease Registry and DNA Bank diagnosed with clinical features overlapping between Leber congenital amaurosis and retinitis pigmentosa. Methods: DNA from the proband was sequenced using a gene panel for inherited retinal disorders, and a single nucleotide polymorphism (SNP) array was conducted to detect the presence of deletions and uniparental disomy. Results: We identified a novel homozygous variant (c.524dupC, p.(Pro176ThrfsTer7)) in TULP1 resulting from maternal uniparental isodisomy of chromosome 6. The patient had clinical features consistent with biallelic pathogenic variants in TULP1, including congenital nystagmus, night blindness, non-recordable electroretinogram, mild myopia, and mild peripheral pigmentary changes in the fundus. Conclusions: This is the first report of uniparental disomy 6 and a homozygous variant in TULP1 associated with a rod-cone dystrophy. Molecular diagnosis of inherited retinal dystrophies is essential to inform the mode of transmission and clinical management, and to identify potential candidates for future gene-specific therapies.

Duplication events downstream of IRX1 cause North Carolina macular dystrophy at the MCDR3 locus.

Autosomal dominant North Carolina macular dystrophy (NCMD) is believed to represent a failure of macular development. The disorder has been linked to two loci, MCDR1 (chromosome 6q16) and MCDR3 (chromosome 5p15-p13). Recently, non-coding variants upstream of PRDM13 (MCDR1) and a duplication including IRX1 (MCDR3) have been identified. However, the underlying disease-causing mechanism remains uncertain. Through a combination of sequencing studies on eighteen NCMD families, we report two novel overlapping duplications at the MCDR3 locus, in a gene desert downstream of IRX1 and upstream of ADAMTS16. One duplication of 43 kb was identified in nine families (with evidence for a shared ancestral haplotype), and another one of 45 kb was found in a single family. Three families carry the previously reported V2 variant (MCDR1), while five remain unsolved. The MCDR3 locus is thus refined to a shared region of 39 kb that contains DNAse hypersensitive sites active at a restricted time window during retinal development. Publicly available data confirmed expression of IRX1 and ADAMTS16 in human fetal retina, with IRX1 preferentially expressed in fetal macula. These findings represent a major advance in our understanding of the molecular genetics of NCMD and provide insights into the genetic pathways involved in human macular development.

FOXO3 longevity interactome on chromosome 6.

FOXO3 has been implicated in longevity in multiple populations. By DNA sequencing in long-lived individuals, we identified all single nucleotide polymorphisms (SNPs) in FOXO3 and showed 41 were associated with longevity. Thirteen of these had predicted alterations in transcription factor binding sites. Those SNPs appeared to be in physical contact, via RNA polymerase II binding chromatin looping, with sites in the FOXO3 promoter, and likely function together as a cis-regulatory unit. The SNPs exhibited a high degree of LD in the Asian population, in which they define a specific longevity haplotype that is relatively common. The haplotype was less frequent in whites and virtually nonexistent in Africans. We identified distant contact points between FOXO3 and 46 neighboring genes, through long-range physical contacts via CCCTC-binding factor zinc finger protein (CTCF) binding sites, over a 7.3 Mb distance on chromosome 6q21. When activated by cellular stress, we visualized movement of FOXO3 toward neighboring genes. FOXO3 resides at the center of this early-replicating and highly conserved syntenic region of chromosome 6. Thus, in addition to its role as a transcription factor regulating gene expression genomewide, FOXO3 may function at the genomic level to help regulate neighboring genes by virtue of its central location in chromatin conformation via topologically associated domains. We believe that the FOXO3 'interactome' on chromosome 6 is a chromatin domain that defines an aging hub. A more thorough understanding of the functions of these neighboring genes may help elucidate the mechanisms through which FOXO3 variants promote longevity and healthy aging.

Rare genetic variants in SMAP1, B3GAT2, and RIMS1 contribute to pediatric venous thromboembolism.

Recent genome-wide association studies (GWAS) have confirmed known risk mutations for venous thromboembolism (VTE) and identified a number of novel susceptibility loci in adults. Here we present a GWAS in 212 nuclear families with pediatric VTE followed by targeted next-generation sequencing (NGS) to identify causative mutations contributing to the association. Three single nucleotide polymorphisms (SNPs) exceeded the threshold for genome-wide significance as determined by permutation testing using 100 000 bootstrap permutations (P < 10-5). These SNPs reside in a region on chromosome 6q13 comprising the genes small ARF GAP1 (SMAP1), an ARF6 guanosine triphosphatase-activating protein that functions in clathrin-dependent endocytosis, and ß-1,3-glucoronyltransferase 2 (B3GAT2), a member of the human natural killer 1 carbohydrate pathway. Rs1304029 and rs2748331 are associated with pediatric VTE with unpermuted/permuted values of P = 1.42 × 10-6/2.0 × 10-6 and P = 6.11 × 10-6/1.8 × 10-5, respectively. Rs2748331 was replicated (P = .00719) in an independent study sample coming from our GWAS on pediatric thromboembolic stroke (combined P = 7.88 × 10-7). Subsequent targeted NGS in 24 discordant sibling pairs identified 17 nonsynonymous coding variants, of which 1 located in SMAP1 and 3 in RIMS1, a member of the RIM family of active zone proteins, are predicted as damaging by Protein Variation Effect Analyzer and/or sorting intolerant from tolerant scores. Three SNPs curtly missed statistical significance in the transmission-disequilibrium test in the full cohort (rs112439957: P = .08326, SMAP1; rs767118962: P = .08326, RIMS1; and rs41265501: P = .05778, RIMS1). In conjunction, our data provide compelling evidence for SMAP1, B3GAT2, and RIMS1 as novel susceptibility loci for pediatric VTE and warrant future functional studies to unravel the underlying molecular mechanisms leading to VTE.

Analysis of Plasminogen Genetic Variants in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a prevalent neurological disease of complex etiology. Here, we describe the characterization of a multi-incident MS family that nominated a rare missense variant (p.G420D) in plasminogen (PLG) as a putative genetic risk factor for MS. Genotyping of PLG p.G420D (rs139071351) in 2160 MS patients, and 886 controls from Canada, identified 10 additional probands, two sporadic patients and one control with the variant. Segregation in families harboring the rs139071351 variant, identified p.G420D in 26 out of 30 family members diagnosed with MS, 14 unaffected parents, and 12 out of 30 family members not diagnosed with disease. Despite considerably reduced penetrance, linkage analysis supports cosegregation of PLG p.G420D and disease. Genotyping of PLG p.G420D in 14446 patients, and 8797 controls from Canada, France, Spain, Germany, Belgium, and Austria failed to identify significant association with disease (P = 0.117), despite an overall higher prevalence in patients (OR = 1.32; 95% CI = 0.93-1.87). To assess whether additional rare variants have an effect on MS risk, we sequenced PLG in 293 probands, and genotyped all rare variants in cases and controls. This analysis identified nine rare missense variants, and although three of them were exclusively observed in MS patients, segregation does not support pathogenicity. PLG is a plausible biological candidate for MS owing to its involvement in immune system response, blood-brain barrier permeability, and myelin degradation. Moreover, components of its activation cascade have been shown to present increased activity or expression in MS patients compared to controls; further studies are needed to clarify whether PLG is involved in MS susceptibility.

PLAGL1 epimutation and bladder exstrophy: Coincidence or concurrent etiology?

BACKGROUND: The bladder exstrophy-epispadias complex (BEEC) is characterized by a spectrum of genitourinary malformations. Both classical bladder exstrophy and the most severe phenotype, exstrophy of the cloaca, display omphaloceles, a cardinal anomaly of some disorders caused by altered imprinting. Therefore, we hypothesized that BEEC in some patients could occur on the basis of an undiagnosed imprinting disorder. Such altered imprinting is associated with changes in the parent-of-origin-specific DNA methylation. METHODS: We analyzed the DNA methylation of 54 imprinted loci in 23 selected patients with different BEEC subtypes (epispadias n = 1, classical bladder exstrophy n = 10, exstrophy of the cloaca n = 12) using the Infinium HumanMethylation450 BeadChip. A total of 471,722 not imprinted autosomal CpG loci and 891 imprinted CpG loci were investigated. Findings were corroborated by methylation-specific-multiplex ligation-dependent probe amplification (MS-MLPA) and microsatellite analysis. RESULTS: No significant differences in the DNA methylation of the not imprinted and imprinted CpG were observed depending on subtype of BEEC. Nevertheless, in 1 of the 23 patients who displayed a classical bladder exstrophy, we detected hypomethylation of the imprinted PLAGL1 locus in chromosome 6q24. We verified this hypomethylation by MS-MLPA and showed further the methylation loss to be caused most likely by a mosaic epimutation. CONCLUSION: Considering that it is highly unlikely to detect a PLAGL1 epimutation among 23 individuals given the low incidence of this alteration in the population, our observations further support a link between BEEC and imprinting disorders. Birth Defects Research (Part A) 106:724-728, 2016. © 2016 Wiley Periodicals, Inc.

Genome-wide association study of serum iron phenotypes in premenopausal women of European descent.

A genome-wide association study was performed on 1130 premenopausal women to detect common variants associated with three serum iron-related phenotypes. Total iron binding capacity was strongly associated (p=10(-14)) with variants in and near the TF gene (transferrin), the serum iron transporting protein, and with variants in HFE (p=4×10(-7)), which encodes the human hemochromatosis gene. Association was also detected between percent iron saturation (p=10(-8)) and variants in the chromosome 6 region containing both HFE and SLC17A2, which encodes a phosphate transport protein. No significant associations were detected with serum iron, but variants in HFE were suggestive (p=10(-6)). Our results corroborate prior studies in older subjects and demonstrate that the association of these genetic variants with iron phenotypes can be detected in premenopausal women.

High-resolution characterization of CPD hotspot formation in human fibroblasts.

Repair of DNA lesions must occur within the chromatin landscape and is associated with alterations in histone modifications and nucleosome rearrangement. To directly associate these chromatin features with DNA damage and repair, it is necessary to be able to map DNA adducts. We have developed a cyclobutane pyrimidine dimer (CPD)-specific immunoprecipitation method and mapped ultraviolet damage hotspots across human chromosomes 1 and 6. CPD hotspots occur almost equally in genic and intergenic regions. However, these hotspots are significantly more prevalent adjacent to repeat elements, especially Alu repeats. Nucleosome mapping studies indicate that nucleosomes are consistently positioned at Alu elements where CPD hotspots form, but by 2 h post-irradiation, these same regions are significantly depleted of nucleosomes. These results indicate that nucleosomes associated with hotspots of CPD formation are readily rearranged, potentially making them accessible to DNA repair machinery. Our results represent the first chromosome scale map of ultraviolet-induced DNA lesions in the human genome, and reveal the sequence features and dynamic chromatin changes associated with CPD hotspots.

Association of genes on chromosome 6, GRIK2 , TMEM217 and TMEM63B (linked to MRPL14 ) with diabetic retinopathy.

PURPOSE: Diabetic retinopathy (DR) is one of the most common complications of diabetes mellitus (DM). The susceptibility genes responsible for increasing the risk for DR in type 2 diabetes (T2D) were sought in this study. METHODS: A case-control study was carried out, comprising 749 unrelated T2D individuals with (n = 174) and without (n = 575) DR. Genotypic distributions of single nucleotide polymorphisms (SNPs) were determined for subjects with and without DR. RESULTS: Eight chromosome 6 SNPs, having the most significant differences, were delineated: rs10499298, rs10499299, rs17827966, rs1224329, rs1150790, rs713050, rs2518344 and rs487083; all were associated with genes TMEM217, MRPL14 and GRIK2. After adjusting for the duration of DM and levels of hemoglobin A(1c), the TT genotype of rs713050, and the AG + AA genotypes of rs2518344 and rs10499298, differed significantly between those with and without DR. Haplotype analysis revealed haplotype C-A-C, residing in rs10499299, rs10499298 and rs17827966, to have significant linkage disequilibrium. CONCLUSIONS: We identified new loci on chromosome 6 associated to DR; all loci showed high levels of linkage disequilibrium.

Description of two new MICA alleles: MICA*058 and MICA*002:03.

We describe two novel alleles, MICA*058 and MICA* 002:03.

Infantile esotropia could be oligogenic and allelic with Duane retraction syndrome.

PURPOSE: To describe phenotyping and linkage analysis results for available members from a consanguineous nuclear family with hereditary congenital strabismus. METHODS: Both parents and all 12 children underwent clinical examination. Available affected and several unaffected family members had venous blood sampling for DNA extraction and 10K single nucleotide polymorphism (SNP) genotyping (Affymetrix Gene Chip® Human). Multipoint logarithm of the odds (LOD) score calculations were performed assuming an autosomal recessive mode of inheritance with 100% penetrance and disease allele frequency of 0.01%. RESULTS: Three children had non-syndromic large-angle infantile esotropia without significant hyperopia. A fourth child had left esotropic Duane retraction syndrome. A fifth child who had esotropia in the setting of prematurity and childhood poliomyelitis was excluded from the analysis. A sixth child had keratoconus and was excluded. Both parents and the remaining 6 children had no significant orthoptic or ophthalmic findings. Using linkage analysis including the 4 esotropic children, disease loci were mapped to regions on chromosomes 3p26.3-26.2 and 6q24.2-25.1 using multipoint linkage analysis with LOD scores of 3.18 and 3.25 respectively. Linkage to these regions persisted when the esotropic Duane retraction syndrome patient was excluded from the linkage analysis (LOD scores of 2.00 and 2.32, respectively). CONCLUSIONS: Non-syndromic infantile esotropia could be related to susceptibility loci on chromosomal regions 3p26.3-26.2 and 6q24.2-25.1 and may share alleles that underlie Duane retraction syndrome.

Multiple independent variants in 6q21-22 associated with susceptibility to celiac disease in the Dutch, Finnish and Hungarian populations.

Celiac disease is an inflammatory enteropathy caused by intolerance to gluten. Previous linkage studies in the Dutch, Finnish and Hungarian populations have revealed a locus on chromosome 6q21-22 conferring susceptibility to celiac disease. This locus has previously been implicated in susceptibility to other autoimmune diseases such as Crohn's disease and type 1 diabetes. We performed fine mapping on 446 independent individuals with celiac disease and 641 controls of Dutch origin, testing 872 tagging SNPs in a 22 Mb region of chromosome 6. The 12 most promising SNPs were followed up in 2071 individuals from 284 Finnish and 357 Hungarian celiac disease families to identify risk variants in this region. Multiple markers in the region were significantly associated with celiac disease in the Dutch material. Two SNPs, rs9391227 and rs4946111, were significantly associated with celiac disease in the Finnish population. The association to rs9391227 represents the strongest association signal found in the Finnish (P = 0.003, OR 0.66) as well as the combined Dutch, Finnish and Hungarian populations (P = 3.6 × 10(-5), OR 0.76). The rs9391227 is situated downstream of the HECT domain and ankyrin repeat containing, E3 ubiquitin protein ligase 1 (HACE1) gene and is contained within a region of strong linkage disequilibrium enclosing HACE1. Two additional, independent, susceptibility variants in the 6q21-22 region were also found in a meta-analysis of the three populations. The 6q21-22 region was confirmed as a celiac disease susceptibility locus and harbors multiple independent associations, some of which may implicate ubiquitin-pathways in celiac disease susceptibility.

The HBS1L-MYB intergenic region on chromosome 6q23.3 influences erythrocyte, platelet, and monocyte counts in humans.

Common sequence variants situated between the HBS1L and MYB genes on chromosome 6q23.3 (HMIP) influence the proportion of F cells (erythrocytes that carry measurable amounts of fetal hemoglobin). Since the physiological processes underlying the F-cell variability are thought to be linked to kinetics of erythrocyte maturation and differentiation, we have investigated the influence of the HMIP locus on other hematologic parameters. Here we show a significant impact of HMIP variability on several types of peripheral blood cells: erythrocyte, platelet, and monocyte counts as well as erythrocyte volume and hemoglobin content in healthy individuals of European ancestry. These results support the notion that changes of F-cell abundance can be an indicator of more general shifts in hematopoietic patterns in humans.

Generic features of tertiary chromatin structure as detected in natural chromosomes.

Knowledge of tertiary chromatin structure in mammalian interphase chromosomes is largely derived from artificial tandem arrays. In these model systems, light microscope images reveal fibers or beaded fibers after high-density targeting of transactivators to insertional domains spanning several megabases. These images of fibers have lent support to chromonema fiber models of tertiary structure. To assess the relevance of these studies to natural mammalian chromatin, we identified two different approximately 400-kb regions on human chromosomes 6 and 22 and then examined light microscope images of interphase tertiary chromatin structure when the regions were transcriptionally active and inactive. When transcriptionally active, these natural chromosomal regions elongated, yielding images characterized by a series of adjacent puncta or "beads", referred to hereafter as beaded images. These elongated structures required transcription for their maintenance. Thus, despite marked differences in the density and the mode of transactivation, the natural and artificial systems showed similarities, suggesting that beaded images are generic features of transcriptionally active tertiary chromatin. We show here, however, that these images do not necessarily favor chromonema fiber models but can also be explained by a radial-loop model or even a simple nucleosome affinity, random-chain model. Thus, light microscope images of tertiary structure cannot distinguish among competing models, although they do impose key constraints: chromatin must be clustered to yield beaded images and then packaged within each cluster to enable decondensation into adjacent clusters.

Multiple excitation confocal analysis of targets in nuclei of cytogenetic preparations.

OBJECTIVE: To demonstrate that cellular preparations requiring color analysis of different domains stained by molecular cytogenetic methods (fluorescence in situ hybridization) can be processed by spectral analysis of fluorescent emissions by either factor analysis of medical image sequences (FAMIS) or a META confocal configuration to isolate fluorescent probes. STUDY DESIGN: Three-dimensional sequences of images obtained by spectral analysis in a META confocal microscope (Carl Zeiss SAS, Jena, Germany) were analyzed by META processing and the FAMIS algorithm, which provides factor curves. META and factor images were then the result of image-processing methods that cover emission spectra. RESULTS: Factor curves and factor or META images can help to analyze targets inside nuclei. CONCLUSION: It is possible to process preparations containing numerous spots on different colors to differentiate stained targets and to improve visualization and detection.

Skew of mononucleotide frequencies, relative abundance of dinucleotides, and DNA strand asymmetry.

Based on 152 mitochondrial genomes and 36 bacterial chromosomes that have been completely sequenced, as well as three long contigs for human chromosomes 6, 21, and 22, we examined skews of mononucleotide frequencies and the relative abundance of dinucleotides in one DNA strand. Each group of these genomes has its own characteristics. Regarding mitochondrial genomes, both CpG and GpT are underrepresented, while either GpG or CpC or both are overrepresented. The relative frequency of nucleotide T vs A and of nucleotide G vs C is strongly skewed, due presumably to strand asymmetry in replication errors and unidirectional DNA replication from single origins. Exceptions are found in the plant and yeast mitochondrial genomes, each of which may replicate from multiple origins. Regarding bacterial genomes, the "universal" rule of CpG deficiency is restricted to archaebacteria and some eubacteria. In other eubacteria, the most underrepresented dinucleotide is either TpA or GpT. In general, there are significant T vs A and G vs C skews in each half of the bacterial genome, although these are almost exactly canceled out over the whole genome. Regarding human chromosomes 6, 21, and 22, dinucleotide CpG tends to be avoided. The relative frequency of mononucleotides exhibits conspicuous local skews, suggesting that each of these chromosomal segments contains more than one DNA replication origin. It is concluded that, when there are several replicons in a genomic region, not only the number of DNA replication origins but also the directionality is important and that the observed patterns of nucleotide frequencies in the genome strongly support the hypothesis of strand asymmetry in replication errors.

Construction of a high-resolution 2.5-Mb transcript map of the human 6p21.2-6p21.3 region immediately centromeric of the major histocompatibility complex.

We have constructed a 2.5-Mb physical and transcription map that spans the human 6p21.2-6p21.3 region and includes the centromeric end of the MHC, using a combination of techniques. In total 88 transcription units including exons, cDNAs, and cDNA contigs were characterized and 60 were confidently positioned on the physical map. These include a number of genes encoding nuclear and splicing factors (Ndr kinase, HSU09564, HSRP20); cell cycle, DNA packaging, and apoptosis related [p21, HMGI(Y), BAK]; immune response (CSBP, SAPK4); transcription activators and zinc finger-containing genes (TEF-5, ZNF76); embryogenesis related (Csa-19); cell signaling (DIPP); structural (HSET), and other genes (TULP1, HSPRARD, DEF-6, EO6811, cyclophilin), as well as a number of RP genes and pseudogenes (RPS10, RPS12-like, RPL12-like, RPL35-like). Furthermore, several novel genes (a Br140-like, a G2S-like, a FBN2-like, a ZNF-like, and B1/KIAA0229) have been identified, as well as cDNAs and cDNA contigs. The detailed map of the gene content of this chromosomal segment provides a number of candidate genes, which may be involved in several biological processes that have been associated with this region, such as spermatogenesis, development, embryogenesis, and neoplasia. The data provide useful tools for synteny studies between mice and humans, for genome structure analysis, gene density comparisons, and studies of nucleotide composition, of different isochores and Giemsa light and Giemsa dark bands.

Direct mapping of the human TATA box-binding protein (TBP) gene to 6q27 by fluorescence in situ hybridization.

TBP (TATA box-binding protein) participates in the expression of eukaryotic genes transcribed by RNA polymerases I, II, and III. Molecular cloning of human TBP revealed that the N-terminal region contains a polymorphic (CAG)n repeat. We report here the direct localization of human TBP gene to chromosome 6q2705-->qter region by fluorescence in situ hybridization, using the cDNA clone with or without the (CAG)n repeat as a probe.

Fine-scale deletion mapping of the distal long arm of chromosome 6 in 70 human ovarian cancers.

To define a small region on chromosome 6q containing a putative tumor suppressor gene for ovarian cancer, we examined loss of heterozygosity in 70 ovarian tumors of three histological types with nine restriction fragment length polymorphism markers located at 6q24-27. Among 33 cancers of serous type that were informative at one or more loci, 17 showed allelic loss at a few or all loci examined, whereas only 1 of 15 mucinous-type tumors and 2 of 12 clear-cell tumors revealed loss of heterozygosity. This result supported our earlier suggestion that alteration of a gene on chromosome 6q may play an important role during development of serous ovarian tumors (Sato et al., Cancer Res., 51: 5118-5122, 1991). Frequent losses were observed between loci defined by CI6-119 (D6S195) at 6q26 and CI6-49 (D6S161) at 6q27. A detailed deletion map indicated a commonly deleted region between loci defined by CI6-111 (D6S193) and CI6-24 (D6S149); these two markers are estimated to be 1.9 cM apart on the basis of linkage analysis. Our results further define a region containing a tumor suppressor gene involved in ovarian carcinoma within an approximately 2-megabase-long segment of chromosome 6q.