The Cell Type-Specific 5hmC Landscape and Dynamics of Healthy Human Hematopoiesis and TET2-Mutant Preleukemia.

The conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) is a key step in DNA demethylation that is mediated by ten-eleven translocation (TET) enzymes, which require ascorbate/vitamin C. Here, we report the 5hmC landscape of normal hematopoiesis and identify cell type-specific 5hmC profiles associated with active transcription and chromatin accessibility of key hematopoietic regulators. We utilized CRISPR/Cas9 to model TET2 loss-of-function mutations in primary human hematopoietic stem and progenitor cells (HSPC). Disrupted cells exhibited increased colonies in serial replating, defective erythroid/megakaryocytic differentiation, and in vivo competitive advantage and myeloid skewing coupled with reduction of 5hmC at erythroid-associated gene loci. Azacitidine and ascorbate restored 5hmC abundance and slowed or reverted the expansion of TET2-mutant clones in vivo. These results demonstrate the key role of 5hmC in normal hematopoiesis and TET2-mutant phenotypes and raise the possibility of utilizing these agents to further our understanding of preleukemia and clonal hematopoiesis. SIGNIFICANCE: We show that 5-hydroxymethylation profiles are cell type-specific and associated with transcriptional abundance and chromatin accessibility across human hematopoiesis. TET2 loss caused aberrant growth and differentiation phenotypes and disrupted 5hmC and transcriptional landscapes. Treatment of TET2 KO HSPCs with ascorbate or azacitidine reverted 5hmC profiles and restored aberrant phenotypes. This article is highlighted in the In This Issue feature, p. 265.

Genome-wide Sequencing of Cell-free DNA Enables Detection of Copy-number Alterations in Patients with Cancer Where Tissue Biopsy is Not Feasible.

When tissue biopsy is not medically prudent or tissue is insufficient for molecular testing, alternative methods are needed. Because cell-free DNA (cfDNA) has been shown to provide a representative surrogate for tumor tissue, we sought to evaluate its utility in this clinical scenario. cfDNA was isolated from the plasma of patients and assayed with low-coverage (∼0.3×), genome-wide sequencing. Copy-number alterations (CNA) were identified and characterized using analytic methods originally developed for noninvasive prenatal testing (NIPT) and quantified using the genomic instability number (GIN), a metric that reflects the quantity and magnitude of CNAs across the genome. The technical variability of the GIN was first evaluated in an independent cohort comprising genome-wide sequencing results from 27,754 women who consented to have their samples used for research and whose NIPT results yielded no detected CNAs to establish a detection threshold. Subsequently, cfDNA sequencing data from 96 patients with known cancers but for whom a tissue biopsy could not be obtained are presented. An elevated GIN was detected in 35% of patients and detection rates varied by tumor origin. Collectively, CNAs covered 96.6% of all autosomes. Survival was significantly reduced in patients with an elevated GIN relative to those without. Overall, these data provide a proof of concept for the use of low-coverage, genome-wide sequencing of cfDNA from patients with cancer to obtain relevant molecular information in instances where tissue is difficult to access. These data may ultimately serve as an informative complement to other molecular tests.

A human tissue map of 5-hydroxymethylcytosines exhibits tissue specificity through gene and enhancer modulation.

DNA 5-hydroxymethylcytosine (5hmC) modification is known to be associated with gene transcription and frequently used as a mark to investigate dynamic DNA methylation conversion during mammalian development and in human diseases. However, the lack of genome-wide 5hmC profiles in different human tissue types impedes drawing generalized conclusions about how 5hmC is implicated in transcription activity and tissue specificity. To meet this need, we describe the development of a 5hmC tissue map by characterizing the genomic distributions of 5hmC in 19 human tissues derived from ten organ systems. Subsequent sequencing results enabled the identification of genome-wide 5hmC distributions that uniquely separates samples by tissue type. Further comparison of the 5hmC profiles with transcriptomes and histone modifications revealed that 5hmC is preferentially enriched on tissue-specific gene bodies and enhancers. Taken together, the results provide an extensive 5hmC map across diverse human tissue types that suggests a potential role of 5hmC in tissue-specific development; as well as a resource to facilitate future studies of DNA demethylation in pathogenesis and the development of 5hmC as biomarkers.

Detection of early stage pancreatic cancer using 5-hydroxymethylcytosine signatures in circulating cell free DNA.

Pancreatic cancer is often detected late, when curative therapies are no longer possible. Here, we present non-invasive detection of pancreatic ductal adenocarcinoma (PDAC) by 5-hydroxymethylcytosine (5hmC) changes in circulating cell free DNA from a PDAC cohort (n = 64) in comparison with a non-cancer cohort (n = 243). Differential hydroxymethylation is found in thousands of genes, most significantly in genes related to pancreas development or function (GATA4, GATA6, PROX1, ONECUT1, MEIS2), and cancer pathogenesis (YAP1, TEAD1, PROX1, IGF1). cfDNA hydroxymethylome in PDAC cohort is differentially enriched for genes that are commonly de-regulated in PDAC tumors upon activation of KRAS and inactivation of TP53. Regularized regression models built using 5hmC densities in genes perform with AUC of 0.92 (discovery dataset, n = 79) and 0.92-0.94 (two independent test sets, n = 228). Furthermore, tissue-derived 5hmC features can be used to classify PDAC cfDNA (AUC = 0.88). These findings suggest that 5hmC changes enable classification of PDAC even during early stage disease.

Genome-Wide Sequencing of Cell-Free DNA Identifies Copy-Number Alterations That Can Be Used for Monitoring Response to Immunotherapy in Cancer Patients.

Inhibitors of the PD-1/PD-L1/CTLA-4 immune checkpoint pathway have revolutionized cancer treatment. Indeed, some patients with advanced, refractory malignancies achieve durable responses; however, only a subset of patients benefit, necessitating new biomarkers to predict outcome. Interrogating cell-free DNA (cfDNA) isolated from plasma (liquid biopsy) provides a promising method for monitoring response. We describe the use of low-coverage, genome-wide sequencing of cfDNA, validated extensively for noninvasive prenatal testing, to detect tumor-specific copy-number alterations, and the development of a new metric-the genome instability number (GIN)-to monitor response to these drugs. We demonstrate how the GIN can be used to discriminate clinical response from progression, differentiate progression from pseudoprogression, and identify hyperprogressive disease. Finally, we provide evidence for delayed kinetics in responses to checkpoint inhibitors relative to molecularly targeted therapies. Overall, these data demonstrate a proof of concept for using this method for monitoring treatment outcome in patients with cancer receiving immunotherapy.

Metastatic basal cell carcinoma with amplification of PD-L1: exceptional response to anti-PD1 therapy.

Metastatic basal cell carcinomas are rare malignancies harbouring Hedgehog pathway alterations targetable by SMO antagonists (vismodegib/sonidegib). We describe, for the first time, the molecular genetics and response of a patient with Hedgehog inhibitor-resistant metastatic basal cell carcinoma who achieved rapid tumour regression (ongoing near complete remission at 4 months) with nivolumab (anti-PD1 antibody). He had multiple hallmarks of anti-PD1 responsiveness including high mutational burden (> 50 mutations per megabase; 19 functional alterations in tissue next-generation sequencing (NGS; 315 genes)) as well as PDL1/PDL2/JAK2 amplification (as determined by both tissue NGS and by analysis of plasma-derived cell-free DNA). The latter was performed using technology originally developed for the genome-wide detection of sub-chromosomal copy-number alterations (CNAs) in noninvasive prenatal testing and showed numerous CNAs including amplification of the 9p24.3-9p22.2 region containing PD-L1, PD-L2 and JAK2. Of interest, PD-L1, PD-L2 and JAK2 amplification is a characteristic of Hodgkin lymphoma, which is exquisitely sensitive to nivolumab. In conclusion, selected SMO antagonist-resistant metastatic basal cell carcinomas may respond to nivolumab based on underlying molecular genetic mechanisms that include PD-L1 amplification and high tumour mutational burden.

Using Targeted Sequencing of Paralogous Sequences for Noninvasive Detection of Selected Fetal Aneuploidies.

BACKGROUND: Current methods for noninvasive prenatal testing (NIPT) ascertain fetal aneuploidies using either direct counting measures of DNA fragments from specific genomic regions or relative measures of single nucleotide polymorphism frequencies. Alternatively, the ratios of paralogous sequence pairs were predicted to reflect fetal aneuploidy. We developed a NIPT assay that uses paralog sequences to enable noninvasive detection of fetal trisomy 21 (T21) and trisomy 18 (T18) using cell-free DNA (cfDNA) from maternal plasma. METHODS: A total of 1060 primer pairs were designed to determine fetal aneuploidy status, fetal sex, and fetal fraction. Each library was prepared from cfDNA by coamplifying all 1060 target pairs together in a single reaction well. Products were measured using massively parallel sequencing and deviations from expected paralog ratios were determined based on the read depth from each paralog. RESULTS: We evaluated this assay in a blinded set of 480 cfDNA samples with fetal aneuploidy status determined by the MaterniT21® PLUS assay. Samples were sequenced (mean = 2.3 million reads) with 432 samples returning a result. Using the MaterniT21 PLUS assay for paired plasma aliquots from the same individuals as a reference, all 385 euploid samples, all 31 T21 samples, and 14 of 16 T18 samples were detected with no false positive results observed. CONCLUSIONS: This study introduces a novel NIPT aneuploidy detection approach using targeted sequencing of paralog motifs and establishes proof-of-concept for a potentially low-cost, highly scalable method for the identification of selected fetal aneuploidies with performance and nonreportable rate similar to other published methods.

Additive genetic architecture underlying a rapidly evolving sexual signaling phenotype in the Hawaiian cricket genus Laupala.

Complex, quantitative traits are often the function of the coordinated action of many physically independent genetic factors. Interactive properties of multilocus genotypes, such as epistasis, are thought to be pervasive components of the genetic architecture of complex phenotypes. Here, we utilize a panel of interspecific backcross introgression lines to evaluate the genetic architecture of song variation, a quantitative sexual signaling phenotype, in the Hawaiian swordtail cricket genus Laupala. Allelic effects across five quantitative trait loci are consistent with a purely additive model of gene action, where alleles at multiple loci are found to have fully independent and discrete effects with respect to the sexual signaling phenotype. Whereas a more complex genetic architecture featuring non-additive dominance and epistasis components may constrain potential evolutionary trajectories and reduce the rate of evolutionary change, the polygenic, additive genetic architecture observed for sexual signaling in Laupala should respond rapidly to directional selection pressures and freely move throughout phenotypic space. This classic type I genetic architecture may facilitate the explosive radiation of song variation observed across the Laupala genus.

Widespread genetic linkage of mating signals and preferences in the Hawaiian cricket Laupala.

The evolution of novel sexual communication systems is integral to the process of speciation, as it discourages gene flow between incipient species. Physical linkage between genes underlying male-female communication (i.e. sexual signals and preferences for them) facilitates both rapid and coordinated divergence of sexual communication systems between populations and reduces recombination in the face of occasional hybridization between diverging populations. Despite these ramifications of the genetic architecture of sexual communication for sexual selection and speciation, few studies have examined this relationship empirically. Previous studies of the closely related Hawaiian crickets Laupala paranigra and Laupala kohalensis have indirectly suggested that many of the genes underlying the difference in pulse rate of male song are physically linked with genes underlying the difference in female preference for pulse rate. Using marker-assisted introgression, we moved 'slow pulse rate' alleles from L. paranigra at five known quantitative trait loci (QTL) underlying male pulse rate into the 'fast pulse rate' genetic background of L. kohalensis and assessed the effect of these loci on female preference. An astounding four out of five song QTL predicted the preferences of female fourth-generation backcrosses, providing direct evidence for the extensive genetic linkage of song and preference in one of the fastest diversifying genera currently known.

Mining non-model genomic libraries for microsatellites: BAC versus EST libraries and the generation of allelic richness.

BACKGROUND: Simple sequence repeats (SSRs) are tandemly repeated sequence motifs common in genomic nucleotide sequence that often harbor significant variation in repeat number. Frequently used as molecular markers, SSRs are increasingly identified via in silico approaches. Two common classes of genomic resources that can be mined are bacterial artificial chromosome (BAC) libraries and expressed sequence tag (EST) libraries. RESULTS: 288 SSR loci were screened in the rapidly radiating Hawaiian swordtail cricket genus Laupala. SSRs were more densely distributed and contained longer repeat structures in BAC library-derived sequence than in EST library-derived sequence, although neither repeat density nor length was exceptionally elevated despite the relatively large genome size of Laupala. A non-random distribution favoring AT-rich SSRs was observed. Allelic diversity of SSRs was positively correlated with repeat length and was generally higher in AT-rich repeat motifs. CONCLUSION: The first large-scale survey of Orthopteran SSR allelic diversity is presented. Selection contributes more strongly to the size and density distributions of SSR loci derived from EST library sequence than from BAC library sequence, although all SSRs likely are subject to similar physical and structural constraints, such as slippage of DNA replication machinery, that may generate increased allelic diversity in AT-rich sequence motifs. Although in silico approaches work well for SSR locus identification in both EST and BAC libraries, BAC library sequence and AT-rich repeat motifs are generally superior SSR development resources for most applications.

Genotype-dependent variation of mitochondrial transcriptional profiles in interpopulation hybrids.

Hybridization between populations can disrupt gene expression, frequently resulting in deleterious hybrid phenotypes. Reduced fitness in interpopulation hybrids of the marine copepod Tigriopus californicus has been traced to interactions between the nuclear and mitochondrial genomes. Here, we determine transcript levels of four to six genes involved in the mitochondrial oxidative phosphorylation pathway for a series of parental and inbred hybrid lines using RT-qPCR. Both nuclear and mitochondrial-encoded genes are included in the analysis. Although all genes studied are up-regulated under salinity stress, only expression of genes located on the mtDNA differed among lines. Because mitochondrial genes are transcribed by a dedicated RNA polymerase encoded in the nuclear genome, we compare transcript levels among hybrid lines with different combinations of mitochondrial RNA polymerase and mtDNA genotypes. Lines bearing certain mtDNA-mitochondrial RNA polymerase genotypic combinations show a diminished capacity to up-regulate mitochondrial genes in response to hypoosmotic stress. Effects on the transcriptional profile depend on the specific interpopulation cross and are correlated with viability effects. We hypothesize that disruption of the mitochondrial transcriptional system in F(2) hybrids may play a central role in hybrid breakdown.

Interpopulation hybrid breakdown maps to the mitochondrial genome.

Hybrid breakdown, or outbreeding depression, is the loss of fitness observed in crosses between genetically divergent populations. The role of maternally inherited mitochondrial genomes in hybrid breakdown has not been widely examined. Using laboratory crosses of the marine copepod Tigriopus californicus, we report that the low fitness of F(3) hybrids is completely restored in the offspring of maternal backcrosses, where parental mitochondrial and nuclear genomic combinations are reassembled. Paternal backcrosses, which result in mismatched mitochondrial and nuclear genomes, fail to restore hybrid fitness. These results suggest that fitness loss in T. californicus hybrids is completely attributable to nuclear-mitochondrial genomic interactions. Analyses of ATP synthetic capacity in isolated mitochondria from hybrid and backcross animals found that reduced ATP synthesis in hybrids was also largely restored in backcrosses, again with maternal backcrosses outperforming paternal backcrosses. The strong fitness consequences of nuclear-mitochondrial interactions have important, and often overlooked, implications for evolutionary and conservation biology.

Disruption of mitochondrial function in interpopulation hybrids of Tigriopus californicus.

Electron transport system (ETS) function in mitochondria is essential for the aerobic production of energy. Because ETS function requires extensive interactions between mitochondrial and nuclear gene products, coadaptation between mitochondrial and nuclear genomes may evolve within populations. Hybridization between allopatric populations may then expose functional incompatibilities between genomes that have not coevolved. The intertidal copepod Tigriopus californicus has high levels of nucleotide divergence among populations at mitochondrial loci and suffers F2 hybrid breakdown in interpopulation hybrids. We hypothesize that hybridization results in incompatibilities among subunits in ETS enzyme complexes and that these incompatibilities result in diminished mitochondrial function and fitness. To test this hypothesis, we measured fitness, mitochondrial function, and ETS enzyme activity in inbred recombinant hybrid lines of Tigriopus californicus. We found that (1) both fitness and mitochondrial function are reduced in hybrid lines, (2) only those ETS enzymes with both nuclear and mitochondrial subunits show a loss of activity in hybrid lines, and (3) positive relationships exist between ETS enzyme activity and mitochondrial function and between mitochondrial function and fitness. We also present evidence that hybrid lines harboring mitochondrial DNA (mtDNA) and mitochondrial RNA polymerase (mtRPOL) from the same parental source population have higher fitness than those with mtDNA and mtRPOL from different populations, suggesting that mitochondrial gene regulation may play a role in disruption of mitochondrial performance and fitness of hybrids. These results suggest that disruption of coadaptation between nuclear and mitochondrial genes contributes to the phenomenon of hybrid breakdown.