Genomic Footprinting Research Articles

Uncovering the genetic signature of quantitative trait evolution with replicated time series data.

The genetic architecture of adaptation in natural populations has not yet been resolved: it is not clear to what extent the spread of beneficial mutations (selective sweeps) or the response of many quantitative trait loci drive adaptation to environmental changes. Although much attention has been given to the genomic footprint of selective sweeps, the importance of selection on quantitative traits is still not well studied, as the associated genomic signature is extremely difficult to detect. We propose 'Evolve and Resequence' as a promising tool, to study polygenic adaptation of quantitative traits in evolving populations. Simulating replicated time series data we show that adaptation to a new intermediate trait optimum has three characteristic phases that are reflected on the genomic level: (1) directional frequency changes towards the new trait optimum, (2) plateauing of allele frequencies when the new trait optimum has been reached and (3) subsequent divergence between replicated trajectories ultimately leading to the loss or fixation of alleles while the trait value does not change. We explore these 3 phase characteristics for relevant population genetic parameters to provide expectations for various experimental evolution designs. Remarkably, over a broad range of parameters the trajectories of selected alleles display a pattern across replicates, which differs both from neutrality and directional selection. We conclude that replicated time series data from experimental evolution studies provide a promising framework to study polygenic adaptation from whole-genome population genetics data.

Genomic analysis reveals selection in Chinese native black pig.

Identification of genomic signatures that help reveal mechanisms underlying desirable traits in domesticated pigs is of significant biological, agricultural and medical importance. To identify the genomic footprints left by selection during domestication of the Enshi black pig, a typical native and meat-lard breed in China, we generated about 72-fold coverage of the pig genome using pools of genomic DNA representing three different populations of Enshi black pigs from three different locations. Combining this data with the available whole genomes of 13 Chinese wild boars, we identified 417 protein-coding genes embedded in the selected regions of Enshi black pigs. These genes are mainly involved in developmental and metabolic processes, response to stimulus, and other biological processes. Signatures of selection were detected in genes involved in body size and immunity (RPS10 and VASN), lipid metabolism (GSK3), male fertility (INSL6) and developmental processes (TBX19). These findings provide a window into the potential genetic mechanism underlying development of desirable phenotypes in Enshi black pigs during domestication and subsequent artificial selection. Thus, our results illustrate how domestication has shaped patterns of genetic variation in Enshi black pigs and provide valuable genetic resources that enable effective use of pigs in agricultural production.

Genome scan reveals selection acting on genes linked to stress response in wild pearl millet.

Uncovering genomic regions involved in adaption is a major goal in evolutionary biology. High-throughput sequencing now makes it possible to tackle this challenge in nonmodel species. Yet, despite the increasing number of methods targeted to specifically detect genomic footprints of selection, the complex demography of natural populations often causes high rates of false positive in gene discoveries. The aim of this study was to identify climate adaptations in wild pearl millet populations, Cenchrus americanus ssp. monodii. We focused on two climate gradients, one in Mali and one in Niger. We used a two-step strategy to limit false-positive outliers. First, we considered gradients as biological replicates and performed RNA sequencing of four populations at the extremities. We combined four methods-three based on differentiation among populations and one based on diversity patterns within populations-to identify outlier SNPs from a set of 87218 high-quality SNPs. Among 11155 contigs of pearl millet reference transcriptome, 540 exhibited selection signals as evidenced by at least one of the four methods. In a second step, we genotyped 762 samples in 11 additional populations distributed along the gradients using SNPs from the detected contigs and random SNPs as control. We further assessed selection on this large data set using a differentiation-based method and a method based on correlations with environmental variables based. Four contigs displayed consistent signatures between the four extreme and 11 additional populations, two of which were linked to abiotic and biotic stress responses.

Adjuvant therapy for resected biliary tract cancer: a review

Surgical resection remains the most important treatment for patients with biliary tract cancer (BTC), but despite radical surgical techniques many patients ultimately develop recurrent disease. BTC encompasses several distinct disease entities-intrahepatic, perihilar and distal bile duct cholangiocarcinoma as well as gallbladder cancer. These tumors are histologically similar, but have different etiologies and recent information regarding their genomic footprint has questioned their biological similarity. Surgical approaches are also necessarily varied based on the site of the tumor. Due to the poor survival rates seen in this disease, there has been significant effort to investigate chemotherapy and radiotherapy as adjuvants in patients whose disease has been successfully resected. The majority of the published evidence supporting this treatment relies on retrospective series or limited prospective studies, making interpretation difficult and complicating treatment decisions. This review summarizes the data regarding these adjunctive therapies.

Chromatin Kinases Act on Transcription Factors and Histone Tails in Regulation of Inducible Transcription

The inflammatory response requires coordinated activation of both transcription factors and chromatin to induce transcription for defense against pathogens and environmental insults. We sought to elucidate the connections between inflammatory signaling pathways and chromatin through genomic footprinting of kinase activity and unbiased identification of prominent histone phosphorylation events.We identified H3 serine 28 phosphorylation (H3S28ph) as the principal stimulation-dependent histone modification and observed its enrichment at induced genes in mouse macrophages stimulated with bacterial lipopolysaccharide. Using pharmacological and genetic approaches, we identified mitogen- and stress-activated protein kinases (MSKs) as primary mediators of H3S28ph in macrophages. Cell-free transcription assays demonstrated that H3S28ph directly promotes p300/CBP-dependent transcription. Further, MSKs can activate both signal-responsive transcription factors and the chromatin template with additive effects on transcription. Specific inhibition of MSKs in macrophages selectively reduced transcription of stimulation-induced genes. Our results suggest that MSKs incorporate upstream signaling inputs and control multiple downstream regulators of inducible transcription.

Preservation of epithelial cell barrier function and muted inflammation in resistance to allergic rhinoconjunctivitis from house dust mite challenge

An emerging paradigm holds that resistance to the development of allergic diseases, including allergic rhinoconjunctivitis, relates to an intact epithelial/epidermal barrier during early childhood. Conceivably, the immunologic and genomic footprint of this resistance is preserved in nonatopic, nonallergic adults and is unmasked during exposure to an aeroallergen. The aim of this study was to obtain direct support of the epithelial/epidermal barrier model for allergic rhinoconjunctivitis. Twenty-three adults allergic to house dust mites (HDMs) (M+) and 15 nonsensitive, nonallergic (M-) participants completed 3-hour exposures to aerosolized HDM (Dermatophagoides pteronyssinus) powder on 4 consecutive days in an allergen challenge chamber. We analyzed: (1) peripheral blood leukocyte levels and immune responses; and (2) RNA sequencing-derived expression profiles of nasal cells, before and after HDM exposure. On HDM challenge: (1) only M+ persons developed allergic rhinoconjunctivitis symptoms; and (2) peripheral blood leukocyte levels/responses and gene expression patterns in nasal cells were largely concordant between M+ and M- participants; gross differences in these parameters were not observed at baseline (pre-exposure). Two key differences were observed. First, peripheral blood CD4+ and CD8+ T-cell activation levels initially decreased in M- participants versus increased in M+ participants. Second, in M- compared with M+ participants, genes that promoted epidermal/epithelial barrier function (eg, filament-aggregating protein [filaggrin]) versus inflammation (eg, chemokines) and innate immunity (interferon) were upregulated versus muted, respectively. An imprint of resistance to HDM challenge in nonatopic, nonallergic adults was muted T-cell activation in the peripheral blood and inflammatory response in the nasal compartment, coupled with upregulation of genes that promote epidermal/epithelial cell barrier function.

Pioneer factors and ATP-dependent chromatin remodeling factors interact dynamically: A new perspective: Multiple transcription factors can effect chromatin pioneer functions through dynamic interactions with ATP-dependent chromatin remodeling factors.

Transcription factor (TF) signaling regulates gene transcription and requires a complex network of proteins. This network includes co-activators, co-repressors, multiple TFs, histone-modifying complexes, and the basal transcription machinery. It has been widely appreciated that pioneer factors, such as FoxA1 and GATA1, play an important role in opening closed chromatin regions, thereby allowing binding of a secondary factor. In this review we will focus on a newly proposed model wherein multiple TFs, such as steroid receptors (SRs), can function in a pioneering role. This model, termed dynamic assisted loading, integrates data from widely divergent methodologies, including genome wide ChIP-Seq, digital genomic footprinting, DHS-Seq, live cell protein dynamics, and biochemical studies of ATP-dependent remodeling complexes, to present a real time view of TF chromatin interactions. Under this view, many TFs can act as initiating factors for chromatin landscape programming. Furthermore, enhancer and promoter states are more accurately described as energy-dependent, non-equilibrium steady states.

Abstract 1987: Recurrent chromosomal translocations mediate reprogramming of Myb dependent regulatory networks in adenoid cystic carcinoma

Abstract Adenoid Cystic Carcinoma (ACC), the second most common malignancy of salivary glands, is a rare tumor with bleak prognosis for which therapeutic targets are unavailable. Recently, we used RNA-sequencing (RNA-Seq) to analyze low-quality RNA from archival, formaldehyde-fixed, paraffin-embedded samples (1). In addition to detecting the most common ACC translocation, t(6;9) fusing the MYB proto-oncogene to NFIB, we also detected previously unknown t(8;9) and t(8;14) translocations fusing the MYBL1 gene to the NFIB and RAD51B genes, respectively. Mutated versions of c-Myb are potent oncoproteins, and Myb proteins are over- or highly-expressed in a wide variety of human tumors. Interestingly, ACC tumors with MYB and MYBL1 translocations displayed similar gene expression profiles, and the combined MYB and MYBL1 expression correlated with outcome, suggesting that the related Myb proteins are interchangeable oncogenic drivers in ACC. While RNA-Seq provides information about gene fusions, alternative RNA splicing and gene expression signatures, this type of data does not address how translocations alter transcription factor function. Previous studies in our lab, and others, have demonstrated that Myb proteins contain C-terminal regulatory domains and are known to interact with several proteins via C-terminal motifs (2). These regulatory and interaction sites are lost upon chromosomal translocation in ACC; however, the downstream affect of this loss is not understood. Do the translocations merely relieve Myb transcription factors of repressive regulatory signals, thereby globally increasing Myb-dependent gene expression or, does the altered protein allow novel protein-protein interactions which re-target Myb transcription factors? Furthermore, the most frequent translocations involve the relocation of a large enhancer regions. What role does the recruitment of novel enhancer regions play in up-regulating Myb protein expression? Using a combination of ChIP-Seq and digital genomic footprinting, we are attempting to understand the mechanism of regulatory rewiring observed in ACC tumors.

Abstract 506: Post-surgical resection monitoring in early stage colorectal carcinoma patients using a circulating cell-free DNA assay with ultra-high accuracy and specificity

Abstract Analysis of cell-free circulating tumor DNA (ctDNA) by next-generation sequencing (NGS) allows non-invasive real-time profiling of actionable genomic alterations. Liquid biopsy provides an option for disease monitoring in early stage cancer patients post surgical resection, with a potential to aid in adjuvant decision making. However, to be applicable, tests must cover a broad enough genomic footprint to not require a priori knowledge of mutations, have high specificity, and sensitivity higher than conventional methods. NGS is necessary, since inactivating mutations are the most common alteration type in many common cancer types such as colorectal carcinomas (CRC). Here we present a highly efficient and specific NGS assay for detection of ctDNA in early stage cancer patients, capable of detecting single molecule variants across a 12 kb gene panel with an analytical sensitivity of >0.02% for single nucleotide variants (SNVs) and indels. This panel was applied to a clinical study involving 14 early stage (II/III) CRC patients with both pre- and post-op blood draws (up to 7 days post surgery). A subset (6 patients) also had tumor samples collected at the time of the surgical resection of the tumor. Overall, the detection rate of ctDNA in pre-op blood draws was 93%. In the post-op blood draws ctDNA was detectable in 43% of cases. The estimated average minor allele frequency (MAF) is 0.58% (± 0.82%) in pre-op, 0.18% (±0.21%) in post-op, and 40% (±18%) in tumor samples. When tumor tissue was available and used as a reference, the clinical sensitivity, specificity, and accuracy in pre-op blood samples were 83%, 99.995%, and 99.99%, respectively. SNVs with MAF as low as 0.04% were confirmed in tissue data. The clinical specificity of variants detected in post-op blood samples using pre-op samples as the reference is 99.996%. Cohort expansion to 50 patients and follow-up for clinical recurrence in both cohorts is ongoing. In conclusion, we have developed an assay with ultra-high accuracy and specificity, for the detection of ctDNA in early stage CRC patients that is capable of detecting alterations present in the tumor post-surgical resection. This technology allows for a promising non-invasive route for molecular monitoring of residual disease post surgery and for early detection of relapse compared to traditional methodologies. Citation Format: Stefanie A. Mortimer, Katharine Dilger, Stephen Fairclough, Diana Abdueva, Darya Chudova, Ankit Sarin, Jim Leng, Jeeyun Lee, Helmy Eltoukhy, AmirAli Talasaz. Post-surgical resection monitoring in early stage colorectal carcinoma patients using a circulating cell-free DNA assay with ultra-high accuracy and specificity. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 506.

8-Oxo-7,8-dihydro-2'-deoxyguanosine and other lesions along the coding strand of the exon 5 of the tumour suppressor gene P53 in a breast cancer case-control study.

The next-generation sequencing studies of breast cancer have reported that the tumour suppressor P53 (TP53) gene is mutated in more than 40% of the tumours. We studied the levels of oxidative lesions, including 8-oxo-7,8-dihydro-2′-deoxyguanosine (8-oxodG), along the coding strand of the exon 5 in breast cancer patients as well as in a reactive oxygen species (ROS)-attacked breast cancer cell line using the ligation-mediated polymerase chain reaction technique. We detected a significant ‘in vitro’ generation of 8-oxodG between the codons 163 and 175, corresponding to a TP53 region with high mutation prevalence, after treatment with xanthine plus xanthine oxidase, a ROS-generating system. Then, we evaluated the occurrence of oxidative lesions in the DNA-binding domain of the TP53 in the core needle biopsies of 113 of women undergoing breast investigation for diagnostic purpose. An increment of oxidative damage at the −G− residues into the codons 163 and 175 was found in the cancer cases as compared to the controls. We found significant associations with the pathological stage and the histological grade of tumours. As the major news of this study, this largest analysis of genomic footprinting of oxidative lesions at the TP53 sequence level to date provided a first roadmap describing the signatures of oxidative lesions in human breast cancer. Our results provide evidence that the generation of oxidative lesions at single nucleotide resolution is not an event highly stochastic, but causes a characteristic pattern of DNA lesions at the site of mutations in the TP53, suggesting causal relationship between oxidative DNA adducts and breast cancer.

Genomic Footprints in Selected and Unselected Beef Cattle Breeds in Korea.

Korean Hanwoo cattle have been subjected to intensive artificial selection over the past four decades to improve meat production traits. Another three cattle varieties very closely related to Hanwoo reside in Korea (Jeju Black and Brindle) and in China (Yanbian). These breeds have not been part of a breeding scheme to improve production traits. Here, we compare the selected Hanwoo against these similar but presumed to be unselected populations to identify genomic regions that have been under recent selection pressure due to the breeding program. Rsb statistics were used to contrast the genomes of Hanwoo versus a pooled sample of the three unselected population (UN). We identified 37 significant SNPs (FDR corrected) in the HW/UN comparison and 21 known protein coding genes were within 1 MB to the identified SNPs. These genes were previously reported to affect traits important for meat production (14 genes), reproduction including mammary gland development (3 genes), coat color (2 genes), and genes affecting behavioral traits in a broader sense (2 genes). We subsequently sequenced (Illumina HiSeq 2000 platform) 10 individuals of the brown Hanwoo and the Chinese Yanbian to identify SNPs within the candidate genomic regions. Based on allele frequency differences, haplotype structures, and literature research, we singled out one non-synonymous SNP in the APP gene (APP: c.569C>T, Ala199Val) and predicted the mutational effect on the protein structure. We found that protein-protein interactions might be impaired due to increased exposed hydrophobic surfaces of the mutated protein. The APP gene has also been reported to affect meat tenderness in pigs and obesity in humans. Meat tenderness has been linked to intramuscular fat content, which is one of the main breeding goals for brown Hanwoo, potentially supporting a causal influence of the herein described nsSNP in the APP gene.

Genomic and functional approaches reveal a case of adaptive introgression from Populus balsamifera (balsam poplar) in P.trichocarpa (black cottonwood).

Natural hybrid zones in forest trees provide systems to study the transfer of adaptive genetic variation by introgression. Previous landscape genomic studies in Populus trichocarpa, a keystone tree species, indicated genomic footprints of admixture with its sister species Populus balsamifera and identified candidate genes for local adaptation. Here, we explored the patterns of introgression and signals of local adaptation in P.trichocarpa and P.balsamifera, employing genome resequencing data from three chromosomes in pure species and admixed individuals from wild populations. Local ancestry analysis in admixed P.trichocarpa revealed a telomeric region in chromosome 15 with P.balsamifera ancestry, containing several candidate genes for local adaptation. Genomic analyses revealed signals of selection in certain genes in this region (e.g. PRR5, COMT1), and functional analyses based on gene expression variation and correlations with adaptive phenotypes suggest distinct functions of the introgressed alleles. In contrast, a block of genes in chromosome 12 paralogous to the introgressed region showed no signs of introgression or signatures of selection. We hypothesize that the introgressed region in chromosome 15 has introduced modular or cassette-like variation into P.trichocarpa. These linked adaptive mutations are associated with a block of genes in chromosome 15 that appear to have undergone neo- or subfunctionalization relative to paralogs in a duplicated region on chromosome 12 that show no signatures of adaptive variation. The association between P.balsamifera introgressed alleles with the expression of adaptive traits in P.trichocarpa supports the hypothesis that this is a case of adaptive introgression in an ecologically important foundation species.

Genomic footprinting.

The advent of DNA footprinting with DNase I more than 35 years ago enabled the systematic analysis of protein-DNA interactions, and the technique has been instrumental in the decoding of cis-regulatory elements and the identification and characterization of transcription factors and other DNA-binding proteins. The ability to analyze millions of individual genomic cleavage events via massively parallel sequencing has enabled in vivo DNase I footprinting on a genomic scale, offering the potential for global analysis of transcription factor occupancy in a single experiment. Genomic footprinting has opened unique vistas on the organization, function and evolution of regulatory DNA; however, the technology is still nascent. Here we discuss both prospects and challenges of genomic footprinting, as well as considerations for its application to complex genomes.

Genome-wide footprinting: ready for prime time?

High-throughput sequencing technologies have allowed many gene locus-level molecular biology assays to become genome-wide profiling methods. DNA-cleaving enzymes such as DNase I have been used to probe accessible chromatin. The accessible regions contain functional regulatory sites, including promoters, insulators and enhancers. Deep sequencing of DNase-seq libraries and computational analysis of the cut profiles have been used to infer protein occupancy in the genome at the nucleotide level, a method introduced as 'digital genomic footprinting'. The approach has been proposed as an attractive alternative to the analysis of transcription factors (TFs) by chromatin immunoprecipitation followed by sequencing (ChIP-seq), and in theory it should overcome antibody issues, poor resolution and batch effects. Recent reports point to limitations of the DNase-based genomic footprinting approach and call into question the scope of detectable protein occupancy, especially for TFs with short-lived chromatin binding. The genomics community is grappling with issues concerning the utility of genomic footprinting and is reassessing the proposed approaches in terms of robust deliverables. Here we summarize the consensus as well as different views emerging from recent reports, and we describe the remaining issues and hurdles for genomic footprinting.

Genome-Scale Analysis of Cell-Specific Regulatory Codes Using Nuclear Enzymes.

High-throughput sequencing technologies have made it possible for biologists to generate genome-wide profiles of chromatin features at the nucleotide resolution. Enzymes such as nucleases or transposes have been instrumental as a chromatin-probing agent due to their ability to target accessible chromatin for cleavage or insertion. On the scale of a few hundred base pairs, preferential action of the nuclear enzymes on accessible chromatin allows mapping of cell state-specific accessibility in vivo. Such accessible regions contain functionally important regulatory sites, including promoters and enhancers, which undergo active remodeling for cells adapting in a dynamic environment. DNase-seq and the more recent ATAC-seq are two assays that are gaining popularity. Deep sequencing of DNA libraries from these assays, termed genomic footprinting, has been proposed to enable the comprehensive construction of protein occupancy profiles over the genome at the nucleotide level. Recent studies have discovered limitations of genomic footprinting which reduce the scope of detectable proteins. In addition, the identification of putative factors that bind to the observed footprints remains challenging. Despite these caveats, the methodology still presents significant advantages over alternative techniques such as ChIP-seq or FAIRE-seq. Here we describe computational approaches and tools for analysis of chromatin accessibility and genomic footprinting. Proper experimental design and assay-specific data analysis ensure the detection sensitivity and maximize retrievable information. The enzyme-based chromatin profiling approaches represent a powerful and evolving methodology which facilitates our understanding of how the genome is regulated.

Population genomic footprints of fine-scale differentiation between habitats in Mediterranean blue tits.

Linking population genetic variation to the spatial heterogeneity of the environment is of fundamental interest to evolutionary biology and ecology, in particular when phenotypic differences between populations are observed at biologically small spatial scales. Here, we applied restriction-site associated DNA sequencing (RAD-Seq) to test whether phenotypically differentiated populations of wild blue tits (Cyanistes caeruleus) breeding in a highly heterogeneous environment exhibit genetic structure related to habitat type. Using 12106 SNPs in 197 individuals from deciduous and evergreen oak woodlands, we applied complementary population genomic analyses, which revealed that genetic variation is influenced by both geographical distance and habitat type. A fine-scale genetic differentiation supported by genome- and transcriptome-wide analyses was found within Corsica, between two adjacent habitats where blue tits exhibit marked differences in breeding time while nesting <6km apart. Using redundancy analysis (RDA), we show that genomic variation remains associated with habitat type when controlling for spatial and temporal effects. Finally, our results suggest that the observed patterns of genomic differentiation were not driven by a small proportion of highly differentiated loci, but rather emerged through a process such as habitat choice, which reduces gene flow between habitats across the entire genome. The pattern of genomic isolation-by-environment closely matches differentiation observed at the phenotypic level, thereby offering significant potential for future inference of phenotype-genotype associations in a heterogeneous environment.

The impact of selection, gene flow and demographic history on heterogeneous genomic divergence: three-spine sticklebacks in divergent environments.

Heterogeneous genomic divergence between populations may reflect selection, but should also be seen in conjunction with gene flow and drift, particularly population bottlenecks. Marine and freshwater three-spine stickleback (Gasterosteus aculeatus) populations often exhibit different lateral armour plate morphs. Moreover, strikingly parallel genomic footprints across different marine-freshwater population pairs are interpreted as parallel evolution and gene reuse. Nevertheless, in some geographic regions like the North Sea and Baltic Sea, different patterns are observed. Freshwater populations in coastal regions are often dominated by marine morphs, suggesting that gene flow overwhelms selection, and genomic parallelism may also be less pronounced. We used RAD sequencing for analysing 28 888 SNPs in two marine and seven freshwater populations in Denmark, Europe. Freshwater populations represented a variety of environments: river populations accessible to gene flow from marine sticklebacks and large and small isolated lakes with and without fish predators. Sticklebacks in an accessible river environment showed minimal morphological and genomewide divergence from marine populations, supporting the hypothesis of gene flow overriding selection. Allele frequency spectra suggested bottlenecks in all freshwater populations, and particularly two small lake populations. However, genomic footprints ascribed to selection could nevertheless be identified. No genomic regions were consistent freshwater-marine outliers, and parallelism was much lower than in other comparable studies. Two genomic regions previously described to be under divergent selection in freshwater and marine populations were outliers between different freshwater populations. We ascribe these patterns to stronger environmental heterogeneity among freshwater populations in our study as compared to most other studies, although the demographic history involving bottlenecks should also be considered in the interpretation of results.

Diverse Gene Expression in Human Regulatory T Cell Subsets Uncovers Connection between Regulatory T Cell Genes and Suppressive Function.

Regulatory T (Treg) cells have a critical role in the control of immunity, and their diverse subpopulations may allow adaptation to different types of immune responses. In this study, we analyzed human Treg cell subpopulations in the peripheral blood by performing genome-wide expression profiling of 40 Treg cell subsets from healthy donors. We found that the human peripheral blood Treg cell population is comprised of five major genomic subgroups, represented by 16 tractable subsets with a particular cell surface phenotype. These subsets possess a range of suppressive function and cytokine secretion and can exert a genomic footprint on target effector T (Teff) cells. Correlation analysis of variability in gene expression in the subsets identified several cell surface molecules associated with Treg suppressive function, and pharmacological interrogation revealed a set of genes having causative effect. The five genomic subgroups of Treg cells imposed a preserved pattern of gene expression on Teff cells, with a varying degree of genes being suppressed or induced. Notably, there was a cluster of genes induced by Treg cells that bolstered an autoinhibitory effect in Teff cells, and this induction appears to be governed by a different set of genes than ones involved in counteracting Teff activation. Our work shows an example of exploiting the diversity within human Treg cell subpopulations to dissect Treg cell biology.

DNase I hypersensitivity mapping, genomic footprinting, and transcription factor networks in plants

Understanding gene regulatory networks in plants requires knowledge of cis-regulatory DNA, trans-acting factors, and their dynamics across development and in response to stimuli. Active cis-regulatory elements are hypersensitive to cleavage by the endonuclease DNase I. Motifs within DNase I hypersensitive sites indicate potential trans-acting factor occupancy and, when combined with DNase I cleavage data, can be used to construct provisional regulatory networks. Several recent studies have applied genome-wide DNase I hypersensitivity mapping to Arabidopsis thaliana and rice, generating chromatin accessibility landscapes for an array of tissues, cell types, and treatments. Here, we discuss these studies, with an emphasis on building regulatory networks and possible directions for improvements.

Abstract 4744: Solving genomic assay trade-offs with an optimized, extended cancer gene panel for research and clinical applications

Abstract Choosing a genomic assay for cancer research is complicated by trade-offs. Cancer gene panels are a common choice, but they target mutational hotspots in a relatively small number of genes, often for cancers that are most commonly tested and that have common genetic etiologies. A common alternative is exome sequencing, which includes all the coding genes but, due to its larger genomic footprint, cannot reach the same depths as panels and therefore is less able to deal with low tumor purity and heterogeneity. Whole genome sequencing trades off very shallow depth and coverage over vast regions of uninterpretable genomic sequence in exchange for the identification of intergenic variants and structural variant breakpoints. All of these assays can be supplemented with RNA sequencing in order to capture gene fusions, allelic expression, splice isoforms, and gene expression. RNAseq comes with its own costs: the need to extract RNA from the same tissue, the need to perform a second assay, and the need to analyze a very different type of data from DNA sequencing. The trade-offs generally come down to three major issues: depth of sequencing, specific genes targeted, and cost. To solve these, we designed an extended, optimized cancer gene panel facilitating high depth sequencing at low cost. We started by identifying a comprehensive list of over 1,300 cancer genes. These genes were chosen through exhaustive cancer gene database and literature curation, and include genes from all major cancer pathways and from the Cancer Gene Census. We then took this gene list and applied an augmented targeting design strategy that we have previously used to create an augmented exome enrichment platform which fills in gaps that standard technical exomes miss. To validate the panel and analysis, we identified test samples including well-described cancer cell lines, cell line mixtures with engineered cancer variations, and formalin-fixed neoplastic tissues. We then performed a series of tests with these samples to measure the panel's small variant sensitivity and specificity, gauge its limits of detection, validate the detection of gene fusions, and demonstrate its ability to identify copy number alterations and loss of heterozygosity. In engineered cell lines, we detected 100% of small variants down to 5% allele frequency. We also mixed the cancer cell lines in various ratios and found similarly high sensitivity as well as very high specificity for small variant detection. We further compared our structural variation calls in the DNA and our fusion calls in the RNA with known data and found that we had very high concordance with known variations. These studies demonstrate that an extended, augmented cancer gene panel strategy solves many genomic assay trade-offs and leads to high accuracy and variant yield for cancer research applications. Citation Format: Michael J. Clark, Sean M. Boyle, Elena Helman, Shujun Luo, Gabor Bartha, Massimo Morra, Anil Patwardhan, Christian Haudenschild, Mirian Karbelashvili, Parin Sripakdeevong, Jason Harris, Deanna Church, Stephen Chervitz, John West, Richard Chen. Solving genomic assay trade-offs with an optimized, extended cancer gene panel for research and clinical applications. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4744. doi:10.1158/1538-7445.AM2015-4744