The IL-4Rα Q576R polymorphism is associated with increased severity of atopic dermatitis and exaggerates allergic skin inflammation in mice

BACKGROUND

Atopic dermatitis (AD) is characterized by TH2-dominated skin inflammation and systemic response to cutaneously encountered antigens. The TH2 cytokines IL-4 and IL-13 play a critical role in the pathogenesis of AD. The Q576->R576 polymorphism in the IL-4 receptor alpha (IL-4Rα) chain common to IL-4 and IL-13 receptors alters IL-4 signaling and is associated with asthma severity.

OBJECTIVE

We sought to investigate whether the IL-4Rα R576 polymorphism is associated with AD severity and exaggerates allergic skin inflammation in mice.

METHODS

Nighttime itching interfering with sleep, Rajka-Langeland, and Eczema Area and Severity Index scores were used to assess AD severity. Allergic skin inflammation following epicutaneous sensitization of mice 1 or 2 IL-4Rα R576 alleles (QR and RR) and IL-4Rα Q576 (QQ) controls was assessed by flow cytometric analysis of cells and quantitative RT-PCR analysis of cytokines in skin.

RESULTS

The frequency of nighttime itching in 190 asthmatic inner-city children with AD, as well as Rajka-Langeland and Eczema Area and Severity Index scores in 1116 White patients with AD enrolled in the Atopic Dermatitis Research Network, was higher in subjects with the IL-4Rα R576 polymorphism compared with those without, with statistical significance for the Rajka-Langeland score. Following epicutaneous sensitization of mice with ovalbumin or house dust mite, skin infiltration by CD4+ cells and eosinophils, cutaneous expression of Il4 and Il13, transepidermal water loss, antigen-specific IgE antibody levels, and IL-13 secretion by antigen-stimulated splenocytes were significantly higher in RR and QR mice compared with QQ controls. Bone marrow radiation chimeras demonstrated that both hematopoietic cells and stromal cells contribute to the mutants' exaggerated allergic skin inflammation.

CONCLUSIONS

The IL-4Rα R576 polymorphism predisposes to more severe AD and increases allergic skin inflammation in mice.

Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

Hunter-gatherer genomes reveal diverse demographic trajectories during the rise of farming in Eastern Africa

The fate of hunting and gathering populations following the rise of agriculture and pastoralism remains a topic of debate in the study of human prehistory. Studies of ancient and modern genomes have found that autochthonous groups were largely replaced by expanding farmer populations with varying levels of gene flow, a characterization that is influenced by the almost universal focus on the European Neolithic.^1-5 We sought to understand the demographic impact of an ongoing cultural transition to farming in Southwest Ethiopia, one of the last regions in Africa to experience such shifts.⁶ Importantly, Southwest Ethiopia is home to several of the world's remaining hunter-gatherer groups, including the Chabu people, who are currently transitioning away from their traditional mode of subsistence.⁷ We generated genome-wide data from the Chabu and four neighboring populations, the Majang, Shekkacho, Bench, and Sheko, to characterize their genetic ancestry and estimate their effective population sizes over the last 60 generations. We show that the Chabu are a distinct population closely related to ancient people who occupied Southwest Ethiopia >4,500 years ago. Furthermore, the Chabu are undergoing a severe population bottleneck, which began approximately 1,400 years ago. By analyzing eleven Eastern African populations, we find evidence for divergent demographic trajectories among hunter-gatherer-descendant groups. Our results illustrate that although foragers respond to encroaching agriculture and pastoralism with multiple strategies, including cultural adoption of agropastoralism, gene flow, and economic specialization, they often face population decline.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Genome-wide association study of asthma, total IgE, and lung function in a cohort of Peruvian children

BACKGROUND

Genetic ancestry plays a role in asthma health disparities.

OBJECTIVE

Our aim was to evaluate the impact of ancestry on and identify genetic variants associated with asthma, total serum IgE level, and lung function.

METHODS

A total of 436 Peruvian children (aged 9-19 years) with asthma and 291 without asthma were genotyped by using the Illumina Multi-Ethnic Global Array. Genome-wide proportions of indigenous ancestry populations from continental America (NAT) and European ancestry from the Iberian populations in Spain (IBS) were estimated by using ADMIXTURE. We assessed the relationship between ancestry and the phenotypes and performed a genome-wide association study.

RESULTS

The mean ancestry proportions were 84.7% NAT (case patients, 84.2%; controls, 85.4%) and 15.3% IBS (15.8%; 14.6%). With adjustment for asthma, NAT was associated with higher total serum IgE levels (P < .001) and IBS was associated with lower total serum IgE levels (P < .001). NAT was associated with higher FEV1 percent predicted values (P < .001), whereas IBS was associated with lower FEV1 values in the controls but not in the case patients. The HLA-DR/DQ region on chromosome 6 (Chr6) was strongly associated with total serum IgE (rs3135348; P = 3.438 × 10-10) and was independent of an association with the haplotype HLA-DQA1∼HLA-DQB1:04.01∼04.02 (P = 1.55 × 10-05). For lung function, we identified a locus (rs4410198; P = 5.536 × 10-11) mapping to Chr19, near a cluster of zinc finger interacting genes that colocalizes to the long noncoding RNA CTD-2537I9.5. This novel locus was replicated in an independent sample of pediatric case patients with asthma with similar admixture from Brazil (P = .005).

CONCLUSION

This study confirms the role of HLA in atopy, and identifies a novel locus mapping to a long noncoding RNA for lung function that may be specific to children with NAT.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma

BACKGROUND

Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent.

METHODS

We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 μg in children, 100 μg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV₁. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 μg twice daily in children and 500 μg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 μg twice daily in children, 250 μg twice daily in adolescents and adults), and 5 × ICS versus 100 μg fluticasone plus a LABA (salmeterol 50 μg twice daily). We used a genome-wide significance threshold of p<1·6 × 10^-4, and tested for replication using independent cohorts of individuals of African descent with asthma.

FINDINGS

We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 μg fluticasone plus salmeterol on chromosome 12 (odds ratio [OR_{local African}] 3·95, 95% CI 2·02-7·72, p=6·1 × 10^-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, OR_{allele dose} 0·17, 95% CI 0·07-0·42, p=8·4 × 10^-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (OR_{local African} 3·35, 1·98-5·67, p=6·8 × 10^-6) containing a locus adjacent to TPST2 (rs5752429, OR_{allele dose} 0·21, 0·09-0·52, p=5·7 × 10^-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts.

INTERPRETATION

BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma.

FUNDING

National Institutes of Health, National Heart, Lung, and Blood Institute.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Whole genome sequencing identifies novel genetic mutations in patients with eczema herpeticum

BACKGROUND

Eczema herpeticum (EH) is a rare complication of atopic dermatitis (AD) caused by disseminated herpes simplex virus (HSV) infection. The role of rare and/or deleterious genetic variants in disease etiology is largely unknown. This study aimed to identify genes that harbor damaging genetic variants associated with HSV infection in AD with a history of recurrent eczema herpeticum (ADEH+).

METHODS

Whole genome sequencing (WGS) was performed on 49 recurrent ADEH+ (≥3 EH episodes), 491 AD without a history of eczema herpeticum (ADEH-) and 237 non-atopic control (NA) subjects. Variants were annotated, and a gene-based approach (SKAT-O) was used to identify genes harboring damaging genetic variants associated with ADEH+. Genes identified through WGS were studied for effects on HSV responses and keratinocyte differentiation.

RESULTS

Eight genes were identified in the comparison of recurrent ADEH+to ADEH-and NA subjects: SIDT2, CLEC7A, GSTZ1, TPSG1, SP110, RBBP8NL, TRIM15, and FRMD3. Silencing SIDT2 and RBBP8NL in normal human primary keratinocytes (NHPKs) led to significantly increased HSV-1 replication. SIDT2-silenced NHPKs had decreased gene expression of IFNk and IL1b in response to HSV-1 infection. RBBP8NL-silenced NHPKs had decreased gene expression of IFNk, but increased IL1b. Additionally, silencing SIDT2 and RBBP8NL also inhibited gene expression of keratinocyte differentiation markers keratin 10 (KRT10) and loricrin (LOR).

CONCLUSION

SIDT2 and RBBP8NL participate in keratinocyte's response to HSV-1 infection. SIDT2 and RBBP8NL also regulate expression of keratinocyte differentiation genes of KRT10 and LOR.

Zika Virus Congenital Syndrome and MTOR gene variants: insights from a family of dizygotic twins

Congenital Zika virus syndrome (CZS) is associated with damage to neural progenitor cells by ZIKA virus infection. There are no accurate statistics on the percentage of pregnant mothers who have had babies affected by the syndrome. Few cases of discordant twins have been described in the literature and, therefore, we hypothesize that the genetic background of the progeny and/or mother may play a role in the fate of the syndrome. We performed a complete exome sequencing in a set of dizygotic individuals and their parents. After that, we selected discordant variants on the MTOR gene between the affected and unaffected twin and we observed a mutation (rs2295079), placed in a region restricted to proximal 5'-UTR, as a strong possible causal variant. In addition, in most brain tissues (including fetal brain) evaluated for expression quantitative trait loci (eQTL), this locus is strongly correlated with post-translational modifications of histones (promoter and enhancer marks) and hypersensitivity to DNAse I (open chromatin mark). Taken together, our data suggest that changes in the MTOR gene may be related to CZS. Additional functional studies should be carried out to prove how and why a MTOR mutation can predispose the fetus to the syndrome.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.¹ Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations leading to clonal expansion was associated with both hematologic cancer ^2–4 and coronary heart disease⁵, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).⁶ Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer resulting in increased hematopoietic stem cell self-renewal. Overall, we observe that germline genetic variation shapes hematopoietic stem cell function leading to CHIP through mechanisms that are both specific to clonal hematopoiesis and shared mechanisms leading to somatic mutations across tissues.

Whole genome sequence analysis of pulmonary function and COPD in 19,996 multi-ethnic participants

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry.

Easy-HLA: a validated web application suite to reveal the full details of HLA typing

Motivation

The HLA system plays a pivotal role in both clinical applications and immunology research. Typing HLA genes in patient and donor is indeed required in hematopoietic stem cell and solid-organ transplantation, and the histocompatibility complex region exhibits countless genetic associations with immune-related pathologies. Since the discovery of HLA antigens, the HLA system nomenclature and typing methods have constantly evolved, which leads to difficulties in using data generated with older methodologies.

Results

Here, we present Easy-HLA, a web-based software suite designed to facilitate analysis and gain knowledge from HLA typing, regardless of nomenclature or typing method. Easy-HLA implements a computational and statistical method of HLA haplotypes inference based on published reference populations containing over 600 000 haplotypes to upgrade missing or partial HLA information: ‘HLA-Upgrade’ tool infers high-resolution HLA typing and ‘HLA-2-Haplo’ imputes haplotype pairs and provides additional functional annotations (e.g. amino acids and KIR ligands). We validated both tools using two independent cohorts (total n = 2500). For HLA-Upgrade, we reached a prediction accuracy of 92% from low- to high-resolution of European genotypes. We observed a 96% call rate and 76% accuracy with HLA-2-Haplo European haplotype pairs prediction. In conclusion, Easy-HLA tools facilitate large-scale immunogenetic analysis and promotes the multi-faceted HLA expertise beyond allelic associations by providing new functional immunogenomics parameters.

Availability and implementation

Easy-HLA is a web application freely available (free account) at: https://hla.univ-nantes.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

A regulatory variant in the C1Q gene cluster is associated with tuberculosis susceptibility and C1qA plasma levels in a South African population

Several genetic studies have implicated genes that encode for components of the innate immune response in tuberculosis (TB) susceptibility. The complement system is an early player in the innate immune response and provides the host with initial protection by promoting phagocytosis of apoptotic or necrotic cells. The C1q molecule is the first component of the classical pathway that leads to the activation of complement by binding to immune complexes and is encoded by the C1Q gene cluster. We investigated variants in this region to determine its association with TB susceptibility. Five single nucleotide polymorphisms (SNPs) (rs12033074, rs631090, rs172378, rs587585, and rs665691) were genotyped using TaqMan® SNP assays in 456 TB cases and 448 healthy controls and analysed by logistic regression models. The rs587585 variant showed a significant additive allelic association where the minor G allele was found more frequently in TB cases than in controls in both the discovery (p = 0.023; OR = 1.30; 95% CI, 1.04-1.64) and validation cohort (p = 0.038; OR = 1.31; 95% CI, 1.22-1.40). In addition, we detected increased C1qA expression when comparing cases and controls (p = 0.037) and linked this to a dosage effect of the G allele, which increased C1qA expression in TB cases. This is the first study to report the association of C1Q gene polymorphisms with progression to tuberculosis.

clustifyr: an R package for automated single-cell RNA sequencing cluster classification

Assignment of cell types from single-cell RNA sequencing (scRNA-seq) data remains a time-consuming and error-prone process. Current packages for identity assignment use limited types of reference data and often have rigid data structure requirements. We developed the clustifyr R package to leverage several external data types, including gene expression profiles to assign likely cell types using data from scRNA-seq, bulk RNA-seq, microarray expression data, or signature gene lists. We benchmark various parameters of a correlation-based approach and implement gene list enrichment methods. clustifyr is a lightweight and effective cell-type assignment tool developed for compatibility with various scRNA-seq analysis workflows. clustifyr is publicly available at https://github.com/rnabioco/clustifyr

clustifyr: an R package for automated single-cell RNA sequencing cluster classification

Assignment of cell types from single-cell RNA sequencing (scRNA-seq) data remains a time-consuming and error-prone process. Current packages for identity assignment use limited types of reference data and often have rigid data structure requirements. We developed the clustifyr R package to leverage several external data types, including gene expression profiles to assign likely cell types using data from scRNA-seq, bulk RNA-seq, microarray expression data, or signature gene lists. We benchmark various parameters of a correlation-based approach and implement gene list enrichment methods. clustifyr is a lightweight and effective cell-type assignment tool developed for compatibility with various scRNA-seq analysis workflows. clustifyr is publicly available at https://github.com/rnabioco/clustifyr.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Whole Genome Sequencing Identifies CRISPLD2 as a Lung Function Gene in Children With Asthma

BACKGROUND

Asthma is a common respiratory disorder with a highly heterogeneous nature that remains poorly understood. The objective was to use whole genome sequencing (WGS) data to identify regions of common genetic variation contributing to lung function in individuals with a diagnosis of asthma.

METHODS

WGS data were generated for 1,053 individuals from trios and extended pedigrees participating in the family-based Genetic Epidemiology of Asthma in Costa Rica study. Asthma affection status was defined through a physician's diagnosis of asthma, and most participants with asthma also had airway hyperresponsiveness (AHR) to methacholine. Family-based association tests for single variants were performed to assess the associations with lung function phenotypes.

RESULTS

A genome-wide significant association was identified between baseline FEV1/FVC ratio and a single-nucleotide polymorphism in the top hit cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) (rs12051168; P = 3.6 × 10-8 in the unadjusted model) that retained suggestive significance in the covariate-adjusted model (P = 5.6 × 10-6). Rs12051168 was also nominally associated with other related phenotypes: baseline FEV1 (P = 3.3 × 10-3), postbronchodilator (PB) FEV1 (7.3 × 10-3), and PB FEV1/FVC ratio (P = 2.7 × 10-3). The identified baseline FEV1/FVC ratio and rs12051168 association was meta-analyzed and replicated in three independent cohorts in which most participants with asthma also had confirmed AHR (combined weighted z-score P = .015) but not in cohorts without information about AHR.

CONCLUSIONS

These findings suggest that using specific asthma characteristics, such as AHR, can help identify more genetically homogeneous asthma subgroups with genotype-phenotype associations that may not be observed in all children with asthma. CRISPLD2 also may be important for baseline lung function in individuals with asthma who also may have AHR.

African American ancestry contribution to asthma and atopic dermatitis

OBJECTIVE

Asthma and atopic dermatitis (AD) are complex diseases with striking disparities across racial and ethnic groups, which may be partly attributable to genetic factors. Here we summarize current knowledge from asthma and AD genome-wide association studies (GWAS) and pharmacogenetic studies in African ancestry populations.

DATA SOURCES

GWAS catalog; PUBMed.

STUDY SELECTIONS

GWAS catalog studies with trait annotations "asthma" and "atopic eczema" and African ancestry individuals in the discovery dataset; the recent CAAPA asthma GWAS; reports on pharmacogenetic studies in asthma and AD.

RESULTS

Although GWASs have revolutionized gene discovery for multiple complex traits, African Americans continue to be severely underrepresented in sufficiently powered genetics studies. Indeed, of the 16 asthma and 21 AD loci that reached genomewide significance in Europeans, very few have replicated in African ancestry populations. Challenges in comparing results from European vs African ancestry cohorts include modest sample size, differences in risk allele frequency, effect size, correlation between genetic variants, and environmental exposure in evolutionary history. African Americans also constitute a small percentage of dermatological and respiratory-focused clinical trials. Pharmacogenetic studies have similarly been focused largely on non-Hispanic whites, despite compelling evidence that genetic variation from different ancestral backgrounds may alter therapeutic efficacy of asthma and AD drugs.

CONCLUSION

Large-scale genetic studies of asthma and AD in African Americans are essential to reduce research and health disparities and empower scientific discoveries.

Author Correction: Assembly of a pan-genome from deep sequencing of 910 humans of African descent

In the version of this article initially published, the statement "there are no pan-genomes for any other animal or plant species" was incorrect. The statement has been corrected to "there are no reported pan-genomes for any other animal species, to our knowledge." We thank David Edwards for bringing this error to our attention. The error has been corrected in the HTML and PDF versions of the article.

Surgery for giant tumours of the breast: a 15 year review

BACKGROUND

Giant tumours of the breast tend to occur in the adolescent age group. Racial predilection has been noted in the literature. The mass often occupies most of the breast, leading to its distortion. Many authors have advocated a mastectomy for benign tumours that severely distort the breast. Giant benign tumours when treated by simple excision risk persistence of asymmetry. To avoid this asymmetry, some authors have resorted to excision and immediate reduction mammaplasty. The aim of this retrospective study was a report on giant tumours of the breast presenting to a plastic surgery unit and to analyse demographic factors, clinical presentations, tumour pathology, management, complications, as well as patient and breast outcomes.

METHODS

Medical records of patients with giant tumours were retrospectively analysed for assessing demographic factors, clinical presentation, tumour pathology, the technique of surgery performed and patient and breast outcomes in a single hospital setting. Breast outcomes were rated by panel of 4 experienced plastic surgeons using the 4 Point Likert scale. Their ratings were statistically analysed for inter-rater agreement.

RESULTS

Twenty-three subjects were identified to have giant tumours of the breast. Of these South African patients, 19 were black, 3 were Indian and 1 was of mixed ethnicity. The age range was 12-49 years (y) with an average of 19y. All masses were palpable. The final pathological diagnosis was fibrocystic disease in 3, giant fibroadenoma in 14, phyllodes tumour in 4, and hamartoma in 2. The size range was 10-45 cm with a median size of 18 cm. All but one patient had simple excision followed by immediate reduction mammaplasty. Twenty patients were assessed after operation. A minimum of 1 to a maximum of 4 patients per reviewer showed unsatisfactory outcomes and a minimum of 18 to a maximum of 21 patients per reviewer showed satisfactory to excellent outcomes. The overall agreement between assessors for this was 84%.

CONCLUSION

Benign giant tumours (> 10 cm) of the breast are suitably managed by excision of the mass and a reduction mammoplasty technique of reconstruction.

A Functional Toll-Interacting Protein Variant Is Associated with Bacillus Calmette-Guérin-Specific Immune Responses and Tuberculosis

RATIONALE

The molecular mechanisms that regulate tuberculosis susceptibility and bacillus Calmette-Guérin (BCG)-induced immunity are mostly unknown. However, induction of the adaptive immune response is a critical step in host control of Mycobacterium tuberculosis. Toll-interacting protein (TOLLIP) is a ubiquitin-binding protein that regulates innate immune responses, including Toll-like receptor signaling, which initiate adaptive immunity. TOLLIP variation is associated with susceptibility to tuberculosis, but the mechanism by which it regulates tuberculosis immunity is poorly understood.

OBJECTIVES

To identify functional TOLLIP variants and evaluate the role of TOLLIP variation on innate and adaptive immune responses to mycobacteria and susceptibility to tuberculosis.

METHODS

We used human cellular immunology approaches to characterize the role of a functional TOLLIP variant on monocyte mRNA expression and M. tuberculosis-induced monocyte immune functions. We also examined the association of TOLLIP variation with BCG-induced T-cell responses and susceptibility to latent tuberculosis infection.

MEASUREMENTS AND MAIN RESULTS

We identified a functional TOLLIP promoter region single-nucleotide polymorphism, rs5743854, which was associated with decreased TOLLIP mRNA expression in infant monocytes. After M. tuberculosis infection, TOLLIP-deficient monocytes demonstrated increased IL-6, increased nitrite, and decreased bacterial replication. The TOLLIP-deficiency G/G genotype was associated with decreased BCG-specific IL-2⁺ CD4⁺ T-cell frequency and proliferation. This genotype was also associated with increased susceptibility to latent tuberculosis infection.

CONCLUSIONS

TOLLIP deficiency is associated with decreased BCG-specific T-cell responses and increased susceptibility to tuberculosis. We hypothesize that the heightened antibacterial monocyte responses after vaccination of TOLLIP-deficient infants are responsible for decreased BCG-specific T-cell responses. Activating TOLLIP may provide a novel adjuvant strategy for BCG vaccination.

The role of ST2 and ST2 genetic variants in schistosomiasis

BACKGROUND

Chronic schistosomiasis and its severe complication, periportal fibrosis, are characterized by a predominant T_h2 response. To date, specific single nucleotide polymorphisms in ST2 have been some of the most consistently associated genetic variants for asthma.

OBJECTIVE

We investigated the role of ST2 (a receptor for the T_h2 cytokine IL-33) in chronic and late-stage schistosomiasis caused by Schistosoma japonicum and the potential effect of ST2 genetic variants on stage of disease and ST2 expression.

METHODS

We recruited 947 adult participants (339 with end-stage schistosomiasis and liver cirrhosis, 307 with chronic infections without liver fibrosis, and 301 health controls) from a S japonicum-endemic area (Hubei, China). Six ST2 single nucleotide polymorphisms were genotyped. Serum soluble ST2 (sST2) was measured by ELISA, and ST2 expression in normal liver tissues, Hepatitis B virus-induced fibrotic liver tissues, and S japonicum-induced fibrotic liver tissues was measured by immunohistochemistry.

RESULTS

We found sST2 levels were significantly higher in the end-stage group (36.04 [95% CI, 33.85-38.37]) compared with chronic cases and controls (22.7 [95% CI, 22.0-23.4], P < 1E-10). In addition, S japonicum-induced fibrotic liver tissues showed increased ST2 staining compared with normal liver tissues (P = .0001). Markers rs12712135, rs1420101, and rs6543119 were strongly associated with sST2 levels (P = 2E-10, 5E-05, and 6E-05, respectively), and these results were replicated in an independent cohort from Brazil living in a S mansoni endemic region.

CONCLUSIONS

We demonstrate for the first time that end-stage schistosomiasis is associated with elevated sST2 levels and show that ST2 genetic variants are associated with sST2 levels in patients with schistosomiasis.

TOLLIP Deficiency is Associated with Increased Mycobacterium tuberculosis-specific Anti-Microbial Monocyte Responses and Protection from Pediatric Tuberculosis in South Africa

Pediatric tuberculosis (TB) is a significant cause of death in children and has a different pathogenesis than adult TB, but BCG-induced immunity is incomplete. Although genetic factors regulate BCG-induced immunity and TB susceptibility, these genes are mostly unknown and have largely been unexamined in pediatric TB. We hypothesized that TOLLIP, a TLR suppressor, orchestrates immune responses to M. tuberculosis (MTb) and susceptibility to pediatric TB. We discovered that a polymorphism of TOLLIP, rs5743854, regulated signaling in a promoter assay and was associated with TOLLIP mRNA deficiency in infant peripheral blood monocytes (p = 0.008, Kruskal-Wallis). The rs5743854 G allele was associated with protection from pediatric TB in a South African cohort (N=226 cases and 626 controls, dominant model, p = 0.002, odds ratio (OR) 0.63 unadjusted analysis; p = 0.011, OR = 0.58 adjusted for ethnicity by ancestry informative markers). The G allele was also associated with decreased frequency of BCG-specific CD4+IL-2 responses (p = 0.003, Kruskal-Wallis test) and CD4+ T-cell proliferation (p = 0.03, generalized linear model) in 174 infants vaccinated with BCG at birth and re-stimulated at 10 weeks of age. To examine mechanisms of protection, we created a TOLLIP-deficient THP-1 cell line. TOLLIP-deficient monocytes had increased MTb-induced TNF and reactive nitrogen secretion and reduced MTb replication compared to control cells. Together, these data suggest that TOLLIP deficiency is associated with enhanced anti-microbial monocyte responses that control MTb replication and provide protection against pediatric TB disease. Suppressing TOLLIP levels by host-directed therapeutics may provide adjunctive therapy for pediatric TB.

Activating KIRs alter susceptibility to pulmonary tuberculosis in a South African population

We investigate the role of killer immunoglobulin-like receptor (KIR) genes and human leukocyte antigen class-I (HLA) variants in susceptibility to tuberculosis in a South African population. In a sample set comprising 408 TB cases and 351 healthy controls, we show that the KIR3DS1 gene and KIR genotypes with five or more activating KIRs, and the presence of 3DS1, protect against developing active TB in the South African Coloured population. Several HLA class-I alleles were identified as susceptibility factors for TB disease. However, none of the KIR-HLA compound genotypes were found to be associated with TB. Our data suggests that the KIR genes may play an important role in TB disease.

TLR1, 2, 4, 6 and 9 Variants Associated with Tuberculosis Susceptibility: A Systematic Review and Meta-Analysis

BACKGROUND

Studies investigating the influence of toll-like receptor (TLR) polymorphisms and tuberculosis susceptibility have yielded varying and often contradictory results in different ethnic groups. A meta-analysis was conducted to investigate the relationship between TLR variants and susceptibility to tuberculosis, both across and within specific ethnic groups.

METHODS

An extensive database search was performed for studies investigating the relationship between TLR and tuberculosis (TB) susceptibility. Data was subsequently extracted from included studies and statistically analysed.

RESULTS

32 articles involving 18907 individuals were included in this meta-analysis, and data was extracted for 14 TLR polymorphisms. Various genetic models were employed. An increased risk of TB was found for individuals with the TLR2 rs3804100 CC and the TLR9 rs352139 GA and GG genotypes, while decreased risk was identified for those with the AG genotype of TLR1 rs4833095. The T allele of TLR6 rs5743810 conferred protection across all ethnic groups. TLR2 rs5743708 subgroup analysis identified the A allele to increase susceptibility to TB in the Asian ethnic group, while conferring protection in the Hispanic group. The T allele of TLR4 rs4986791 was also found to increase the risk of TB in the Asian subgroup. All other TLR gene variants investigated were not found to be associated with TB in this meta-analysis.

DISCUSSION

Although general associations were identified, most TLR variants showed no significant association with TB, indicating that additional studies investigating a wider range of pattern recognition receptors is required to gain a better understanding of this complex disease.

Association of toll-like receptors with susceptibility to tuberculosis suggests sex-specific effects of TLR8 polymorphisms

BACKGROUND

Toll-like receptors (TLRs) are involved in the recognition of conserved microbial structures, leading to activation of an inflammatory response and formation of an adaptive immune response.

METHODS

Twenty-three polymorphisms in five TLR genes were genotyped in 729 tuberculosis cases and 487 healthy controls in a population-based case-control association study in a South African population.

RESULTS

We detected sex-specific associations for TLR8 polymorphisms, with rs3761624 (OR=1.54, p<0.001), rs3764879 (OR=1.41, p=0.011) and rs3764880 (OR=1.42, p=0.011) associated in females and rs3764879 (OR=0.72, p=0.013) and rs3764880 (OR=0.75, p=0.036) associated in males. Epistatic interactions between the TLR genes were investigated and the TLR1_rs4833095 polymorphism was shown to interact with TLR2_rs3804100 and (GT)n microsatellite (p=0.002) and alter susceptibility to TB. We also studied the role of TLRs in disease caused by different Mycobacterium tuberculosis genotypes in 257 tuberculosis cases, and identified associations between specific TLR polymorphisms and disease caused by specific strains.

CONCLUSION

This study provides further evidence that the TLRs play an important role in the outcome of tuberculosis disease, and suggests a partial explanation for the male bias in tuberculosis ratios.

Investigating the Role of Gene-Gene Interactions in TB Susceptibility

Tuberculosis (TB) is the second leading cause of mortality from infectious disease worldwide. One of the factors involved in developing disease is the genetics of the host, yet the field of TB susceptibility genetics has not yielded the answers that were expected. A commonly posited explanation for the missing heritability of complex disease is gene-gene interactions, also referred to as epistasis. In this study we investigate the role of gene-gene interactions in genetic susceptibility to TB using a cohort recruited from a high TB incidence community from Cape Town, South Africa. Our discovery data set incorporates genotypes from a large a number of candidate gene studies as well as genome-wide data. After limiting our search space to pairs of putative TB susceptibility genes, as well as pairs of genes that have been curated in online databases as potential interactors, we use statistical modelling to identify pairs of interacting SNPs. We attempt to validate the top models identified in our discovery data set using an independent genome-wide TB case-control data set from The Gambia. A number of models were successfully validated, indicating that interplay between the NRG1 - NRG3, GRIK1 - GRIK3 and IL23R - ATG4C gene pairs may modify susceptibility to TB. Gene pairs involved in the NF-κB pathway were also identified in the discovery data set (SFTPD - NOD2, ISG15 - TLR8 and NLRC5 - IL12RB1), but could not be tested in the Gambian study group due to lack of overlapping data.

Using multi-way admixture mapping to elucidate TB susceptibility in the South African Coloured population

Background

The admixed South African Coloured population is ideally suited to the discovery of tuberculosis susceptibility genetic variants and their probable ethnic origins, but previous attempts at finding such variants using genome-wide admixture mapping were hampered by the inaccuracy of local ancestry inference. In this study, we infer local ancestry using the novel algorithm implemented in RFMix, with the emphasis on identifying regions of excess San or Bantu ancestry, which we hypothesize may harbour TB susceptibility genes.

Results

Using simulated data, we demonstrate reasonable accuracy of local ancestry inference by RFMix, with a tendency towards miss-calling San ancestry as Bantu. Regions with either excess San ancestry or excess African (San or Bantu) ancestry are less likely to be affected by this bias, and we therefore proceeded to identify such regions, found in cases but not in controls (642 cases and 91 controls). A number of promising regions were found (overall p-values of 7.19×10^-5 for San ancestry and <2.00×10^-16 for African ancestry), including chromosomes 15q15 and 17q22, which are close to genomic regions previously implicated in TB. Promising immune-related susceptibility genes such as the GADD45A, OSM and B7-H5 genes are also harboured in the identified regions.

Conclusion

Admixture mapping is feasible in the South African Coloured population and a number of novel TB susceptibility genomic regions were uncovered.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1021) contains supplementary material, which is available to authorized users.

A panel of ancestry informative markers for the complex five-way admixed South African coloured population

Admixture is a well known confounder in genetic association studies. If genome-wide data is not available, as would be the case for candidate gene studies, ancestry informative markers (AIMs) are required in order to adjust for admixture. The predominant population group in the Western Cape, South Africa, is the admixed group known as the South African Coloured (SAC). A small set of AIMs that is optimized to distinguish between the five source populations of this population (African San, African non-San, European, South Asian, and East Asian) will enable researchers to cost-effectively reduce false-positive findings resulting from ignoring admixture in genetic association studies of the population. Using genome-wide data to find SNPs with large allele frequency differences between the source populations of the SAC, as quantified by Rosenberg et. al's -statistic, we developed a panel of AIMs by experimenting with various selection strategies. Subsets of different sizes were evaluated by measuring the correlation between ancestry proportions estimated by each AIM subset with ancestry proportions estimated using genome-wide data. We show that a panel of 96 AIMs can be used to assess ancestry proportions and to adjust for the confounding effect of the complex five-way admixture that occurred in the South African Coloured population.

Genome-wide association study of ancestry-specific TB risk in the South African Coloured population

The worldwide burden of tuberculosis (TB) remains an enormous problem, and is particularly severe in the admixed South African Coloured (SAC) population residing in the Western Cape. Despite evidence from twin studies suggesting a strong genetic component to TB resistance, only a few loci have been identified to date. In this work, we conduct a genome-wide association study (GWAS), meta-analysis and trans-ethnic fine mapping to attempt the replication of previously identified TB susceptibility loci. Our GWAS results confirm the WT1 chr11 susceptibility locus (rs2057178: odds ratio = 0.62, P = 2.71e(-06)) previously identified by Thye et al., but fail to replicate previously identified polymorphisms in the TLR8 gene and locus 18q11.2. Our study demonstrates that the genetic contribution to TB risk varies between continental populations, and illustrates the value of including admixed populations in studies of TB risk and other complex phenotypes. Our evaluation of local ancestry based on the real and simulated data demonstrates that case-only admixture mapping is currently impractical in multi-way admixed populations, such as the SAC, due to spurious deviations in average local ancestry generated by current local ancestry inference methods. This study provides insights into identifying disease genes and ancestry-specific disease risk in multi-way admixed populations.

Determining ancestry proportions in complex admixture scenarios in South Africa using a novel proxy ancestry selection method

Admixed populations can make an important contribution to the discovery of disease susceptibility genes if the parental populations exhibit substantial variation in susceptibility. Admixture mapping has been used successfully, but is not designed to cope with populations that have more than two or three ancestral populations. The inference of admixture proportions and local ancestry and the imputation of missing genotypes in admixed populations are crucial in both understanding variation in disease and identifying novel disease loci. These inferences make use of reference populations, and accuracy depends on the choice of ancestral populations. Using an insufficient or inaccurate ancestral panel can result in erroneously inferred ancestry and affect the detection power of GWAS and meta-analysis when using imputation. Current algorithms are inadequate for multi-way admixed populations. To address these challenges we developed PROXYANC, an approach to select the best proxy ancestral populations. From the simulation of a multi-way admixed population we demonstrate the capability and accuracy of PROXYANC and illustrate the importance of the choice of ancestry in both estimating admixture proportions and imputing missing genotypes. We applied this approach to a complex, uniquely admixed South African population. Using genome-wide SNP data from over 764 individuals, we accurately estimate the genetic contributions from the best ancestral populations: isiXhosa [Formula: see text], ‡Khomani SAN [Formula: see text], European [Formula: see text], Indian [Formula: see text], and Chinese [Formula: see text]. We also demonstrate that the ancestral allele frequency differences correlate with increased linkage disequilibrium in the South African population, which originates from admixture events rather than population bottlenecks.

NOMENCLATURE

The collective term for people of mixed ancestry in southern Africa is "Coloured," and this is officially recognized in South Africa as a census term, and for self-classification. Whilst we acknowledge that some cultures may use this term in a derogatory manner, these connotations are not present in South Africa, and are certainly not intended here.

Associations between human leukocyte antigen class I variants and the Mycobacterium tuberculosis subtypes causing disease

BACKGROUND

 The development of active tuberculosis disease has been shown to be multifactorial. Interactions between host and bacterial genotype may influence disease outcome, with some studies indicating the adaptation of M. tuberculosis strains to specific human populations. Here we investigate the role of the human leukocyte antigen (HLA) class I genes in this biological process.

METHODS

 Three hundred patients with tuberculosis from South Africa were typed for their HLA class I alleles by direct sequencing. Mycobacterium tuberculosis genotype classification was done by IS6110 restriction fragment length polymorphism genotyping and spoligotyping.

RESULTS

 We showed that Beijing strain occurred more frequently in individuals with multiple disease episodes (P < .001) with the HLA-B27 allele lowering the odds of having an additional episode (odds ratio, 0.21; P = .006). Associations were also identified for specific HLA types and disease caused by the Beijing, LAM, LCC, and Quebec strains. HLA types were also associated with disease caused by strains from the Euro-American or East Asian lineages, and the frequencies of these alleles in their sympatric human populations identified potential coevolutionary events between host and pathogen.

CONCLUSIONS

 This is the first report of the association of human HLA types and M. tuberculosis strain genotype, highlighting that both host and pathogen genetics need to be taken into consideration when studying tuberculosis disease development.

Advanced Marjolin's ulcer of the scalp in a 13-year-old boy treated by excision and free tissue transfer: Case report and review of literature

Marjolin's ulcer originates in chronic scars and wounds of diverse origin. This relatively rare tumour is most commonly a squamous cell carcinoma. The reason for malignant transformation is not well understood. The burn scar is a common predilection. We present the youngest patient in the literature, a 13-year-old male with a locally advanced squamous cell carcinoma of the scalp with intracranial extension following an unhealed burn injury at the age of three. Bilateral cervical lymphadenopathy was also noted. The tumour was excised and the large defect overlying the brain was covered by free latissimus dorsi musculocutaneous flap. At four weeks a therapeutic bilateral selective neck dissection was done. Adjuvant chemotherapy was administered. This report reiterates the importance of early diagnosis. Free tissue transfer further enhances our ability to cover complex defects associated with excision of advanced lesions.