Searching thousands of genomes to classify somatic and novel structural variants using STIX.

Structural variants are associated with cancers and developmental disorders, but challenges with estimating population frequency remain a barrier to prioritizing mutations over inherited variants. In particular, variability in variant calling heuristics and filtering limits the use of current structural variant catalogs. We present STIX, a method that, instead of relying on variant calls, indexes and searches the raw alignments from thousands of samples to enable more comprehensive allele frequency estimation.

CaBagE: A Cas9-based Background Elimination strategy for targeted, long-read DNA sequencing.

A substantial fraction of the human genome is difficult to interrogate with short-read DNA sequencing technologies due to paralogy, complex haplotype structures, or tandem repeats. Long-read sequencing technologies, such as Oxford Nanopore's MinION, enable direct measurement of complex loci without introducing many of the biases inherent to short-read methods, though they suffer from relatively lower throughput. This limitation has motivated recent efforts to develop amplification-free strategies to target and enrich loci of interest for subsequent sequencing with long reads. Here, we present CaBagE, a method for target enrichment that is efficient and useful for sequencing large, structurally complex targets. The CaBagE method leverages the stable binding of Cas9 to its DNA target to protect desired fragments from digestion with exonuclease. Enriched DNA fragments are then sequenced with Oxford Nanopore's MinION long-read sequencing technology. Enrichment with CaBagE resulted in a median of 116X coverage (range 39-416) of target loci when tested on five genomic targets ranging from 4-20kb in length using healthy donor DNA. Four cancer gene targets were enriched in a single reaction and multiplexed on a single MinION flow cell. We further demonstrate the utility of CaBagE in two ALS patients with C9orf72 short tandem repeat expansions to produce genotype estimates commensurate with genotypes derived from repeat-primed PCR for each individual. With CaBagE there is a physical enrichment of on-target DNA in a given sample prior to sequencing. This feature allows adaptability across sequencing platforms and potential use as an enrichment strategy for applications beyond sequencing. CaBagE is a rapid enrichment method that can illuminate regions of the 'hidden genome' underlying human disease.

Somalier: rapid relatedness estimation for cancer and germline studies using efficient genome sketches.

BACKGROUND: When interpreting sequencing data from multiple spatial or longitudinal biopsies, detecting sample mix-ups is essential, yet more difficult than in studies of germline variation. In most genomic studies of tumors, genetic variation is detected through pairwise comparisons of the tumor and a matched normal tissue from the sample donor. In many cases, only somatic variants are reported, which hinders the use of existing tools that detect sample swaps solely based on genotypes of inherited variants. To address this problem, we have developed Somalier, a tool that operates directly on alignments and does not require jointly called germline variants. Instead, Somalier extracts a small sketch of informative genetic variation for each sample. Sketches from hundreds of germline or somatic samples can then be compared in under a second, making Somalier a useful tool for measuring relatedness in large cohorts. Somalier produces both text output and an interactive visual report that facilitates the detection and correction of sample swaps using multiple relatedness metrics. RESULTS: We introduce the tool and demonstrate its utility on a cohort of five glioma samples each with a normal, tumor, and cell-free DNA sample. Applying Somalier to high-coverage sequence data from the 1000 Genomes Project also identifies several related samples. We also demonstrate that it can distinguish pairs of whole-genome and RNA-seq samples from the same individuals in the Genotype-Tissue Expression (GTEx) project. CONCLUSIONS: Somalier is a tool that can rapidly evaluate relatedness from sequencing data. It can be applied to diverse sequencing data types and genome builds and is available under an MIT license at github.com/brentp/somalier .

Publisher Correction: Combating subclonal evolution of resistant cancer phenotypes.

The originally published version of this Article contained an error in Figure 4. In panel a, grey boxes surrounding the subclones associated with patients #2 and #4 obscured adjacent portions of the heatmap. This error has now been corrected in both the PDF and HTML versions of the Article.

GIGGLE: a search engine for large-scale integrated genome analysis.

GIGGLE is a genomics search engine that identifies and ranks the significance of genomic loci shared between query features and thousands of genome interval files. GIGGLE (https://github.com/ryanlayer/giggle) scales to billions of intervals and is over three orders of magnitude faster than existing methods. Its speed extends the accessibility and utility of resources such as ENCODE, Roadmap Epigenomics, and GTEx by facilitating data integration and hypothesis generation.

Methylene-tetrahydrofolate reductase contributes to allergic airway disease.

RATIONALE: Environmental exposures strongly influence the development and progression of asthma. We have previously demonstrated that mice exposed to a diet enriched with methyl donors during vulnerable periods of fetal development can enhance the heritable risk of allergic airway disease through epigenetic changes. There is conflicting evidence on the role of folate (one of the primary methyl donors) in modifying allergic airway disease. OBJECTIVES: We hypothesized that blocking folate metabolism through the loss of methylene-tetrahydrofolate reductase (Mthfr) activity would reduce the allergic airway disease phenotype through epigenetic mechanisms. METHODS: Allergic airway disease was induced in C57BL/6 and C57BL/6Mthfr-/- mice through house dust mite (HDM) exposure. Airway inflammation and airway hyperresponsiveness (AHR) were measured between the two groups. Gene expression and methylation profiles were generated for whole lung tissue. Disease and molecular outcomes were evaluated in C57BL/6 and C57BL/6Mthfr-/- mice supplemented with betaine. MEASUREMENTS AND MAIN RESULTS: Loss of Mthfr alters single carbon metabolite levels in the lung and serum including elevated homocysteine and cystathionine and reduced methionine. HDM-treated C57BL/6Mthfr-/- mice demonstrated significantly less airway hyperreactivity (AHR) compared to HDM-treated C57BL/6 mice. Furthermore, HDM-treated C57BL/6Mthfr-/- mice compared to HDM-treated C57BL/6 mice have reduced whole lung lavage (WLL) cellularity, eosinophilia, and Il-4/Il-5 cytokine concentrations. Betaine supplementation reversed parts of the HDM-induced allergic airway disease that are modified by Mthfr loss. 737 genes are differentially expressed and 146 regions are differentially methylated in lung tissue from HDM-treated C57BL/6Mthfr-/- mice and HDM-treated C57BL/6 mice. Additionally, analysis of methylation/expression relationships identified 503 significant correlations. CONCLUSION: Collectively, these findings indicate that the loss of folate as a methyl donor is a modifier of allergic airway disease, and that epigenetic and expression changes correlate with this modification. Further investigation into the mechanisms that drive this observation is warranted.

Combating subclonal evolution of resistant cancer phenotypes.

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves.

Regulation of MUC5B Expression in Idiopathic Pulmonary Fibrosis.

The gain-of-function mucin 5B (MUC5B) promoter variant, rs35705950, confers the largest risk, genetic or otherwise, for the development of idiopathic pulmonary fibrosis; however, the mechanisms underlying the regulation of MUC5B expression have yet to be elucidated. Here, we identify a critical regulatory domain that contains the MUC5B promoter variant and has a highly conserved forkhead box protein A2 (FOXA2) binding motif. This region is differentially methylated in association with idiopathic pulmonary fibrosis, MUC5B expression, and rs35705950. In addition, we show that this locus binds FOXA2 dynamically, and that binding of FOXA2 is necessary for enhanced expression of MUC5B. In aggregate, our findings identify novel targets to regulate the expression of MUC5B.

Novel Innate Immune Genes Regulating the Macrophage Response to Gram Positive Bacteria.

Host variation in Toll-like receptors and other innate immune signaling molecules alters infection susceptibility. However, only a portion of the variability observed in the innate immune response is accounted for by known genes in these pathways. Thus, the identification of additional genes that regulate the response to Gram positive bacteria is warranted. Bone marrow-derived macrophages (BMMs) from 43 inbred mouse strains were stimulated with lipotechoic acid (LTA), a major component of the Gram positive bacterial cell wall. Concentrations of the proinflammatory cytokines IL-6, IL-12, and TNF-α were measured. In silico whole genome association (WGA) mapping was performed using cytokine responses followed by network analysis to prioritize candidate genes. To determine which candidate genes could be responsible for regulating the LTA response, candidate genes were inhibited using RNA interference (RNAi) and were overexpressed in RAW264.7 macrophages. BMMs from Bdkrb1-deficient mice were used to assess the effect of Bdkrb1 gene deletion on the response to LTA, heat-killed Streptococcus pneumoniae, and heat-killed Staphylococcus aureus WGA mapping identified 117 loci: IL-6 analysis yielded 20 loci (average locus size = 0.133 Mb; 18 genes), IL-12 analysis produced 5 loci (0.201 Mb average; 7 genes), and TNF-α analysis yielded 92 loci (0.464 Mb average; 186 genes of which 46 were prioritized by network analysis). The follow-up small interfering RNA screen of 71 target genes identified four genes (Bdkrb1, Blnk, Fbxo17, and Nkx6-1) whose inhibition resulted in significantly reduced cytokine production following LTA stimulation. Overexpression of these four genes resulted in significantly increased cytokine production in response to LTA. Bdkrb1-deficient macrophages were less responsive to LTA and heat-killed S. aureus, validating the genetic and RNAi approach to identify novel regulators of the response to LTA. We have identified four innate immune response genes that may contribute to Gram positive bacterial susceptibility.

Genome-wide imputation study identifies novel HLA locus for pulmonary fibrosis and potential role for auto-immunity in fibrotic idiopathic interstitial pneumonia.

BACKGROUND: Fibrotic idiopathic interstitial pneumonias (fIIP) are a group of fatal lung diseases with largely unknown etiology and without definitive treatment other than lung transplant to prolong life. There is strong evidence for the importance of both rare and common genetic risk alleles in familial and sporadic disease. We have previously used genome-wide single nucleotide polymorphism data to identify 10 risk loci for fIIP. Here we extend that work to imputed genome-wide genotypes and conduct new RNA sequencing studies of lung tissue to identify and characterize new fIIP risk loci. RESULTS: We performed genome-wide genotype imputation association analyses in 1616 non-Hispanic white (NHW) cases and 4683 NHW controls followed by validation and replication (878 cases, 2017 controls) genotyping and targeted gene expression in lung tissue. Following meta-analysis of the discovery and replication populations, we identified a novel fIIP locus in the HLA region of chromosome 6 (rs7887 P meta = 3.7 × 10(-09)). Imputation of classic HLA alleles identified two in high linkage disequilibrium that are associated with fIIP (DRB1*15:01 P = 1.3 × 10(-7) and DQB1*06:02 P = 6.1 × 10(-8)). Targeted RNA-sequencing of the HLA locus identified 21 genes differentially expressed between fibrotic and control lung tissue (Q < 0.001), many of which are involved in immune and inflammatory response regulation. In addition, the putative risk alleles, DRB1*15:01 and DQB1*06:02, are associated with expression of the DQB1 gene among fIIP cases (Q < 1 × 10(-16)). CONCLUSIONS: We have identified a genome-wide significant association between the HLA region and fIIP. Two HLA alleles are associated with fIIP and affect expression of HLA genes in lung tissue, indicating that the potential genetic risk due to HLA alleles may involve gene regulation in addition to altered protein structure. These studies reveal the importance of the HLA region for risk of fIIP and a basis for the potential etiologic role of auto-immunity in fIIP.

Desmoplakin Variants Are Associated with Idiopathic Pulmonary Fibrosis.

RATIONALE: Sequence variation, methylation differences, and transcriptional changes in desmoplakin (DSP) have been observed in patients with idiopathic pulmonary fibrosis (IPF). OBJECTIVES: To identify novel variants in DSP associated with IPF and to characterize the relationship of these IPF sequence variants with DSP gene expression in human lung. METHODS: A chromosome 6 locus (7,370,061-7,606,946) was sequenced in 230 subjects with IPF and 228 control subjects. Validation genotyping of disease-associated variants was conducted in 936 subjects with IPF and 936 control subjects. DSP gene expression was measured in lung tissue from 334 subjects with IPF and 201 control subjects. MEASUREMENTS AND MAIN RESULTS: We identified 23 sequence variants in the chromosome 6 locus associated with IPF. Genotyping of selected variants in our validation cohort revealed that noncoding intron 1 variant rs2744371 (odds ratio = 0.77, 95% confidence interval [CI] = 0.66-0.91, P = 0.002) is protective for IPF, and a previously described IPF-associated intron 5 variant (rs2076295) is associated with increased risk of IPF (odds ratio = 1.36, 95% CI = 1.19-1.56, P < 0.001) after controlling for sex and age. DSP expression is 2.3-fold increased (95% CI = 1.91-2.71) in IPF lung tissue (P < 0.0001). Only the minor allele at rs2076295 is associated with decreased DSP expression (P = 0.001). Staining of fibrotic and normal human lung tissue localized DSP to airway epithelia. CONCLUSIONS: Sequence variants in DSP are associated with IPF, and rs2076295 genotype is associated with differential expression of DSP in the lung. DSP expression is increased in IPF lung and concentrated in the airway epithelia, suggesting a potential role for DSP in the pathogenesis of IPF.

In Utero Cigarette Smoke Affects Allergic Airway Disease But Does Not Alter the Lung Methylome.

Prenatal and postnatal cigarette smoke exposure enhances the risk of developing asthma. Despite this as well as other smoking related risks, 11% of women still smoke during pregnancy. We hypothesized that cigarette smoke exposure during prenatal development generates long lasting differential methylation altering transcriptional activity that correlates with disease. In a house dust mite (HDM) model of allergic airway disease, we measured airway hyperresponsiveness (AHR) and airway inflammation between mice exposed prenatally to cigarette smoke (CS) or filtered air (FA). DNA methylation and gene expression were then measured in lung tissue. We demonstrate that HDM-treated CS mice develop a more severe allergic airway disease compared to HDM-treated FA mice including increased AHR and airway inflammation. While DNA methylation changes between the two HDM-treated groups failed to reach genome-wide significance, 99 DMRs had an uncorrected p-value < 0.001. 6 of these 99 DMRs were selected for validation, based on the immune function of adjacent genes, and only 2 of the 6 DMRs confirmed the bisulfite sequencing data. Additionally, genes near these 6 DMRs (Lif, Il27ra, Tle4, Ptk7, Nfatc2, and Runx3) are differentially expressed between HDM-treated CS mice and HDM-treated FA mice. Our findings confirm that prenatal exposure to cigarette smoke is sufficient to modify allergic airway disease; however, it is unlikely that specific methylation changes account for the exposure-response relationship. These findings highlight the important role in utero cigarette smoke exposure plays in the development of allergic airway disease.

Signatures of accelerated somatic evolution in gene promoters in multiple cancer types.

Cancer-associated somatic mutations outside protein-coding regions remain largely unexplored. Analyses of the TERT locus have indicated that non-coding regulatory mutations can be more frequent than previously suspected and play important roles in oncogenesis. Using a computational method called SASE-hunter, developed here, we identified a novel signature of accelerated somatic evolution (SASE) marked by a significant excess of somatic mutations localized in a genomic locus, and prioritized those loci that carried the signature in multiple cancer patients. Interestingly, even when an affected locus carried the signature in multiple individuals, the mutations contributing to SASE themselves were rarely recurrent at the base-pair resolution. In a pan-cancer analysis of 906 samples from 12 tumor types, we detected SASE in the promoters of several genes, including known cancer genes such as MYC, BCL2, RBM5 and WWOX. Nucleotide substitution patterns consistent with oxidative DNA damage and local somatic hypermutation appeared to contribute to this signature in selected gene promoters (e.g. MYC). SASEs in selected cancer gene promoters were associated with over-expression, and also correlated with the age of onset of cancer, aggressiveness of the disease and survival. Taken together, our work detects a hitherto under-appreciated and clinically important class of regulatory changes in cancer genomes.

DNA methylation and childhood asthma in the inner city.

BACKGROUND: Epigenetic marks are heritable, influenced by the environment, direct the maturation of T lymphocytes, and in mice enhance the development of allergic airway disease. Thus it is important to define epigenetic alterations in asthmatic populations. OBJECTIVE: We hypothesize that epigenetic alterations in circulating PBMCs are associated with allergic asthma. METHODS: We compared DNA methylation patterns and gene expression in inner-city children with persistent atopic asthma versus healthy control subjects by using DNA and RNA from PBMCs. Results were validated in an independent population of asthmatic patients. RESULTS: Comparing asthmatic patients (n = 97) with control subjects (n = 97), we identified 81 regions that were differentially methylated. Several immune genes were hypomethylated in asthma, including IL13, RUNX3, and specific genes relevant to T lymphocytes (TIGIT). Among asthmatic patients, 11 differentially methylated regions were associated with higher serum IgE concentrations, and 16 were associated with percent predicted FEV1. Hypomethylated and hypermethylated regions were associated with increased and decreased gene expression, respectively (P < 6 × 10(-12) for asthma and P < .01 for IgE). We further explored the relationship between DNA methylation and gene expression using an integrative analysis and identified additional candidates relevant to asthma (IL4 and ST2). Methylation marks involved in T-cell maturation (RUNX3), TH2 immunity (IL4), and oxidative stress (catalase) were validated in an independent asthmatic cohort of children living in the inner city. CONCLUSIONS: Our results demonstrate that DNA methylation marks in specific gene loci are associated with asthma and suggest that epigenetic changes might play a role in establishing the immune phenotype associated with asthma.

Relationship of DNA methylation and gene expression in idiopathic pulmonary fibrosis.

RATIONALE: Idiopathic pulmonary fibrosis (IPF) is an untreatable and often fatal lung disease that is increasing in prevalence and is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control gene expression and are likely to regulate the IPF transcriptome. OBJECTIVES: To identify methylation marks that modify gene expression in IPF lung. METHODS: We assessed DNA methylation (comprehensive high-throughput arrays for relative methylation arrays [CHARM]) and gene expression (Agilent gene expression arrays) in 94 patients with IPF and 67 control subjects, and performed integrative genomic analyses to define methylation-gene expression relationships in IPF lung. We validated methylation changes by a targeted analysis (Epityper), and performed functional validation of one of the genes identified by our analysis. MEASUREMENTS AND MAIN RESULTS: We identified 2,130 differentially methylated regions (DMRs; <5% false discovery rate), of which 738 are associated with significant changes in gene expression and enriched for expected inverse relationship between methylation and expression (P < 2.2 × 10(-16)). We validated 13/15 DMRs by targeted analysis of methylation. Methylation-expression quantitative trait loci (methyl-eQTL) identified methylation marks that control cis and trans gene expression, with an enrichment for cis relationships (P < 2.2 × 10(-16)). We found five trans methyl-eQTLs where a methylation change at a single DMR is associated with transcriptional changes in a substantial number of genes; four of these DMRs are near transcription factors (castor zinc finger 1 [CASZ1], FOXC1, MXD4, and ZDHHC4). We studied the in vitro effects of change in CASZ1 expression and validated its role in regulation of target genes in the methyl-eQTL. CONCLUSIONS: These results suggest that DNA methylation may be involved in the pathogenesis of IPF.

Loss of heterozygosity preferentially occurs in early replicating regions in cancer genomes.

Erroneous repair of DNA double-strand breaks by homologous recombination (HR) leads to loss of heterozygosity (LOH). Analysing 22,392 and 74,415 LOH events in 363 glioblastoma and 513 ovarian cancer samples, respectively, and using three different metrics, we report that LOH selectively occurs in early replicating regions; this pattern differs from the trends for point mutations and somatic deletions, which are biased toward late replicating regions. Our results are independent of BRCA1 and BRCA2 mutation status. The LOH events are significantly clustered near RNA polII-bound transcription start sites, consistent with the reports that slow replication near paused RNA polII might initiate HR-mediated repair. The frequency of LOH events is higher in the chromosomes with shorter inter-homolog distance inside the nucleus. We propose that during early replication, HR-mediated rescue of replication near paused RNA polII using homologous chromosomes as template leads to LOH. The difference in the preference for replication timing between different classes of genomic alterations in cancer genomes also provokes a testable hypothesis that replicating cells show changing preference between various DNA repair pathways, which have different levels of efficiency and fidelity, as the replication progresses.

Copy neutral loss of heterozygosity is more frequent in older ovarian cancer patients.

Loss of heterozygosity (LOH) is a common type of genomic alterations in ovarian cancer. Analyzing 74,415 copy neutral LOH events in 513 serous ovarian adenocarcinomas samples from the Cancer Genome Atlas, we report that the frequency of LOH events increases with age. Similar trend is observed for LOH involving chromosome 17, which is frequently implicated in ovarian cancer. The results are consistent when we analyze data from the Boston high-grade serous cancer cohort. We further show that germ line and somatic mutations in BRCA1 (in chromosome 17) and BRCA2 (in chromosome 13) loci are not necessary to establish the pattern. We also report significant age-related changes in expression patterns for several genes in the homologous recombination (HR) pathway, such as BRCA1, RAD50, RAD52, XRCC2, XRCC3, and MRE11A in these patient samples. Furthermore, we develop a metric for pathway-level imbalance, and show that increased imbalance in the HR pathway, i.e., increase in expression of some HR genes and decrease in expression of others, is common and correlates significantly with the frequency of LOH events in the patient samples. Taken together, it is highly likely that aging and deregulation of HR pathway contribute to the increased incidence of copy-neutral LOH in ovarian cancer patients.

Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis.

We performed a genome-wide association study of non-Hispanic, white individuals with fibrotic idiopathic interstitial pneumonias (IIPs; n = 1,616) and controls (n = 4,683), with follow-up replication analyses in 876 cases and 1,890 controls. We confirmed association with TERT at 5p15, MUC5B at 11p15 and the 3q26 region near TERC, and we identified seven newly associated loci (Pmeta = 2.4 × 10(-8) to 1.1 × 10(-19)), including FAM13A (4q22), DSP (6p24), OBFC1 (10q24), ATP11A (13q34), DPP9 (19p13) and chromosomal regions 7q22 and 15q14-15. Our results suggest that genes involved in host defense, cell-cell adhesion and DNA repair contribute to risk of fibrotic IIPs.

Altered patterns of fractionation and exon deletions in Brassica rapa support a two-step model of paleohexaploidy.

The genome sequence of the paleohexaploid Brassica rapa shows that fractionation is biased among the three subgenomes and that the least fractionated subgenome has approximately twice as many orthologs as its close (and relatively unduplicated) relative Arabidopsis than had either of the other two subgenomes. One evolutionary scenario is that the two subgenomes with heavy gene losses (I and II) were in the same nucleus for a longer period of time than the third subgenome (III) with the fewest gene losses. This "two-step" hypothesis is essentially the same as that proposed previously for the eudicot paleohexaploidy; however, the more recent nature of the B. rapa paleohexaploidy makes this model more testable. We found that subgenome II suffered recent small deletions within exons more frequently than subgenome I, as would be expected if the genes in subgenome I had already been near maximally fractionated before subgenome III was introduced. We observed that some sequences, before these deletions, were flanked by short direct repeats, a unique signature of intrachromosomal illegitimate recombination. We also found, through simulations, that short--single or two-gene--deletions appear to dominate the fractionation patterns in B. rapa. We conclude that the observed patterns of the triplicated regions in the Brassica genome are best explained by a two-step fractionation model. The triplication and subsequent mode of fractionation could influence the potential to generate morphological diversity--a hallmark of the Brassica genus.

Dose-sensitivity, conserved non-coding sequences, and duplicate gene retention through multiple tetraploidies in the grasses.

Whole genome duplications, or tetraploidies, are an important source of increased gene content. Following whole genome duplication, duplicate copies of many genes are lost from the genome. This loss of genes is biased both in the classes of genes deleted and the subgenome from which they are lost. Many or all classes are genes preferentially retained as duplicate copies are engaged in dose sensitive protein-protein interactions, such that deletion of any one duplicate upsets the status quo of subunit concentrations, and presumably lowers fitness as a result. Transcription factors are also preferentially retained following every whole genome duplications studied. This has been explained as a consequence of protein-protein interactions, just as for other highly retained classes of genes. We show that the quantity of conserved noncoding sequences (CNSs) associated with genes predicts the likelihood of their retention as duplicate pairs following whole genome duplication. As many CNSs likely represent binding sites for transcriptional regulators, we propose that the likelihood of gene retention following tetraploidy may also be influenced by dose-sensitive protein-DNA interactions between the regulatory regions of CNS-rich genes - nicknamed bigfoot genes - and the proteins that bind to them. Using grass genomes, we show that differential loss of CNSs from one member of a pair following the pre-grass tetraploidy reduces its chance of retention in the subsequent maize lineage tetraploidy.