Lira Mamanova

Target-enrichment strategies for next-generation sequencing

We have not yet reached a point at which routine sequencing of large numbers of whole eukaryotic genomes is feasible, and so it is often necessary to select genomic regions of interest and to enrich these regions before sequencing. There are several enrichment approaches, each with unique advantages and disadvantages. Here we describe our experiences with the leading target-enrichment technologies, the optimizations that we have performed and typical results that can be obtained using each. We also provide detailed protocols for each technology so that end users can find the best compromise between sensitivity, specificity and uniformity for their particular project.

Mutation spectrum revealed by breakpoint sequencing of human germline CNVs

Precisely characterizing the breakpoints of copy number variants (CNVs) is crucial for assessing their functional impact. However, fewer than 10% of known germline CNVs have been mapped to the single-nucleotide level. We characterized the sequence breakpoints from a dataset of all CNVs detected in three unrelated individuals in previous array-based CNV discovery experiments. We used targeted hybridization-based DNA capture and 454 sequencing to sequence 324 CNV breakpoints, including 315 deletions. We observed two major breakpoint signatures: 70% of the deletion breakpoints have 1−30 bp of microhomology, whereas 33% of deletion breakpoints contain 1−367 bp of inserted sequence. The co-occurrence of microhomology and inserted sequence is low (10%), suggesting that there are at least two different mutational mechanisms. Approximately 5% of the breakpoints represent more complex rearrangements, including local microinversions, suggesting a replication-based strand switching mechanism. Despite a rich literature on DNA repair processes, reconstruction of the molecular events generating each of these mutations is not yet possible.

FRT-seq: amplification-free, strand-specific transcriptome sequencing

We report an alternative approach to transcriptome sequencing for the Illumina Genome Analyzer, in which the reverse transcription reaction takes place on the flowcell. No amplification is performed during the library preparation, so PCR biases and duplicates are avoided, and because the template is poly(A)+ RNA rather than cDNA, the resulting sequences are necessarily strand-specific. The method is compatible with paired- or single-end sequencing.

Simvastatin and sildenafil combine to attenuate pulmonary hypertension

Statins have been proposed to be a potential treatment for pulmonary arterial hypertension. If introduced into clinical practice, the statin would have to be used in conjunction with established therapy. We investigated the effects of combining simvastatin with a phosphodiesterase type-5 inhibitor, sildenafil, in the rat model of hypoxia-induced pulmonary hypertension. Rats were allocated to either: 1) a prevention protocol, to receive simvastatin 20 mg·kg−1·day−1 by intraperitoneal injection or sildenafil 75 mg·kg−1·day−1 orally or the combination (or vehicle) for 2 weeks beginning at the start of exposure to hypoxia (10% inspired oxygen); or 2) a treatment protocol, where the same agents were administered in the last 2 weeks of a 4-week period of hypoxia. In both protocols, the combination of sildenafil and simvastatin lowered pulmonary artery pressure and produced a significantly greater reduction in right ventricular hypertrophy and pulmonary vascular muscularisation than either drug alone. Moreover, the combination augmented significantly endothelial nitric oxide synthase expression and cGMP levels in the lung and right ventricle above that produced by either drug independently and resulted in greater inhibition of RhoA activity. These data suggest that simvastatin can be usefully combined with sildenafil in the treatment of pulmonary arterial hypertension to achieve greater therapeutic benefit.

A comparison of key aspects of gene regulation in Streptomyces coelicolor and Escherichia coli using nucleotide-resolution transcription maps produced in parallel by global and differential RNA sequencing

Streptomyces coelicolor is a model for studying bacteria renowned as the foremost source of natural products used clinically. Post‐genomic studies have revealed complex patterns of gene expression and links to growth, morphological development and individual genes. However, the underlying regulation remains largely obscure, but undoubtedly involves steps after transcription initiation. Here we identify sites involved in RNA processing and degradation as well as transcription within a nucleotide‐resolution map of the transcriptional landscape. This was achieved by combining RNA‐sequencing approaches suited to the analysis of GC‐rich organisms. Escherichia coli was analysed in parallel to validate the methodology and allow comparison. Previously, sites of RNA processing and degradation had not been mapped on a transcriptome‐wide scale for E. coli. Through examples, we show the value of our approach and data sets. This includes the identification of new layers of transcriptional complexity associated with several key regulators of secondary metabolism and morphological development in S. coelicolor and the identification of host‐encoded leaderless mRNA and rRNA processing associated with the generation of specialized ribosomes in E. coli. New regulatory small RNAs were identified for both organisms. Overall the results illustrate the diversity in mechanisms used by different bacterial groups to facilitate and regulate gene expression.

A combination of improved differential and global RNA-seq reveals pervasive transcription initiation and events in all stages of the life-cycle of functional RNAs in Propionibacterium acnes, a major contributor to wide-spread human disease.

BackgroundSequencing of the genome of Propionibacterium acnes produced a catalogue of genes many of which enable this organism to colonise skin and survive exposure to the elements. Despite this platform, there was little understanding of the gene regulation that gives rise to an organism that has a major impact on human health and wellbeing and causes infections beyond the skin. To address this situation, we have undertaken a genome–wide study of gene regulation using a combination of improved differential and global RNA-sequencing and an analytical approach that takes into account the inherent noise within the data.ResultsWe have produced nucleotide-resolution transcriptome maps that identify and differentiate sites of transcription initiation from sites of stable RNA processing and mRNA cleavage. Moreover, analysis of these maps provides strong evidence for ‘pervasive’ transcription and shows that contrary to initial indications it is not biased towards the production of antisense RNAs. In addition, the maps reveal an extensive array of riboswitches, leaderless mRNAs and small non-protein-coding RNAs alongside vegetative promoters and post-transcriptional events, which includes unusual tRNA processing. The identification of such features will inform models of complex gene regulation, as illustrated here for ribonucleotide reductases and a potential quorum-sensing, two-component system.ConclusionsThe approach described here, which is transferable to any bacterial species, has produced a step increase in whole-cell knowledge of gene regulation in P. acnes. Continued expansion of our maps to include transcription associated with different growth conditions and genetic backgrounds will provide a new platform from which to computationally model the gene expression that determines the physiology of P. acnes and its role in human disease.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer

Summary Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.

Low-bias, strand-specific transcriptome Illumina sequencing by on-flowcell reverse transcription (FRT-seq)

The unifying feature of second-generation sequencing technologies is that single template strands are amplified clonally onto a solid surface prior to the sequencing reaction. To convert template strands into a compatible state for attachment to this surface, a multistep library preparation is required, which typically culminates in amplification by the PCR. PCR is an inherently biased process, which decreases the efficiency of data acquisition. Flowcell reverse transcription sequencing is a method of transcriptome sequencing for Illumina sequencers in which the reverse transcription reaction is performed on the flowcell by using unamplified, adapter-ligated mRNA as a template. This approach removes PCR biases and duplicates, generates strand-specific paired-end data and is highly reproducible. The procedure can be performed quickly, taking 2 d to generate clusters from mRNA.

Genome-Wide Scan for Premature Hypertension Supports Linkage to Chromosome 2 in a Large Kyrgyz Family

We report a genome-wide scan for susceptibility loci to hypertension in a single Kyrgyz family where 10 of the affected relatives developed hypertension before the age of 35 years, and some members have suffered stroke. The early onset of disease and the geographic isolation of the Kyrgyz population are both expected to select for an increased influence of genetic factors in hypertension. We genotyped 44 individuals from this Krygyz family with 374 microsatellite markers, covering a 10-centimorgan map. Nonparametric analysis suggests that affected status is linked to loci in the chromosome 2q23 to q37 genomic interval, whereas 2-point parametric analysis returned a logarithm of odds score of 2.67 for marker D2S2330 (2q24.3). Multipoint linkage analysis substantiated the evidence for a hypertension susceptibility allele in the chromosome 2q23 to q36 region. Fine mapping and haplotype analysis implicate that the genetic lesion resides between markers D2S2380 (166.5 cM) and D2S335 (175.9 cM). This finding supports other recent studies of early onset hypertension suggesting that the region 2q24.3 to q31.1 encompasses a novel locus for premature hypertension.

A novel innate lymphoid cell delineates childhood autoimmune arthritis

Inflammation in autoimmune disease is mediated by a complex network of interacting cells. Their identity and cross-talk are encoded in messenger RNA (mRNA). Juvenile idiopathic arthritis (JIA), a chronic autoimmune arthritis of childhood, is characterised by synovial inflammation with infiltration of both innate and adaptive immune cells1. Activated T cells play a role in disease2 but the cell types that drive the recruitment and activation of immune cells within the synovium are not known. Here, we utilised droplet-based and full length single cell mRNA sequencing to obtain a quantitative map of the cellular landscape of JIA. We studied 45,715 cells from the synovial fluid of inflamed knee joints and peripheral blood. We identified a population of synovial innate lymphoid cells (ILCs), shared across patients, that exhibited a unique transcriptional profile in comparison to canonical ILC subtypes. Validation at protein-level across a spectrum of autoimmune arthritides revealed that these ILCs are pathologically expanded in a particular type of JIA. Using statistical tools to assess cellular interactions in synovial fluid, ILCs emerged as a central node of communication, expressing the full repertoire of genes required to orchestrate and maintain the inflammatory milieu. Several ILC-mediated signalling pathways may lend themselves as novel therapeutic targets. Together our findings demonstrate a distinct ILC subtype associated with a tissue-specific childhood autoimmune disease.