Genome-wide homozygosity signature and risk of Hodgkin lymphoma

Aberrant Expression and Oncogenic Activity of SPP1 in Hodgkin Lymphoma

Background: Hodgkin lymphoma (HL) is a B-cell-derived malignancy and one of the most frequent types of lymphoma. The tumour cells typically exhibit multiple genomic alterations together with aberrantly activated signalling pathways, driven by paracrine and/or autocrine modes. SPP1 (alias osteopontin) is a cytokine acting as a signalling activator and has been connected with relapse in HL patients. To understand its pathogenic role, here, we investigated the mechanisms and function of deregulated SPP1 in HL. Methods: We screened public patient datasets and cell lines for aberrant SPP1 expression. HL cell lines were stimulated with SPP1 and subjected to siRNA-mediated knockdown. Gene and protein activities were analyzed by RQ-PCR, ELISA, Western blot, and immuno-cytology. Results: SPP1 expression was detected in 8.3% of classic HL patients and in HL cell line SUP-HD1, chosen to serve as an experimental model. The gene encoding SPP1 is located at chromosomal position 4q22 and is genomically amplified in SUP-HD1. Transcription factor binding site analysis revealed TALE and HOX factors as potential regulators. Consistent with this finding, we showed that aberrantly expressed PBX1 and HOXB9 mediate the transcriptional activation of SPP1. RNA-seq data and knockdown experiments indicated that SPP1 signals via integrin ITGB1 in SUP-HD1. Accordingly, SPP1 activated NFkB in addition to MAPK/ERK which in turn mediated the nuclear import of ETS2, activating oncogenic JUNB expression. Conclusions: SPP1 is aberrantly activated in HL cell line SUP-HD1 via genomic copy number gain and by homeodomain transcription factors PBX1 and HOXB9. SPP1-activated NFkB and MAPK merit further investigation as potential therapeutic targets in affected HL patients.

Loss of smarcad1a accelerates tumorigenesis of malignant peripheral nerve sheath tumors in zebrafish

Malignant peripheral nerve sheath tumors (MPNSTs) are a type of sarcoma that generally originates from Schwann cells. The prognosis for this type of malignancy is relatively poor due to complicated genetic alterations and the lack of specific targeted therapy. Chromosome fragment 4q22-23 is frequently deleted in MPNSTs and other human tumors, suggesting tumor suppressor genes may reside in this region. Here, we provide evidence that SMARCAD1, a known chromatin remodeler, is a novel tumor suppressor gene located in 4q22-23. We identified two human homologous smarcad1 genes (smarcad1a and smarcad1b) in zebrafish, and both genes share overlapping expression patterns during embryonic development. We demonstrated that two smarcad1a loss-of-function mutants, sa1299 and p403, can accelerate MPNST tumorigenesis in the tp53 mutant background, suggesting smarcad1a is a bona fide tumor suppressor gene for MPNSTs. Moreover, we found that DNA double-strand break (DSB) repair might be compromised in both mutants compared to wildtype zebrafish, as indicated by pH2AX, a DNA DSB marker. In addition, both SMARCAD1 gene knockdown and overexpression in human cells were able to inhibit tumor growth and displayed similar DSB repair responses, suggesting proper SMARCAD1 gene expression level or gene dosage is critical for cell growth. Given that mutations of SMARCAD1 sensitize cells to poly ADP ribose polymerase inhibitors in yeast and the human U2OS osteosarcoma cell line, the identification of SMARCAD1 as a novel tumor suppressor gene might contribute to the development of new cancer therapies for MPNSTs.

Inbreeding, native American ancestry and child mortality: Linking human selection and paediatric medicine

The children of related parents show increased risk of early mortality. The Native American genome typically exhibits long stretches of homozygosity, and Latin Americans are highly heterogeneous regarding the individual burden of homozygosity, the proportion, and the type of Native American ancestry. We analysed nationwide mortality and genome-wide genotype data from admixed Chileans to investigate the relationship between common causes of child mortality, homozygosity and Native American ancestry. Results from two-stage linear-Poisson regression revealed a strong association between the sum length of runs of homozygosity (SROH) above 1.5 Megabases (Mb) in each genome and mortality due to intracranial non-traumatic haemorrhage of foetus and new-born (5% increased risk of death per Mb in SROH, P = 1 × 10-3) and disorders related to short gestation and low birth weight (P = 3 × 10-4). The major indigenous populations in Chile are Aymara-Quechua in the north of the country, and the Mapuche-Huilliche in the south. The individual proportion of Aymara-Quechua ancestry was associated with an increased risk of death due to anencephaly and similar malformations (P = 4 × 10-5), and the risk of death due to Edwards and Patau trisomy syndromes decreased 4% per 1% Aymara-Quechua ancestry proportion (P = 4 × 10-4) and 5% per 1% Mapuche-Huilliche ancestry proportion (P = 2 × 10-3). The present results suggest that short gestation, low birth weight and intracranial non-traumatic haemorrhage mediate the negative effect of inbreeding on human selection. Independent validation of the identified associations between common causes of child death, homozygosity and fine-scale ancestry proportions may inform paediatric medicine.

Genome-wide homozygosity and risk of four non-Hodgkin lymphoma subtypes

AIM

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

METHODS

We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.

RESULTS

We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10-6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10-5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.

CONCLUSION

Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

SMAD1 as a biomarker and potential therapeutic target in drug-resistant multiple myeloma

BACKGROUND

SMAD1, a central mediator in TGF-β signaling, is involved in a broad range of biological activities including cell growth, apoptosis, development and immune response, and is implicated in diverse type of malignancies. Whether SMAD1 plays an important role in multiple myeloma (MM) pathogenesis and can serve as a therapeutic target are largely unknown.

METHODS

Myeloma cell lines and primary MM samples were used. Cell culture, cytotoxicity and apoptosis assay, siRNA transfection, Western blot, RT-PCR, Soft-agar colony formation, and migration assay, Chromatin immunoprecipitation (Chip), animal xenograft model studies and statistical analysis were applied in this study.

RESULTS

We demonstrate that SMAD1 is highly expressed in myeloma cells of MM patients with advanced stages or relapsed disease, and is associated with significantly shorter progression-free and overall survivals. Mechanistically, we show that SMAD1 is required for TGFβ-mediated proliferation in MM via an ID1/p21/p27 pathway. TGF-β also enhanced TNFα-Induced protein 8 (TNFAIP8) expression and inhibited apoptosis through SMAD1-mediated induction of NF-κB1. Accordingly, depletion of SMAD1 led to downregulation of NF-κB1 and TNFAIP8, resulting in caspase-8-induced apoptosis. In turn, inhibition of NF-κB1 suppressed SMAD1 and ID1 expression uncovering an autoregulatory loop. Dorsomorphin (DM), a SMAD1 inhibitor, exerted a dose-dependent cytotoxic effect on drug-resistant MM cells with minimal cytotoxicity to normal hematopoietic cells, and further synergized with the proteasomal-inhibitor bortezomib to effectively kill drug-resistant MM cells in vitro and in a myeloma xenograft model.

CONCLUSIONS

This study identifies SMAD1 regulation of NF-κB1/TNFAIP8 and ID1-p21/p27 as critical axes of MM drug resistance and provides a potentially new therapeutic strategy to treat drug resistance MM through targeted inhibition of SMAD1.

Genome-wide association study of multiple yield traits in a diversity panel of polyploid sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) is an important economic crop, contributing up to 80% of sugar and approximately 60% of biofuel globally. To meet the increased demand for sugar and biofuel supplies, it is critical to breed sugarcane cultivars with robust performance in yield traits. Therefore, dissection of causal DNA sequence variants is of great importance, as it provides genetic resources and fundamental information for crop improvement. In this study, we analyzed nine yield traits in a sugarcane diversity panel consisting of 308 accessions primarily selected from the World Collection of Sugarcane and Related Grasses. By genotyping the diversity panel via target enrichment sequencing, we identified a large number of sequence variants. Genome-wide association studies between the markers and traits were conducted, taking dosages and gene actions into consideration. In total, 217 nonredundant markers and 225 candidate genes were identified to be significantly associated with the yield traits, which can serve as a comprehensive genetic resource database for future gene identification, characterization, and selection for sugarcane improvement. We further investigated runs of homozygosity (ROH) in the sugarcane diversity panel. We characterized 282 ROHs and found that the occurrence of ROHs in the genome were nonrandom and probably under selection. The ROHs were associated with total weight and dry weight, and high ROHs resulted in a decrease in the two traits. This study suggests that genomic inbreeding has led to negative impacts on sugarcane yield.

Runs of homozygosity and testicular cancer risk

BACKGROUND

Testicular germ cell tumour (TGCT) is highly heritable but > 50% of the genetic risk remains unexplained. Epidemiological observation of greater relative risk to brothers of men with TGCT compared to sons has long alluded to recessively acting TGCT genetic susceptibility factors, but to date none have been reported. Runs of homozygosity (RoH) are a signature indicating underlying recessively acting alleles and have been associated with increased risk of other cancer types.

OBJECTIVE

To examine whether RoH are associated with TGCT risk.

METHODS

We performed a genome-wide RoH analysis using GWAS data from 3206 TGCT cases and 7422 controls uniformly genotyped using the OncoArray platform.

RESULTS

Global measures of homozygosity were not significantly different between cases and controls, and the frequency of individual consensus RoH was not significantly different between cases and controls, after correction for multiple testing. RoH at three regions, 11p13-11p14.3, 5q14.1-5q22.3 and 13q14.11-13q.14.13, were, however, nominally statistically significant at p < 0.01. Intriguingly, RoH200 at 11p13-11p14.3 encompasses Wilms tumour 1 (WT1), a recognized cancer susceptibility gene with roles in sex determination and developmental transcriptional regulation, processes repeatedly implicated in TGCT aetiology.

DISCUSSION AND CONCLUSION

Overall, our data do not support a major role in the risk of TGCT for recessively acting alleles acting through homozygosity, as measured by RoH in outbred populations of cases and controls.

GARLIC: Genomic Autozygosity Regions Likelihood-based Inference and Classification

Summary

Runs of homozygosity (ROH) are important genomic features that manifest when identical-by-descent haplotypes are inherited from parents. Their length distributions and genomic locations are informative about population history and they are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. Here, we present software implementing a model-based method ( Pemberton et al., 2012 ) for inferring ROH in genome-wide SNP datasets that incorporates population-specific parameters and a genotyping error rate as well as provides a length-based classification module to identify biologically interesting classes of ROH. Using simulations, we evaluate the performance of this method.

Availability and Implementation

GARLIC is written in C ++. Source code and pre-compiled binaries (Windows, OSX and Linux) are hosted on GitHub ( https://github.com/szpiech/garlic ) under the GNU General Public License version 3.

Contact

zachary.szpiech@ucsf.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.

Breast cancer protection by genomic imprinting in close kin families

Human inbreeding generally reduces breast cancer risk (BCR). When the parents are biologically related, their infants have a lower birth weight due to smaller body organs. The undersized breasts, because of fewer mammary stem cells, have a lower likelihood of malignant conversion. Fetal growth is regulated by genomically imprinted genes which are in conflict; they promote growth when derived from the father and suppress growth when derived from the mother. The kinship theory explicates that the intensity of conflict between these genes affects growth and therefore the size of the newborn. In descendants of closely related parents, this gene clash is less resulting in a smaller infant. In this review, we elucidate the different mechanisms by which human inbreeding affects BCR, and why this risk is dissimilar in different inbred populations.

Evidence for genetic association between chromosome 1q loci and predisposition to colorectal neoplasia

BACKGROUND

A substantial fraction of familial colorectal cancer (CRC) and polyposis heritability remains unexplained. This study aimed to identify predisposing loci in patients with these disorders.

METHODS

Homozygosity mapping was performed using 222 563 SNPs in 302 index patients with various colorectal neoplasms and 3367 controls. Linkage analysis, exome and whole-genome sequencing were performed in a family affected by microsatellite stable CRCs. Candidate variants were genotyped in 10 554 cases and 21 480 controls. Gene expression was assessed at the mRNA and protein level.

RESULTS

Homozygosity mapping revealed a disease-associated region at 1q32.3 which was part of the linkage region 1q32.2-42.2 identified in the CRC family. This includes a region previously associated with risk of CRC. Sequencing identified the p.Asp1432Glu variant in the MIA3 gene (known as TANGO1 or TANGO) and 472 additional rare, shared variants within the linkage region. In both cases and controls the population frequency was 0.02% for this MIA3 variant. The MIA3 mutant allele showed predominant mRNA expression in normal, cancer and precancerous tissues. Furthermore, immunohistochemistry revealed increased expression of MIA3 in adenomatous tissues.

CONCLUSIONS

Taken together, our two independent strategies associate genetic variations in chromosome 1q loci and predisposition to familial CRC and polyps, which warrants further investigation.

Evidence of Inbreeding in Hodgkin Lymphoma

Genome-wide association studies (GWASs) have identified several, mainly co-dominantly acting, single-nucleotide polymorphisms (SNPs) associated with Hodgkin lymphoma (HL). We searched for recessively acting disease loci by performing an analysis of runs of homozygosity (ROH) based on windows of homozygous SNP-blocks and by calculating genomic inbreeding coefficients on a SNP-wise basis. We used data from a previous GWAS with 906 cases and 1217 controls from a population with a long history of no matings between relatives. Ten recurrent ROHs were identified among 25 055 ROHs across all individuals but their association with HL was not genome-wide significant. All recurrent ROHs showed significant evidence for natural selection. As a novel finding genomic inbreeding among cases was significantly higher than among controls (P = 2.11*10-14) even after correcting for covariates. Higher inbreeding among the cases was mainly based on a group of individuals with a higher average length of ROHs per person. This result suggests a correlation of higher levels of inbreeding with higher cancer incidence and might reflect the existence of recessive alleles causing HL. Genomic inbreeding may result in a higher expression of deleterious recessive genes within a population.