Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Nucleus-Translocated ACSS2 Promotes Gene Transcription for Lysosomal Biogenesis and Autophagy.

Overcoming metabolic stress is a critical step in tumor growth. Acetyl coenzyme A (acetyl-CoA) generated from glucose and acetate uptake is important for histone acetylation and gene expression. However, how acetyl-CoA is produced under nutritional stress is unclear. We demonstrate here that glucose deprivation results in AMP-activated protein kinase (AMPK)-mediated acetyl-CoA synthetase 2 (ACSS2) phosphorylation at S659, which exposed the nuclear localization signal of ACSS2 for importin α5 binding and nuclear translocation. In the nucleus, ACSS2 binds to transcription factor EB and translocates to lysosomal and autophagy gene promoter regions, where ACSS2 incorporates acetate generated from histone acetylation turnover to locally produce acetyl-CoA for histone H3 acetylation in these regions and promote lysosomal biogenesis, autophagy, cell survival, and brain tumorigenesis. In addition, ACSS2 S659 phosphorylation positively correlates with AMPK activity in glioma specimens and grades of glioma malignancy. These results underscore the significance of nuclear ACSS2-mediated histone acetylation in maintaining cell homeostasis and tumor development.

Recurrent mutations in topoisomerase IIα cause a previously undescribed mutator phenotype in human cancers.

Topoisomerases nick and reseal DNA to relieve torsional stress associated with transcription and replication and to resolve structures such as knots and catenanes. Stabilization of the yeast Top2 cleavage intermediates is mutagenic in yeast, but whether this extends to higher eukaryotes is less clear. Chemotherapeutic topoisomerase poisons also elevate cleavage, resulting in mutagenesis. Here, we describe p.K743N mutations in human topoisomerase hTOP2α and link them to a previously undescribed mutator phenotype in cancer. Overexpression of the orthologous mutant protein in yeast generated a characteristic pattern of 2- to 4-base pair (bp) duplications resembling those in tumors with p.K743N. Using mutant strains and biochemical analysis, we determined the genetic requirements of this mutagenic process and showed that it results from trapping of the mutant yeast yTop2 cleavage complex. In addition to 2- to 4-bp duplications, hTOP2α p.K743N is also associated with deletions that are absent in yeast. We call the combined pattern of duplications and deletions ID_TOP2α. All seven tumors carrying the hTOP2α p.K743N mutation showed ID_TOP2α, while it was absent from all other tumors examined (n = 12,269). Each tumor with the ID_TOP2α signature had indels in several known cancer genes, which included frameshift mutations in tumor suppressors PTEN and TP53 and an activating insertion in BRAF. Sequence motifs found at ID_TOP2α mutations were present at 80% of indels in cancer-driver genes, suggesting that ID_TOP2α mutagenesis may contribute to tumorigenesis. The results reported here shed further light on the role of topoisomerase II in genome instability.

Characterization of colibactin-associated mutational signature in an Asian oral squamous cell carcinoma and in other mucosal tumor types.

Mutational signatures can reveal the history of mutagenic processes that cells were exposed to before and during tumorigenesis. We expect that as-yet-undiscovered mutational processes will shed further light on mutagenesis leading to carcinogenesis. With this in mind, we analyzed the mutational spectra of 36 Asian oral squamous cell carcinomas. The mutational spectra of two samples from patients who presented with oral bacterial infections showed novel mutational signatures. One of these novel signatures, SBS_AnT, is characterized by a preponderance of thymine mutations, strong transcriptional strand bias, and enrichment for adenines in the 4 bp 5' of mutation sites. The mutational signature described in this manuscript was shown to be caused by colibactin, a bacterial mutagen produced by E. coli carrying the pks-island. Examination of publicly available sequencing data revealed SBS_AnT in 25 tumors from several mucosal tissue types, expanding the list of tissues in which this mutational signature is observed.

A tumor-associated splice-isoform of MAP2K7 drives dedifferentiation in MBNL1-low cancers via JNK activation.

Master splicing regulator MBNL1 shapes large transcriptomic changes that drive cellular differentiation during development. Here we demonstrate that MBNL1 is a suppressor of tumor dedifferentiation. We surveyed MBNL1 expression in matched tumor/normal pairs across The Cancer Genome Atlas and found that MBNL1 was down-regulated in several common cancers. Down-regulation of MBNL1 predicted poor overall survival in breast, lung, and stomach adenocarcinomas and increased relapse and distant metastasis in triple-negative breast cancer. Down-regulation of MBNL1 led to increased tumorigenic and stem/progenitor-like properties in vitro and in vivo. A discrete set of alternative splicing events (ASEs) are shared between MBNL1-low cancers and embryonic stem cells including a MAP2K7∆exon2 splice variant that leads to increased stem/progenitor-like properties via JNK activation. Accordingly, JNK inhibition is capable of reversing MAP2K7∆exon2-driven tumor dedifferentiation in MBNL1-low cancer cells. Our work elucidates an alternative-splicing mechanism that drives tumor dedifferentiation and identifies biomarkers that predict enhanced susceptibility to JNK inhibition.

Experimental and pan-cancer genome analyses reveal widespread contribution of acrylamide exposure to carcinogenesis in humans.

Humans are frequently exposed to acrylamide, a probable human carcinogen found in commonplace sources such as most heated starchy foods or tobacco smoke. Prior evidence has shown that acrylamide causes cancer in rodents, yet epidemiological studies conducted to date are limited and, thus far, have yielded inconclusive data on association of human cancers with acrylamide exposure. In this study, we experimentally identify a novel and unique mutational signature imprinted by acrylamide through the effects of its reactive metabolite glycidamide. We next show that the glycidamide mutational signature is found in a full one-third of approximately 1600 tumor genomes corresponding to 19 human tumor types from 14 organs. The highest enrichment of the glycidamide signature was observed in the cancers of the lung (88% of the interrogated tumors), liver (73%), kidney (>70%), bile duct (57%), cervix (50%), and, to a lesser extent, additional cancer types. Overall, our study reveals an unexpectedly extensive contribution of acrylamide-associated mutagenesis to human cancers.

The mutational landscape of early- and typical-onset oral tongue squamous cell carcinoma.

BACKGROUND: The incidence of oral tongue squamous cell carcinoma (OTSCC) is increasing among younger birth cohorts. The etiology of early-onset OTSCC (diagnosed before the age of 50 years) and cancer driver genes remain largely unknown. METHODS: The Sequencing Consortium of Oral Tongue Cancer was established through the pooling of somatic mutation data of oral tongue cancer specimens (n = 227 [107 early-onset cases]) from 7 studies and The Cancer Genome Atlas. Somatic mutations at microsatellite loci and Catalog of Somatic Mutations in Cancer mutation signatures were identified. Cancer driver genes were identified with the MutSigCV and WITER algorithms. Mutation comparisons between early- and typical-onset OTSCC were evaluated via linear regression with adjustments for patient-related factors. RESULTS: Two novel driver genes (ATXN1 and CDC42EP1) and 5 previously reported driver genes (TP53, CDKN2A, CASP8, NOTCH1, and FAT1) were identified. Six recurrent mutations were identified, with 4 occurring in TP53. Early-onset OTSCC had significantly fewer nonsilent mutations even after adjustments for tobacco use. No associations of microsatellite locus mutations and mutation signatures with the age of OTSCC onset were observed. CONCLUSIONS: This international, multicenter consortium is the largest study to characterize the somatic mutational landscape of OTSCC and the first to suggest differences by age of onset. This study validates multiple previously identified OTSCC driver genes and proposes 2 novel cancer driver genes. In analyses by age, early-onset OTSCC had a significantly smaller somatic mutational burden that was not explained by differences in tobacco use. LAY SUMMARY: This study identifies 7 specific areas in the human genetic code that could be responsible for promoting the development of tongue cancer. Tongue cancer in young patients (under the age of 50 years) has fewer overall changes to the genetic code in comparison with tongue cancer in older patients, but the authors do not think that this is due to differences in smoking rates between the 2 groups. The cause of increasing cases of tongue cancer in young patients remains unclear.

In-depth characterization of the cisplatin mutational signature in human cell lines and in esophageal and liver tumors.

Cisplatin reacts with DNA and thereby likely generates a characteristic pattern of somatic mutations, called a mutational signature. Despite widespread use of cisplatin in cancer treatment and its role in contributing to secondary malignancies, its mutational signature has not been delineated. We hypothesize that cisplatin's mutational signature can serve as a biomarker to identify cisplatin mutagenesis in suspected secondary malignancies. Knowledge of which tissues are at risk of developing cisplatin-induced secondary malignancies could lead to guidelines for noninvasive monitoring for secondary malignancies after cisplatin chemotherapy. We performed whole genome sequencing of 10 independent clones of cisplatin-exposed MCF-10A and HepG2 cells and delineated the patterns of single and dinucleotide mutations in terms of flanking sequence, transcription strand bias, and other characteristics. We used the mSigAct signature presence test and nonnegative matrix factorization to search for cisplatin mutagenesis in hepatocellular carcinomas and esophageal adenocarcinomas. All clones showed highly consistent patterns of single and dinucleotide substitutions. The proportion of dinucleotide substitutions was high: 8.1% of single nucleotide substitutions were part of dinucleotide substitutions, presumably due to cisplatin's propensity to form intra- and interstrand crosslinks between purine bases in DNA. We identified likely cisplatin exposure in nine hepatocellular carcinomas and three esophageal adenocarcinomas. All hepatocellular carcinomas for which clinical data were available and all esophageal cancers indeed had histories of cisplatin treatment. We experimentally delineated the single and dinucleotide mutational signature of cisplatin. This signature enabled us to detect previous cisplatin exposure in human hepatocellular carcinomas and esophageal adenocarcinomas with high confidence.

Beyond fitness tracking: The use of consumer-grade wearable data from normal volunteers in cardiovascular and lipidomics research.

The use of consumer-grade wearables for purposes beyond fitness tracking has not been comprehensively explored. We generated and analyzed multidimensional data from 233 normal volunteers, integrating wearable data, lifestyle questionnaires, cardiac imaging, sphingolipid profiling, and multiple clinical-grade cardiovascular and metabolic disease markers. We show that subjects can be stratified into distinct clusters based on daily activity patterns and that these clusters are marked by distinct demographic and behavioral patterns. While resting heart rates (RHRs) performed better than step counts in being associated with cardiovascular and metabolic disease markers, step counts identified relationships between physical activity and cardiac remodeling, suggesting that wearable data may play a role in reducing overdiagnosis of cardiac hypertrophy or dilatation in active individuals. Wearable-derived activity levels can be used to identify known and novel activity-modulated sphingolipids that are in turn associated with insulin sensitivity. Our findings demonstrate the potential for wearables in biomedical research and personalized health.

Integrated paired-end enhancer profiling and whole-genome sequencing reveals recurrent CCNE1 and IGF2 enhancer hijacking in primary gastric adenocarcinoma.

OBJECTIVE: Genomic structural variations (SVs) causing rewiring of cis-regulatory elements remain largely unexplored in gastric cancer (GC). To identify SVs affecting enhancer elements in GC (enhancer-based SVs), we integrated epigenomic enhancer profiles revealed by paired-end H3K27ac ChIP-sequencing from primary GCs with tumour whole-genome sequencing (WGS) data (PeNChIP-seq/WGS). DESIGN: We applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of >200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs. RESULTS: PeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targeting CCNE1, a key driver of therapy resistance, occurred in 8% of patients frequently juxtaposing diverse distal enhancers to CCNE1 proximal regions. CCNE1-rearranged GCs were associated with high CCNE1 expression, disrupted CCNE1 topologically associating domain (TAD) boundaries, and novel TAD interactions in CCNE1-rearranged primary tumours. We also observed IGF2 enhancer-based SVs, previously noted in colorectal cancer, highlighting a common non-coding genetic driver alteration in gastric and colorectal malignancies. CONCLUSION: Integrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.

Genome-scale mutational signatures of aflatoxin in cells, mice, and human tumors.

Aflatoxin B1 (AFB1) is a mutagen and IARC (International Agency for Research on Cancer) Group 1 carcinogen that causes hepatocellular carcinoma (HCC). Here, we present the first whole-genome data on the mutational signatures of AFB1 exposure from a total of >40,000 mutations in four experimental systems: two different human cell lines, in liver tumors in wild-type mice, and in mice that carried a hepatitis B surface antigen transgene-this to model the multiplicative effects of aflatoxin exposure and hepatitis B in causing HCC. AFB1 mutational signatures from all four experimental systems were remarkably similar. We integrated the experimental mutational signatures with data from newly sequenced HCCs from Qidong County, China, a region of well-studied aflatoxin exposure. This indicated that COSMIC mutational signature 24, previously hypothesized to stem from aflatoxin exposure, indeed likely represents AFB1 exposure, possibly combined with other exposures. Among published somatic mutation data, we found evidence of AFB1 exposure in 0.7% of HCCs treated in North America, 1% of HCCs from Japan, but 16% of HCCs from Hong Kong. Thus, aflatoxin exposure apparently remains a substantial public health issue in some areas. This aspect of our study exemplifies the promise of future widespread resequencing of tumor genomes in providing new insights into the contribution of mutagenic exposures to cancer incidence.

SigProfilerMatrixGenerator: a tool for visualizing and exploring patterns of small mutational events.

BACKGROUND: Cancer genomes are peppered with somatic mutations imprinted by different mutational processes. The mutational pattern of a cancer genome can be used to identify and understand the etiology of the underlying mutational processes. A plethora of prior research has focused on examining mutational signatures and mutational patterns from single base substitutions and their immediate sequencing context. We recently demonstrated that further classification of small mutational events (including substitutions, insertions, deletions, and doublet substitutions) can be used to provide a deeper understanding of the mutational processes that have molded a cancer genome. However, there has been no standard tool that allows fast, accurate, and comprehensive classification for all types of small mutational events. RESULTS: Here, we present SigProfilerMatrixGenerator, a computational tool designed for optimized exploration and visualization of mutational patterns for all types of small mutational events. SigProfilerMatrixGenerator is written in Python with an R wrapper package provided for users that prefer working in an R environment. SigProfilerMatrixGenerator produces fourteen distinct matrices by considering transcriptional strand bias of individual events and by incorporating distinct classifications for single base substitutions, doublet base substitutions, and small insertions and deletions. While the tool provides a comprehensive classification of mutations, SigProfilerMatrixGenerator is also faster and more memory efficient than existing tools that generate only a single matrix. CONCLUSIONS: SigProfilerMatrixGenerator provides a standardized method for classifying small mutational events that is both efficient and scalable to large datasets. In addition to extending the classification of single base substitutions, the tool is the first to provide support for classifying doublet base substitutions and small insertions and deletions. SigProfilerMatrixGenerator is freely available at https://github.com/AlexandrovLab/SigProfilerMatrixGenerator with an extensive documentation at https://osf.io/s93d5/wiki/home/ .

Wnt proteins synergize to activate ß-catenin signaling.

Wnt ligands are involved in diverse signaling pathways that are active during development, maintenance of tissue homeostasis and in various disease states. While signaling regulated by individual Wnts has been extensively studied, Wnts are rarely expressed alone, and the consequences of Wnt gene co-expression are not well understood. Here, we studied the effect of co-expression of Wnts on the ß-catenin signaling pathway. While some Wnts are deemed 'non-canonical' due to their limited ability to activate ß-catenin when expressed alone, unexpectedly, we find that multiple Wnt combinations can synergistically activate ß-catenin signaling in multiple cell types. WNT1- and WNT7B-mediated synergistic Wnt signaling requires FZD5, FZD8 and LRP6, as well as the WNT7B co-receptors GPR124 (also known as ADGRA2) and RECK. Unexpectedly, this synergistic signaling occurs downstream of ß-catenin stabilization, and is correlated with increased lysine acetylation of ß-catenin. Wnt synergy provides a general mechanism to confer increased combinatorial control over this important regulatory pathway.

Experimental Delineation of Mutational Signatures Is an Essential Tool in Cancer Epidemiology and Prevention.

Advances in experimental modeling of the mutational signatures of environmental exposures and endogenous mutagenic processes will elucidate the role of mutagenesis in cancer, facilitate carcinogen classification, and enable new molecular cancer epidemiology studies.

Functional genomics identifies specific vulnerabilities in PTEN-deficient breast cancer.

BACKGROUND: Phosphatase and tensin homolog (PTEN) is one of the most frequently inactivated tumor suppressors in breast cancer. While PTEN itself is not considered a druggable target, PTEN synthetic-sick or synthetic-lethal (PTEN-SSL) genes are potential drug targets in PTEN-deficient breast cancers. Therefore, with the aim of identifying potential targets for precision breast cancer therapy, we sought to discover PTEN-SSL genes present in a broad spectrum of breast cancers. METHODS: To discover broad-spectrum PTEN-SSL genes in breast cancer, we used a multi-step approach that started with (1) a genome-wide short interfering RNA (siRNA) screen of ~ 21,000 genes in a pair of isogenic human mammary epithelial cell lines, followed by (2) a short hairpin RNA (shRNA) screen of ~ 1200 genes focused on hits from the first screen in a panel of 11 breast cancer cell lines; we then determined reproducibility of hits by (3) identification of overlaps between our results and reanalyzed data from 3 independent gene-essentiality screens, and finally, for selected candidate PTEN-SSL genes we (4) confirmed PTEN-SSL activity using either drug sensitivity experiments in a panel of 19 cell lines or mutual exclusivity analysis of publicly available pan-cancer somatic mutation data. RESULTS: The screens (steps 1 and 2) and the reproducibility analysis (step 3) identified six candidate broad-spectrum PTEN-SSL genes (PIK3CB, ADAMTS20, AP1M2, HMMR, STK11, and NUAK1). PIK3CB was previously identified as PTEN-SSL, while the other five genes represent novel PTEN-SSL candidates. Confirmation studies (step 4) provided additional evidence that NUAK1 and STK11 have PTEN-SSL patterns of activity. Consistent with PTEN-SSL status, inhibition of the NUAK1 protein kinase by the small molecule drug HTH-01-015 selectively impaired viability in multiple PTEN-deficient breast cancer cell lines, while mutations affecting STK11 and PTEN were largely mutually exclusive across large pan-cancer data sets. CONCLUSIONS: Six genes showed PTEN-SSL patterns of activity in a large proportion of PTEN-deficient breast cancer cell lines and are potential specific vulnerabilities in PTEN-deficient breast cancer. Furthermore, the NUAK1 PTEN-SSL vulnerability identified by RNA interference techniques can be recapitulated and exploited using the small molecule kinase inhibitor HTH-01-015. Thus, NUAK1 inhibition may be an effective strategy for precision treatment of PTEN-deficient breast tumors.

Abundant copy-number loss of CYCLOPS and STOP genes in gastric adenocarcinoma.

BACKGROUND: Gastric cancer, a leading cause of cancer death worldwide, has been little studied compared with other cancers that impose similar health burdens. Our goal is to assess genomic copy-number loss and the possible functional consequences and therapeutic implications thereof across a large series of gastric adenocarcinomas. METHODS: We used high-density single-nucleotide polymorphism microarrays to determine patterns of copy-number loss and allelic imbalance in 74 gastric adenocarcinomas. We investigated whether suppressor of tumorigenesis and/or proliferation (STOP) genes are associated with genomic copy-number loss. We also analyzed the extent to which copy-number loss affects Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS (CYCLOPS) genes-genes that may be attractive targets for therapeutic inhibition when partially deleted. RESULTS: The proportion of the genome subject to copy-number loss varies considerably from tumor to tumor, with a median of 5.5 %, and a mean of 12 % (range 0-58.5 %). On average, 91 STOP genes were subject to copy-number loss per tumor (median 35, range 0-452), and STOP genes tended to have lower copy-number compared with the rest of the genes. Furthermore, on average, 1.6 CYCLOPS genes per tumor were both subject to copy-number loss and downregulated, and 51.4 % of the tumors had at least one such gene. CONCLUSIONS: The enrichment of STOP genes in regions of copy-number loss indicates that their deletion may contribute to gastric carcinogenesis. Furthermore, the presence of several deleted and downregulated CYCLOPS genes in some tumors suggests potential therapeutic targets in these tumors.

Extensive diversity in circadian regulation of plasma lipids and evidence for different circadian metabolic phenotypes in humans.

The circadian system regulates daily rhythms in lipid metabolism and adipose tissue function. Although disruption of circadian clock function is associated with negative cardiometabolic end points, very little is known about interindividual variation in circadian-regulated metabolic pathways. Here, we used targeted lipidomics-based approaches to profile the time course of 263 lipids in blood plasma in 20 healthy individuals. Over a span of 28 h, blood was collected every 4 h and plasma lipids were analyzed by HPLC/MS. Across subjects, about 13% of lipid metabolites showed circadian variation. Rhythmicity spanned all metabolite classes examined, suggesting widespread circadian control of lipid-mediated energy storage, transport, and signaling. Intersubject agreement for lipids identified as rhythmic was only about 20%, however, and the timing of lipid rhythms ranged up to 12 h apart between individuals. Healthy subjects therefore showed substantial variation in the timing and strength of rhythms across different lipid species. Strong interindividual differences were also observed for rhythms of blood glucose and insulin, but not cortisol. Using consensus clustering with iterative feature selection, subjects clustered into different groups based on strength of rhythmicity for a subset of triglycerides and phosphatidylcholines, suggesting that there are different circadian metabolic phenotypes in the general population. These results have potential implications for lipid metabolism disorders linked to circadian clock disruption.

AZFc deletions and spermatogenic failure: a population-based survey of 20,000 Y chromosomes.

Deletions involving the Y chromosome's AZFc region are the most common known genetic cause of severe spermatogenic failure (SSF). Six recurrent interstitial deletions affecting the region have been reported, but their population genetics are largely unexplored. We assessed the deletions' prevalence in 20,884 men in five populations and found four of the six deletions (presented here in descending order of prevalence): gr/gr, b2/b3, b1/b3, and b2/b4. One of every 27 men carried one of these four deletions. The 1.6 Mb gr/gr deletion, found in one of every 41 men, almost doubles the risk of SSF and accounts for ∼2% of SSF, although <2% of men with the deletion are affected. The 1.8 Mb b2/b3 deletion, found in one of every 90 men, does not appear to be a risk factor for SSF. The 1.6 Mb b1/b3 deletion, found in one of every 994 men, appears to increase the risk of SSF by a factor of 2.5, although <2% of men with the deletion are affected, and it accounts for only 0.15% of SSF. The 3.5 Mb b2/b4 deletion, found in one of every 2,320 men, increases the risk of SSF 145 times and accounts for ∼6% of SSF; the observed prevalence should approximate the rate at which the deletion arises anew in each generation. We conclude that a single rare variant of major effect (the b2/b4 deletion) and a single common variant of modest effect (the gr/gr deletion) are largely responsible for the AZFc region's contribution to SSF in the population.