RNA-RNA Interactions Enable Specific Targeting of Noncoding RNAs to Nascent Pre-mRNAs and Chromatin Sites

Intermolecular RNA-RNA interactions are used by many noncoding RNAs (ncRNAs) to achieve their diverse functions. To identify these contacts, we developed a method based on RNA antisense purification to systematically map RNA-RNA interactions (RAP-RNA) and applied it to investigate two ncRNAs implicated in RNA processing: U1 small nuclear RNA, a component of the spliceosome, and Malat1, a large ncRNA that localizes to nuclear speckles. U1 and Malat1 interact with nascent transcripts through distinct targeting mechanisms. Using differential crosslinking, we confirmed that U1 directly hybridizes to 5′ splice sites and 5′ splice site motifs throughout introns and found that Malat1 interacts with pre-mRNAs indirectly through protein intermediates. Interactions with nascent pre-mRNAs cause U1 and Malat1 to localize proximally to chromatin at active genes, demonstrating that ncRNAs can use RNA-RNA interactions to target specific pre-mRNAs and genomic sites. RAP-RNA is sensitive to lower abundance RNAs as well, making it generally applicable for investigating ncRNAs

Mining data from 1000 genomes to identify the causal variant in regions under positive selection

The human genome contains hundreds of regions in which the patterns of genetic variation indicate recent positive natural selection, yet for most of these the underlying gene and the advantageous mutation remain unknown. We recently reported the development of a method, Composite of Multiple Signals (CMS), that combines tests for multiple signals of natural selection and increases resolution by up to 100-fold

Longitudinal assessment of neuronal 3D genomes in mouse prefrontal cortex

Neuronal epigenomes, including chromosomal loopings moving distal cis-regulatory elements into proximity of target genes, could serve as molecular proxy linking present-day-behaviour to past exposures. However, longitudinal assessment of chromatin state is challenging, because conventional chromosome conformation capture assays essentially provide single snapshots at a given time point, thus reflecting genome organization at the time of brain harvest and therefore are non-informative about the past. Here we introduce ‘NeuroDam’ to assess epigenome status retrospectively. Short-term expression of the bacterial DNA adenine methyltransferase Dam, tethered to the Gad1 gene promoter in mouse prefrontal cortex neurons, results in stable G[superscriptmethyl]ATC tags at Gad1-bound chromosomal contacts. We show by NeuroDam that mice with defective cognition 4 months after pharmacological NMDA receptor blockade already were affected by disrupted chromosomal conformations shortly after drug exposure. Retrospective profiling of neuronal epigenomes is likely to illuminate epigenetic determinants of normal and diseased brain development in longitudinal context.United States. National Institutes of Healt

Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation

As researchers begin probing deep coverage sequencing data for increasingly rare mutations and subclonal events, the fidelity of next generation sequencing (NGS) laboratory methods will become increasingly critical. Although error rates for sequencing and polymerase chain reaction (PCR) are well documented, the effects that DNA extraction and other library preparation steps could have on downstream sequence integrity have not been thoroughly evaluated. Here, we describe the discovery of novel C > A/G > T transversion artifacts found at low allelic fractions in targeted capture data. Characteristics such as sequencer read orientation and presence in both tumor and normal samples strongly indicated a non-biological mechanism. We identified the source as oxidation of DNA during acoustic shearing in samples containing reactive contaminants from the extraction process. We show generation of 8-oxoguanine (8-oxoG) lesions during DNA shearing, present analysis tools to detect oxidation in sequencing data and suggest methods to reduce DNA oxidation through the introduction of antioxidants. Further, informatics methods are presented to confidently filter these artifacts from sequencing data sets. Though only seen in a low percentage of reads in affected samples, such artifacts could have profoundly deleterious effects on the ability to confidently call rare mutations, and eliminating other possible sources of artifacts should become a priority for the research community.National Human Genome Research Institute (U.S.) (HG03067-05

The state of the Martian climate

60°N was +2.0°C, relative to the 1981–2010 average value (Fig. 5.1). This marks a new high for the record. The average annual surface air temperature (SAT) anomaly for 2016 for land stations north of starting in 1900, and is a significant increase over the previous highest value of +1.2°C, which was observed in 2007, 2011, and 2015. Average global annual temperatures also showed record values in 2015 and 2016. Currently, the Arctic is warming at more than twice the rate of lower latitudes

Identification and Functional Validation of the Novel Antimalarial Resistance Locus PF10_0355 in Plasmodium falciparum

The Plasmodium falciparum parasite's ability to adapt to environmental pressures, such as the human immune system and antimalarial drugs, makes malaria an enduring burden to public health. Understanding the genetic basis of these adaptations is critical to intervening successfully against malaria. To that end, we created a high-density genotyping array that assays over 17,000 single nucleotide polymorphisms (~1 SNP/kb), and applied it to 57 culture-adapted parasites from three continents. We characterized genome-wide genetic diversity within and between populations and identified numerous loci with signals of natural selection, suggesting their role in recent adaptation. In addition, we performed a genome-wide association study (GWAS), searching for loci correlated with resistance to thirteen antimalarials; we detected both known and novel resistance loci, including a new halofantrine resistance locus, PF10_0355. Through functional testing we demonstrated that PF10_0355 overexpression decreases sensitivity to halofantrine, mefloquine, and lumefantrine, but not to structurally unrelated antimalarials, and that increased gene copy number mediates resistance. Our GWAS and follow-on functional validation demonstrate the potential of genome-wide studies to elucidate functionally important loci in the malaria parasite genome.Bill & Melinda Gates FoundationEllison Medical FoundationExxon Mobil FoundationFogarty International CenterNational Institute of Allergy and Infectious Diseases (U.S.)Burroughs Wellcome FundDavid & Lucile Packard FoundationNational Science Foundation (U.S.). Graduate Research Fellowship Progra

Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel

A major use of the 1000 Genomes Project (1000GP) data is genotype imputation in genome-wide association studies (GWAS). Here we develop a method to estimate haplotypes from low-coverage sequencing data that can take advantage of single-nucleotide polymorphism (SNP) microarray genotypes on the same samples. First the SNP array data are phased to build a backbone (or \u27scaffold\u27) of haplotypes across each chromosome. We then phase the sequence data \u27onto\u27 this haplotype scaffold. This approach can take advantage of relatedness between sequenced and non-sequenced samples to improve accuracy. We use this method to create a new 1000GP haplotype reference set for use by the human genetic community. Using a set of validation genotypes at SNP and bi-allelic indels we show that these haplotypes have lower genotype discordance and improved imputation performance into downstream GWAS samples, especially at low-frequency variants. © 2014 Macmillan Publishers Limited. All rights reserved

Investigating the complex genetic architecture of ankle-brachial index, a measure of peripheral arterial disease, in non-Hispanic whites

<p>Abstract</p> <p>Background</p> <p>Atherosclerotic peripheral arterial disease (PAD) affects 8–10 million people in the United States and is associated with a marked impairment in quality of life and an increased risk of cardiovascular events. Noninvasive assessment of PAD is performed by measuring the ankle-brachial index (ABI). Complex traits, such as ABI, are influenced by a large array of genetic and environmental factors and their interactions. We attempted to characterize the genetic architecture of ABI by examining the main and interactive effects of individual single nucleotide polymorphisms (SNPs) and conventional risk factors.</p> <p>Methods</p> <p>We applied linear regression analysis to investigate the association of 435 SNPs in 112 positional and biological candidate genes with ABI and related physiological and biochemical traits in 1046 non-Hispanic white, hypertensive participants from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. The main effects of each SNP, as well as SNP-covariate and SNP-SNP interactions, were assessed to investigate how they contribute to the inter-individual variation in ABI. Multivariable linear regression models were then used to assess the joint contributions of the top SNP associations and interactions to ABI after adjustment for covariates. We reduced the chance of false positives by 1) correcting for multiple testing using the false discovery rate, 2) internal replication, and 3) four-fold cross-validation.</p> <p>Results</p> <p>When the results from these three procedures were combined, only two SNP main effects in <it>NOS3</it>, three SNP-covariate interactions (<it>ADRB2 </it>Gly 16 – lipoprotein(a) and <it>SLC4A5 </it>– diabetes interactions), and 25 SNP-SNP interactions (involving SNPs from 29 different genes) were significant, replicated, and cross-validated. Combining the top SNPs, risk factors, and their interactions into a model explained nearly 18% of variation in ABI in the sample. SNPs in six genes (<it>ADD2, ATP6V1B1, PRKAR2B, SLC17A2, SLC22A3, and TGFB3</it>) were also influencing triglycerides, C-reactive protein, homocysteine, and lipoprotein(a) levels.</p> <p>Conclusion</p> <p>We found that candidate gene SNP main effects, SNP-covariate and SNP-SNP interactions contribute to the inter-individual variation in ABI, a marker of PAD. Our findings underscore the importance of conducting systematic investigations that consider context-dependent frameworks for developing a deeper understanding of the multidimensional genetic and environmental factors that contribute to complex diseases.</p

The genetic landscape of high-risk neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%1. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 cases using a combination of whole exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per megabase (0.48 non-silent), and remarkably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, an additional 7.1% had focal deletions), MYCN (1.7%, a recurrent p.Pro44Leu alteration), and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1, and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies reliant upon frequently altered oncogenic drivers

Development of a High-Density Linkage Map and Tagging Leaf Spot Resistance in Pearl Millet Using Genotyping-by-Sequencing Markers

Pearl millet [Pennisetum glaucum (L.) R. Br; also Cenchrus americanus (L.) Morrone] is an important crop throughout the world but better genomic resources for this species are needed to facilitate crop improvement. Genome mapping studies are a prerequisite for tagging agronomically important traits. Genotyping-by-sequencing (GBS) markers can be used to build high-density linkage maps, even in species lacking a reference genome. A recombinant inbred line (RIL) mapping population was developed from a cross between the lines ‘Tift 99D2B1’ and ‘Tift 454’. DNA from 186 RILs, the parents, and the F1 was used for 96-plex ApeKI GBS library development, which was further used for sequencing. The sequencing results showed that the average number of good reads per individual was 2.2 million, the pass filter rate was 88%, and the CV was 43%. High-quality GBS markers were developed with stringent filtering on sequence data from 179 RILs. The reference genetic map developed using 150 RILs contained 16,650 single-nucleotide polymorphisms (SNPs) and 333,567 sequence tags spread across all seven chromosomes. The overall average density of SNP markers was 23.23 SNP/cM in the final map and 1.66 unique linkage bins per cM covering a total genetic distance of 716.7 cM. The linkage map was further validated for its utility by using it in mapping quantitative trait loci (QTLs) for flowering time and resistance to Pyricularia leaf spot [Pyricularia grisea (Cke.) Sacc.]. This map is the densest yet reported for this crop and will be a valuable resource for the pearl millet community