Altered splicing machinery in lung carcinoids unveils NOVA1, PRPF8 and SRSF10 as novel candidates to understand tumor biology and expand biomarker discovery.

BACKGROUND: Lung neuroendocrine neoplasms (LungNENs) comprise a heterogeneous group of tumors ranging from indolent lesions with good prognosis to highly aggressive cancers. Carcinoids are the rarest LungNENs, display low to intermediate malignancy and may be surgically managed, but show resistance to radiotherapy/chemotherapy in case of metastasis. Molecular profiling is providing new information to understand lung carcinoids, but its clinical value is still limited. Altered alternative splicing is emerging as a novel cancer hallmark unveiling a highly informative layer. METHODS: We primarily examined the status of the splicing machinery in lung carcinoids, by assessing the expression profile of the core spliceosome components and selected splicing factors in a cohort of 25 carcinoids using a microfluidic array. Results were validated in an external set of 51 samples. Dysregulation of splicing variants was further explored in silico in a separate set of 18 atypical carcinoids. Selected altered factors were tested by immunohistochemistry, their associations with clinical features were assessed and their putative functional roles were evaluated in vitro in two lung carcinoid-derived cell lines. RESULTS: The expression profile of the splicing machinery was profoundly dysregulated. Clustering and classification analyses highlighted five splicing factors: NOVA1, SRSF1, SRSF10, SRSF9 and PRPF8. Anatomopathological analysis showed protein differences in the presence of NOVA1, PRPF8 and SRSF10 in tumor versus non-tumor tissue. Expression levels of each of these factors were differentially related to distinct number and profiles of splicing events, and were associated to both common and disparate functional pathways. Accordingly, modulating the expression of NOVA1, PRPF8 and SRSF10 in vitro predictably influenced cell proliferation and colony formation, supporting their functional relevance and potential as actionable targets. CONCLUSIONS: These results provide primary evidence for dysregulation of the splicing machinery in lung carcinoids and suggest a plausible functional role and therapeutic targetability of NOVA1, PRPF8 and SRSF10.

Differential Orthopedia Homeobox expression in pulmonary carcinoids is associated with changes in DNA methylation.

Limited number of tumor types have been examined for Orthopedia Homeobox (OTP) expression. In pulmonary carcinoids, loss of expression is a strong indicator of poor prognosis. Here, we investigated OTP expression in 37 different tumor types, and the association between OTP expression and DNA methylation levels in lung neuroendocrine neoplasms. We analyzed publicly available multi-omics data (whole-exome-, whole-genome-, RNA sequencing and Epic 850K-methylation array) of 58 typical carcinoids, 27 atypical carcinoids, 69 large cell neuroendocrine carcinoma and 51 small cell lung cancer patients and TCGA (The Cancer Genome Atlas) data of 33 tumor types. 850K-methylation analysis was cross-validated using targeted pyrosequencing on 35 carcinoids. We report bimodality of OTP expression in carcinoids (OTPhigh vs OTPlow group, likelihood-ratio test P = 1.5 × 10-2 ), with the OTPhigh group specific to pulmonary carcinoids while absent from all other cohorts analyzed. Significantly different DNA methylation levels were observed between OTPhigh and OTPlow carcinoids in 12/34 OTP infinium probes (FDR < 0.05 and ß-value effect size > .2). OTPlow carcinoids harbor high DNA methylation levels as compared to OTPhigh carcinoids. OTPlow carcinoids showed a significantly worse overall survival (log-rank test P = .0052). Gene set enrichment analysis for somatically mutated genes associated with hallmarks of cancer showed robust enrichment of three hallmarks in the OTPlow group, that is, sustaining proliferative signaling, evading growth suppressor and genome instability and mutation. Together our data suggest that high OTP expression is a unique feature of pulmonary carcinoids with a favorable prognosis and that in poor prognostic patients, OTP expression is lost, most likely due to changes in DNA methylation levels.

TP53 Targeted Deep Sequencing of Cell-Free DNA in Esophageal Squamous Cell Carcinoma Using Low-Quality Serum: Concordance with Tumor Mutation.

Circulating cell-free DNA (cfDNA) is emerging as a potential tumor biomarker. CfDNA-based biomarkers may be applicable in tumors without an available non-invasive screening method among at-risk populations. Esophageal squamous cell carcinoma (ESCC) and residents of the Asian cancer belt are examples of those malignancies and populations. Previous epidemiological studies using cfDNA have pointed to the need for high volumes of good quality plasma (i.e., >1 mL plasma with 0 or 1 cycles of freeze-thaw) rather than archival serum, which is often the main available source of cfDNA in retrospective studies. Here, we have investigated the concordance of TP53 mutations in tumor tissue and cfDNA extracted from archival serum left-over from 42 cases and 39 matched controls (age, gender, residence) in a high-risk area of Northern Iran (Golestan). Deep sequencing of TP53 coding regions was complemented with a specialized variant caller (Needlestack). Overall, 23% to 31% of mutations were concordantly detected in tumor and serum cfDNA (based on two false discovery rate thresholds). Concordance was positively correlated with high cfDNA concentration, smoking history (p-value = 0.02) and mutations with a high potential of neoantigen formation (OR; 95%CI = 1.9 (1.11-3.29)), suggesting that tumor DNA release in the bloodstream might reflect the effects of immune and inflammatory context on tumor cell turnover. We identified TP53 mutations in five controls, one of whom was subsequently diagnosed with ESCC. Overall, the results showed that cfDNA mutations can be reliably identified by deep sequencing of archival serum, with a rate of success comparable to plasma. Nonetheless, 70% non-identifiable mutations among cancer patients and 12% mutation detection in controls are the main challenges in applying cfDNA to detect tumor-related variants when blindly targeting whole coding regions of the TP53 gene in ESCC.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow.

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

The influence of obesity-related factors in the etiology of renal cell carcinoma-A mendelian randomization study.

BACKGROUND: Several obesity-related factors have been associated with renal cell carcinoma (RCC), but it is unclear which individual factors directly influence risk. We addressed this question using genetic markers as proxies for putative risk factors and evaluated their relation to RCC risk in a mendelian randomization (MR) framework. This methodology limits bias due to confounding and is not affected by reverse causation. METHODS AND FINDINGS: Genetic markers associated with obesity measures, blood pressure, lipids, type 2 diabetes, insulin, and glucose were initially identified as instrumental variables, and their association with RCC risk was subsequently evaluated in a genome-wide association study (GWAS) of 10,784 RCC patients and 20,406 control participants in a 2-sample MR framework. The effect on RCC risk was estimated by calculating odds ratios (ORSD) for a standard deviation (SD) increment in each risk factor. The MR analysis indicated that higher body mass index increases the risk of RCC (ORSD: 1.56, 95% confidence interval [CI] 1.44-1.70), with comparable results for waist-to-hip ratio (ORSD: 1.63, 95% CI 1.40-1.90) and body fat percentage (ORSD: 1.66, 95% CI 1.44-1.90). This analysis further indicated that higher fasting insulin (ORSD: 1.82, 95% CI 1.30-2.55) and diastolic blood pressure (DBP; ORSD: 1.28, 95% CI 1.11-1.47), but not systolic blood pressure (ORSD: 0.98, 95% CI 0.84-1.14), increase the risk for RCC. No association with RCC risk was seen for lipids, overall type 2 diabetes, or fasting glucose. CONCLUSIONS: This study provides novel evidence for an etiological role of insulin in RCC, as well as confirmatory evidence that obesity and DBP influence RCC risk.

Thinking too positive? Revisiting current methods of population genetic selection inference.

In the age of next-generation sequencing, the availability of increasing amounts and improved quality of data at decreasing cost ought to allow for a better understanding of how natural selection is shaping the genome than ever before. However, alternative forces, such as demography and background selection (BGS), obscure the footprints of positive selection that we would like to identify. In this review, we illustrate recent developments in this area, and outline a roadmap for improved selection inference. We argue (i) that the development and obligatory use of advanced simulation tools is necessary for improved identification of selected loci, (ii) that genomic information from multiple time points will enhance the power of inference, and (iii) that results from experimental evolution should be utilized to better inform population genomic studies.

Widespread signals of convergent adaptation to high altitude in Asia and america.

Living at high altitude is one of the most difficult challenges that humans had to cope with during their evolution. Whereas several genomic studies have revealed some of the genetic bases of adaptations in Tibetan, Andean, and Ethiopian populations, relatively little evidence of convergent evolution to altitude in different continents has accumulated. This lack of evidence can be due to truly different evolutionary responses, but it can also be due to the low power of former studies that have mainly focused on populations from a single geographical region or performed separate analyses on multiple pairs of populations to avoid problems linked to shared histories between some populations. We introduce here a hierarchical Bayesian method to detect local adaptation that can deal with complex demographic histories. Our method can identify selection occurring at different scales, as well as convergent adaptation in different regions. We apply our approach to the analysis of a large SNP data set from low- and high-altitude human populations from America and Asia. The simultaneous analysis of these two geographic areas allows us to identify several candidate genome regions for altitudinal selection, and we show that convergent evolution among continents has been quite common. In addition to identifying several genes and biological processes involved in high-altitude adaptation, we identify two specific biological pathways that could have evolved in both continents to counter toxic effects induced by hypoxia.

Adaptive, convergent origins of the pygmy phenotype in African rainforest hunter-gatherers.

The evolutionary history of the human pygmy phenotype (small body size), a characteristic of African and Southeast Asian rainforest hunter-gatherers, is largely unknown. Here we use a genome-wide admixture mapping analysis to identify 16 genomic regions that are significantly associated with the pygmy phenotype in the Batwa, a rainforest hunter-gatherer population from Uganda (east central Africa). The identified genomic regions have multiple attributes that provide supporting evidence of genuine association with the pygmy phenotype, including enrichments for SNPs previously associated with stature variation in Europeans and for genes with growth hormone receptor and regulation functions. To test adaptive evolutionary hypotheses, we computed the haplotype-based integrated haplotype score (iHS) statistic and the level of population differentiation (FST) between the Batwa and their agricultural neighbors, the Bakiga, for each genomic SNP. Both |iHS| and FST values were significantly higher for SNPs within the Batwa pygmy phenotype-associated regions than the remainder of the genome, a signature of polygenic adaptation. In contrast, when we expanded our analysis to include Baka rainforest hunter-gatherers from Cameroon and Gabon (west central Africa) and Nzebi and Nzime neighboring agriculturalists, we did not observe elevated |iHS| or FST values in these genomic regions. Together, these results suggest adaptive and at least partially convergent origins of the pygmy phenotype even within Africa, supporting the hypothesis that small body size confers a selective advantage for tropical rainforest hunter-gatherers but raising questions about the antiquity of this behavior.

Influenza virus drug resistance: a time-sampled population genetics perspective.

The challenge of distinguishing genetic drift from selection remains a central focus of population genetics. Time-sampled data may provide a powerful tool for distinguishing these processes, and we here propose approximate Bayesian, maximum likelihood, and analytical methods for the inference of demography and selection from time course data. Utilizing these novel statistical and computational tools, we evaluate whole-genome datasets of an influenza A H1N1 strain in the presence and absence of oseltamivir (an inhibitor of neuraminidase) collected at thirteen time points. Results reveal a striking consistency amongst the three estimation procedures developed, showing strongly increased selection pressure in the presence of drug treatment. Importantly, these approaches re-identify the known oseltamivir resistance site, successfully validating the approaches used. Enticingly, a number of previously unknown variants have also been identified as being positively selected. Results are interpreted in the light of Fisher's Geometric Model, allowing for a quantification of the increased distance to optimum exerted by the presence of drug, and theoretical predictions regarding the distribution of beneficial fitness effects of contending mutations are empirically tested. Further, given the fit to expectations of the Geometric Model, results suggest the ability to predict certain aspects of viral evolution in response to changing host environments and novel selective pressures.

Positive Selection Drives Preferred Segment Combinations during Influenza Virus Reassortment.

Influenza A virus (IAV) has a segmented genome that allows for the exchange of genome segments between different strains. This reassortment accelerates evolution by breaking linkage, helping IAV cross species barriers to potentially create highly virulent strains. Challenges associated with monitoring the process of reassortment in molecular detail have limited our understanding of its evolutionary implications. We applied a novel deep sequencing approach with quantitative analysis to assess the in vitro temporal evolution of genomic reassortment in IAV. The combination of H1N1 and H3N2 strains reproducibly generated a new H1N2 strain with the hemagglutinin and nucleoprotein segments originating from H1N1 and the remaining six segments from H3N2. By deep sequencing the entire viral genome, we monitored the evolution of reassortment, quantifying the relative abundance of all IAV genome segments from the two parent strains over time and measuring the selection coefficients of the reassorting segments. Additionally, we observed several mutations coemerging with reassortment that were not found during passaging of pure parental IAV strains. Our results demonstrate how reassortment of the segmented genome can accelerate viral evolution in IAV, potentially enabled by the emergence of a small number of individual mutations.

Inferring the age of a fixed beneficial allele.

Estimating the age and strength of beneficial alleles is central to understanding how adaptation proceeds in response to changing environmental conditions. Several haplotype-based estimators exist for inferring the age of segregating beneficial mutations. Here, we develop an approximate Bayesian-based approach that rather estimates these parameters for fixed beneficial mutations in single populations. We integrate a range of existing diversity, site frequency spectrum, haplotype- and linkage disequilibrium-based summary statistics. We show that for strong selective sweeps on de novo mutations the method can estimate allele age and selection strength even in nonequilibrium demographic scenarios. We extend our approach to models of selection on standing variation, and co-infer the frequency at which selection began to act upon the mutation. Finally, we apply our method to estimate the age and selection strength of a previously identified mutation underpinning cryptic colour adaptation in a wild deer mouse population, and compare our findings with previously published estimates as well as with geological data pertaining to the presumed shift in selective pressure.

Robust demographic inference from genomic and SNP data.

We introduce a flexible and robust simulation-based framework to infer demographic parameters from the site frequency spectrum (SFS) computed on large genomic datasets. We show that our composite-likelihood approach allows one to study evolutionary models of arbitrary complexity, which cannot be tackled by other current likelihood-based methods. For simple scenarios, our approach compares favorably in terms of accuracy and speed with ∂a∂i, the current reference in the field, while showing better convergence properties for complex models. We first apply our methodology to non-coding genomic SNP data from four human populations. To infer their demographic history, we compare neutral evolutionary models of increasing complexity, including unsampled populations. We further show the versatility of our framework by extending it to the inference of demographic parameters from SNP chips with known ascertainment, such as that recently released by Affymetrix to study human origins. Whereas previous ways of handling ascertained SNPs were either restricted to a single population or only allowed the inference of divergence time between a pair of populations, our framework can correctly infer parameters of more complex models including the divergence of several populations, bottlenecks and migration. We apply this approach to the reconstruction of African demography using two distinct ascertained human SNP panels studied under two evolutionary models. The two SNP panels lead to globally very similar estimates and confidence intervals, and suggest an ancient divergence (>110 Ky) between Yoruba and San populations. Our methodology appears well suited to the study of complex scenarios from large genomic data sets.

Genomic data reveal a complex making of humans.

In the last few years, two paradigms underlying human evolution have crumbled. Modern humans have not totally replaced previous hominins without any admixture, and the expected signatures of adaptations to new environments are surprisingly lacking at the genomic level. Here we review current evidence about archaic admixture and lack of strong selective sweeps in humans. We underline the need to properly model differential admixture in various populations to correctly reconstruct past demography. We also stress the importance of taking into account the spatial dimension of human evolution, which proceeded by a series of range expansions that could have promoted both the introgression of archaic genes and background selection.

Approximate Bayesian computation.

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

Multi-omic dataset of patient-derived tumor organoids of neuroendocrine neoplasms.

BACKGROUND: Organoids are 3-dimensional experimental models that summarize the anatomical and functional structure of an organ. Although a promising experimental model for precision medicine, patient-derived tumor organoids (PDTOs) have currently been developed only for a fraction of tumor types. RESULTS: We have generated the first multi-omic dataset (whole-genome sequencing [WGS] and RNA-sequencing [RNA-seq]) of PDTOs from the rare and understudied pulmonary neuroendocrine tumors (n = 12; 6 grade 1, 6 grade 2) and provide data from other rare neuroendocrine neoplasms: small intestine (ileal) neuroendocrine tumors (n = 6; 2 grade 1 and 4 grade 2) and large-cell neuroendocrine carcinoma (n = 5; 1 pancreatic and 4 pulmonary). This dataset includes a matched sample from the parental sample (primary tumor or metastasis) for a majority of samples (21/23) and longitudinal sampling of the PDTOs (1 to 2 time points), for a total of n = 47 RNA-seq and n = 33 WGS. We here provide quality control for each technique and the raw and processed data as well as all scripts for genomic analyses to ensure an optimal reuse of the data. In addition, we report gene expression data and somatic small variant calls and describe how they were generated, in particular how we used WGS somatic calls to train a random forest classifier to detect variants in tumor-only RNA-seq. We also report all histopathological images used for medical diagnosis: hematoxylin and eosin-stained slides, brightfield images, and immunohistochemistry images of protein markers of clinical relevance. CONCLUSIONS: This dataset will be critical to future studies relying on this PDTO biobank, such as drug screens for novel therapies and experiments investigating the mechanisms of carcinogenesis in these understudied diseases.

Clinical management of typical and atypical carcinoids/neuroendocrine tumors in ENETS centres of excellence (CoE): Survey from the ENETS lung NET task force.

Lung carcinoid tumours are neuroendocrine neoplasms originating from the bronchopulmonary tract's neuroendocrine cells, accounting for only 1%-3% of all lung cancers but 30% of all neuroendocrine tumours. The incidence of lung carcinoids, both typical and atypical, has been increasing over the years due to improved diagnostic methods and increased awareness among clinicians and pathologists. The most recent WHO classification includes a subgroup of lung carcinoids with atypical morphology and higher mitotic count and/or Ki67 labelling index. Despite appropriate surgery, the 5-year survival rate for atypical carcinoids barely exceeds 50%-70%. The role of adjuvant therapy in lung carcinoids is not well-defined, and clinical decisions are generally based on the presence of high-risk features. Long-term follow-up is essential to monitor for recurrence, although the optimal follow-up protocol remains unclear. To address the lack of consensus in clinical management decisions, the European Neuroendocrine Tumor Society (ENETS) initiated a survey among 20 expert centres. The survey identified varied opinions on approaches to imaging, surgery, use of adjuvant therapy, and follow-up protocols. Notably, the absence of dedicated multidisciplinary lung neuroendocrine tumour boards in some centres was evident. Experts agreed on the need for a prospective adjuvant trial in high-risk patients, emphasizing the feasibility of such a study. In conclusion, the study highlights the need for a more uniform adoption of existing guidelines in the management of lung carcinoid tumours and emphasizes the importance of international collaboration to advance research and patient care. Close collaboration between healthcare providers and patients is vital for effective long-term surveillance and management of these rare tumours.

Spotlight on Small-Cell Lung Cancer and Other Lung Neuroendocrine Neoplasms.

Lung neuroendocrine neoplasms (NENs) encompass a spectrum of neoplasms that are subdivided into the well-differentiated neuroendocrine tumors comprising the low- and intermediate-grade typical and atypical carcinoids, respectively, and the poorly differentiated, high-grade neuroendocrine carcinomas including large-cell neuroendocrine carcinomas and small-cell lung carcinoma (SCLC). Here, we review the current morphological and molecular classifications of the NENs on the basis of the updated WHO Classification of Thoracic Tumors and discuss the emerging subclassifications on the basis of molecular profiling and the potential therapeutic implications. We focus on the efforts in subtyping SCLC, a particularly aggressive tumor with few treatment options, and the recent advances in therapy with the adoption of immune checkpoint inhibitors in the frontline setting for patients with extensive-stage SCLC. We further highlight the promising immunotherapy strategies in SCLC that are currently under investigation.

Druggable growth dependencies and tumor evolution analysis in patient-derived organoids of neuroendocrine neoplasms from multiple body sites.

Neuroendocrine neoplasms (NENs) comprise well-differentiated neuroendocrine tumors (NETs) and poorly differentiated neuroendocrine carcinomas (NECs). Treatment options for patients with NENs are limited, in part due to lack of accurate models. We establish patient-derived tumor organoids (PDTOs) from pulmonary NETs and derive PDTOs from an understudied subtype of NEC, large cell neuroendocrine carcinoma (LCNEC), arising from multiple body sites. PDTOs maintain the gene expression patterns, intra-tumoral heterogeneity, and evolutionary processes of parental tumors. Through hypothesis-driven drug sensitivity analyses, we identify ASCL1 as a potential biomarker for response of LCNEC to treatment with BCL-2 inhibitors. Additionally, we discover a dependency on EGF in pulmonary NET PDTOs. Consistent with these findings, we find that, in an independent cohort, approximately 50% of pulmonary NETs express EGFR. This study identifies an actionable vulnerability for a subset of pulmonary NETs, emphasizing the utility of these PDTO models.

Evolutionary forces shaping genomic islands of population differentiation in humans.

BACKGROUND: Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. METHODS: Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. RESULTS: Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. CONCLUSIONS: Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

fastsimcoal: a continuous-time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios.

MOTIVATION: Genetic studies focus on increasingly larger genomic regions of both extant and ancient DNA, and there is a need for simulation software to match these technological advances. We present here a new coalescent-based simulation program fastsimcoal, which is able to quickly simulate a variety of genetic markers scattered over very long genomic regions with arbitrary recombination patterns under complex evolutionary scenarios. AVAILABILITY AND IMPLEMENTATION: fastsimcoal is a C++ program compiled for Windows, MacOsX and Linux platforms. It is freely available at cmpg.unibe.ch/software/fastsimcoal/, together with its detailed user manual and example input files.