Evidence of Inbreeding in Hodgkin Lymphoma

A genomic tale of inbreeding in western Mediterranean human populations

Consanguineous marriages are common in many worldwide human populations, and the biological consequences for offspring can be relevant at the biomedical level. The current genomic revolution displayed through genome-wide studies is challenging the paradigm in the analysis of consanguinity. Here, we analyzed genomic inbreeding patterns in human populations located at the western edge of the Mediterranean region (Iberia and Morocco). Runs of Homozygosity (ROH) (autozygosity fragments) were identified in 139 autochthonous individuals originating from southern Iberia and Morocco via microarray data. All individuals analyzed carried at least one ROH in their genomes. The genomic inbreeding coefficient (FROH) and the presence of ROH islands (ROHi) revealed interesting patterns in the target populations as well as in the rest of the Mediterranean basin. Moroccan Berbers presented signals of recent inbreeding, relying on high coverage of long ROH (> 5 Mb) and FROH. The location and structure of ROHi among people in the western Mediterranean could be interpreted as a signature of common genetic links across the Strait of Gibraltar. We found a significant enrichment of some relevant biological functions in the estimated ROHi hotspots associated with the immune system and chemosensation. Genomic inbreeding approaches allow us to understand past population histories and can be used as a proxy to scan the genome in search of selection signals.

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.

Regions of homozygosity as risk factors for multiple myeloma

Genomic regions of homozygosity (ROH), detectable in outbred populations, have been implicated as determinants of inherited risk. To examine whether ROH is associated with risk of multiple myeloma (MM), we performed whole-genome homozygosity analysis using single-nucleotide polymorphism genotype data from 2,282 MM cases and 5,197 controls, with replication in an additional 878 MM cases and 7,083 controls. Globally, the distribution of ROH between cases and controls was not significantly different. However, one ROH at chromosome 9q21, harboring the B-cell transcription factor gene KLF9, showed evidence of a consistent association and may therefore warrant further investigation as a candidate risk factor for MM. Overall, our analysis provides little support for a homozygosity signature being a significant factor in MM risk.

Genetic diversity, inbreeding and cancer

Genetic diversity is essential for adaptive capacities, providing organisms with the potential of successfully responding to intrinsic and extrinsic challenges. Although a clear reciprocal link between genetic diversity and resistance to parasites and pathogens has been established across taxa, the impact of loss of genetic diversity by inbreeding on the emergence and progression of non-communicable diseases, such as cancer, has been overlooked. Here we provide an overview of such associations and show that low genetic diversity and inbreeding associate with an increased risk of cancer in both humans and animals. Cancer being a multifaceted disease, loss of genetic diversity can directly (via accumulation of oncogenic homozygous mutations) and indirectly (via increased susceptibility to oncogenic pathogens) impact abnormal cell emergence and escape of immune surveillance. The observed link between reduced genetic diversity and cancer in wildlife may further imperil the long-term survival of numerous endangered species, highlighting the need to consider the impact of cancer in conservation biology. Finally, the somewhat incongruent data originating from human studies suggest that the association between genetic diversity and cancer development is multifactorial and may be tumour specific. Further studies are therefore crucial in order to elucidate the underpinnings of the interactions between genetic diversity, inbreeding and cancer.

A meta-analysis of consanguinity and breast cancer

BACKGROUND

There have been various publications stating that consanguinity both increases and decreases the risk of breast cancer.

AIMS

The objective of this study was to determine the impact of consanguinity upon breast cancer. We conducted a systematic review of the literature and meta-analysis.

METHODS

Eligible studies were identified on Medline and EMBASE updated to the 19 of September 2017. Studies with sufficient comparative data were included in a meta-analysis. Analyses were carried out using RevMan software.

RESULTS

Three comparative studies with a total of 317 individuals with breast cancer and 1459 controls. Reviewing the literature demonstrated conflicting conclusions of the influence of consanguinity upon breast cancer. The meta-analysis showed that there were no statistically significant associations between consanguinity and breast cancer though there was a trend protection from a history of consanguinity.

CONCLUSION

Though there is limited literature published on the effects of parental consanguinity, the available data does not demonstrate that it is a risk factor for breast cancer.

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.

Evidence of Recent Intricate Adaptation in Human Populations

Recent human adaptations have shaped population differentiation in genomic regions containing putative functional variants, mostly located in predicted regulatory elements. However, their actual functionalities and the underlying mechanism of recent adaptation remain poorly understood. In the current study, regions of genes and repeats were investigated for functionality depending on the degree of population differentiation, F_ST or ΔDAF (a difference in derived allele frequency). The high F_ST in the 5´ or 3´ untranslated regions (UTRs), in particular, confirmed that population differences arose mainly from differences in regulation. Expression quantitative trait loci (eQTL) analyses using lymphoblastoid cell lines indicated that the majority of the highly population-specific regions represented cis- and/or trans-eQTL. However, groups having the highest ΔDAFs did not necessarily have higher proportions of eQTL variants; in these groups, the patterns were complex, indicating recent intricate adaptations. The results indicated that East Asian (EAS) and European populations (EUR) experienced mutual selection pressures. The mean derived allele frequency of the high ΔDAF groups suggested that EAS and EUR underwent strong adaptation; however, the African population in Africa (AFR) experienced slight, yet broad, adaptation. The DAF distributions of variants in the gene regions showed clear selective pressure in each population, which implies the existence of more recent regulatory adaptations in cells other than lymphoblastoid cell lines. In-depth analysis of population-differentiated regions indicated that the coding gene, RNF135, represented a trans-regulation hotspot via cis-regulation by the population-specific variants in the region of selective sweep. Together, the results provide strong evidence of actual intricate adaptation of human populations via regulatory manipulation.