High altitude adaptation in Daghestani populations from the Caucasus

We have surveyed 15 high-altitude adaptation candidate genes for signals of positive selection in North Caucasian highlanders using targeted re-sequencing. A total of 49 unrelated Daghestani from three ethnic groups (Avars, Kubachians, and Laks) living in ancient villages located at around 2,000 m above sea level were chosen as the study population. Caucasian (Adygei living at sea level, N = 20) and CEU (CEPH Utah residents with ancestry from northern and western Europe; N = 20) were used as controls. Candidate genes were compared with 20 putatively neutral control regions resequenced in the same individuals. The regions of interest were amplified by long-PCR, pooled according to individual, indexed by adding an eight-nucleotide tag, and sequenced using the Illumina GAII platform. 1,066 SNPs were called using false discovery and false negative thresholds of ~6%. The neutral regions provided an empirical null distribution to compare with the candidate genes for signals of selection. Two genes stood out. In Laks, a non-synonymous variant within HIF1A already known to be associated with improvement in oxygen metabolism was rediscovered, and in Kubachians a cluster of 13 SNPs located in a conserved intronic region within EGLN1 showing high population differentiation was found. These variants illustrate both the common pathways of adaptation to high altitude in different populations and features specific to the Daghestani populations, showing how even a mildly hypoxic environment can lead to genetic adaptation.     

A modified method for whole exome resequencing from minimal amounts of starting DNA

Next generation DNA sequencing (NGS) technologies have revolutionized the pace at which whole genome and exome sequences can be generated. However, despite these advances, many of the methods for targeted resequencing, such as the generation of high-depth exome sequences, are somewhat limited by the relatively large amounts of starting DNA that are normally required. In the case of tumour analysis this is particularly pertinent as many tumour biopsies often return submicrogram quantities of DNA, especially when tumours are microdissected prior to analysis. Here, we present a method for exome capture and resequencing using as little as 50 ng of starting DNA. The sequencing libraries generated by this minimal starting amount (MSA-Cap) method generate datasets that are comparable to standard amount (SA) whole exome libraries that use three micrograms of starting DNA. This method, which can be performed in most laboratories using commonly available reagents, has the potential to enhance large scale profiling efforts such as the resequencing of tumour exomes.     

A Genetic Screen Using the PiggyBac Transposon in Haploid Cells Identifies Parp1 as a Mediator of Olaparib Toxicity

Genetic perturbation screens have the potential to dissect a wide range of cellular phenotypes. Such screens have historically been difficult in diploid mammalian cells. The recent derivation of haploid embryonic stem cells provides an opportunity to cause loss of function mutants with a random mutagen in a mammalian cell with a normal genetic background. We describe an approach to genetic screens that exploits the highly active piggyBac transposon in haploid mammalian cells. As an example of haploid transposon (HTP) screening, we apply this approach to identifying determinants of cancer drug toxicity and resistance. In a screen for 6-thioguanine resistance we recovered components of the DNA mismatch repair pathway, a known requirement for toxicity. In a further screen for resistance to the clinical poly(ADP-ribose) polymerase (PARP) inhibitor olaparib we recovered multiple Parp1 mutants. Our results show that olaparib toxicity to normal cells is mediated predominantly via Parp1, and suggest that the clinical side effects of olaparib may be on target. The transposon mutant libraries are stable and can be readily reused to screen other drugs. The screening protocol described has several advantages over other methods such as RNA interference: it is rapid and low cost, and mutations can be easily reverted to establish causality.     

Whole genome sequencing highlights genetic changes associated with laboratory domestication of C. elegans

Defining the mutational landscape when individuals of a species grow separately and diverge over many generations can provide insights into trait evolution. A specific example of this involves studying changes associated with domestication where different lines of the same wild stock have been cultivated independently in different standard environments. Whole genome sequence comparison of such lines permits estimation of mutation rates, inference of genes' ancestral states and ancestry of existing strains, and correction of sequencing errors in genome databases. Here we study domestication of the C. elegans Bristol strain as a model, and report the genome sequence of LSJ1 (Bristol), a sibling of the standard C. elegans reference wild type N2 (Bristol). The LSJ1 and N2 lines were cultivated separately from shortly after the Bristol strain was isolated until methods to freeze C. elegans were developed. We find that during this time the two strains have accumulated 1208 genetic differences. We describe phenotypic variation between N2 and LSJ1 in the rate at which embryos develop, the rate of production of eggs, the maturity of eggs at laying, and feeding behavior, all the result of post-isolation changes. We infer the ancestral alleles in the original Bristol isolate and highlight 2038 likely sequencing errors in the original N2 reference genome sequence. Many of these changes modify genome annotation. Our study provides a starting point to further investigate genotype-phenotype association and offers insights into the process of selection as a result of laboratory domestication.     

Alteration of PHYA expression change circadian rhythms and timing of bud set in Populus

In many temperate woody species, dormancy is induced by short photoperiods. Earlier studies have shown that the photoreceptor phytochrome A (phyA) promotes growth. Specifically, Populus plants that over-express the oat PHYA gene (oatPHYAox) show daylength-independent growth and do not become dormant. However, we show that oatPHYAox plants could be induced to set bud and become cold hardy by exposure to a shorter, non-24 h diurnal cycle that significantly alters the relative position between endogenous rhythms and perceived light/dark cycles. Furthermore, we describe studies in which the expression of endogenous Populus tremula × P. tremuloides PHYTOCHROME A (PttPHYA) was reduced in Populus trees by antisense inhibition. The antisense plants showed altered photoperiodic requirements, resulting in earlier growth cessation and bud formation in response to daylength shortening, an effect that was explained by an altered innate period that leads to phase changes of clock-associated genes such as PttCO2. Moreover, gene expression studies following far-red light pulses show a phyA-mediated repression of PttLHY1 and an induction of PttFKF1 and PttFT. We conclude that the level of PttPHYA expression strongly influences seasonally regulated growth in Populus and is central to co-ordination between internal clock-regulated rhythms and external light/dark cycles through its dual effect on the pace of clock rhythms and in light signaling.     

Comprehensive genomic analysis of a BRCA2 deficient human pancreatic cancer

Capan-1 is a well-characterised BRCA2-deficient human cell line isolated from a liver metastasis of a pancreatic adenocarcinoma. Here we report a genome-wide assessment of structural variations and high-depth exome characterization of single nucleotide variants and small insertion/deletions in Capan-1. To identify potential somatic and tumour-associated variations in the absence of a matched-normal cell line, we devised a novel method based on the analysis of HapMap samples. We demonstrate that Capan-1 has one of the most rearranged genomes sequenced to date. Furthermore, small insertions and deletions are detected more frequently in the context of short sequence repeats than in other genomes. We also identify a number of novel mutations that may represent genetic changes that have contributed to tumour progression. These data provide insight into the genomic effects of loss of BRCA2 function.     

Genomic characterisation of acral melanoma cell lines

Acral melanoma is a rare melanoma subtype with distinct epidemiological, clinical and genetic features. To determine if acral melanoma cell lines are representative of this melanoma subtype, six lines were analysed by whole-exome sequencing and array comparative genomic hybridisation. We demonstrate that the cell lines display a mutation rate that is comparable to that of published primary and metastatic acral melanomas and observe a mutational signature suggestive of UV-induced mutagenesis in two of the cell lines. Mutations were identified in oncogenes and tumour suppressors previously linked to melanoma including BRAF, NRAS, KIT, PTEN and TP53, in cancer genes not previously linked to melanoma and in genes linked to DNA repair such as BRCA1 and BRCA2. Our findings provide strong circumstantial evidence to suggest that acral melanoma cell lines and acral tumours share genetic features in common and that these cells are therefore valuable tools to investigate the biology of this aggressive melanoma subtype. Data are available at: http://rock.icr.ac.uk/collaborations/Furney_et_al_2012/.     

Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach

Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogene-driven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCR™ ddPCR™ Mutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma.     

High-throughput RNA interference screening using pooled shRNA libraries and next generation sequencing

RNA interference (RNAi) screening is a state-of-the-art technology that enables the dissection of biological processes and disease-related phenotypes. The commercial availability of genome-wide, short hairpin RNA (shRNA) libraries has fueled interest in this area but the generation and analysis of these complex data remain a challenge. Here, we describe complete experimental protocols and novel open source computational methodologies, shALIGN and shRNAseq, that allow RNAi screens to be rapidly deconvoluted using next generation sequencing. Our computational pipeline offers efficient screen analysis and the flexibility and scalability to quickly incorporate future developments in shRNA library technology.