Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Clonal dynamics of haematopoiesis across the human lifespan.

Age-related change in human haematopoiesis causes reduced regenerative capacity1, cytopenias2, immune dysfunction3 and increased risk of blood cancer4-6, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at constant rates throughout life, affecting more genes than identified in blood cancers. Loss of the Y chromosome conferred selective benefits in males. Simulations of haematopoiesis, with constant stem cell population size and constant acquisition of driver mutations conferring moderate fitness benefits, entirely explained the abrupt change in clonal structure in the elderly. Rapidly decreasing clonal diversity is a universal feature of haematopoiesis in aged humans, underpinned by pervasive positive selection acting on many more genes than currently identified.

Somatic mutation landscapes at single-molecule resolution.

Somatic mutations drive the development of cancer and may contribute to ageing and other diseases1,2. Despite their importance, the difficulty of detecting mutations that are only present in single cells or small clones has limited our knowledge of somatic mutagenesis to a minority of tissues. Here, to overcome these limitations, we developed nanorate sequencing (NanoSeq), a duplex sequencing protocol with error rates of less than five errors per billion base pairs in single DNA molecules from cell populations. This rate is two orders of magnitude lower than typical somatic mutation loads, enabling the study of somatic mutations in any tissue independently of clonality. We used this single-molecule sensitivity to study somatic mutations in non-dividing cells across several tissues, comparing stem cells to differentiated cells and studying mutagenesis in the absence of cell division. Differentiated cells in blood and colon displayed remarkably similar mutation loads and signatures to their corresponding stem cells, despite mature blood cells having undergone considerably more divisions. We then characterized the mutational landscape of post-mitotic neurons and polyclonal smooth muscle, confirming that neurons accumulate somatic mutations at a constant rate throughout life without cell division, with similar rates to mitotically active tissues. Together, our results suggest that mutational processes that are independent of cell division are important contributors to somatic mutagenesis. We anticipate that the ability to reliably detect mutations in single DNA molecules could transform our understanding of somatic mutagenesis and enable non-invasive studies on large-scale cohorts.

Breaking through to the Other Side: Microenvironment Contributions to DCIS Initiation and Progression.

Refinements in early detection, surgical and radiation therapy, and hormone receptor-targeted treatments have improved the survival rates for breast cancer patients. However, the ability to reliably identify which non-invasive lesions and localized tumors have the ability to progress and/or metastasize remains a major unmet need in the field. The current diagnostic and therapeutic strategies focus on intrinsic alterations within carcinoma cells that are closely associated with proliferation. However, substantial accumulating evidence has indicated that permissive changes in the stromal tissues surrounding the carcinoma play an integral role in breast cancer tumor initiation and progression. Numerous studies have suggested that the stromal environment surrounding ductal carcinoma in situ (DCIS) lesions actively contributes to enhancing tumor cell invasion and immune escape. This review will describe the current state of knowledge regarding the mechanisms through which the microenvironment interacts with DCIS lesions focusing on recent studies that describe the contributions of myoepithelial cells, fibroblasts and immune cells to invasion and subsequent progression. These mechanisms will be considered in the context of developing biomarkers for identifying lesions that will progress to invasive carcinoma and/or developing approaches for therapeutic intervention.

Modeling Human Ductal Carcinoma In Situ in the Mouse.

Breast cancer development is a multi-step process in which genetic and molecular heterogeneity occurs at multiple stages. Ductal carcinoma arises from pre-invasive lesions such as atypical ductal hyperplasia (ADH) and ductal carcinoma in situ (DCIS), which progress to invasive and metastatic cancer. The feasibility of obtaining tissue samples from all stages of progression from the same patient is low, and thus molecular studies dissecting the mechanisms that mediate the transition from pre-invasive DCIS to invasive carcinoma have been hampered. In the past 25 years, numerous mouse models have been developed that partly recapitulate the histological and biological properties of early stage lesions. In this review, we discuss in vivo model systems of breast cancer progression from syngeneic mouse models to human xenografts, with particular focus on how accurately these models mimic human disease.

Gene expression signatures of mating system evolution.

The diversity of mating systems among animals is astounding. Importantly, similar mating systems have evolved even across distantly related taxa. However, our understanding of the mechanisms underlying these convergently evolved phenotypes is limited. Here, we examine on a genomic scale the neuromolecular basis of social organization in cichlids of the tribe Ectodini from Lake Tanganyika. Using field-collected males and females of four closely related species representing two independent evolutionary transitions from polygyny to monogamy, we take a comparative transcriptomic approach to test the hypothesis that these independent transitions have recruited similar gene sets. Our results demonstrate that while lineage and species exert a strong influence on neural gene expression profiles, social phenotype can also drive gene expression evolution. Specifically, 331 genes (∼6% of those assayed) were associated with monogamous mating systems independent of species or sex. Among these genes, we find a strong bias (4:1 ratio) toward genes with increased expression in monogamous individuals. A highly conserved nonapeptide system known to be involved in the regulation of social behavior across animals was not associated with mating system in our analysis. Overall, our findings suggest deep molecular homologies underlying the convergent or parallel evolution of monogamy in different cichlid lineages of Ectodini.

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster.

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.

Cancer Stem Cells and Macrophages: Implications in Tumor Biology and Therapeutic Strategies.

Cancer stem cells (CSCs) are a unique subset of cells within tumors with stemlike properties that have been proposed to be key drivers of tumor initiation and progression. CSCs are functionally defined by their unlimited self-renewal capacity and their ability to initiate tumor formation in vivo. Like normal stem cells, CSCs exist in a cellular niche comprised of numerous cell types including tumor-associated macrophages (TAMs) which provides a unique microenvironment to protect and promote CSC functions. TAMs provide pivotal signals to promote CSC survival, self-renewal, maintenance, and migratory ability, and in turn, CSCs deliver tumor-promoting cues to TAMs that further enhance tumorigenesis. Studies in the last decade have aimed to understand the molecular mediators of CSCs and TAMs, and recent advances have begun to elucidate the complex cross talk that occurs between these two cell types. In this review, we discuss the molecular interactions that define CSC-TAM cross talk at each stage of tumor progression and examine the clinical implications of targeting these interactions.

Drosophila suzukii: the genetic footprint of a recent, worldwide invasion.

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.

Gene duplication in an African cichlid adaptive radiation.

BACKGROUND: Gene duplication is a source of evolutionary innovation and can contribute to the divergence of lineages; however, the relative importance of this process remains to be determined. The explosive divergence of the African cichlid adaptive radiations provides both a model for studying the general role of gene duplication in the divergence of lineages and also an exciting foray into the identification of genomic features that underlie the dramatic phenotypic and ecological diversification in this particular lineage. We present the first genome-wide study of gene duplication in African cichlid fishes, identifying gene duplicates in three species belonging to the Lake Malawi adaptive radiation (Metriaclima estherae, Protomelas similis, Rhamphochromis "chilingali") and one closely related species from a non-radiated riverine lineage (Astatotilapia tweddlei). RESULTS: Using Astatotilapia burtoni as reference, microarray comparative genomic hybridization analysis of 5689 genes reveals 134 duplicated genes among the four cichlid species tested. Between 51 and 55 genes were identified as duplicated in each of the three species from the Lake Malawi radiation, representing a 38%-49% increase in number of duplicated genes relative to the non-radiated lineage (37 genes). Duplicated genes include several that are involved in immune response, ATP metabolism and detoxification. CONCLUSIONS: These results contribute to our understanding of the abundance and type of gene duplicates present in cichlid fish lineages. The duplicated genes identified in this study provide candidates for the analysis of functional relevance with regard to phenotype and divergence. Comparative sequence analysis of gene duplicates can address the role of positive selection and adaptive evolution by gene duplication, while further study across the phylogenetic range of cichlid radiations (and more generally in other adaptive radiations) will determine whether the patterns of gene duplication seen in this study consistently accompany rapid radiation.

Emerging functions of C/EBPß in breast cancer.

Breast tumorigenesis relies on complex interactions between tumor cells and their surrounding microenvironment, orchestrated by tightly regulated transcriptional networks. C/EBPß is a key transcription factor that regulates the proliferation and differentiation of multiple cell types and modulates a variety of biological processes such as tissue homeostasis and the immune response. In addition, C/EBPß has well-established roles in mammary gland development, is overexpressed in breast cancer, and has tumor-promoting functions. In this review, we discuss context-specific roles of C/EBPß during breast tumorigenesis, isoform-specific gene regulation, and regulation of the tumor immune response. We present challenges in C/EBPß biology and discuss the importance of C/EBPß isoform-specific gene regulation in devising new therapeutic strategies.

Convergent somatic evolution commences in utero in a germline ribosomopathy.

Clonal tracking of cells using somatic mutations permits exploration of clonal dynamics in human disease. Here, we perform whole genome sequencing of 323 haematopoietic colonies from 10 individuals with the inherited ribosomopathy Shwachman-Diamond syndrome to reconstruct haematopoietic phylogenies. In ~30% of colonies, we identify mutually exclusive mutations in TP53, EIF6, RPL5, RPL22, PRPF8, plus chromosome 7 and 15 aberrations that increase SBDS and EFL1 gene dosage, respectively. Target gene mutations commence in utero, resulting in a profusion of clonal expansions, with only a few haematopoietic stem cell lineages (mean 8, range 1-24) contributing ~50% of haematopoietic colonies across 8 individuals (range 4-100% clonality) by young adulthood. Rapid clonal expansion during disease transformation is associated with biallelic TP53 mutations and increased mutation burden. Our study highlights how convergent somatic mutation of the p53-dependent nucleolar surveillance pathway offsets the deleterious effects of germline ribosomopathy but increases opportunity for TP53-mutated cancer evolution.

Molecular Targets of Triple-Negative Breast Cancer: Where Do We Stand?

Triple-negative breast cancer (TNBC) is a highly aggressive form of breast cancer. Due to its heterogeneity and lack of hormone receptor expression, this subtype is more likely to metastasize and resist treatment attempts than are other forms of breast cancer. Due to the absence of targetable receptors, chemotherapy and breast conserving surgery have been the predominant treatment options for patients. However, resistance to chemotherapy and local recurrence of the tumors is frequent. Emerging immunotherapies have begun to change treatment plans for patients diagnosed with TNBC. In this review, we discuss the various immune pathways identified in TNBC and the role they play as targets for new potential treatment choices. Various therapeutic options that inhibit key pathways in cellular growth cycles, DNA repair mechanisms, epithelial mesenchymal transition, and immunosuppression have been shown to improve survival in patients with this disease. With promising results thus far, continued studies of immunotherapy and neoadjuvant therapy options for TNBC are likely to alter the treatment course for these diagnoses in the future.

FOXO3 Expression in Macrophages Is Lowered by a High-Fat Diet and Regulates Colonic Inflammation and Tumorigenesis.

Obesity, characterized by augmented inflammation and tumorigenesis, is linked to genetic predispositions, such as FOXO3 polymorphisms. As obesity is associated with aberrant macrophages infiltrating different tissues, including the colon, we aimed to identify FOXO3-dependent transcriptomic changes in macrophages that drive obesity-mediated colonic inflammation and tumorigenesis. We found that in mouse colon, high-fat-diet-(HFD)-related obesity led to diminished FOXO3 levels and increased macrophages. Transcriptomic analysis of mouse peritoneal FOXO3-deficient macrophages showed significant differentially expressed genes (DEGs; FDR < 0.05) similar to HFD obese colons. These DEG-related pathways, linked to mouse colonic inflammation and tumorigenesis, were similar to those in inflammatory bowel disease (IBD) and human colon cancer. Additionally, we identified a specific transcriptional signature for the macrophage-FOXO3 axis (MAC-FOXO382), which separated the transcriptome of affected tissue from control in both IBD (p = 5.2 × 10−8 and colon cancer (p = 1.9 × 10−11), revealing its significance in human colonic pathobiologies. Further, we identified (heatmap) and validated (qPCR) DEGs specific to FOXO3-deficient macrophages with established roles both in IBD and colon cancer (IL-1B, CXCR2, S100A8, S100A9, and TREM1) and those with unexamined roles in these colonic pathobiologies (STRA6, SERPINH1, LAMB1, NFE2L3, OLR1, DNAJC28 and VSIG10). These findings establish an important understanding of how HFD obesity and related metabolites promote colonic pathobiologies.

Breast cancer cells survive chemotherapy by activating targetable immune-modulatory programs characterized by PD-L1 or CD80.

Breast cancer cells must avoid intrinsic and extrinsic cell death to relapse following chemotherapy. Entering senescence enables survival from mitotic catastrophe, apoptosis and nutrient deprivation, but mechanisms of immune evasion are poorly understood. Here we show that breast tumors surviving chemotherapy activate complex programs of immune modulation. Characterization of residual disease revealed distinct tumor cell populations. The first population was characterized by interferon response genes, typified by Cd274, whose expression required chemotherapy to enhance chromatin accessibility, enabling recruitment of IRF1 transcription factor. A second population was characterized by p53 signaling, typified by CD80 expression. Treating mammary tumors with chemotherapy followed by targeting the PD-L1 and/or CD80 axes resulted in marked accumulation of T cells and improved response; however, even combination strategies failed to fully eradicate tumors in the majority of cases. Our findings reveal the challenge of eliminating residual disease populated by senescent cells expressing redundant immune inhibitory pathways and highlight the need for rational immune targeting strategies.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila.

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization.

BACKGROUND: Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates. RESULTS: Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species. CONCLUSIONS: Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila.

BACKGROUND: Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence. RESULTS: We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) approximately 84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (>or= 97% identity), only 13% of genes were identified as diverged. While approximately 40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation. CONCLUSIONS: Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.

Pervasive Strong Selection at the Level of Codon Usage Bias in Drosophila melanogaster.

Codon usage bias (CUB), where certain codons are used more frequently than expected by chance, is a ubiquitous phenomenon and occurs across the tree of life. The dominant paradigm is that the proportion of preferred codons is set by weak selection. While experimental changes in codon usage have at times shown large phenotypic effects in contrast to this paradigm, genome-wide population genetic estimates have supported the weak selection model. Here we use deep genomic population sequencing of two Drosophila melanogaster populations to measure selection on synonymous sites in a way that allowed us to estimate the prevalence of both weak and strong purifying selection. We find that selection in favor of preferred codons ranges from weak (|Nes| ∼ 1) to strong (|Nes| > 10), with strong selection acting on 10-20% of synonymous sites in preferred codons. While previous studies indicated that selection at synonymous sites could be strong, this is the first study to detect and quantify strong selection specifically at the level of CUB. Further, we find that CUB-associated polymorphism accounts for the majority of strong selection on synonymous sites, with secondary contributions of splicing (selection on alternatively spliced genes, splice junctions, and spliceosome-bound sites) and transcription factor binding. Our findings support a new model of CUB and indicate that the functional importance of CUB, as well as synonymous sites in general, have been underestimated.

Gas6 expression is reduced in advanced breast cancers.

Growth arrest-specific gene 6 (Gas6) is a cytokine that binds to receptor tyrosine kinases Tyro3, Axl, and Mer. Numerous studies have suggested that macrophage-derived Gas6 interacts with Axl to promote cancer progression, and Axl has been associated with poor clinical outcome. However, the expression and relevance of Gas6 in human breast cancer patients has not been studied. Analysis of tissue microarrays showed that Gas6 was highly expressed in ductal carcinoma in situ (DCIS) but markedly decreased in invasive breast cancer. Gas6 and Axl were weakly correlated, suggesting that their functions may not exclusively rely on each other. Analyses of publicly available databases showed significantly improved overall and relapse-free survival in patients with high Gas6 mRNA, particularly in luminal A breast cancers. These findings indicate that tumor-derived Gas6 is not overexpressed in invasive breast cancer, and may not be a negative prognostic factor in human breast cancer.