Diagnostic applications of next generation sequencing: working towards quality standards

Over the past 6 years, next generation sequencing (NGS) has been established as a valuable high-throughput method for research in molecular genetics and has successfully been employed in the identification of rare and common genetic variations. All major NGS technology companies providing commercially available instruments (Roche 454, Illumina, Life Technologies) have recently marketed bench top sequencing instruments with lower throughput and shorter run times, thereby broadening the applications of NGS and opening the technology to the potential use for clinical diagnostics. Although the high expectations regarding the discovery of new diagnostic targets and an overall reduction of cost have been achieved, technological challenges in instrument handling, robustness of the chemistry and data analysis need to be overcome. To facilitate the implementation of NGS as a routine method in molecular diagnostics, consistent quality standards need to be developed. Here the authors give an overview of the current standards in protocols and workflows and discuss possible approaches to define quality criteria for NGS in molecular genetic diagnostics

Identification of gene pathways implicated in Alzheimer's disease using longitudinal imaging phenotypes with sparse regression

We present a new method for the detection of gene pathways associated with a multivariate quantitative trait, and use it to identify causal pathways associated with an imaging endophenotype characteristic of longitudinal structural change in the brains of patients with Alzheimer's disease (AD). Our method, known as pathways sparse reduced-rank regression (PsRRR), uses group lasso penalised regression to jointly model the effects of genome-wide single nucleotide polymorphisms (SNPs), grouped into functional pathways using prior knowledge of gene-gene interactions. Pathways are ranked in order of importance using a resampling strategy that exploits finite sample variability. Our application study uses whole genome scans and MR images from 464 subjects in the Alzheimer's Disease Neuroimaging Initiative (ADNI) database. 66,182 SNPs are mapped to 185 gene pathways from the KEGG pathways database. Voxel-wise imaging signatures characteristic of AD are obtained by analysing 3D patterns of structural change at 6, 12 and 24 months relative to baseline. High-ranking, AD endophenotype-associated pathways in our study include those describing chemokine, Jak-stat and insulin signalling pathways, and tight junction interactions. All of these have been previously implicated in AD biology. In a secondary analysis, we investigate SNPs and genes that may be driving pathway selection, and identify a number of previously validated AD genes including CR1, APOE and TOMM40

Unique transcriptomic landscapes identified in idiopathic spontaneous and infection related preterm births compared to normal term births.

Preterm birth (PTB) is leading contributor to infant death in the United States and globally, yet the underlying mechanistic causes are not well understood. Histopathological studies of preterm birth suggest advanced villous maturity may have a role in idiopathic spontaneous preterm birth (isPTB). To better understand pathological and molecular basis of isPTB, we compared placental villous transcriptomes from carefully phenotyped cohorts of PTB due to infection or isPTB between 28-36 weeks gestation and healthy term placentas. Transcriptomic analyses revealed a unique expression signature for isPTB distinct from the age-matched controls that were delivered prematurely due to infection. This signature included the upregulation of three IGF binding proteins (IGFBP1, IGFBP2, and IGFBP6), supporting a role for aberrant IGF signaling in isPTB. However, within the isPTB expression signature, we detected secondary signature of inflammatory markers including TNC, C3, CFH, and C1R, which have been associated with placental maturity. In contrast, the expression signature of the gestational age-matched infected samples included upregulation of proliferative genes along with cell cycling and mitosis pathways. Together, these data suggest an isPTB molecular signature of placental hypermaturity, likely contributing to the premature activation of inflammatory pathways associated with birth and providing a molecular basis for idiopathic spontaneous birth

A novel integrative risk index of papillary thyroid cancer progression combining genomic alterations and clinical factors.

Although the majority of papillary thyroid cancer (PTC) is indolent, a subset of PTC behaves aggressively despite the best available treatment. A major clinical challenge is to reliably distinguish early on between those patients who need aggressive treatment from those who do not. Using a large cohort of PTC samples obtained from The Cancer Genome Atlas (TCGA), we analyzed the association between disease progression and multiple forms of genomic data, such as transcriptome, somatic mutations, and somatic copy number alterations, and found that genes related to FOXM1 signaling pathway were significantly associated with PTC progression. Integrative genomic modeling was performed, controlling for demographic and clinical characteristics, which included patient age, gender, TNM stages, histological subtypes, and history of other malignancy, using a leave-one-out elastic net model and 10-fold cross validation. For each subject, the model from the remaining subjects was used to determine the risk index, defined as a linear combination of the clinical and genomic variables from the elastic net model, and the stability of the risk index distribution was assessed through 2,000 bootstrap resampling. We developed a novel approach to combine genomic alterations and patient-related clinical factors that delineates the subset of patients who have more aggressive disease from those whose tumors are indolent and likely will require less aggressive treatment and surveillance (p = 4.62 × 10-10, log-rank test). Our results suggest that risk index modeling that combines genomic alterations with current staging systems provides an opportunity for more effective anticipation of disease prognosis and therefore enhanced precision management of PTC

Differential expression analysis for sequence count data

*Motivation:* High-throughput nucleotide sequencing provides quantitative readouts in assays for RNA expression (RNA-Seq), protein-DNA binding (ChIP-Seq) or cell counting (barcode sequencing). Statistical inference of differential signal in such data requires estimation of their variability throughout the dynamic range. When the number of replicates is small, error modelling is needed to achieve statistical power.

*Results:* We propose an error model that uses the negative binomial distribution, with variance and mean linked by local regression, to model the null distribution of the count data. The method controls type-I error and provides good detection power. 

*Availability:* A free open-source R software package, _DESeq_, is available from the Bioconductor project and from "http://www-huber.embl.de/users/anders/DESeq":http://www-huber.embl.de/users/anders/DESeq

Minority Stress and Leukocyte Gene Expression In Sexual Minority Men Living With Treated HIV Infection

Sexual minority (i.e., non-heterosexual) individuals experience poorer mental and physical health, accounted for in part by the additional burden of sexual minority stress occurring from being situated in a culture favoring heteronormativity. Informed by previous research, the purpose of this study was to identify the relationship between sexual minority stress and leukocyte gene expression related to inflammation, cancer, immune function, and cardiovascular function. Sexual minority men living with HIV who were on anti-retroviral medication, had viral load \u3c 200 copies/mL, and had biologically confirmed, recent methamphetamine use completed minority stress measures and submitted blood samples for RNA sequencing on leukocytes. Differential gene expression and pathway analyses were conducted comparing those with clinically elevated minority stress (n = 18) and those who did not meet the clinical cutoff (n = 20), covarying reactive urine toxicology results for very recent stimulant use. In total, 90 differentially expressed genes and 138 gene set pathways evidencing 2-directional perturbation were observed at false discovery rate (FDR) \u3c 0.10. Of these, 41 of the differentially expressed genes and 35 of the 2-directionally perturbed pathways were identified as functionally related to hypothesized mechanisms of inflammation, cancer, immune function, and cardiovascular function. The neuroactive-ligand receptor pathway (implicated in cancer development) was identified using signaling pathway impact analysis. Our results suggest several potential biological pathways for future work investigating the relationship between sexual minority stress and health