Question: 1,000 genomes project enables mapping of genetic sequence variation consisting of rare variants with larger penetrance effects than common variants.

Title: Rare Variants Create Synthetic Genome-Wide Associations Content: Genome-wide association studies (GWAS) have now identified at least 2,000 common variants that appear associated with common diseases or related traits (http://www.genome.gov/gwastudies), hundreds of which have been convincingly replicated. It is generally thought that the associated markers reflect the effect of a nearby common (minor allele frequency >0.05) causal site, which is associated with the marker, leading to extensive resequencing efforts to find causal sites. We propose as an alternative explanation that variants much less common than the associated one may create "synthetic associations" by occurring, stochastically, more often in association with one of the alleles at the common site versus the other allele. Although synthetic associations are an obvious theoretical possibility, they have never been systematically explored as a possible explanation for GWAS findings. Here, we use simple computer simulations to show the conditions under which such synthetic associations will arise and how they may be recognized. We show that they are not only possible, but inevitable, and that under simple but reasonable genetic models, they are likely to account for or contribute to many of the recently identified signals reported in genome-wide association studies. We also illustrate the behavior of synthetic associations in real datasets by showing that rare causal mutations responsible for both hearing loss and sickle cell anemia create genome-wide significant synthetic associations, in the latter case extending over a 2.5-Mb interval encompassing scores of "blocks" of associated variants. In conclusion, uncommon or rare genetic variants can easily create synthetic associations that are credited to common variants, and this possibility requires careful consideration in the interpretation and follow up of GWAS signals.

Title: Rare and common variants: twenty arguments Content: Genome-wide association studies have greatly improved our understanding of the genetic basis of disease risk. The fact that they tend not to identify more than a fraction of the specific causal loci has led to divergence of opinion over whether most of the variance is hidden as numerous rare variants of large effect or as common variants of very small effect. Here I review 20 arguments for and against each of these models of the genetic basis of complex traits and conclude that both classes of effect can be readily reconciled.

Title: A fine-scale map of recombination rates and hotspots across the human genome. Content: Genetic maps, which document the way in which recombination rates vary over a genome, are an essential tool for many genetic analyses. We present a high-resolution genetic map of the human genome, based on statistical analyses of genetic variation data, and identify more than 25,000 recombination hotspots, together with motifs and sequence contexts that play a role in hotspot activity. Differences between the behavior of recombination rates over large (megabase) and small (kilobase) scales lead us to suggest a two-stage model for recombination in which hotspots are stochastic features, within a framework in which large-scale rates are constrained.

Title: Large-Scale Profiling Reveals the Influence of Genetic Variation on Gene Expression in Human Induced Pluripotent Stem Cells. Content: In this study, we used whole-genome sequencing and gene expression profiling of 215 human induced pluripotent stem cell (iPSC) lines from different donors to identify genetic variants associated with RNA expression for 5,746 genes. We were able to predict causal variants for these expression quantitative trait loci (eQTLs) that disrupt transcription factor binding and validated a subset of them experimentally. We also identified copy-number variant (CNV) eQTLs, including some that appear to affect gene expression by altering the copy number of intergenic regulatory regions. In addition, we were able to identify effects on gene expression of rare genic CNVs and regulatory single-nucleotide variants and found that reactivation of gene expression on the X chromosome depends on gene chromosomal position. Our work highlights the value of iPSCs for genetic association analyses and provides a unique resource for investigating the genetic regulation of gene expression in pluripotent cells.

Title: Common and rare variants in multifactorial susceptibility to common diseases Content: Here, we give a historical overview of the search for genetic variants that influence the susceptibility of an individual to a chronic disease, from RA Fisher's seminal work to the current excitement of whole-genome association studies (WGAS). We then discuss the concepts behind the identification of common variants as disease causal factors and contrast them to the basic ideas that underlie the rare variant hypothesis. The identification of rare variants involves the careful selection of candidate genes to examine, the availability of highly efficient resequencing techniques and the appropriate assessment of the functional consequences of the implicated variant. We believe that this strategy can be successfully applied at present in order to unravel the contribution of rare variants to the multifactorial inheritance of common diseases, which could lead to the implementation of much needed preventative screening schemes.

Title: Linkage Disequilibrium Mapping of CHEK2: Common Variation and Breast Cancer Risk Content: Background Checkpoint kinase 2 (CHEK2) averts cancer development by promoting cell cycle arrest and activating DNA repair in genetically damaged cells. Previous investigation has established a role for the CHEK2 gene in breast cancer aetiology, but studies have largely been limited to the rare 1100delC mutation. Whether common polymorphisms in this gene influence breast cancer risk remains unknown. In this study, we aimed to assess the importance of common CHEK2 variants on population risk for breast cancer by capturing the majority of diversity in the gene using haplotype tagging single nucleotide polymorphisms (tagSNPs). Methods and Findings We analyzed 14 common SNPs spanning 52 kilobases (kb) of the CHEK2 gene in 92 Swedish women. Coverage evaluation indicated that these typed SNPs would efficiently convey association signal also from untyped SNPs in the same region. Six of the 14 SNPs predicted well both the haplotypic and single SNP variations within CHEK2. We genotyped these six tagSNPs in 1,577 postmenopausal breast cancer cases and 1,513 population controls, but found no convincing association between any common CHEK2 haplotype and breast cancer risk. The 1100delC mutation was rare in our Swedish population—0.7% in cases and 0.4% in controls— with a corresponding odds ratio for carriers versus noncarriers of 2.26 (95% confidence interval, 0.99–5.15). Estimates of the population frequency and the odds ratio of 1100delC indicate that our sample is representative of a Northern European population.

Title: Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer Content: Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case–control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10−5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10−4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10−9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Title: A common variant in BRCA2 is associated with both breast cancer risk and prenatal viability Content: Inherited mutations in the gene BRCA2 predispose carriers to early onset breast cancer, but such mutations account for fewer than 2% of all cases in East Anglia. It is likely that low penetrance alleles explain the greater part of inherited susceptibility to breast cancer; polymorphic variants in strongly predisposing genes, such as BRCA2, are candidates for this role. BRCA2 is thought to be involved in DNA double strand break-repair. Few mice in which Brca2 is truncated survive to birth; of those that do, most are male, smaller than their normal littermates and have high cancer incidence. Here we show that a common human polymorphism (N372H) in exon 10 of BRCA2 confers an increased risk of breast cancer: the HH homozygotes have a 1.31-fold (95% CI, 1.07–1.61) greater risk than the NN group. Moreover, in normal female controls of all ages there is a significant deficiency of homozygotes compared with that expected from Hardy-Weinberg equilibrium, whereas in males there is an excess of homozygotes: the HH group has an estimated fitness of 0.82 in females and 1.38 in males. Therefore, this variant of BRCA2 appears also to affect fetal survival in a sex-dependent manner.

Title: The transcriptional landscape of the mammalian genome. Content: This study describes comprehensive polling of transcription start and termination sites and analysis of previously unidentified full-length complementary DNAs derived from the mouse genome. We identify the 5' and 3' boundaries of 181,047 transcripts with extensive variation in transcripts arising from alternative promoter usage, splicing, and polyadenylation. There are 16,247 new mouse protein-coding transcripts, including 5154 encoding previously unidentified proteins. Genomic mapping of the transcriptome reveals transcriptional forests, with overlapping transcription on both strands, separated by deserts in which few transcripts are observed. The data provide a comprehensive platform for the comparative analysis of mammalian transcriptional regulation in differentiation and development.

Title: High-throughput sequencing reveals extensive variation in human-specific L1 content in individual human genomes. Content: Using high-throughput sequencing, we devised a technique to determine the insertion sites of virtually all members of the human-specific L1 retrotransposon family in any human genome. Using diagnostic nucleotides, we were able to locate the approximately 800 L1Hs copies corresponding specifically to the pre-Ta, Ta-0, and Ta-1 L1Hs subfamilies, with over 90% of sequenced reads corresponding to human-specific elements. We find that any two individual genomes differ at an average of 285 sites with respect to L1 insertion presence or absence. In total, we assayed 25 individuals, 15 of which are unrelated, at 1139 sites, including 772 shared with the reference genome and 367 nonreference L1 insertions. We show that L1Hs profiles recapitulate genetic ancestry, and determine the chromosomal distribution of these elements. Using these data, we estimate that the rate of L1 retrotransposition in humans is between 1/95 and 1/270 births, and the number of dimorphic L1 elements in the human population with gene frequencies greater than 0.05 is between 3000 and 10,000.

Title: Understanding mechanisms underlying human gene expression variation with RNA sequencing Content: Understanding the genetic mechanisms underlying natural variation in gene expression is a central goal of both medical and evolutionary genetics, and studies of expression quantitative trait loci (eQTLs) have become an important tool for achieving this goal. Although all eQTL studies so far have assayed messenger RNA levels using expression microarrays, recent advances in RNA sequencing enable the analysis of transcript variation at unprecedented resolution. We sequenced RNA from 69 lymphoblastoid cell lines derived from unrelated Nigerian individuals that have been extensively genotyped by the International HapMap Project. By pooling data from all individuals, we generated a map of the transcriptional landscape of these cells, identifying extensive use of unannotated untranslated regions and more than 100 new putative protein-coding exons. Using the genotypes from the HapMap project, we identified more than a thousand genes at which genetic variation influences overall expression levels or splicing. We demonstrate that eQTLs near genes generally act by a mechanism involving allele-specific expression, and that variation that influences the inclusion of an exon is enriched within and near the consensus splice sites. Our results illustrate the power of high-throughput sequencing for the joint analysis of variation in transcription, splicing and allele-specific expression across individuals.

Title: The Sequence Alignment/Map format and SAMtools Content: SUMMARY The Sequence Alignment/Map (SAM) format is a generic alignment format for storing read alignments against reference sequences, supporting short and long reads (up to 128 Mbp) produced by different sequencing platforms. It is flexible in style, compact in size, efficient in random access and is the format in which alignments from the 1000 Genomes Project are released. SAMtools implements various utilities for post-processing alignments in the SAM format, such as indexing, variant caller and alignment viewer, and thus provides universal tools for processing read alignments. AVAILABILITY http://samtools.sourceforge.net.

Title: Meta-analysis identifies common variants associated with body mass index in East Asians Content: Multiple genetic loci associated with obesity or body mass index (BMI) have been identified through genome-wide association studies conducted predominantly in populations of European ancestry. We performed a meta-analysis of associations between BMI and approximately 2.4 million SNPs in 27,715 east Asians, which was followed by in silico and de novo replication studies in 37,691 and 17,642 additional east Asians, respectively. We identified ten BMI-associated loci at genome-wide significance (P < 5.0 × 10(-8)), including seven previously identified loci (FTO, SEC16B, MC4R, GIPR-QPCTL, ADCY3-DNAJC27, BDNF and MAP2K5) and three novel loci in or near the CDKAL1, PCSK1 and GP2 genes. Three additional loci nearly reached the genome-wide significance threshold, including two previously identified loci in the GNPDA2 and TFAP2B genes and a newly identified signal near PAX6, all of which were associated with BMI with P < 5.0 × 10(-7). Findings from this study may shed light on new pathways involved in obesity and demonstrate the value of conducting genetic studies in non-European populations.

Title: The essence of SNPs. Content: Single nucleotide polymorphisms (SNPs) are an abundant form of genome variation, distinguished from rare variations by a requirement for the least abundant allele to have a frequency of 1% or more. A wide range of genetics disciplines stand to benefit greatly from the study and use of SNPs. The recent surge of interest in SNPs stems from, and continues to depend upon, the merging and coincident maturation of several research areas, i.e. (i) large-scale genome analysis and related technologies, (ii) bio-informatics and computing, (iii) genetic analysis of simple and complex disease states, and (iv) global human population genetics. These fields will now be propelled forward, often into uncharted territories, by ongoing discovery efforts that promise to yield hundreds of thousands of human SNPs in the next few years. Major questions are now being asked, experimentally, theoretically and ethically, about the most effective ways to unlock the full potential of the upcoming SNP revolution.

Title: Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes Content: Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Title: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells Content: We report the application of single-molecule-based sequencing technology for high-throughput profiling of histone modifications in mammalian cells. By obtaining over four billion bases of sequence from chromatin immunoprecipitated DNA, we generated genome-wide chromatin-state maps of mouse embryonic stem cells, neural progenitor cells and embryonic fibroblasts. We find that lysine 4 and lysine 27 trimethylation effectively discriminates genes that are expressed, poised for expression, or stably repressed, and therefore reflect cell state and lineage potential. Lysine 36 trimethylation marks primary coding and non-coding transcripts, facilitating gene annotation. Trimethylation of lysine 9 and lysine 20 is detected at satellite, telomeric and active long-terminal repeats, and can spread into proximal unique sequences. Lysine 4 and lysine 9 trimethylation marks imprinting control regions. Finally, we show that chromatin state can be read in an allele-specific manner by using single nucleotide polymorphisms. This study provides a framework for the application of comprehensive chromatin profiling towards characterization of diverse mammalian cell populations.

Title: Improved variant discovery through local re-alignment of short-read next-generation sequencing data using SRMA Content: A primary component of next-generation sequencing analysis is to align short reads to a reference genome, with each read aligned independently. However, reads that observe the same non-reference DNA sequence are highly correlated and can be used to better model the true variation in the target genome. A novel short-read micro realigner, SRMA, that leverages this correlation to better resolve a consensus of the underlying DNA sequence of the targeted genome is described here.

Title: Automated genome sequence analysis and annotation. Content: MOTIVATION Large-scale genome projects generate a rapidly increasing number of sequences, most of them biochemically uncharacterized. Research in bioinformatics contributes to the development of methods for the computational characterization of these sequences. However, the installation and application of these methods require experience and are time consuming. RESULTS We present here an automatic system for preliminary functional annotation of protein sequences that has been applied to the analysis of sets of sequences from complete genomes, both to refine overall performance and to make new discoveries comparable to those made by human experts. The GeneQuiz system includes a Web-based browser that allows examination of the evidence leading to an automatic annotation and offers additional information, views of the results, and links to biological databases that complement the automatic analysis. System structure and operating principles concerning the use of multiple sequence databases, underlying sequence analysis tools, lexical analyses of database annotations and decision criteria for functional assignments are detailed. The system makes automatic quality assessments of results based on prior experience with the underlying sequence analysis tools; overall error rates in functional assignment are estimated at 2.5-5% for cases annotated with highest reliability ('clear' cases). Sources of over-interpretation of results are discussed with proposals for improvement. A conservative definition for reporting 'new findings' that takes account of database maturity is presented along with examples of possible kinds of discoveries (new function, family and superfamily) made by the system. System performance in relation to sequence database coverage, database dynamics and database search methods is analysed, demonstrating the inherent advantages of an integrated automatic approach using multiple databases and search methods applied in an objective and repeatable manner. AVAILABILITY The GeneQuiz system is publicly available for analysis of protein sequences through a Web server at http://www.sander.ebi.ac. uk/gqsrv/submit

Title: Common variants in DGKK are strongly associated with risk of hypospadias Content: Hypospadias is a common congenital malformation of the male external genitalia. We performed a genome-wide association study using pooled DNA from 436 individuals with hypospadias (cases) and 494 controls of European descent and selected the highest ranked SNPs for individual genotyping in the discovery sample, an additional Dutch sample of 133 cases and their parents, and a Swedish series of 266 cases and 402 controls. Individual genotyping of two SNPs (rs1934179 and rs7063116) in DGKK, encoding diacylglycerol kinase κ, produced compelling evidence for association with hypospadias in the discovery sample (allele-specific odds ratio (OR) = 2.5, P = 2.5 × 10−11 and OR = 2.3, P = 2.9 × 10−9, respectively) and in the Dutch (OR = 3.9, P = 2.4 × 10−5 and OR = 3.8, P = 3.4 × 10−5) and Swedish (OR = 2.5, P = 2.6 × 10−8 and OR = 2.2, P = 2.7 × 10−6) replication samples. Expression studies showed expression of DGKK in preputial tissue of cases and controls, which was lower in carriers of the risk allele of rs1934179 (P = 0.047). We propose DGKK as a major risk gene for hypospadias.

Title: Ribosomal DNA copy number loss and sequence variation in cancer Content: Ribosomal DNA is one of the most variable regions in the human genome with respect to copy number. Despite the importance of rDNA for cellular function, we know virtually nothing about what governs its copy number, stability, and sequence in the mammalian genome due to challenges associated with mapping and analysis. We applied computational and droplet digital PCR approaches to measure rDNA copy number in normal and cancer states in human and mouse genomes. We find that copy number and sequence can change in cancer genomes. Counterintuitively, human cancer genomes show a loss of copies, accompanied by global copy number co-variation. The sequence can also be more variable in the cancer genome. Cancer genomes with lower copies have mutational evidence of mTOR hyperactivity. The PTEN phosphatase is a tumor suppressor that is critical for genome stability and a negative regulator of the mTOR kinase pathway. Surprisingly, but consistent with the human cancer genomes, hematopoietic cancer stem cells from a Pten-/- mouse model for leukemia have lower rDNA copy number than normal tissue, despite increased proliferation, rRNA production, and protein synthesis. Loss of copies occurs early and is associated with hypersensitivity to DNA damage. Therefore, copy loss is a recurrent feature in cancers associated with mTOR activation. Ribosomal DNA copy number may be a simple and useful indicator of whether a cancer will be sensitive to DNA damaging treatments.

Title: Gene–environment interactions in 7610 women with breast cancer: prospective evidence from the Million Women Study Content: BACKGROUND Information is scarce about the combined effects on breast cancer incidence of low-penetrance genetic susceptibility polymorphisms and environmental factors (reproductive, behavioural, and anthropometric risk factors for breast cancer). To test for evidence of gene-environment interactions, we compared genotypic relative risks for breast cancer across the other risk factors in a large UK prospective study. METHODS We tested gene-environment interactions in 7610 women who developed breast cancer and 10 196 controls without the disease, studying the effects of 12 polymorphisms (FGFR2-rs2981582, TNRC9-rs3803662, 2q35-rs13387042, MAP3K1-rs889312, 8q24-rs13281615, 2p-rs4666451, 5p12-rs981782, CASP8-rs1045485, LSP1-rs3817198, 5q-rs30099, TGFB1-rs1982073, and ATM-rs1800054) in relation to prospectively collected information about ten established environmental risk factors (age at menarche, parity, age at first birth, breastfeeding, menopausal status, age at menopause, use of hormone replacement therapy, body-mass index, height, and alcohol consumption). FINDINGS After allowance for multiple testing none of the 120 comparisons yielded significant evidence of a gene-environment interaction. By contrast with previous suggestions, there was little evidence that the genotypic relative risks were affected by use of hormone replacement therapy, either overall or for oestrogen-receptor-positive disease. Only one of the 12 polymorphisms was correlated with any of the ten other risk factors: carriers of the high-risk C allele of MAP3K1-rs889312 were significantly shorter than non-carriers (mean height 162.4 cm [95% CI 162.1-162.7] vs 163.1 cm [162.9-163.2]; p=0.01 after allowance for multiple testing). INTERPRETATION Risks of breast cancer associated with low-penetrance susceptibility polymorphisms do not vary significantly with these ten established environmental risk factors. FUNDING Cancer Research UK and the UK Medical Research Council.

Title: Circos: an information aesthetic for comparative genomics. Content: We created a visualization tool called Circos to facilitate the identification and analysis of similarities and differences arising from comparisons of genomes. Our tool is effective in displaying variation in genome structure and, generally, any other kind of positional relationships between genomic intervals. Such data are routinely produced by sequence alignments, hybridization arrays, genome mapping, and genotyping studies. Circos uses a circular ideogram layout to facilitate the display of relationships between pairs of positions by the use of ribbons, which encode the position, size, and orientation of related genomic elements. Circos is capable of displaying data as scatter, line, and histogram plots, heat maps, tiles, connectors, and text. Bitmap or vector images can be created from GFF-style data inputs and hierarchical configuration files, which can be easily generated by automated tools, making Circos suitable for rapid deployment in data analysis and reporting pipelines.

Title: Natural variation in the parameters of innate immune cells is preferentially driven by genetic factors Content: The quantification and characterization of circulating immune cells provide key indicators of human health and disease. To identify the relative effects of environmental and genetic factors on variation in the parameters of innate and adaptive immune cells in homeostatic conditions, we combined standardized flow cytometry of blood leukocytes and genome-wide DNA genotyping of 1,000 healthy, unrelated people of Western European ancestry. We found that smoking, together with age, sex and latent infection with cytomegalovirus, were the main non-genetic factors that affected variation in parameters of human immune cells. Genome-wide association studies of 166 immunophenotypes identified 15 loci that showed enrichment for disease-associated variants. Finally, we demonstrated that the parameters of innate cells were more strongly controlled by genetic variation than were those of adaptive cells, which were driven by mainly environmental exposure. Our data establish a resource that will generate new hypotheses in immunology and highlight the role of innate immunity in susceptibility to common autoimmune diseases. Both environmental factors and genetic factors influence human immunity. Albert and colleagues leverage data from the Milieu Intérieur Consortium to comprehensively describe the effects of lifestyle, environment and genetics on human innate and adaptive immunity.

Title: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer Content: COSMIC (http://www.sanger.ac.uk/cosmic) curates comprehensive information on somatic mutations in human cancer. Release v48 (July 2010) describes over 136,000 coding mutations in almost 542,000 tumour samples; of the 18,490 genes documented, 4803 (26%) have one or more mutations. Full scientific literature curations are available on 83 major cancer genes and 49 fusion gene pairs (19 new cancer genes and 30 new fusion pairs this year) and this number is continually increasing. Key amongst these is TP53, now available through a collaboration with the IARC p53 database. In addition to data from the Cancer Genome Project (CGP) at the Sanger Institute, UK, and The Cancer Genome Atlas project (TCGA), large systematic screens are also now curated. Major website upgrades now make these data much more mineable, with many new selection filters and graphics. A Biomart is now available allowing more automated data mining and integration with other biological databases. Annotation of genomic features has become a significant focus; COSMIC has begun curating full-genome resequencing experiments, developing new web pages, export formats and graphics styles. With all genomic information recently updated to GRCh37, COSMIC integrates many diverse types of mutation information and is making much closer links with Ensembl and other data resources.

Title: Whole genomes redefine the mutational landscape of pancreatic cancer Content: Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Title: Landscape of genetic lesions in 944 patients with myelodysplastic syndromes Content: High-throughput DNA sequencing significantly contributed to diagnosis and prognostication in patients with myelodysplastic syndromes (MDS). We determined the biological and prognostic significance of genetic aberrations in MDS. In total, 944 patients with various MDS subtypes were screened for known/putative mutations/deletions in 104 genes using targeted deep sequencing and array-based genomic hybridization. In total, 845/944 patients (89.5%) harbored at least one mutation (median, 3 per patient; range, 0-12). Forty-seven genes were significantly mutated with TET2, SF3B1, ASXL1, SRSF2, DNMT3A, and RUNX1 mutated in >10% of cases. Many mutations were associated with higher risk groups and/or blast elevation. Survival was investigated in 875 patients. By univariate analysis, 25/48 genes (resulting from 47 genes tested significantly plus PRPF8) affected survival (P<0.05). The status of 14 genes combined with conventional factors revealed a novel prognostic model ('Model-1') separating patients into four risk groups ('low', 'intermediate', 'high', 'very high risk') with 3-year survival of 95.2, 69.3, 32.8, and 5.3% (P<0.001). Subsequently, a 'gene-only model' ('Model-2') was constructed based on 14 genes also yielding four significant risk groups (P<0.001). Both models were reproducible in the validation cohort (n=175 patients; P<0.001 each). Thus, large-scale genetic and molecular profiling of multiple target genes is invaluable for subclassification and prognostication in MDS patients.

Title: Mitochondrial genome variation and the origin of modern humans. Content: The analysis of mitochondrial DNA (mtDNA) has been a potent tool in our understanding of human evolution, owing to characteristics such as high copy number, apparent lack of recombination, high substitution rate and maternal mode of inheritance. However, almost all studies of human evolution based on mtDNA sequencing have been confined to the control region, which constitutes less than 7% of the mitochondrial genome. These studies are complicated by the extreme variation in substitution rate between sites, and the consequence of parallel mutations causing difficulties in the estimation of genetic distance and making phylogenetic inferences questionable. Most comprehensive studies of the human mitochondrial molecule have been carried out through restriction-fragment length polymorphism analysis, providing data that are ill suited to estimations of mutation rate and therefore the timing of evolutionary events. Here, to improve the information obtained from the mitochondrial molecule for studies of human evolution, we describe the global mtDNA diversity in humans based on analyses of the complete mtDNA sequence of 53 humans of diverse origins. Our mtDNA data, in comparison with those of a parallel study of the Xq13.3 region in the same individuals, provide a concurrent view on human evolution with respect to the age of modern humans.

Title: Comprehensive genomic characterization defines human glioblastoma genes and core pathways Content: Human cancer cells typically harbour multiple chromosomal aberrations, nucleotide substitutions and epigenetic modifications that drive malignant transformation. The Cancer Genome Atlas ( TCGA) pilot project aims to assess the value of large- scale multi- dimensional analysis of these molecular characteristics in human cancer and to provide the data rapidly to the research community. Here we report the interim integrative analysis of DNA copy number, gene expression and DNA methylation aberrations in 206 glioblastomas - the most common type of primary adult brain cancer - and nucleotide sequence aberrations in 91 of the 206 glioblastomas. This analysis provides new insights into the roles of ERBB2, NF1 and TP53, uncovers frequent mutations of the phosphatidylinositol- 3- OH kinase regulatory subunit gene PIK3R1, and provides a network view of the pathways altered in the development of glioblastoma. Furthermore, integration of mutation, DNA methylation and clinical treatment data reveals a link between MGMT promoter methylation and a hypermutator phenotype consequent to mismatch repair deficiency in treated glioblastomas, an observation with potential clinical implications. Together, these findings establish the feasibility and power of TCGA, demonstrating that it can rapidly expand knowledge of the molecular basis of cancer.

Title: Large-scale genetic fine mapping and genotype-phenotype associations implicate polymorphism in the IL2RA region in type 1 diabetes Content: Genome-wide association studies are now identifying disease-associated chromosome regions. However, even after convincing replication, the localization of the causal variant(s) requires comprehensive resequencing, extensive genotyping and statistical analyses in large sample sets leading to targeted functional studies. Here, we have localized the type 1 diabetes (T1D) association in the interleukin 2 receptor alpha (IL2RA) gene region to two independent groups of SNPs, spanning overlapping regions of 14 and 40 kb, encompassing IL2RA intron 1 and the 5′ regions of IL2RA and RBM17 (odds ratio = 2.04, 95% confidence interval = 1.70–2.45; P = 1.92 × 10−28; control frequency = 0.635). Furthermore, we have associated IL2RA T1D susceptibility genotypes with lower circulating levels of the biomarker, soluble IL-2RA (P = 6.28 × 10−28), suggesting that an inherited lower immune responsiveness predisposes to T1D.

Title: Common BRCA1 variants and susceptibility to breast and ovarian cancer in the general population. Content: Most multiple case families of young onset breast cancer and ovarian cancer are thought to be due to highly penetrant mutations in the predisposing genes BRCA1 and BRCA2. However, these mutations are uncommon in the population and they probably account for only a few percent of all breast cancer incidence. A much larger fraction of breast cancer might, in principle, be due to common variants which confer more modest individual risks. There are several common polymorphisms in the BRCA1 gene which generate amino acid substitutions. We have examined the frequency of four of these polymorphisms: Gln356Arg, Pro871Leu, Glu1038Gly and Ser1613Gly in large series of breast and ovarian cancer cases and matched controls. Due to strong linkage disequilibrium, these four sites generate only three haplotypes with a frequency > 1.3%. The most common haplotypes, defined by the alleles Gln356Pro871Glu1038Ser1613 and Gln356Leu871Gly1038Gly1613, have frequencies of 0.57 and 0.32 respectively, and these frequencies do not differ significantly between patient and control groups. Thus the most common polymorphisms of the BRCA1 gene do not make a significant contribution to breast or ovarian cancer risk. However, our data suggest that the Arg356 allele may have a different genotype distribution in breast cancer patients from that in controls (Arg356 homozygotes are more frequent in the control groups, P = 0.01), indicating that it may be protective against breast cancer. If this finding can be confirmed, it may provide an insight into the structural features of the BRCA1 protein that are important for its function.

Title: LINE-1 Retrotransposition Activity in Human Genomes Content: Highly active (i.e., "hot") long interspersed element-1 (LINE-1 or L1) sequences comprise the bulk of retrotransposition activity in the human genome; however, the abundance of hot L1s in the human population remains largely unexplored. Here, we used a fosmid-based, paired-end DNA sequencing strategy to identify 68 full-length L1s that are differentially present among individuals but are absent from the human genome reference sequence. The majority of these L1s were highly active in a cultured cell retrotransposition assay. Genotyping 26 elements revealed that two L1s are only found in Africa and that two more are absent from the H952 subset of the Human Genome Diversity Panel. Therefore, these results suggest that hot L1s are more abundant in the human population than previously appreciated, and that ongoing L1 retrotransposition continues to be a major source of interindividual genetic variation.

Title: Identification of Gene Positioning Factors Using High-Throughput Imaging Mapping Content: Genomes are arranged non-randomly in the 3D space of the cell nucleus. Here, we have developed HIPMap, a high-precision, high-throughput, automated fluorescent in situ hybridization imaging pipeline, for mapping of the spatial location of genome regions at large scale. High-throughput imaging position mapping (HIPMap) enabled an unbiased siRNA screen for factors involved in genome organization in human cells. We identify 50 cellular factors required for proper positioning of a set of functionally diverse genomic loci. Positioning factors include chromatin remodelers, histone modifiers, and nuclear envelope and pore proteins. Components of the replication and post-replication chromatin re-assembly machinery are prominently represented among positioning factors, and timely progression of cells through replication, but not mitosis, is required for correct gene positioning. Our results establish a method for the large-scale mapping of genome locations and have led to the identification of a compendium of cellular factors involved in spatial genome organization.

Title: Haplotype-based variant detection from short-read sequencing Content: The direct detection of haplotypes from short-read DNA sequencing data requires changes to existing small-variant detection methods. Here, we develop a Bayesian statistical framework which is capable of modeling multiallelic loci in sets of individuals with non-uniform copy number. We then describe our implementation of this framework in a haplotype-based variant detector, FreeBayes.

Title: High-resolution mapping of open chromatin in the rice genome. Content: Gene expression is controlled by the complex interaction of transcription factors binding to promoters and other regulatory DNA elements. One common characteristic of the genomic regions associated with regulatory proteins is a pronounced sensitivity to DNase I digestion. We generated genome-wide high-resolution maps of DNase I hypersensitive (DH) sites from both seedling and callus tissues of rice (Oryza sativa). Approximately 25% of the DH sites from both tissues were found in putative promoters, indicating that the vast majority of the gene regulatory elements in rice are not located in promoter regions. We found 58% more DH sites in the callus than in the seedling. For DH sites detected in both the seedling and callus, 31% displayed significantly different levels of DNase I sensitivity within the two tissues. Genes that are differentially expressed in the seedling and callus were frequently associated with DH sites in both tissues. The DNA sequences contained within the DH sites were hypomethylated, consistent with what is known about active gene regulatory elements. Interestingly, tissue-specific DH sites located in the promoters showed a higher level of DNA methylation than the average DNA methylation level of all the DH sites located in the promoters. A distinct elevation of H3K27me3 was associated with intergenic DH sites. These results suggest that epigenetic modifications play a role in the dynamic changes of the numbers and DNase I sensitivity of DH sites during development.

Title: Genetic Variants at Newly Identified Lipid Loci Are Associated with Coronary Heart Disease in a Chinese Han Population Content: BACKGROUND Recent genome-wide association studies (GWAS) have mapped several novel loci influencing blood lipid levels in Caucasians. We sought to explore whether the genetic variants at newly identified lipid-associated loci were associated with CHD susceptibility in a Chinese Han population. METHODOLOGY/PRINCIPAL FINDINGS We conducted a two-stage case-control study in a Chinese Han population. The first-stage, consisting of 1,376 CHD cases and 1,376 sex and age- frequency matched controls, examined 5 novel lipid-associated single-nucleotide polymorphisms (SNPs) identified from GWAS among Caucasians in relation to CHD risk in Chinese. We then validated significant SNPs in the second-stage, consisting of 1,269 cases and 2,745 controls. We also tested associations between SNPs within the five novel loci and blood lipid levels in 4,121 controls. We identified two novel SNPs (rs599839 in CELSR2-PSRC1-SORT1 and rs16996148 in NCAN-CILP2) that were significantly associated with reduced CHD risk in Chinese (odds ratios (95% confidence intervals) in the dominant model 0.76 (0.61-0.90; P = 0.001), 0.67 (0.57-0.77; P = 3.4×10(-8)), respectively). Multiple linear regression analyses using dominant model showed that rs599839 was significantly associated with decreased LDL levels (P = 0.022) and rs16996148 was significantly associated with increased LDL and HDL levels (P = 2.9×10(-4) and 0.001, respectively). CONCLUSIONS/SIGNIFICANCE We identified two novel SNPs (rs599839 and rs16996148) at newly identified lipid-associated loci that were significantly associated with CHD susceptibility in a Chinese Han population.

Title: Integrated maps in quail (Coturnix japonica) confirm the high degree of synteny conservation with chicken (Gallus gallus) despite 35 million years of divergence Content: BackgroundBy comparing the quail genome with that of chicken, chromosome rearrangements that have occurred in these two galliform species over 35 million years of evolution can be detected. From a more practical point of view, the definition of conserved syntenies helps to predict the position of genes in quail, based on information taken from the chicken sequence, thus enhancing the utility of this species in biological studies through a better knowledge of its genome structure. A microsatellite and an Amplified Fragment Length Polymorphism (AFLP) genetic map were previously published for quail, as well as comparative cytogenetic data with chicken for macrochromosomes. Quail genomics will benefit from the extension and the integration of these maps. ResultsThe integrated linkage map presented here is based on segregation analysis of both anonymous markers and functional gene loci in 1,050 quail from three independent F2 populations. Ninety-two loci are resolved into 14 autosomal linkage groups and a Z chromosome-specific linkage group, aligned with the quail AFLP map. The size of linkage groups ranges from 7.8 cM to 274.8 cM. The total map distance covers 904.3 cM with an average spacing of 9.7 cM between loci. The coverage is not complete, as macrochromosome CJA08, the gonosome CJAW and 23 microchromosomes have no marker assigned yet. Significant sequence identities of quail markers with chicken enabled the alignment of the quail linkage groups on the chicken genome sequence assembly. This, together with interspecific Fluorescence In Situ Hybridization (FISH), revealed very high similarities in marker order between the two species for the eight macrochromosomes and the 14 microchromosomes studied. ConclusionIntegrating the two microsatellite and the AFLP quail genetic maps greatly enhances the quality of the resulting information and will thus facilitate the identification of Quantitative Trait Loci (QTL). The alignment with the chicken chromosomes confirms the high conservation of gene order that was expected between the two species for macrochromosomes. By extending the comparative study to the microchromosomes, we suggest that a wealth of information can be mined in chicken, to be used for genome analyses in quail.

Title: Inter- and Intra-Host Viral Diversity in a Large Seasonal DENV2 Outbreak Content: BACKGROUND High genetic diversity at both inter- and intra-host level are hallmarks of RNA viruses due to the error-prone nature of their genome replication. Several groups have evaluated the extent of viral variability using different RNA virus deep sequencing methods. Although much of this effort has been dedicated to pathogens that cause chronic infections in humans, few studies investigated arthropod-borne, acute viral infections. METHODS AND PRINCIPAL FINDINGS We deep sequenced the complete genome of ten DENV2 isolates from representative classical and severe cases sampled in a large outbreak in Brazil using two different approaches. Analysis of the consensus genomes confirmed the larger extent of the 2010 epidemic in comparison to a previous epidemic caused by the same viruses in another city two years before (genetic distance = 0.002 and 0.0008 respectively). Analysis of viral populations within the host revealed a high level of conservation. After excluding homopolymer regions of 454/Roche generated sequences, we found 10 to 44 variable sites per genome population at a frequency of >1%, resulting in very low intra-host genetic diversity. While up to 60% of all variable sites at intra-host level were non-synonymous changes, only 10% of inter-host variability resulted from non-synonymous mutations, indicative of purifying selection at the population level. CONCLUSIONS AND SIGNIFICANCE Despite the error-prone nature of RNA-dependent RNA-polymerase, dengue viruses maintain low levels of intra-host variability.

Title: The role of regulatory variation in complex traits and disease Content: We are in a phase of unprecedented progress in identifying genetic loci that cause variation in traits ranging from growth and fitness in simple organisms to disease in humans. However, a mechanistic understanding of how these loci influence traits is lacking for the majority of loci. Studies of the genetics of gene expression have emerged as a key tool for linking DNA sequence variation to phenotypes. Here, we review recent insights into the molecular nature of regulatory variants and describe their influence on the transcriptome and the proteome. We discuss conceptual advances from studies in model organisms and present examples of complete chains of causality that link individual polymorphisms to changes in gene expression, which in turn result in physiological changes and, ultimately, disease risk.

Title: Genetics and Beyond – The Transcriptome of Human Monocytes and Disease Susceptibility Content: BACKGROUND Variability of gene expression in human may link gene sequence variability and phenotypes; however, non-genetic variations, alone or in combination with genetics, may also influence expression traits and have a critical role in physiological and disease processes. METHODOLOGY/PRINCIPAL FINDINGS To get better insight into the overall variability of gene expression, we assessed the transcriptome of circulating monocytes, a key cell involved in immunity-related diseases and atherosclerosis, in 1,490 unrelated individuals and investigated its association with >675,000 SNPs and 10 common cardiovascular risk factors. Out of 12,808 expressed genes, 2,745 expression quantitative trait loci were detected (P<5.78x10(-12)), most of them (90%) being cis-modulated. Extensive analyses showed that associations identified by genome-wide association studies of lipids, body mass index or blood pressure were rarely compatible with a mediation by monocyte expression level at the locus. At a study-wide level (P<3.9x10(-7)), 1,662 expression traits (13.0%) were significantly associated with at least one risk factor. Genome-wide interaction analyses suggested that genetic variability and risk factors mostly acted additively on gene expression. Because of the structure of correlation among expression traits, the variability of risk factors could be characterized by a limited set of independent gene expressions which may have biological and clinical relevance. For example expression traits associated with cigarette smoking were more strongly associated with carotid atherosclerosis than smoking itself. CONCLUSIONS/SIGNIFICANCE This study demonstrates that the monocyte transcriptome is a potent integrator of genetic and non-genetic influences of relevance for disease pathophysiology and risk assessment.

Title: Haplotype estimation using sequencing reads. Content: High-throughput sequencing technologies produce short sequence reads that can contain phase information if they span two or more heterozygote genotypes. This information is not routinely used by current methods that infer haplotypes from genotype data. We have extended the SHAPEIT2 method to use phase-informative sequencing reads to improve phasing accuracy. Our model incorporates the read information in a probabilistic model through base quality scores within each read. The method is primarily designed for high-coverage sequence data or data sets that already have genotypes called. One important application is phasing of single samples sequenced at high coverage for use in medical sequencing and studies of rare diseases. Our method can also use existing panels of reference haplotypes. We tested the method by using a mother-father-child trio sequenced at high-coverage by Illumina together with the low-coverage sequence data from the 1000 Genomes Project (1000GP). We found that use of phase-informative reads increases the mean distance between switch errors by 22% from 274.4 kb to 328.6 kb. We also used male chromosome X haplotypes from the 1000GP samples to simulate sequencing reads with varying insert size, read length, and base error rate. When using short 100 bp paired-end reads, we found that using mixtures of insert sizes produced the best results. When using longer reads with high error rates (5-20 kb read with 4%-15% error per base), phasing performance was substantially improved.

Title: High density mapping of the MHC identifies a shared role for HLA-DRB1*01:03 in inflammatory bowel diseases and heterozygous advantage in ulcerative colitis Content: Genome-wide association studies of the related chronic inflammatory bowel diseases (IBD) known as Crohn's disease and ulcerative colitis have shown strong evidence of association to the major histocompatibility complex (MHC). This region encodes a large number of immunological candidates, including the antigen-presenting classical human leukocyte antigen (HLA) molecules. Studies in IBD have indicated that multiple independent associations exist at HLA and non-HLA genes, but they have lacked the statistical power to define the architecture of association and causal alleles. To address this, we performed high-density SNP typing of the MHC in >32,000 individuals with IBD, implicating multiple HLA alleles, with a primary role for HLA-DRB1*01:03 in both Crohn's disease and ulcerative colitis. Noteworthy differences were observed between these diseases, including a predominant role for class II HLA variants and heterozygous advantage observed in ulcerative colitis, suggesting an important role of the adaptive immune response in the colonic environment in the pathogenesis of IBD.

Title: Genetic variation in West Nile virus from naturally infected mosquitoes and birds suggests quasispecies structure and strong purifying selection. Content: Intrahost genetic diversity was analysed in naturally infected mosquitoes and birds to determine whether West Nile virus (WNV) exists in nature as a quasispecies and to quantify selective pressures within and between hosts. WNV was sampled from ten infected birds and ten infected mosquito pools collected on Long Island, NY, USA, during the peak of the 2003 WNV transmission season. A 1938 nt fragment comprising the 3' 1159 nt of the WNV envelope (E) coding region and the 5' 779 nt of the non-structural protein 1 (NS1) coding region was amplified and cloned and 20 clones per specimen were sequenced. Results from this analysis demonstrate that WNV infections are derived from a genetically diverse population of genomes in nature. The mean nucleotide diversity was 0.016 % within individual specimens and the mean percentage of clones that differed from the consensus sequence was 19.5 %. WNV sequences in mosquitoes were significantly more genetically diverse than WNV in birds. No host-dependent bias for particular types of mutations was observed and estimates of genetic diversity did not differ significantly between E and NS1 coding sequences. Non-consensus clones obtained from two avian specimens had highly similar genetic signatures, providing preliminary evidence that WNV genetic diversity may be maintained throughout the enzootic transmission cycle, rather than arising independently during each infection. Evidence of purifying selection was obtained from both intra- and interhost WNV populations. Combined, these data support the observation that WNV populations may be structured as a quasispecies and document strong purifying natural selection in WNV populations.

Title: Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Content: Cataloging the very large number of undescribed species of insects could be greatly accelerated by automated DNA based approaches, but procedures for large-scale species discovery from sequence data are currently lacking. Here, we use mitochondrial DNA variation to delimit species in a poorly known beetle radiation in the genus Rivacindela from arid Australia. Among 468 individuals sampled from 65 sites and multiple morphologically distinguishable types, sequence variation in three mtDNA genes (cytochrome oxidase subunit 1, cytochrome b, 16S ribosomal RNA) was strongly partitioned between 46 or 47 putative species identified with quantitative methods of species recognition based on fixed unique ("diagnostic") characters. The boundaries between groups were also recognizable from a striking increase in branching rate in clock-constrained calibrated trees. Models of stochastic lineage growth (Yule models) were combined with coalescence theory to develop a new likelihood method that determines the point of transition from species-level (speciation and extinction) to population-level (coalescence) evolutionary processes. Fitting the location of the switches from speciation to coalescent nodes on the ultrametric tree of Rivacindela produced a transition in branching rate occurring at 0.43 Mya, leading to an estimate of 48 putative species (confidence interval for the threshold ranging from 47 to 51 clusters within 2 logL units). Entities delimited in this way exhibited biological properties of traditionally defined species, showing coherence of geographic ranges, broad congruence with morphologically recognized species, and levels of sequence divergence typical for closely related species of insects. The finding of discontinuous evolutionary groupings that are readily apparent in patterns of sequence variation permits largely automated species delineation from DNA surveys of local communities as a scaffold for taxonomy in this poorly known insect group.

Title: Genetic variability in a frozen batch of MCF-7 cells invisible in routine authentication affecting cell function Content: Common recommendations for cell line authentication, annotation and quality control fall short addressing genetic heterogeneity. Within the Human Toxome Project, we demonstrate that there can be marked cellular and phenotypic heterogeneity in a single batch of the human breast adenocarcinoma cell line MCF-7 obtained directly from a cell bank that are invisible with the usual cell authentication by short tandem repeat (STR) markers. STR profiling just fulfills the purpose of authentication testing, which is to detect significant cross-contamination and cell line misidentification. Heterogeneity needs to be examined using additional methods. This heterogeneity can have serious consequences for reproducibility of experiments as shown by morphology, estrogenic growth dose-response, whole genome gene expression and untargeted mass-spectroscopy metabolomics for MCF-7 cells. Using Comparative Genomic Hybridization (CGH), differences were traced back to genetic heterogeneity already in the cells from the original frozen vials from the same ATCC lot, however, STR markers did not differ from ATCC reference for any sample. These findings underscore the need for additional quality assurance in Good Cell Culture Practice and cell characterization, especially using other methods such as CGH to reveal possible genomic heterogeneity and genetic drifts within cell lines.

Title: MotifMap: integrative genome-wide maps of regulatory motif sites for model species Content: BACKGROUND A central challenge of biology is to map and understand gene regulation on a genome-wide scale. For any given genome, only a small fraction of the regulatory elements embedded in the DNA sequence have been characterized, and there is great interest in developing computational methods to systematically map all these elements and understand their relationships. Such computational efforts, however, are significantly hindered by the overwhelming size of non-coding regions and the statistical variability and complex spatial organizations of regulatory elements and interactions. Genome-wide catalogs of regulatory elements for all model species simply do not yet exist. RESULTS The MotifMap system uses databases of transcription factor binding motifs, refined genome alignments, and a comparative genomic statistical approach to provide comprehensive maps of candidate regulatory elements encoded in the genomes of model species. The system is used to derive new genome-wide maps for yeast, fly, worm, mouse, and human. The human map contains 519,108 sites for 570 matrices with a False Discovery Rate of 0.1 or less. The new maps are assessed in several ways, for instance using high-throughput experimental ChIP-seq data and AUC statistics, providing strong evidence for their accuracy and coverage. The maps can be usefully integrated with many other kinds of omic data and are available at http://motifmap.igb.uci.edu/. CONCLUSIONS MotifMap and its integration with other data provide a foundation for analyzing gene regulation on a genome-wide scale, and for automatically generating regulatory pathways and hypotheses. The power of this approach is demonstrated and discussed using the P53 apoptotic pathway and the Gli hedgehog pathways as examples.

Title: Inter-individual variability and genetic influences on cytokine responses to bacteria and fungi Content: Little is known about the inter-individual variation of cytokine responses to different pathogens in healthy individuals. To systematically describe cytokine responses elicited by distinct pathogens and to determine the effect of genetic variation on cytokine production, we profiled cytokines produced by peripheral blood mononuclear cells from 197 individuals of European origin from the 200 Functional Genomics (200FG) cohort in the Human Functional Genomics Project (http://www.humanfunctionalgenomics.org), obtained over three different years. We compared bacteria- and fungi-induced cytokine profiles and found that most cytokine responses were organized around a physiological response to specific pathogens, rather than around a particular immune pathway or cytokine. We then correlated genome-wide single-nucleotide polymorphism (SNP) genotypes with cytokine abundance and identified six cytokine quantitative trait loci (QTLs). Among them, a cytokine QTL at the NAA35-GOLM1 locus markedly modulated interleukin (IL)-6 production in response to multiple pathogens and was associated with susceptibility to candidemia. Furthermore, the cytokine QTLs that we identified were enriched among SNPs previously associated with infectious diseases and heart diseases. These data reveal and begin to explain the variability in cytokine production by human immune cells in response to pathogens.

Title: Trans-ancestry mutational landscape of hepatocellular carcinoma genomes Content: Diverse epidemiological factors are associated with hepatocellular carcinoma (HCC) prevalence in different populations. However, the global landscape of the genetic changes in HCC genomes underpinning different epidemiological and ancestral backgrounds still remains uncharted. Here a collection of data from 503 liver cancer genomes from different populations uncovered 30 candidate driver genes and 11 core pathway modules. Furthermore, a collaboration of two large-scale cancer genome projects comparatively analyzed the trans-ancestry substitution signatures in 608 liver cancer cases and identified unique mutational signatures that predominantly contribute to Asian cases. This work elucidates previously unexplored ancestry-associated mutational processes in HCC development. A combination of hotspot TERT promoter mutation, TERT focal amplification and viral genome integration occurs in more than 68% of cases, implicating TERT as a central and ancestry-independent node of hepatocarcinogenesis. Newly identified alterations in genes encoding metabolic enzymes, chromatin remodelers and a high proportion of mTOR pathway activations offer potential therapeutic and diagnostic opportunities.

Title: Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease Content: The gene encoding apolipoprotein E (APOE) on chromosome 19 is the only confirmed susceptibility locus for late-onset Alzheimer's disease. To identify other risk loci, we conducted a large genome-wide association study of 2,032 individuals from France with Alzheimer's disease (cases) and 5,328 controls. Markers outside APOE with suggestive evidence of association (P < 10−5) were examined in collections from Belgium, Finland, Italy and Spain totaling 3,978 Alzheimer's disease cases and 3,297 controls. Two loci gave replicated evidence of association: one within CLU (also called APOJ), encoding clusterin or apolipoprotein J, on chromosome 8 (rs11136000, OR = 0.86, 95% CI 0.81–0.90, P = 7.5 × 10−9 for combined data) and the other within CR1, encoding the complement component (3b/4b) receptor 1, on chromosome 1 (rs6656401, OR = 1.21, 95% CI 1.14–1.29, P = 3.7 × 10−9 for combined data). Previous biological studies support roles of CLU and CR1 in the clearance of β amyloid (Aβ) peptide, the principal constituent of amyloid plaques, which are one of the major brain lesions of individuals with Alzheimer's disease.

Title: Genetic diagnosis and prognosis of Alzheimer's disease: challenges and opportunities. Content: Alzheimer's disease (AD), the most common form of dementia in western societies, is a pathologically and clinically heterogeneous disease with a strong genetic component. The recent advances in high-throughput genome technologies allowing for the rapid analysis of millions of polymorphisms in thousands of subjects has significantly advanced our understanding of the genomic underpinnings of AD susceptibility. During the last 5 years, genome-wide association and whole-exome- and whole-genome sequencing studies have mapped more than 20 disease-associated loci, providing insights into the molecular pathways involved in AD pathogenesis and hinting at potential novel therapeutic targets. This review article summarizes the challenges and opportunities of when using genomic information for the diagnosis and prognosis of AD.

Title: The genetics of Mexico recapitulates Native American substructure and affects biomedical traits Content: Mexico harbors great cultural and ethnic diversity, yet fine-scale patterns of human genome-wide variation from this region remain largely uncharacterized. We studied genomic variation within Mexico from over 1000 individuals representing 20 indigenous and 11 mestizo populations. We found striking genetic stratification among indigenous populations within Mexico at varying degrees of geographic isolation. Some groups were as differentiated as Europeans are from East Asians. Pre-Columbian genetic substructure is recapitulated in the indigenous ancestry of admixed mestizo individuals across the country. Furthermore, two independently phenotyped cohorts of Mexicans and Mexican Americans showed a significant association between subcontinental ancestry and lung function. Thus, accounting for fine-scale ancestry patterns is critical for medical and population genetic studies within Mexico, in Mexican-descent populations, and likely in many other populations worldwide.

Title: Genetic variation and gastric cancer risk: a field synopsis and meta-analysis. Content: BACKGROUND Data on genetic susceptibility to sporadic gastric carcinoma have been published at a growing pace, but to date no comprehensive overview and quantitative summary has been available. METHODS We conducted a systematic review and meta-analysis of the evidence on the association between DNA variation and risk of developing stomach cancer. To assess result credibility, summary evidence was graded according to the Venice criteria and false positive report probability (FPRP) was calculated to further validate result noteworthiness. Meta-analysis was also conducted for subgroups, which were defined by ethnicity (Asian vs Caucasian), tumour histology (intestinal vs diffuse), tumour site (cardia vs non-cardia) and Helicobacter pylori infection status (positive vs negative). RESULTS Literature search identified 824 eligible studies comprising 2 530 706 subjects (cases: 261 386 (10.3%)) and investigating 2841 polymorphisms involving 952 distinct genes. Overall, we performed 456 primary and subgroup meta-analyses on 156 variants involving 101 genes. We identified 11 variants significantly associated with disease risk and assessed to have a high level of summary evidence: MUC1 rs2070803 at 1q22 (diffuse carcinoma subgroup), MTX1 rs2075570 at 1q22 (diffuse), PSCA rs2294008 at 8q24.2 (non-cardia), PRKAA1 rs13361707 5p13 (non-cardia), PLCE1 rs2274223 10q23 (cardia), TGFBR2 rs3087465 3p22 (Asian), PKLR rs3762272 1q22 (diffuse), PSCA rs2976392 (intestinal), GSTP1 rs1695 11q13 (Asian), CASP8 rs3834129 2q33 (mixed) and TNF rs1799724 6p21.3 (mixed), with the first nine variants characterised by a low FPRP. We also identified polymorphisms with lower quality significant associations (n=110). CONCLUSIONS We have identified several high-quality biomarkers of gastric cancer susceptibility. These data will form the backbone of an annually updated online resource that will be integral to the study of gastric carcinoma genetics and may inform future screening programmes.

Title: Global analysis of ATM polymorphism reveals significant functional constraint. Content: ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.

Title: Mapping Meiotic Single-Strand DNA Reveals a New Landscape of DNA Double-Strand Breaks in Saccharomyces cerevisiae Content: DNA double-strand breaks (DSBs), which are formed by the Spo11 protein, initiate meiotic recombination. Previous DSB-mapping studies have used rad50S or sae2Δ mutants, which are defective in break processing, to accumulate Spo11-linked DSBs, and report large (≥ 50 kb) “DSB-hot” regions that are separated by “DSB-cold” domains of similar size. Substantial recombination occurs in some DSB-cold regions, suggesting that DSB patterns are not normal in rad50S or sae2Δ mutants. We therefore developed a novel method to map genome-wide, single-strand DNA (ssDNA)–associated DSBs that accumulate in processing-capable, repair-defective dmc1Δ and dmc1Δ rad51Δ mutants. DSBs were observed at known hot spots, but also in most previously identified “DSB-cold” regions, including near centromeres and telomeres. Although approximately 40% of the genome is DSB-cold in rad50S mutants, analysis of meiotic ssDNA from dmc1Δ shows that most of these regions have substantial DSB activity. Southern blot assays of DSBs in selected regions in dmc1Δ, rad50S, and wild-type cells confirm these findings. Thus, DSBs are distributed much more uniformly than was previously believed. Comparisons of DSB signals in dmc1, dmc1 rad51, and dmc1 spo11 mutant strains identify Dmc1 as a critical strand-exchange activity genome-wide, and confirm previous conclusions that Spo11-induced lesions initiate all meiotic recombination.

Title: Resequencing studies of nonmodel organisms using closely related reference genomes: optimal experimental designs and bioinformatics approaches for population genomics. Content: Decreasing costs of next-generation sequencing (NGS) experiments have made a wide range of genomic questions open for study with nonmodel organisms. However, experimental designs and analysis of NGS data from less well-known species are challenging because of the lack of genomic resources. In this work, we investigate the performance of alternative experimental designs and bioinformatics approaches in estimating variability and neutrality tests based on the site-frequency-spectrum (SFS) from individual resequencing data. We pay particular attention to challenges faced in the study of nonmodel organisms, in particular the absence of a species-specific reference genome, although phylogenetically close genomes are assumed to be available. We compare the performance of three alternative bioinformatics approaches – genotype calling, genotype–haplotype calling and direct estimation without calling genotypes. We find that relying on genotype calls provides biased estimates of population genetic statistics at low to moderate read depth (2–8X). Genotype–haplotype calling returns more accurate estimates irrespective of the divergence to the reference genome, but requires moderate depth (8–20X). Direct estimation without calling genotypes returns the most accurate estimates of variability and of most SFS tests investigated, including at low read depth (2–4X). Studies without species-specific reference genome should thus aim for low read depth and avoid genotype calling whenever individual genotypes are not essential. Otherwise, aiming for moderate to high depth at the expense of number of individuals, and using genotype–haplotype calling, is recommended.

Title: Candidate Causal Regulatory Effects by Integration of Expression QTLs with Complex Trait Genetic Associations Content: The recent success of genome-wide association studies (GWAS) is now followed by the challenge to determine how the reported susceptibility variants mediate complex traits and diseases. Expression quantitative trait loci (eQTLs) have been implicated in disease associations through overlaps between eQTLs and GWAS signals. However, the abundance of eQTLs and the strong correlation structure (LD) in the genome make it likely that some of these overlaps are coincidental and not driven by the same functional variants. In the present study, we propose an empirical methodology, which we call Regulatory Trait Concordance (RTC) that accounts for local LD structure and integrates eQTLs and GWAS results in order to reveal the subset of association signals that are due to cis eQTLs. We simulate genomic regions of various LD patterns with both a single or two causal variants and show that our score outperforms SNP correlation metrics, be they statistical (r(2)) or historical (D'). Following the observation of a significant abundance of regulatory signals among currently published GWAS loci, we apply our method with the goal to prioritize relevant genes for each of the respective complex traits. We detect several potential disease-causing regulatory effects, with a strong enrichment for immunity-related conditions, consistent with the nature of the cell line tested (LCLs). Furthermore, we present an extension of the method in trans, where interrogating the whole genome for downstream effects of the disease variant can be informative regarding its unknown primary biological effect. We conclude that integrating cellular phenotype associations with organismal complex traits will facilitate the biological interpretation of the genetic effects on these traits.

Title: Genetic Variation in the Interleukin-28B Gene Is Associated with Spontaneous Clearance and Progression of Hepatitis C Virus in Moroccan Patients Content: BACKGROUND Genetic variation in the IL28B gene has been strongly associated with treatment outcomes, spontaneous clearance and progression of the hepatitis C virus infection (HCV). The aim of the present study was to investigate the role of polymorphisms at this locus with progression and outcome of HCV infection in a Moroccan population. METHODS We analyzed a cohort of 438 individuals among them 232 patients with persistent HCV infection, of whom 115 patients had mild chronic hepatitis and 117 had advanced liver disease (cirrhosis and hepatocellular carcinoma), 68 individuals who had naturally cleared HCV and 138 healthy subjects. The IL28B SNPs rs12979860 and rs8099917 were genotyped using a TaqMan 5' allelic discrimination assay. RESULTS The protective rs12979860-C and rs8099917-T alleles were more common in subjects with spontaneous clearance (77.9% vs 55.2%; p = 0.00001 and 95.6% vs 83.2%; p = 0.0025, respectively). Individuals with clearance were 4.69 (95% CI, 1.99-11.07) times more likely to have the C/C genotype for rs12979860 polymorphism (p = 0.0017) and 3.55 (95% CI, 0.19-66.89) times more likely to have the T/T genotype at rs8099917. Patients with advanced liver disease carried the rs12979860-T/T genotype more frequently than patients with mild chronic hepatitis C (OR = 1.89; 95% CI, 0.99-3.61; p = 0.0532) and this risk was even more pronounced when we compared them with healthy controls (OR = 4.27; 95% CI, 2.08-8.76; p = 0.0005). The rs8099917-G allele was also associated with advanced liver disease (OR = 2.34; 95% CI, 1.40-3.93; p = 0.0100). CONCLUSIONS In the Moroccan population, polymorphisms near the IL28B gene play a role both in spontaneous clearance and progression of HCV infection.

Title: Large-scale candidate gene scan reveals the role of chemoreceptor genes in host plant specialization and speciation in the pea aphid. Content: Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Title: Mapping of a gene for severe pediatric gastroesophageal reflux to chromosome 13q14. Content: CONTEXT Gastroesophageal reflux (GER) has not previously been widely regarded as a hereditary disease. A few reports have suggested, however, that a genetic component may contribute to the incidence of GER, especially in its severe or chronic forms. OBJECTIVE To identify a genetic locus that cosegregates with a severe pediatric GER phenotype in families with multiple affected members. DESIGN A genome-wide scan of families affected by severe pediatric GER using polymorphic microsatellite markers spaced at an average of 8 centimorgans (cM), followed by haplotyping and by pairwise and multipoint linkage analyses. SETTING General US community, with research performed in a university tertiary care hospital. SUBJECTS Affected and unaffected family members from 5 families having multiple individuals affected by severe pediatric GER, identified through a patient support group. MAIN OUTCOME MEASURES Determination of inheritance patterns and linkage of a genetic locus with the severe pediatric GER phenotype by logarithm-of-odds (lod) score analysis, considering a lod score of 3 or greater as evidence of linkage. RESULTS In these families, severe pediatric GER followed an autosomal dominant hereditary pattern with high penetrance. A gene for severe pediatric GER was mapped to a 13-cM region on chromosome 13q between microsatellite markers D13S171 and D13S263. A maximum multifamily 2-point lod score of 5.58 and a maximum multifamily multipoint lod score of 7.15 were obtained for marker D13S1253 at map position 35 cM when presumptively affected persons were modeled as unknown (a maximum multipoint score of 4.88 was obtained when presumptively affected persons were modeled as unaffected). CONCLUSION These data suggest that a gene for severe pediatric GER maps to chromosome 13q14. JAMA. 2000;284:325-334

Title: Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly Content: Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.

Title: Fine-scale chromatin interaction maps reveal the cis-regulatory landscape of human lincRNA genes Content: High-throughput methods based on chromosome conformation capture have greatly advanced our understanding of the three-dimensional (3D) organization of genomes but are limited in resolution by their reliance on restriction enzymes. Here we describe a method called DNase Hi-C for comprehensively mapping global chromatin contacts. DNase Hi-C uses DNase I for chromatin fragmentation, leading to greatly improved efficiency and resolution over that of Hi-C. Coupling this method with DNA-capture technology provides a high-throughput approach for targeted mapping of fine-scale chromatin architecture. We applied targeted DNase Hi-C to characterize the 3D organization of 998 large intergenic noncoding RNA (lincRNA) promoters in two human cell lines. Our results revealed that expression of lincRNAs is tightly controlled by complex mechanisms involving both super-enhancers and the Polycomb repressive complex. Our results provide the first glimpse of the cell type-specific 3D organization of lincRNA genes.

Title: High-throughput parallel proteogenomics: a bacterial case study. Content: In recent years, a new paradigm for genome annotation has emerged, termed "proteogenomics," that leverages peptide MS to annotate a genome. This is achieved by mapping peptides to a six-frame translation of a genome, including available splice databases, which may suggest refinements to gene models. Using this approach, it is possible to refine gene regions such as exon boundaries, novel genes, gene boundaries, frame shifts, reverse strands, translated UTRs, and novel splice junctions. One of the challenges of proteogenomics is how best to (1) tackle assigning confidence to any resulting annotation and (2) apply these gene model refinements, either through manual annotation or through an automated process via training gene prediction tools. This is not a straightforward process, as many gene prediction tools have their defined suitability for niche genomes (either eukaryotic or prokaryotic) trained on and refined with model organisms such as Arabidopsis thaliana and Escherichia coli, and varying degrees of features that can leverage the use of external evidence. In this study, we outline a suitable approach toward preprocessing mass spectra and optimizing the MS/MS search for a given dataset. We also discuss future challenges, which continue to pose a problem in the field of proteogenomics, and better strategies to successfully tackle them with, using existing tools. We use Bradyrhizobium diazoefficiens (Nitrogen-fixing bacteria), with a 9.1 Mb genome as a case study, utilizing the latest in second-generation proteogenomics tools with multiple gene models for cross-validation of proteogenomics annotations.

Title: Smallest angiosperm genomes found in lentibulariaceae, with chromosomes of bacterial size. Content: Nuclear holoploid genome sizes (C-values) have been estimated to vary about 800-fold in angiosperms, with the smallest established 1C-value of 157 Mbp recorded in Arabidopsis thaliana. In the highly specialized carnivorous family Lentibulariaceae now three taxa have been found that exhibit significantly lower values: Genlisea margaretae with 63 Mbp, G. aurea with 64 Mbp, and Utricularia gibba with 88 Mbp. The smallest mitotic anaphase chromatids in G. aurea have 2.1 Mbp and are thus of bacterial size (NB: E. coli has ca. 4 Mbp). Several Utricularia species range somewhat lower than A. thaliana or are similar in genome size. The highest 1C-value known from species of Lentibulariaceae was found in Genlisea hispidula with 1510 Mbp, and results in about 24-fold variation for Genlisea and the Lentibulariaceae. Taking into account these new measurements, genome size variation in angiosperms is now almost 2000-fold. Genlisea and Utricularia are plants with terminal positions in the phylogeny of the eudicots, so that the findings are relevant for the understanding of genome miniaturization. Moreover, the Genlisea-Utricularia clade exhibits one of the highest mutational rates in several genomic regions in angiosperms, what may be linked to specialized patterns of genome evolution. Ultrasmall genomes have not been found in Pinguicula, which is the sister group of the Genlisea-Utricularia clade, and which does not show accelerated mutational rates. C-values in Pinguicula varied only 1.7-fold from 487 to 829 Mbp.

Title: Etiology of type 1 diabetes. Content: Recent genetic mapping and gene-phenotype studies have revealed the genetic architecture of type 1 diabetes. At least ten genes so far can be singled out as strong causal candidates. The known functions of these genes indicate the primary etiological pathways of this disease, including HLA class II and I molecules binding to preproinsulin peptides and T cell receptors, T and B cell activation, innate pathogen-viral responses, chemokine and cytokine signaling, and T regulatory and antigen-presenting cell functions. This review considers research in the field of type 1 diabetes toward identifying disease mechanisms using genetic approaches. The expression and functions of these pathways, and, therefore, disease susceptibility, will be influenced by epigenetic and environmental factors. Certain inherited immune phenotypes will be early precursors of type 1 diabetes and could be useful in future clinical trials.

Title: Multiplex genome engineering using CRISPR/Cas systems. Content: Functional elucidation of causal genetic variants and elements requires precise genome editing technologies. The type II prokaryotic CRISPR (clustered regularly interspaced short palindromic repeats)/Cas adaptive immune system has been shown to facilitate RNA-guided site-specific DNA cleavage. We engineered two different type II CRISPR/Cas systems and demonstrate that Cas9 nucleases can be directed by short RNAs to induce precise cleavage at endogenous genomic loci in human and mouse cells. Cas9 can also be converted into a nicking enzyme to facilitate homology-directed repair with minimal mutagenic activity. Lastly, multiple guide sequences can be encoded into a single CRISPR array to enable simultaneous editing of several sites within the mammalian genome, demonstrating easy programmability and wide applicability of the RNA-guided nuclease technology.

Title: Genome-wide association analyses in East Asians identify new susceptibility loci for colorectal cancer Content: To identify new genetic factors for colorectal cancer (CRC), we conducted a genome-wide association study in east Asians. By analyzing genome-wide data in 2,098 cases and 5,749 controls, we selected 64 promising SNPs for replication in an independent set of samples, including up to 5,358 cases and 5,922 controls. We identified four SNPs with association P values of 8.58 × 10(-7) to 3.77 × 10(-10) in the combined analysis of all east Asian samples. Three of the four were replicated in a study conducted in 26,060 individuals of European descent, with combined P values of 1.22 × 10(-10) for rs647161 (5q31.1), 6.64 × 10(-9) for rs2423279 (20p12.3) and 3.06 × 10(-8) for rs10774214 (12p13.32 near the CCND2 gene), derived from meta-analysis of data from both east Asian and European-ancestry populations. This study identified three new CRC susceptibility loci and provides additional insight into the genetics and biology of CRC.

Title: Model-based Analysis of ChIP-Seq (MACS) Content: We present Model-based Analysis of ChIP-Seq data, MACS, which analyzes data generated by short read sequencers such as Solexa's Genome Analyzer. MACS empirically models the shift size of ChIP-Seq tags, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome, allowing for more robust predictions. MACS compares favorably to existing ChIP-Seq peak-finding algorithms, and is freely available.

Title: A comprehensive genomic binding map of gene and chromatin regulatory proteins in Saccharomyces. Content: Hundreds of different proteins regulate and implement transcription in Saccharomyces. Yet their interrelationships have not been investigated on a comprehensive scale. Here we determined the genome-wide binding locations of 200 transcription-related proteins, under normal and acute heat-shock conditions. This study distinguishes binding between distal versus proximal promoter regions as well as the 3' ends of genes for nearly all mRNA and tRNA genes. This study reveals (1) a greater diversity and specialization of regulation associated with the SAGA transcription pathway compared to the TFIID pathway, (2) new regulators enriched at tRNA genes, (3) a global co-occupancy network of >20,000 unique regulator combinations that show a high degree of regulatory interconnections among lowly expressed genes, (4) regulators of the SAGA pathway located largely distal to the core promoter and regulators of the TFIID pathway located proximally, and (5) distinct mobilization of SAGA- versus TFIID-linked regulators during acute heat shock.

Title: Polycomb Associates Genome-wide with a Specific RNA Polymerase II Variant, and Regulates Metabolic Genes in ESCs Content: Polycomb repressor complexes (PRCs) are important chromatin modifiers fundamentally implicated in pluripotency and cancer. Polycomb silencing in embryonic stem cells (ESCs) can be accompanied by active chromatin and primed RNA polymerase II (RNAPII), but the relationship between PRCs and RNAPII remains unclear genome-wide. We mapped PRC repression markers and four RNAPII states in ESCs using ChIP-seq, and found that PRC targets exhibit a range of RNAPII variants. First, developmental PRC targets are bound by unproductive RNAPII (S5p(+)S7p(-)S2p(-)) genome-wide. Sequential ChIP, Ring1B depletion, and genome-wide correlations show that PRCs and RNAPII-S5p physically bind to the same chromatin and functionally synergize. Second, we identify a cohort of genes marked by PRC and elongating RNAPII (S5p(+)S7p(+)S2p(+)); they produce mRNA and protein, and their expression increases upon PRC1 knockdown. We show that this group of PRC targets switches between active and PRC-repressed states within the ESC population, and that many have roles in metabolism.

Title: Genome architecture, rearrangements and genomic disorders. Content: An increasing number of human diseases are recognized to result from recurrent DNA rearrangements involving unstable genomic regions. These are termed genomic disorders, in which the clinical phenotype is a consequence of abnormal dosage of gene(s) located within the rearranged genomic fragments. Both inter- and intrachromosomal rearrangements are facilitated by the presence of region-specific low-copy repeats (LCRs) and result from nonallelic homologous recombination (NAHR) between paralogous genomic segments. LCRs usually span approximately 10-400 kb of genomic DNA, share >or= 97% sequence identity, and provide the substrates for homologous recombination, thus predisposing the region to rearrangements. Moreover, it has been suggested that higher order genomic architecture involving LCRs plays a significant role in karyotypic evolution accompanying primate speciation.

Title: Sequence features surrounding the translation initiation sites assigned on the genome sequence of Synechocystis sp. strain PCC6803 by amino-terminal protein sequencing. Content: To characterize the sequence features surrounding the translation initiation sites on the genome of Synechocystis sp. strain 6803, the total proteins extracted from the cell were resolved by two-dimensional electrophoresis, and the amino-terminal sequences of the relatively abundant protein spots were determined. By comparison of the determined amino-terminal sequences with the nucleotide sequence of the entire genome, the translation initiation sites of a total of 72 proteins were successfully assigned on the genome. The sequence features emerged from the nucleotide sequences at and surrounding the translation initiation sites were as follows: (1) In addition to the three initiation codons, ATG, GTG, and TTG, evidence was obtained that ATT was also used as a rare initiation codon; (2) the core sequences (GAGG, GGAG and AGGA) of the Shine-Dalgarno sequence were identified in the appropriate position preceding the 35 initiation sites (48.6%); and (3) the preferential sequence surrounding the initiation codons was formulated as 5'-YY[...]R-3' where Y and R denote pyrimidine and purine nucleotides, respectively, and three dots represent the initiation codons. The result obtained would provide valuable information for improvement of the gene-finding software, and the approach used in this study should be applicable for comprehensive analysis of the expression profiles of cellular proteins.

Title: Genomewide Analysis of PRC1 and PRC2 Occupancy Identifies Two Classes of Bivalent Domains Content: In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes -- the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.

Title: Genome Sequencing and Analysis of the Tasmanian Devil and Its Transmissible Cancer Content: The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

Title: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Content: High-throughput mRNA sequencing (RNA-Seq) promises simultaneous transcript discovery and abundance estimation. However, this would require algorithms that are not restricted by prior gene annotations and that account for alternative transcription and splicing. Here we introduce such algorithms in an open-source software program called Cufflinks. To test Cufflinks, we sequenced and analyzed >430 million paired 75-bp RNA-Seq reads from a mouse myoblast cell line over a differentiation time series. We detected 13,692 known transcripts and 3,724 previously unannotated ones, 62% of which are supported by independent expression data or by homologous genes in other species. Over the time series, 330 genes showed complete switches in the dominant transcription start site (TSS) or splice isoform, and we observed more subtle shifts in 1,304 other genes. These results suggest that Cufflinks can illuminate the substantial regulatory flexibility and complexity in even this well-studied model of muscle development and that it can improve transcriptome-based genome annotation.

Title: Array-based Comparative Genomic Hybridization for Genome-Wide Screening of DNA Copy Number in Bladder Tumors Content: Genome-wide copy number profiles were characterized in 41 primary bladder tumors using array-based comparative genomic hybridization (array CGH). In addition to previously identified alterations in large chromosomal regions, alterations were identified in many small genomic regions, some with high-level amplifications or homozygous deletions. High-level amplifications were detected for 192 genomic clones, most frequently at 6p22.3 (E2F3), 8p12 (FGFR1), 8q22.2 (CMYC), 11q13 (CCND1, EMS1, INT2), and 19q13.1 (CCNE). Homozygous deletions were detected in 51 genomic clones, with four showing deletions in more than one case: two clones mapping to 9p21.3 (CDKN2A/p16, in nine cases), one at 8p23.1 (three cases), and one at 11p13 (two cases). Significant correlations were observed between copy number gain of clones containing CCNE1 and gain of ERBB2, and between gain of CCND1 and deletion of TP53. In addition, there was a significant complementary association between gain of CCND1 and gain of E2F3. Although there was no significant relationship between copy number changes and tumor stage or grade, the linked behavior among genomic loci suggests that array CGH will be increasingly important in understanding pathways critical to bladder tumor biology.

Title: Variant Histone H2A.Z Is Globally Localized to the Promoters of Inactive Yeast Genes and Regulates Nucleosome Positioning Content: H2A.Z is an evolutionary conserved histone variant involved in transcriptional regulation, antisilencing, silencing, and genome stability. The mechanism(s) by which H2A.Z regulates these various biological functions remains poorly defined, in part due to the lack of knowledge regarding its physical location along chromosomes and the bearing it has in regulating chromatin structure. Here we mapped H2A.Z across the yeast genome at an approximately 300-bp resolution, using chromatin immunoprecipitation combined with tiling microarrays. We have identified 4,862 small regions—typically one or two nucleosomes wide—decorated with H2A.Z. Those “Z loci” are predominantly found within specific nucleosomes in the promoter of inactive genes all across the genome. Furthermore, we have shown that H2A.Z can regulate nucleosome positioning at the GAL1 promoter. Within HZAD domains, the regions where H2A.Z shows an antisilencing function, H2A.Z is localized in a wider pattern, suggesting that the variant histone regulates a silencing and transcriptional activation via different mechanisms. Our data suggest that the incorporation of H2A.Z into specific promoter-bound nucleosomes configures chromatin structure to poise genes for transcriptional activation. The relevance of these findings to higher eukaryotes is discussed.

Title: Shared and distinct genetic variants in type 1 diabetes and celiac disease. Content: BACKGROUND Two inflammatory disorders, type 1 diabetes and celiac disease, cosegregate in populations, suggesting a common genetic origin. Since both diseases are associated with the HLA class II genes on chromosome 6p21, we tested whether non-HLA loci are shared. METHODS We evaluated the association between type 1 diabetes and eight loci related to the risk of celiac disease by genotyping and statistical analyses of DNA samples from 8064 patients with type 1 diabetes, 9339 control subjects, and 2828 families providing 3064 parent-child trios (consisting of an affected child and both biologic parents). We also investigated 18 loci associated with type 1 diabetes in 2560 patients with celiac disease and 9339 control subjects. RESULTS Three celiac disease loci--RGS1 on chromosome 1q31, IL18RAP on chromosome 2q12, and TAGAP on chromosome 6q25--were associated with type 1 diabetes (P<1.00x10(-4)). The 32-bp insertion-deletion variant on chromosome 3p21 was newly identified as a type 1 diabetes locus (P=1.81x10(-8)) and was also associated with celiac disease, along with PTPN2 on chromosome 18p11 and CTLA4 on chromosome 2q33, bringing the total number of loci with evidence of a shared association to seven, including SH2B3 on chromosome 12q24. The effects of the IL18RAP and TAGAP alleles confer protection in type 1 diabetes and susceptibility in celiac disease. Loci with distinct effects in the two diseases included INS on chromosome 11p15, IL2RA on chromosome 10p15, and PTPN22 on chromosome 1p13 in type 1 diabetes and IL12A on 3q25 and LPP on 3q28 in celiac disease. CONCLUSIONS A genetic susceptibility to both type 1 diabetes and celiac disease shares common alleles. These data suggest that common biologic mechanisms, such as autoimmunity-related tissue damage and intolerance to dietary antigens, may be etiologic features of both diseases.

Title: Choosing Haplotype-Tagging SNPS Based on Unphased Genotype Data Using a Preliminary Sample of Unrelated Subjects with an Example from the Multiethnic Cohort Study Content: We describe an approach for picking haplotype-tagging single nucleotide polymorphisms (htSNPs) that is presently being taken in two large nested case-control studies within a multiethnic cohort (MEC), which are engaged in a search for associations between risk of prostate and breast cancer and common genetic variations in candidate genes. Based on a preliminary sample of 70 control subjects chosen at random from each of the 5 ethnic groups in the MEC we estimate haplotype frequencies using a variant of the Excoffier-Slatkin E-M algorithm after genotyping a high density of SNPs selected every 3–5 kb in and surrounding a candidate gene. In order to evaluate the performance of a candidate set of htSNPS (which will be genotyped in the much larger case-control sample) we treat the haplotype frequencies estimate above as known, and carry out a formal calculation of the uncertainty of the number of copies of common haplotypes carried by an individual, summarizing this calculation as a coefficient of determination, R2h. A candidate set of htSNPS of a given size is chosen so as to maximize the minimum value of R2h over the common haplotypes, h.

Title: Genetic variants and non-genetic factors predict circulating vitamin D levels in Hispanic and non-Hispanic White women: the Breast Cancer Health Disparities Study. Content: Genome-wide association studies (GWAS) have identified common polymorphisms in or near GC, CYP2R1, CYP24A1, and NADSYN1/DHCR7 genes to be associated with circulating levels of 25-hydroxyvitamin D [25(OH)D] in European populations. To replicate these GWAS findings, we examined six selected polymorphisms from these regions and their relation with circulating 25(OH)D levels in 1,605 Hispanic women (629 U.S. Hispanics and 976 Mexicans) and 354 non-Hispanic White (NHW) women. We also assessed the potential interactions between these variants and known non-genetic predictors of 25(OH)D levels, including body mass index (BMI), sunlight exposure and vitamin D intake from diet and supplements. The minor alleles of the two GC polymorphisms (rs7041 and rs2282679) were significantly associated with lower 25(OH)D levels in both Hispanic and NHW women. The CYP2R1 polymorphism, rs2060793, also was significantly associated with 25(OH)D levels in both groups. We found no significant associations for the polymorphisms in the CYP24A1. In Hispanic controls, 25(OH)D levels were significantly associated with the rs12785878T and rs1790349G haplotype in the NADSYN1/DHCR7 region. Significant interactions between GC rs2282679 and BMI and between rs12785878 and time spent in outdoor activities were observed. These results provide further support for the contribution of common genetic variants to individual variability in circulating 25(OH)D levels. The observed interactions between SNPs and non-genetic factors warrant confirmation.

Title: CHROMOTHRIPSIS FROM DNA DAMAGE IN MICRONUCLEI Content: Genome sequencing has uncovered a new mutational phenomenon in cancer and congenital disorders called chromothripsis. Chromothripsis is characterized by extensive genomic rearrangements and an oscillating pattern of DNA copy number levels, all curiously restricted to one or a few chromosomes. The mechanism for chromothripsis is unknown, but we previously proposed that it could occur through the physical isolation of chromosomes in aberrant nuclear structures called micronuclei. Here, using a combination of live cell imaging and single-cell genome sequencing, we demonstrate that micronucleus formation can indeed generate a spectrum of genomic rearrangements, some of which recapitulate all known features of chromothripsis. These events are restricted to the mis-segregated chromosome and occur within one cell division. We demonstrate that the mechanism for chromothripsis can involve the fragmentation and subsequent reassembly of a single chromatid from a micronucleus. Collectively, these experiments establish a new mutational process of which chromothripsis is one extreme outcome.

Title: A common missense variant in NUDT15 confers susceptibility to thiopurine-induced leukopenia Content: Thiopurine therapy, commonly used in autoimmune conditions, can be complicated by life-threatening leukopenia. This leukopenia is associated with genetic variation in TPMT (encoding thiopurine S-methyltransferase). Despite a lower frequency of TPMT mutations in Asians, the incidence of thiopurine-induced leukopenia is higher in Asians than in individuals of European descent. Here we performed an Immunochip-based 2-stage association study in 978 Korean subjects with Crohn's disease treated with thiopurines. We identified a nonsynonymous SNP in NUDT15 (encoding p. Arg139Cys) that was strongly associated with thiopurine-induced early leukopenia (odds ratio (OR) = 35.6; Pcombined = 4.88 × 10−94). In Koreans, this variant demonstrated sensitivity and specificity of 89.4% and 93.2%, respectively, for thiopurine-induced early leukopenia (in comparison to 12.1% and 97.6% for TPMT variants). Although rare, this SNP was also strongly associated with thiopurine-induced leukopenia in subjects with inflammatory bowel disease of European descent (OR = 9.50; P = 4.64 × 10−4). Thus, NUDT15 is a pharmacogenetic determinant for thiopurine-induced leukopenia in diverse populations.

Title: Update of the FANTOM web resource: high resolution transcriptome of diverse cell types in mammals Content: Upon the first publication of the fifth iteration of the Functional Annotation of Mammalian Genomes collaborative project, FANTOM5, we gathered a series of primary data and database systems into the FANTOM web resource (http://fantom.gsc.riken.jp) to facilitate researchers to explore transcriptional regulation and cellular states. In the course of the collaboration, primary data and analysis results have been expanded, and functionalities of the database systems enhanced. We believe that our data and web systems are invaluable resources, and we think the scientific community will benefit for this recent update to deepen their understanding of mammalian cellular organization. We introduce the contents of FANTOM5 here, report recent updates in the web resource and provide future perspectives.

Title: Genetic and Functional Diversification of Small RNA Pathways in Plants Content: Multicellular eukaryotes produce small RNA molecules (approximately 21–24 nucleotides) of two general types, microRNA (miRNA) and short interfering RNA (siRNA). They collectively function as sequence-specific guides to silence or regulate genes, transposons, and viruses and to modify chromatin and genome structure. Formation or activity of small RNAs requires factors belonging to gene families that encode DICER (or DICER-LIKE [DCL]) and ARGONAUTE proteins and, in the case of some siRNAs, RNA-dependent RNA polymerase (RDR) proteins. Unlike many animals, plants encode multiple DCL and RDR proteins. Using a series of insertion mutants of Arabidopsis thaliana, unique functions for three DCL proteins in miRNA (DCL1), endogenous siRNA (DCL3), and viral siRNA (DCL2) biogenesis were identified. One RDR protein (RDR2) was required for all endogenous siRNAs analyzed. The loss of endogenous siRNA in dcl3 and rdr2 mutants was associated with loss of heterochromatic marks and increased transcript accumulation at some loci. Defects in siRNA-generation activity in response to turnip crinkle virus in dcl2 mutant plants correlated with increased virus susceptibility. We conclude that proliferation and diversification of DCL and RDR genes during evolution of plants contributed to specialization of small RNA-directed pathways for development, chromatin structure, and defense.

Title: LDL-cholesterol concentrations: a genome-wide association study Content: BACKGROUND LDL cholesterol has a causal role in the development of cardiovascular disease. Improved understanding of the biological mechanisms that underlie the metabolism and regulation of LDL cholesterol might help to identify novel therapeutic targets. We therefore did a genome-wide association study of LDL-cholesterol concentrations. METHODS We used genome-wide association data from up to 11,685 participants with measures of circulating LDL-cholesterol concentrations across five studies, including data for 293 461 autosomal single nucleotide polymorphisms (SNPs) with a minor allele frequency of 5% or more that passed our quality control criteria. We also used data from a second genome-wide array in up to 4337 participants from three of these five studies, with data for 290,140 SNPs. We did replication studies in two independent populations consisting of up to 4979 participants. Statistical approaches, including meta-analysis and linkage disequilibrium plots, were used to refine association signals; we analysed pooled data from all seven populations to determine the effect of each SNP on variations in circulating LDL-cholesterol concentrations. FINDINGS In our initial scan, we found two SNPs (rs599839 [p=1.7x10(-15)] and rs4970834 [p=3.0x10(-11)]) that showed genome-wide statistical association with LDL cholesterol at chromosomal locus 1p13.3. The second genome screen found a third statistically associated SNP at the same locus (rs646776 [p=4.3x10(-9)]). Meta-analysis of data from all studies showed an association of SNPs rs599839 (combined p=1.2x10(-33)) and rs646776 (p=4.8x10(-20)) with LDL-cholesterol concentrations. SNPs rs599839 and rs646776 both explained around 1% of the variation in circulating LDL-cholesterol concentrations and were associated with about 15% of an SD change in LDL cholesterol per allele, assuming an SD of 1 mmol/L. INTERPRETATION We found evidence for a novel locus for LDL cholesterol on chromosome 1p13.3. These results potentially provide insight into the biological mechanisms that underlie the regulation of LDL cholesterol and might help in the discovery of novel therapeutic targets for cardiovascular disease.

Title: Effective population size and patterns of molecular evolution and variation Content: The effective size of a population, Ne, determines the rate of change in the composition of a population caused by genetic drift, which is the random sampling of genetic variants in a finite population. Ne is crucial in determining the level of variability in a population, and the effectiveness of selection relative to drift. This article reviews the properties of Ne in a variety of different situations of biological interest, and the factors that influence it. In particular, the action of selection means that Ne varies across the genome, and advances in genomic techniques are giving new insights into how selection shapes Ne.

Title: Tolerance of whole-genome doubling propagates chromosomal instability and accelerates cancer genome evolution. Content: UNLABELLED The contribution of whole-genome doubling to chromosomal instability (CIN) and tumor evolution is unclear. We use long-term culture of isogenic tetraploid cells from a stable diploid colon cancer progenitor to investigate how a genome-doubling event affects genome stability over time. Rare cells that survive genome doubling demonstrate increased tolerance to chromosome aberrations. Tetraploid cells do not exhibit increased frequencies of structural or numerical CIN per chromosome. However, the tolerant phenotype in tetraploid cells, coupled with a doubling of chromosome aberrations per cell, allows chromosome abnormalities to evolve specifically in tetraploids, recapitulating chromosomal changes in genomically complex colorectal tumors. Finally, a genome-doubling event is independently predictive of poor relapse-free survival in early-stage disease in two independent cohorts in multivariate analyses [discovery data: hazard ratio (HR), 4.70, 95% confidence interval (CI), 1.04-21.37; validation data: HR, 1.59, 95% CI, 1.05-2.42]. These data highlight an important role for the tolerance of genome doubling in driving cancer genome evolution. SIGNIFICANCE Our work sheds light on the importance of whole-genome–doubling events in colorectal cancer evolution. We show that tetraploid cells undergo rapid genomic changes and recapitulate the genetic alterations seen in chromosomally unstable tumors. Furthermore, we demonstrate that a genome-doubling event is prognostic of poor relapse-free survival in this disease type.

Title: The Mammalian Epigenome Content: Chemical modifications to DNA and histone proteins form a complex regulatory network that modulates chromatin structure and genome function. The epigenome refers to the complete description of these potentially heritable changes across the genome. The composition of the epigenome within a given cell is a function of genetic determinants, lineage, and environment. With the sequencing of the human genome completed, investigators now seek a comprehensive view of the epigenetic changes that determine how genetic information is made manifest across an incredibly varied background of developmental stages, tissue types, and disease states. Here we review current research efforts, with an emphasis on large-scale studies, emerging technologies, and challenges ahead.

Title: EIF2AK3, encoding translation initiation factor 2-α kinase 3, is mutated in patients with Wolcott-Rallison syndrome Content: Wolcott-Rallison syndrome (WRS) is a rare, autosomal recessive disorder characterized by permanent neonatal or early infancy insulin-dependent diabetes. Epiphyseal dysplasia, osteoporosis and growth retardation occur at a later age. Other frequent multisystemic manifestations include hepatic and renal dysfunction, mental retardation and cardiovascular abnormalities. On the basis of two consanguineous families, we mapped WRS to a region of less than 3 cM on chromosome 2p12, with maximal evidence of linkage and homozygosity at 4 microsatellite markers within an interval of approximately 1 cM. The gene encoding the eukaryotic translation initiation factor 2-α kinase 3 (EIF2AK3) resides in this interval; thus we explored it as a candidate. We identified distinct mutations of EIF2AK3 that segregated with the disorder in each of the families. The first mutation produces a truncated protein in which the entire catalytic domain is missing. The other changes an amino acid, located in the catalytic domain of the protein, that is highly conserved among kinases from the same subfamily. Our results provide evidence for the role of EIF2AK3 in WRS. The identification of this gene may provide insight into the understanding of the more common forms of diabetes and other pathologic manifestations of WRS.

Title: Exome sequencing identifies somatic gain-of-function PPM1D mutations in brainstem gliomas Content: Gliomas arising in the brainstem and thalamus are devastating tumors that are difficult to surgically resect. To determine the genetic and epigenetic landscape of these tumors, we performed exomic sequencing of 14 brainstem gliomas (BSGs) and 12 thalamic gliomas. We also performed targeted mutational analysis of an additional 24 such tumors and genome-wide methylation profiling of 45 gliomas. This study led to the discovery of tumor-specific mutations in PPM1D, encoding wild-type p53-induced protein phosphatase 1D (WIP1), in 37.5% of the BSGs that harbored hallmark H3F3A mutations encoding p. Lys27Met substitutions. PPM1D mutations were mutually exclusive with TP53 mutations in BSG and attenuated p53 activation in vitro. PPM1D mutations were truncating alterations in exon 6 that enhanced the ability of PPM1D to suppress the activation of the DNA damage response checkpoint protein CHK2. These results define PPM1D as a frequent target of somatic mutation and as a potential therapeutic target in brainstem gliomas.

Title: Sperm cross-over activity in regions of the human genome showing extreme breakdown of marker association. Content: Population diversity data have recently provided profound, albeit inferential, insights into meiotic recombination across the human genome, revealing a landscape dominated by thousands of cross-over hotspots. However, very few of these putative hotspots have been directly analyzed for cross-over activity. We now describe a search for very active hotspots, by using extreme breakdown of marker association as a guide for high-resolution sperm cross-over analysis. This strategy has led to the isolation of the most active cross-over hotspots yet described. Their morphology, sequence attributes, and cross-over processes are very similar to those seen at less active hotspots, but their activity in sperm is poorly predicted from population diversity information. Several of these hotspots showed evidence for biased gene conversion accompanying cross-over, in some cases associated with variation between men in cross-over activity and with two hotspots showing complete presence/absence polymorphism in different men. Hotspot polymorphism is very common at less active hotspots but curiously was not seen at any of the most active hotspots. This contrasts with the prediction that extreme hotspots should be the most vulnerable to attenuation by meiotic drive in favor of mutations that suppress recombination and should therefore show rapid rate evolution and thus variation in activity between men. Finally, these very intense hotspots provide a valuable resource for dissecting meiotic recombination processes and pathways in humans.

Title: The miniature genome of a carnivorous plant Genlisea aurea contains a low number of genes and short non-coding sequences Content: BACKGROUND Genlisea aurea (Lentibulariaceae) is a carnivorous plant with unusually small genome size - 63.6 Mb - one of the smallest known among higher plants. Data on the genome sizes and the phylogeny of Genlisea suggest that this is a derived state within the genus. Thus, G. aurea is an excellent model organism for studying evolutionary mechanisms of genome contraction. RESULTS Here we report sequencing and de novo draft assembly of G. aurea genome. The assembly consists of 10,687 contigs of the total length of 43.4 Mb and includes 17,755 complete and partial protein-coding genes. Its comparison with the genome of Mimulus guttatus, another representative of higher core Lamiales clade, reveals striking differences in gene content and length of non-coding regions. CONCLUSIONS Genome contraction was a complex process, which involved gene loss and reduction of lengths of introns and intergenic regions, but not intron loss. The gene loss is more frequent for the genes that belong to multigenic families indicating that genetic redundancy is an important prerequisite for genome size reduction.

Title: A genome-wide siRNA screen reveals diverse cellular processes and pathways that mediate genome stability. Content: Signaling pathways that respond to DNA damage are essential for the maintenance of genome stability and are linked to many diseases, including cancer. Here, a genome-wide siRNA screen was employed to identify additional genes involved in genome stabilization by monitoring phosphorylation of the histone variant H2AX, an early mark of DNA damage. We identified hundreds of genes whose downregulation led to elevated levels of H2AX phosphorylation (gammaH2AX) and revealed links to cellular complexes and to genes with unclassified functions. We demonstrate a widespread role for mRNA-processing factors in preventing DNA damage, which in some cases is caused by aberrant RNA-DNA structures. Furthermore, we connect increased gammaH2AX levels to the neurological disorder Charcot-Marie-Tooth (CMT) syndrome, and we find a role for several CMT proteins in the DNA-damage response. These data indicate that preservation of genome stability is mediated by a larger network of biological processes than previously appreciated.

Title: Quantitative analysis of tumor-derived methylated p16INK4a sequences in plasma, serum, and blood cells of hepatocellular carcinoma patients. Content: PURPOSE AND EXPERIMENTAL DESIGN Using real-time quantitative methylation-specific PCR (RTQ-MSP), we quantified methylated p16INK4a sequences and determined the fractional concentrations of circulating tumor DNA in plasma, serum, and peripheral blood cells collected preoperatively, intraoperatively, and postoperatively from 49 patients with hepatocellular carcinoma (HCC). RESULTS RTQ-MSP was sufficiently sensitive to detect down to 10 genome-equivalents of methylated p16INK4a sequences. Quantitative MSP data were expressed in terms of the methylation index, which was the percentage of bisulfite converted unmethylated and methylated p16INK4a sequences that consisted of methylated p16INK4a sequences. Quantities of methylated p16INK4a sequences were detected in peripheral circulation of 80% (23 of 29) of HCC patients. No significant difference was seen in the detectability and concentrations of methylated p16INK4a sequences (range: 10-4046 genome-equivalents/ml) between preoperative plasma and serum samples from HCC patients. Preoperatively, the p16INK4a methylation indices ranged from 0.2 to 100% and from 0.012 to 0.075% in the patients' plasma and buffy coat samples, respectively. After surgical resection, the median p16INK4a methylation indices in plasma and buffy coat concordantly decreased 12- and 15-fold, respectively. These results demonstrated the clinical usefulness and effectiveness of peripheral blood RTQ-MSP for detecting and monitoring HCC after treatment. Furthermore, none of the intraoperative plasma samples and only two of the intraoperative buffy coat samples were p16INK4a methylation positive. CONCLUSIONS Quantification of epigenetic changes in peripheral blood by RTQ-MSP is useful for the detection and monitoring of HCC.

Title: Large scale replication and meta-analysis of variants on chromosome 4q25 associated with atrial fibrillation. Content: AIMS A recent genome-wide association study identified a haplotype block on chromosome 4q25 associated with atrial fibrillation (AF). We sought to replicate this association in four independent cohorts. METHODS AND RESULTS The Framingham Heart Study and Rotterdam Study are community-based longitudinal studies. The Vanderbilt AF Registry and German AF Network (AFNet) are case-control studies. Participants with AF (n = 3508) were more likely to be male and were older than referent participants (n = 12 173; Framingham 82 +/- 10 vs. 71 +/- 13 years; Rotterdam 73 +/- 8 vs. 69 +/- 9 years; Vanderbilt 54 +/- 14 vs. 57 +/- 14 years; AFNet 62 +/- 12 vs. 49 +/- 14 years). Single nucleotide polymorphism (SNP) rs2200733 was associated with AF in all four cohorts, with odds ratios (ORs) ranging from 1.37 in Rotterdam [95% confidence interval (CI) 1.18-1.59; P = 3.1 x 10(-5)] to 2.52 in AFNet (95% CI 2.22-2.8; P = 1.8 x 10(-49)). There also was a significant association between AF and rs10033464 in Framingham (OR 1.34; 95% CI 1.03-1.75; P = 0.031) and AFNet (OR 1.30; 95% CI 1.13-1.51; P = 0.0002), but not Vanderbilt (OR 1.16; 95% CI 0.86-1.56; P = 0.33). A meta-analysis of the current and prior AF studies revealed an OR of 1.90 (95% CI 1.60-2.26; P = 3.3 x 10(-13)) for rs2200733 and of 1.36 (95% CI 1.26-1.47; P = 6.7 x 10(-15)) for rs10033464. CONCLUSION The non-coding SNPs rs2200733 and rs10033464 are strongly associated with AF in four cohorts of European descent. These results confirm the significant relations between AF and intergenic variants on chromosome 4.

Title: The DNA Methylome of Human Peripheral Blood Mononuclear Cells Content: DNA methylation plays an important role in biological processes in human health and disease. Recent technological advances allow unbiased whole-genome DNA methylation (methylome) analysis to be carried out on human cells. Using whole-genome bisulfite sequencing at 24.7-fold coverage (12.3-fold per strand), we report a comprehensive (92.62%) methylome and analysis of the unique sequences in human peripheral blood mononuclear cells (PBMC) from the same Asian individual whose genome was deciphered in the YH project. PBMC constitute an important source for clinical blood tests world-wide. We found that 68.4% of CpG sites and <0.2% of non-CpG sites were methylated, demonstrating that non-CpG cytosine methylation is minor in human PBMC. Analysis of the PBMC methylome revealed a rich epigenomic landscape for 20 distinct genomic features, including regulatory, protein-coding, non-coding, RNA-coding, and repeat sequences. Integration of our methylome data with the YH genome sequence enabled a first comprehensive assessment of allele-specific methylation (ASM) between the two haploid methylomes of any individual and allowed the identification of 599 haploid differentially methylated regions (hDMRs) covering 287 genes. Of these, 76 genes had hDMRs within 2 kb of their transcriptional start sites of which >80% displayed allele-specific expression (ASE). These data demonstrate that ASM is a recurrent phenomenon and is highly correlated with ASE in human PBMCs. Together with recently reported similar studies, our study provides a comprehensive resource for future epigenomic research and confirms new sequencing technology as a paradigm for large-scale epigenomics studies.

Title: Genome-wide dFOXO targets and topology of the transcriptomic response to stress and insulin signalling Content: FoxO transcription factors, inhibited by insulin/insulin-like growth factor signalling (IIS), are crucial players in numerous organismal processes including lifespan. Using genomic tools, we uncover over 700 direct dFOXO targets in adult female Drosophila. dFOXO is directly required for transcription of several IIS components and interacting pathways, such as TOR, in the wild-type fly. The genomic locations occupied by dFOXO in adults are different from those observed in larvae or cultured cells. These locations remain unchanged upon activation by stresses or reduced IIS, but the binding is increased and additional targets activated upon genetic reduction in IIS. We identify the part of the IIS transcriptional response directly controlled by dFOXO and the indirect effects and show that parts of the transcriptional response to IIS reduction do not require dfoxo. Promoter analyses revealed GATA and other forkhead factors as candidate mediators of the indirect and dfoxo-independent effects. We demonstrate genome-wide evolutionary conservation of dFOXO targets between the fly and the worm Caenorhabditis elegans, enriched for a second tier of regulators including the dHR96/daf-12 nuclear hormone receptor.

Title: Two classes of endogenous small RNAs in Tetrahymena thermophila. Content: Endogenous small RNAs function in RNA interference (RNAi) pathways to guide RNA cleavage, translational repression, or methylation of DNA or chromatin. In Tetrahymena thermophila, developmentally regulated DNA elimination is governed by an RNAi mechanism involving approximately 27-30-nucleotide (nt) RNAs. Here we characterize the sequence features of the approximately 27-30-nt RNAs and a approximately 23-24-nt RNA class representing a second RNAi pathway. The approximately 23-24-nt RNAs accumulate strain-specifically manner and map to the genome in clusters that are antisense to predicted genes. These findings reveal the existence of distinct endogenous RNAi pathways in the unicellular T. thermophila, a complexity previously demonstrated only in multicellular organisms.

Title: Genetic variants rs1994016 and rs3825807 in ADAMTS7 affect its mRNA expression in atherosclerotic occlusive peripheral arterial disease Content: AIM Peripheral artery disease (PAD) is a vascular disease affecting peripheral circulation. Recently, genome-wide association studies revealed a relationship between single nucleotide polymorphisms (SNPs) in ADAMTS7 (a disintegrin and metalloprotease with thrombospondin motif 7) and atherosclerosis. In this study, we aimed to determine ADAMTS7 expression in peripheral blood mononuclear cells (PBMCs) and the frequency of ADAMTS7 rs1994016 and rs3825807 polymorphisms in a sample of Turkish patients with PAD, and to evaluate the association of matrix metalloproteinase (MMP) levels with PAD development. METHODS In this case-control study, ADAMTS7mRNA and protein expression was determined using reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) and western blot, respectively, and rs1994016 and rs3825807 variants in ADAMTS7 were determined by real-time PCR in 115 PAD patients and 116 healthy controls. Plasma levels of nine MMPs were determined using a multiplex immunoassay system. RESULTS ADAMTS7mRNA levels were significantly higher in PAD patients than in controls (t=-2.75, P=.007). There was no significant difference in the frequencies of rs1994016 and rs3825807 between PAD patients and controls (P>.05). In PAD patients, ADAMTS7mRNA levels were significantly increased for the CC genotype of rs1994016 (t=-2.31, P=.026) and TT genotype of rs3825807 (t=-2.23, P=.032). Furthermore, plasma levels of MMP-1, MMP-3, MMP-7, MMP-10, MMP-12, and MMP-13 were significantly higher in PAD patients than in controls (P<.05). CONCLUSION This is the first report of the relationship between PAD and ADAMTS7 expression and the effects of the rs1994016 and rs3825807 variants on PAD development. ADAMTS7 may be associated with PAD development.

Title: Familial Barrett's esophagus associated with adenocarcinoma. Content: A family with six cases of Barrett's esophagus over three generations is described. The Barrett's esophagus affected only males, and there were three associated adenocarcinomas. The mechanism of inheritance is compatible with an autosomal dominant pattern with nearly complete penetrance. Other case reports in the literature are compared with this familial case of Barrett's esophagus with associated adenocarcinomas.

Title: A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4. Content: IMPORTANCE Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain. OBJECTIVE To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture. DESIGN, SETTING, AND PATIENTS In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain. MAIN OUTCOMES AND MEASURES The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples. RESULTS During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered. CONCLUSIONS AND RELEVANCE These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

Title: Functional genomics and proteomics approaches to study the ERBB network in cancer. Content: Substantial progress in functional genomic and proteomic technologies has opened new perspectives in biomedical research. The sequence of the human genome has been mostly determined and opened new visions on its complexity and regulation. New technologies, like RNAi and protein arrays, allow gathering knowledge beyond single gene analysis. Increasingly, biological processes are studied with systems biological approaches, where qualitative and quantitative data of the components are utilized to model the respective processes, to predict effects of perturbations, and to then refine these models after experimental testing. Here, we describe the potential of applying functional genomics and proteomics, taking the ERBB family of growth-factor receptors as an example to study the signaling network and its impact on cancer.