Decoding bioactive signals of the RNA secretome: the cell-free messenger RNA catalogue.

Despite gene-expression profiling being one of the most common methods to evaluate molecular dysregulation in tissues, the utilization of cell-free messenger RNA (cf-mRNA) as a blood-based non-invasive biomarker analyte has been limited compared to other RNA classes. Recent advancements in low-input RNA-sequencing and normalization techniques, however, have enabled characterization as well as accurate quantification of cf-mRNAs allowing direct pathological insights. The molecular profile of the cell-free transcriptome in multiple diseases has subsequently been characterized including, prenatal diseases, neurological disorders, liver diseases and cancers suggesting this biological compartment may serve as a disease agnostic platform. With mRNAs packaged in a myriad of extracellular vesicles and particles, these signals may be used to develop clinically actionable, non-invasive disease biomarkers. Here, we summarize the recent scientific developments of extracellular mRNA, biology of extracellular mRNA carriers, clinical utility of cf-mRNA as disease biomarkers, as well as proposed functions in cell and tissue pathophysiology.

Noninvasive stratification of nonalcoholic fatty liver disease by whole transcriptome cell-free mRNA characterization.

Hepatic fibrosis stage is the most important determinant of outcomes in patients with nonalcoholic fatty liver disease (NAFLD). There is an urgent need for noninvasive tests that can accurately stage fibrosis and determine efficacy of interventions. Here, we describe a novel cell-free (cf)-mRNA sequencing approach that can accurately and reproducibly profile low levels of circulating mRNAs and evaluate the feasibility of developing a cf-mRNA-based NAFLD fibrosis classifier. Using separate discovery and validation cohorts with biopsy-confirmed NAFLD (n = 176 and 59, respectively) and healthy subjects (n = 23), we performed serum cf-mRNA RNA-Seq profiling. Differential expression analysis identified 2,498 dysregulated genes between patients with NAFLD and healthy subjects and 134 fibrosis-associated genes in patients with NAFLD. Comparison between cf-mRNA and liver tissue transcripts revealed significant overlap of fibrosis-associated genes and pathways indicating that the circulating cf-mRNA transcriptome reflects molecular changes in the livers of patients with NAFLD. In particular, metabolic and immune pathways reflective of known underlying steatosis and inflammation were highly dysregulated in the cf-mRNA profile of patients with advanced fibrosis. Finally, we used an elastic net ordinal logistic model to develop a classifier that predicts clinically significant fibrosis (F2-F4). In an independent cohort, the cf-mRNA classifier was able to identify 50% of patients with at least 90% probability of clinically significant fibrosis. We demonstrate a novel and robust cf-mRNA-based RNA-Seq platform for noninvasive identification of diverse hepatic molecular disruptions and for fibrosis staging with promising potential for clinical trials and clinical practice.NEW & NOTEWORTHY This work is the first study, to our knowledge, to utilize circulating cell-free mRNA sequencing to develop an NAFLD diagnostic classifier.

Cell-Free DNA and Active Rejection in Kidney Allografts.

Histologic analysis of the allograft biopsy specimen is the standard method used to differentiate rejection from other injury in kidney transplants. Donor-derived cell-free DNA (dd-cfDNA) is a noninvasive test of allograft injury that may enable more frequent, quantitative, and safer assessment of allograft rejection and injury status. To investigate this possibility, we prospectively collected blood specimens at scheduled intervals and at the time of clinically indicated biopsies. In 102 kidney recipients, we measured plasma levels of dd-cfDNA and correlated the levels with allograft rejection status ascertained by histology in 107 biopsy specimens. The dd-cfDNA level discriminated between biopsy specimens showing any rejection (T cell-mediated rejection or antibody-mediated rejection [ABMR]) and controls (no rejection histologically), P<0.001 (receiver operating characteristic area under the curve [AUC], 0.74; 95% confidence interval [95% CI], 0.61 to 0.86). Positive and negative predictive values for active rejection at a cutoff of 1.0% dd-cfDNA were 61% and 84%, respectively. The AUC for discriminating ABMR from samples without ABMR was 0.87 (95% CI, 0.75 to 0.97). Positive and negative predictive values for ABMR at a cutoff of 1.0% dd-cfDNA were 44% and 96%, respectively. Median dd-cfDNA was 2.9% (ABMR), 1.2% (T cell-mediated types ≥IB), 0.2% (T cell-mediated type IA), and 0.3% in controls (P=0.05 for T cell-mediated rejection types ≥IB versus controls). Thus, dd-cfDNA may be used to assess allograft rejection and injury; dd-cfDNA levels <1% reflect the absence of active rejection (T cell-mediated type ≥IB or ABMR) and levels >1% indicate a probability of active rejection.

Proportionally more deleterious genetic variation in European than in African populations.

Quantifying the number of deleterious mutations per diploid human genome is of crucial concern to both evolutionary and medical geneticists. Here we combine genome-wide polymorphism data from PCR-based exon resequencing, comparative genomic data across mammalian species, and protein structure predictions to estimate the number of functionally consequential single-nucleotide polymorphisms (SNPs) carried by each of 15 African American (AA) and 20 European American (EA) individuals. We find that AAs show significantly higher levels of nucleotide heterozygosity than do EAs for all categories of functional SNPs considered, including synonymous, non-synonymous, predicted 'benign', predicted 'possibly damaging' and predicted 'probably damaging' SNPs. This result is wholly consistent with previous work showing higher overall levels of nucleotide variation in African populations than in Europeans. EA individuals, in contrast, have significantly more genotypes homozygous for the derived allele at synonymous and non-synonymous SNPs and for the damaging allele at 'probably damaging' SNPs than AAs do. For SNPs segregating only in one population or the other, the proportion of non-synonymous SNPs is significantly higher in the EA sample (55.4%) than in the AA sample (47.0%; P < 2.3 x 10(-37)). We observe a similar proportional excess of SNPs that are inferred to be 'probably damaging' (15.9% in EA; 12.1% in AA; P < 3.3 x 10(-11)). Using extensive simulations, we show that this excess proportion of segregating damaging alleles in Europeans is probably a consequence of a bottleneck that Europeans experienced at about the time of the migration out of Africa.

An optimized five-gene multi-platform predictor of hormone receptor negative and triple negative breast cancer metastatic risk.

INTRODUCTION: Outcome predictors in use today are prognostic only for hormone receptor-positive (HRpos) breast cancer. Although microarray-derived multigene predictors of hormone receptor-negative (HRneg) and/or triple negative (Tneg) breast cancer recurrence risk are emerging, to date none have been transferred to clinically suitable assay platforms (for example, RT-PCR) or validated against formalin-fixed paraffin-embedded (FFPE) HRneg/Tneg samples. METHODS: Multiplexed RT-PCR was used to assay two microarray-derived HRneg/Tneg prognostic signatures IR-7 and Buck-4) in a pooled FFPE collection of 139 chemotherapy-naïve HRneg breast cancers. The prognostic value of the RTPCR measured gene signatures were evaluated as continuous and dichotomous variables, and in conditional risk models incorporating clinical parameters. An optimized five-gene index was derived by evaluating gene combinations from both signatures. RESULTS: RT-PCR measured IR-7 and Buck-4 signatures proved prognostic as continuous variables; and conditional risk modeling chose nodal status, the IR-7 signature, and tumor grade as significant predictors of distant recurrence (DR). From the Buck-4 and IR-7 signatures, an optimized five-gene (TNFRSF17, CLIC5, HLA-F, CXCL13, XCL2) predictor was generated, referred to as the Integrated Cytokine Score (ICS) based on its functional pathway linkage through interferon-γ and IL-10. Across all FFPE cases, the ICS was prognostic as either a continuous or dichotomous variable, and conditional risk modeling selected nodal status and ICS as DR predictors. Further dichotomization of node-negative/ICS-low FFPE cases identified a subset of low-grade HRneg tumors with <10% 5-year DR risk. The prognostic value of ICS was reaffirmed in two previously studied microarray assayed cohorts containing 274 node-negative and chemotherapy naive HRneg breast cancers, including 95 Tneg cases where it proved prognostically independent of Tneg molecular subtyping. In additional HRneg/Tneg microarray assayed cohorts, the five-gene ICS also proved prognostic irrespective of primary tumor nodal status and adjuvant chemotherapy intervention. CONCLUSION: We advanced the measurement of two previously reported microarray-derived HRneg/Tneg breast cancer prognostic signatures for use in FFPE samples, and derived an optimized five-gene Integrated Cytokine Score (ICS) with multi-platform capability of predicting metastatic outcome from primary HRneg/Tneg tumors independent of nodal status, adjuvant chemotherapy use, and Tneg molecular subtype.

Circulating cell-free messenger RNA secretome characterization of primary sclerosing cholangitis.

BACKGROUND: Primary sclerosing cholangitis (PSC) is a rare chronic cholestatic liver disease characterized by multifocal bile duct strictures. To date, underlying molecular mechanisms of PSC remain unclear, and therapeutic options are limited. METHODS: We performed cell-free messenger RNA (cf-mRNA) sequencing to characterize the circulating transcriptome of PSC and noninvasively investigate potentially bioactive signals that are associated with PSC. Serum cf-mRNA profiles were compared among 50 individuals with PSC, 20 healthy controls, and 235 individuals with NAFLD. Tissue and cell type-of-origin genes that are dysregulated in subjects with PSC were evaluated. Subsequently, diagnostic classifiers were developed using PSC dysregulated cf-mRNA genes. RESULTS: Differential expression analysis of the cf-mRNA transcriptomes of PSC and healthy controls resulted in identification of 1407 dysregulated genes. Furthermore, differentially expressed genes between PSC and healthy controls or NAFLD shared common genes known to be involved in liver pathophysiology. In particular, genes from liver- and specific cell type-origin, including hepatocyte, HSCs, and KCs, were highly abundant in cf-mRNA of subjects with PSC. Gene cluster analysis revealed that liver-specific genes dysregulated in PSC form a distinct cluster, which corresponded to a subset of the PSC subject population. Finally, we developed a cf-mRNA diagnostic classifier using liver-specific genes that discriminated PSC from healthy control subjects using gene transcripts of liver origin. CONCLUSIONS: Blood-based whole-transcriptome cf-mRNA profiling revealed high abundance of liver-specific genes in sera of subjects with PSC, which may be used to diagnose patients with PSC. We identified several unique cf-mRNA profiles of subjects with PSC. These findings may also have utility for noninvasive molecular stratification of subjects with PSC for pharmacotherapy safety and response studies.

A functional polymorphism in the epidermal growth factor gene is associated with risk for hepatocellular carcinoma.

BACKGROUND & AIMS: A single nucleotide polymorphism 61*G (rs4444903) in the epidermal growth factor (EGF) gene has been associated, in 2 case-control studies, with hepatocellular carcinoma (HCC). We tested associations between demographic, clinical, and genetic data and development of HCC, and developed a simple predictive model in a cohort of patients with chronic hepatitis C and advanced fibrosis. METHODS: Black and white subjects from the Hepatitis C Antiviral Long-term Treatment against Cirrhosis (HALT-C) trial (n=816) were followed up prospectively for development of a definite or presumed case of HCC for a median time period of 6.1 years. We used the Cox proportional hazards regression model to determine the hazard ratio for risk of HCC and to develop prediction models. RESULTS: Subjects with EGF genotype G/G had a higher adjusted risk for HCC than those with genotype A/A (hazard ratio, 2.10; 95% confidence interval, 1.05-4.23; P=.03). After adjusting for EGF genotype, blacks had no increased risk of HCC risk compared with whites. Higher serum levels of EGF were observed among subjects with at least one G allele (P=.08); the subset of subjects with EGF G/G genotype and above-median serum levels of EGF had the highest risk of HCC. We developed a simple prediction model that included the EGF genotype to identify patients at low, intermediate, and high risk for HCC; 6-year cumulative HCC incidences were 2.3%, 10.4%, and 26%, respectively. CONCLUSIONS: We associated the EGF genotype G/G with increased risk for HCC; differences in its frequency among black and white subjects might account for differences in HCC incidence between these groups. We developed a model that incorporates EGF genotype and demographic and clinical variables to identify patients at low, intermediate, and high risk for HCC.

A polymorphism that delays fibrosis in hepatitis C promotes alternative splicing of AZIN1, reducing fibrogenesis.

UNLABELLED: Among several single-nucleotide polymorphisms (SNPs) that correlate with fibrosis progression in chronic HCV, an SNP in the antizyme inhibitor (AzI) gene is most strongly associated with slow fibrosis progression. Our aim was to identify the mechanism(s) underlying this observation by exploring the impact of the AzI SNP on hepatic stellate cell (HSC) activity. Seven novel AZIN1 splice variants ("SV2-8") were cloned by polymerase chain reaction from the LX2 human HSC line. Expression of a minigene in LX2 containing the AZIN1 slow-fibrosis SNP yielded a 1.67-fold increase in AZIN1 splice variant 2 (AZIN1 SV2) messenger RNA (mRNA) (P = 0.05). In healthy human leukocytes, the SNP variant also correlated with significantly increased SV2 mRNA. Cells (293T) transfected with short hairpin RNA (shRNA) complementary to the exonic splicing chaperone SRp40 expressed 30% less SRp40 (P = 0.044) and 43% more AzI SV2 (P = 0.021) than control shRNA-expressing cells, mimicking the effect of the sequence variant. LX2 cells transfected with AZIN1 full-length complementary DNA expressed 35% less collagen I mRNA (P = 0.09) and 18% less α-smooth muscle actin mRNA (P = 0.09). Transient transfection of AZIN1 SV2 complementary DNA into LX2 cells reduced collagen I gene expression by 64% (P = 0.001) and α-smooth muscle actin by 43% (P = 0.005) compared to vector-transfected controls, paralleling changes in protein expression. Both AZIN1 and AZIN-SV2 mRNAs are detectable in normal human liver and reduced in HCV cirrhotic livers. The AZIN1-SV2 acts via a polyamine-independent pathway, as it neither interacts with antizyme nor affects the ability of AZIN1 lacking this variant to neutralize antizyme. CONCLUSION: An SNP variant in the AZIN1 gene leads to enhanced generation of a novel alternative splice form that modifies the fibrogenic potential of HSCs.

Impact of patatin-like phospholipase-3 (rs738409 C>G) polymorphism on fibrosis progression and steatosis in chronic hepatitis C.

UNLABELLED: Only 20% of patients with chronic hepatitis C (CHC) will develop cirrhosis, and fibrosis progression remains highly unpredictable. A recent genome-wide association study identified a genetic variant in the patatin-like phospholipase-3 (PNPLA3) gene (rs738409 C>G) associated with steatosis that was further demonstrated to influence severity of fibrosis in nonalcoholic fatty liver disease. The aim of this study was to assess the impact of this polymorphism on histological liver damage and response to antiviral therapy in CHC. We recruited 537 Caucasian CHC patients from three European centers (Brussels, Belgium [n = 229]; Hannover, Germany [n = 171]; Lyon, France [n = 137]); these patients were centrally genotyped for the PNPLA3 (rs738409 C>G) polymorphism. We studied the influence of rs738409 and other variants in the PNPLA3 region on steatosis and fibrosis assessed both in a cross-sectional and longitudinal manner. Seven other variants previously associated with fibrosis progression were included. Finally, we explored the impact of rs738409 on response to standard antiviral therapy using the interferon lambda 3 (IL28B) [rs12979860 C>T] variant both as a comparator and as a positive control. After adjustment for age, sex, body mass index, alcohol consumption, and diabetes, rs738409 mutant G allele homozygote carriers remained at higher risk for steatosis (odds ratio [OR] 2.55, 95% confidence interval [CI] 1.08-6.03, P = 0.034), fibrosis (OR 3.13, 95% CI 1.50-6.51, P = 0.002), and fibrosis progression (OR 2.64, 95% CI 1.22-5.67, P = 0.013). Conversely, rs738409 was not independently associated with treatment failure (OR 1.07, 95% CI 0.46-2.49, P = 0.875) and did not influence clinical or biological variables. CONCLUSION: The PNPLA3 (rs738409 C>G) polymorphism favors steatosis and fibrosis progression in CHC. This polymorphism may represent a valuable genetic predictor and a potential therapeutic target in CHC liver damage.

Evolutionary processes acting on candidate cis-regulatory regions in humans inferred from patterns of polymorphism and divergence.

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.

Survey of extracellular communication of systemic and organ-specific inflammatory responses through cell free messenger RNA profiling in mice.

BACKGROUND: Inflammatory and immune responses are essential and dynamic biological processes that protect the body against acute and chronic adverse stimuli. While conventional protein markers have been used to evaluate systemic inflammatory response, the immunological response to stimulation is complex and involves modulation of a large set of genes and interacting signalling pathways of innate and adaptive immune systems. There is a need for a non-invasive tool that can comprehensively evaluate and monitor molecular dysregulations associated with inflammatory and immune responses in circulation and in inaccessible solid organs. METHODS: Here we utilized cell-free messenger RNA (cf-mRNA) RNA-Seq whole transcriptome profiling and computational biology to temporally assess lipopolysaccharide (LPS) induced and JAK inhibitor modulated inflammatory and immune responses in mouse plasma samples. FINDINGS: Cf-mRNA profiling displayed a pattern of systemic immune responses elicited by LPS and dysregulation of associated pathways. Moreover, attenuation of several inflammatory pathways, including STAT and interferon pathways, were observed following the treatment of JAK inhibitor. We further identified the dysregulation of liver-specific transcripts in cf-mRNA which reflected changes in the gene-expression pattern in this generally inaccessible biological compartment. INTERPRETATION: Using a preclinical mouse model, we demonstrated the potential of plasma cf-mRNA profiling for systemic and organ-specific characterization of drug-induced molecular alterations that are associated with inflammatory and immune responses. FUNDING: Molecular Stethoscope.

A DDX5 S480A polymorphism is associated with increased transcription of fibrogenic genes in hepatic stellate cells.

We recently identified a missense single nucleotide polymorphism (SNP) in DDX5 (rs1140409, p.S480A) that enhances the risk of developing cirrhosis. DDX5 is an ATP-dependent RNA helicase and transcriptional modulator. We hypothesized that the activity of DDX5 in regulating fibrogenic gene transcription in hepatic stellate cells (HSCs) is altered by the S480A SNP. To test this, we employed two approaches: 1) transient overexpression of DDX5 cDNA or siRNA knockdown of endogenous DDX5, with replacement by either DDX5 wild type (WT) or SNP cDNA, or 2) stable expression of exogenous DDX5 WT and SNP in HSC lines. WT DDX5 mRNA in HSCs was inversely correlated with gene expression for alpha2(I) collagen, tissue inhibitor of metalloproteinase-1, and transforming growth factor-beta1. Stable DDX5 SNP-expressing cells had higher basal and transforming growth factor-beta1-stimulated expression and enhanced promoter activities of fibrogenic genes. DDX5 variant-expressing cells also had higher Smad3 and AP-1-responsive reporter activities. In a one-hybrid GAL4 system, co-expression of the DDX5 SNP variant with chimeras of GAL4 DNA binding domain linked to JunD or Sp1 displayed higher transactivation of a GAL4-responsive reporter than that of DDX5 WT. Increased fibrogenic gene expression in DDX5 SNP-expressing cells was associated with reduced recruitment of DDX5 homodimers to responsive promoters, but there was no difference in the recruitment of the co-repressor HDAC1 (histone deacetylase 1). These data suggest that DDX5 is a repressor of fibrogenic genes in HSCs through interaction with transcriptional complexes. The enhanced fibrogenic activity of the DDX5 risk variant is linked to a reduced repressive function toward these target genes.

Predicting cirrhosis and clinical outcomes in patients with advanced chronic hepatitis C with a panel of genetic markers (CRS7).

OBJECTIVES: Genetic factors may play a role in fibrosis progression in patients with chronic hepatitis C (CHC). A cirrhosis risk score (CRS7) with seven single nucleotide polymorphisms was previously shown to correlate with cirrhosis in patients with CHC. This study aimed to assess the validity of CRS7 as a marker of fibrosis progression and cirrhosis and as a predictor of clinical outcomes in patients with CHC. METHODS: A total of 938 patients (677 Caucasians, 165 African-Americans, and 96 Hispanic/Other) in the Hepatitis C Antiviral Long-term Treatment against Cirrhosis Trial were studied. CRS7 was categorized a priori as high risk (n=440), medium risk (n=310), or low risk (n=188). Patients were assessed for four possible outcomes: fibrosis progression, cirrhosis, clinical outcomes [decompensation or hepatocellular carcinoma (HCC)], or HCC alone. RESULTS: Twenty-nine percent (142/493) developed an increase in fibrosis score by greater than or equal to 2 points on follow-up biopsies, 58% had cirrhosis on one or more biopsies, 35% developed at least one clinical outcome, and 13% developed HCC. CRS7 (trend test) was associated with risk for fibrosis progression (P=0.04) with adjusted hazard ratio of 1.27 (95% confidence interval: 1.01-1.58) and with cirrhosis (P=0.05) with adjusted odds ratio of 1.19 (1.00-1.41). Rates of HCC and clinical outcomes were increased in patients with higher CRS7 scores, but were not statistically significant (P=0.12 clinical outcomes, and P=0.07 HCC). A single nucleotide polymorphism in AZIN1 was significantly associated with fibrosis progression. CONCLUSION: CRS7 was validated as a predictor of fibrosis progression and cirrhosis among Hepatitis C Antiviral Long-term Treatment against Cirrhosis patients, who all had advanced fibrosis. CRS7 was not predictive of clinical outcome.

Role of a cirrhosis risk score for the early prediction of fibrosis progression in hepatitis C patients with minimal liver disease.

BACKGROUND & AIMS: Fibrosis progression in patients with chronic hepatitis C (CHC) is highly variable. A Cirrhosis Risk Score (CRS) based on seven genetic variants has been recently developed for identifying patients at risk for cirrhosis. The objective of this study was to assess the role of the CRS for the early prediction of fibrosis progression in CHC patients with mild liver fibrosis. In addition, we evaluated the potential benefit, for prediction accuracy, of a recently described non-invasive fibrosis staging assay, the Enhanced Liver Fibrosis (ELF) test. METHODS: Two separate cohorts of HCV patients (Brussels, Belgium/Hannover, Germany) were retrospectively analyzed. Only patients with a fibrosis Ishak or METAVIR score of F0-F1 at baseline were included. Patients were classified as progressors if they showed an increase ≥2 fibrosis stages at the second histological evaluation after a follow-up ≥5years. The CRS was calculated locally. Genotyping was performed by PCR and oligonucleotide ligation with the resulting signal detected with a Luminex® 200TM and computer analysis. RESULTS: In Brussels, 12/25 patients progressed (48%); similarly in Hannover, 16/31 (52%) patients progressed. In both sample sets, the CRS was significantly associated with fibrosis progression (p=0.050 in Brussels; p=0.018 in Hannover). The ELF test was only a significant predictor in Hannover (p=0.015). In multivariate analysis the CRS remained the only variable associated with fibrosis progression (odds-ratio=2.23, 95%CI 1.21-4.11 p=0.01). CONCLUSIONS: Although conducted on a limited number of patients, this study in two independent centres confirms that the CRS predicts fibrosis progression in initially mild CHC.

Genetic testing for early detection of individuals at risk of coronary heart disease and monitoring response to therapy: challenges and promises.

Coronary heart disease (CHD) often presents suddenly with little warning. Traditional risk factors are inadequate to identify the asymptomatic high-risk individuals. Early identification of patients with subclinical coronary artery disease using noninvasive imaging modalities would allow the early adoption of aggressive preventative interventions. Currently, it is impractical to screen the entire population with noninvasive coronary imaging tools. The use of relatively simple and inexpensive genetic markers of increased CHD risk can identify a population subgroup in which benefit of atherosclerotic imaging modalities would be increased despite nominal cost and radiation exposure. Additionally, genetic markers are fixed and need only be measured once in a patient's lifetime, can help guide therapy selection, and may be of utility in family counseling.

Natural selection on protein-coding genes in the human genome.

Comparisons of DNA polymorphism within species to divergence between species enables the discovery of molecular adaptation in evolutionarily constrained genes as well as the differentiation of weak from strong purifying selection. The extent to which weak negative and positive darwinian selection have driven the molecular evolution of different species varies greatly, with some species, such as Drosophila melanogaster, showing strong evidence of pervasive positive selection, and others, such as the selfing weed Arabidopsis thaliana, showing an excess of deleterious variation within local populations. Here we contrast patterns of coding sequence polymorphism identified by direct sequencing of 39 humans for over 11,000 genes to divergence between humans and chimpanzees, and find strong evidence that natural selection has shaped the recent molecular evolution of our species. Our analysis discovered 304 (9.0%) out of 3,377 potentially informative loci showing evidence of rapid amino acid evolution. Furthermore, 813 (13.5%) out of 6,033 potentially informative loci show a paucity of amino acid differences between humans and chimpanzees, indicating weak negative selection and/or balancing selection operating on mutations at these loci. We find that the distribution of negatively and positively selected genes varies greatly among biological processes and molecular functions, and that some classes, such as transcription factors, show an excess of rapidly evolving genes, whereas others, such as cytoskeletal proteins, show an excess of genes with extensive amino acid polymorphism within humans and yet little amino acid divergence between humans and chimpanzees.

Assessing the evolutionary impact of amino acid mutations in the human genome.

Quantifying the distribution of fitness effects among newly arising mutations in the human genome is key to resolving important debates in medical and evolutionary genetics. Here, we present a method for inferring this distribution using Single Nucleotide Polymorphism (SNP) data from a population with non-stationary demographic history (such as that of modern humans). Application of our method to 47,576 coding SNPs found by direct resequencing of 11,404 protein coding-genes in 35 individuals (20 European Americans and 15 African Americans) allows us to assess the relative contribution of demographic and selective effects to patterning amino acid variation in the human genome. We find evidence of an ancient population expansion in the sample with African ancestry and a relatively recent bottleneck in the sample with European ancestry. After accounting for these demographic effects, we find strong evidence for great variability in the selective effects of new amino acid replacing mutations. In both populations, the patterns of variation are consistent with a leptokurtic distribution of selection coefficients (e.g., gamma or log-normal) peaked near neutrality. Specifically, we predict 27-29% of amino acid changing (nonsynonymous) mutations are neutral or nearly neutral (|s|<0.01%), 30-42% are moderately deleterious (0.01%<|s|<1%), and nearly all the remainder are highly deleterious or lethal (|s|>1%). Our results are consistent with 10-20% of amino acid differences between humans and chimpanzees having been fixed by positive selection with the remainder of differences being neutral or nearly neutral. Our analysis also predicts that many of the alleles identified via whole-genome association mapping may be selectively neutral or (formerly) positively selected, implying that deleterious genetic variation affecting disease phenotype may be missed by this widely used approach for mapping genes underlying complex traits.

[Polymorphism of ornithine decarboxylase antizyme inhibitor 1 gene is associated with liver cirrhosis in Chinese hepatitis B patients].

A cirrhosis risk score (CRS) comprised of single nucleotide polymorphisms (SNPs) in seven genes that predicts the risk of cirrhosis in Caucasian hepatitis C has been reported. The present study was to evaluate the association of 11 separate but related SNPs and the CRS with cirrhosis risk in Chinese hepatitis B patients. A total of 563 Chinese subjects with persistent HBV infection (349 with evident liver cirrhosis and 214 without cirrhosis clinically or pathologically) were studied. The candidate SNPs were detected with a matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) method. The allele frequency and genotype distribution of each polymorphism as well as the CRS value within the cirrhosis and non-cirrhosis subjects were compared. The rs2679757 polymorphism of the antizyme inhibitor 1 (AZIN1) gene was associated with the risk of cirrhosis (x2 = 6.79, P = 0.03, odds ratio for GG+AG versus AA = 1.63, 95% confidence interval = 1.13-2.35). A gene variant (rs886277) in the transient receptor potential cation channel subfamily M, member 5 gene (TRPM5) was associated with liver cirrhosis, but did not reach statistical significance (x2 = 5.77, P = 0.06). Two SNPs (rs4986791, rs62522600) are not polymorphic in Chinese. Genotype frequencies of other SNPs were not different between the cirrhosis and non-cirrhosis groups. The overall CRS values were not different between the cirrhotic and non-cirrhotic groups (median value 0.57 versus 0.62, Z = -1.05, P = 0.29). SNP rs2679757 in the AZIN1 gene is associated with the risk of HBV-related liver cirrhosis in Chinese. The CRS for Caucasian population has limited applicability for predicting liver cirrhosis in Chinese hepatitis B patients. SNPs associated with cirrhosis prognosis in hepatitis B patients and liver diseases with other etiologies warrant further clinical validation.

A multigene predictor of metastatic outcome in early stage hormone receptor-negative and triple-negative breast cancer.

INTRODUCTION: Various multigene predictors of breast cancer clinical outcome have been commercialized, but proved to be prognostic only for hormone receptor (HR) subsets overexpressing estrogen or progesterone receptors. Hormone receptor negative (HRneg) breast cancers, particularly those lacking HER2/ErbB2 overexpression and known as triple-negative (Tneg) cases, are heterogeneous and generally aggressive breast cancer subsets in need of prognostic subclassification, since most early stage HRneg and Tneg breast cancer patients are cured with conservative treatment yet invariably receive aggressive adjuvant chemotherapy. METHODS: An unbiased search for genes predictive of distant metastatic relapse was undertaken using a training cohort of 199 node-negative, adjuvant treatment naive HRneg (including 154 Tneg) breast cancer cases curated from three public microarray datasets. Prognostic gene candidates were subsequently validated using a different cohort of 75 node-negative, adjuvant naive HRneg cases curated from three additional datasets. The HRneg/Tneg gene signature was prognostically compared with eight other previously reported gene signatures, and evaluated for cancer network associations by two commercial pathway analysis programs. RESULTS: A novel set of 14 prognostic gene candidates were identified as outcome predictors: CXCL13, CLIC5, RGS4, RPS28, RFX7, EXOC7, HAPLN1, ZNF3, SSX3, HRBL, PRRG3, ABO, PRTN3, MATN1. A composite HRneg/Tneg gene signature index proved more accurate than any individual candidate gene or other reported multigene predictors in identifying cases likely to remain free of metastatic relapse. Significant positive correlations between the HRneg/Tneg index and three independent immune-related signatures (STAT1, IFN, and IR) were observed, as were consistent negative associations between the three immune-related signatures and five other proliferation module-containing signatures (MS-14, ONCO-RS, GGI, CSR/wound and NKI-70). Network analysis identified 8 genes within the HRneg/Tneg signature as being functionally linked to immune/inflammatory chemokine regulation. CONCLUSIONS: A multigene HRneg/Tneg signature linked to immune/inflammatory cytokine regulation was identified from pooled expression microarray data and shown to be superior to other reported gene signatures in predicting the metastatic outcome of early stage and conservatively managed HRneg and Tneg breast cancer. Further validation of this prognostic signature may lead to new therapeutic insights and spare many newly diagnosed breast cancer patients the need for aggressive adjuvant chemotherapy.

Functional linkage of cirrhosis-predictive single nucleotide polymorphisms of Toll-like receptor 4 to hepatic stellate cell responses.

UNLABELLED: In a recent study, a single nucleotide polymorphism (SNP) of the Toll-like receptor 4 (TLR4) gene (c.1196C>T [rs4986791, p.T399I]) emerged as conferring protection from fibrosis progression compared to a major, wild-type (WT) CC allele (p.T399). The present study examined the functional linkage of this SNP, along with another common, highly cosegregated TLR4 SNP (c.896A>G [rs4986790, p.D299G]), to hepatic stellate cell (HSC) responses. Both HSCs from TLR4(-/-) mice and a human HSC line (LX-2) reconstituted with either TLR4 D299G and/or T399I complementary DNAs were hyporesponsive to lipopolysaccharide (LPS) stimulation compared to those expressing WT TLR4, as assessed by the expression and secretion of LPS-induced inflammatory and chemotactic cytokines (i.e., monocyte chemoattractant protein-1, interleukin-6), down-regulation of bone morphogenic protein and the activin membrane-bound inhibitor expression (an inhibitory transforming growth factor beta pseudoreceptor), and activation of a nuclear factor kappaB (NF-kappaB)-responsive luciferase reporter. In addition, spontaneous apoptosis, as well as apoptosis induced by pathway inhibitors of NF-kappaB, extracellular signal-regulated kinase (ERK), and phosphatidylinositol 3-kinase were greatly increased in HSCs from either TLR4(-/-) or myeloid differentiation factor 88(-/-) (a TLR adaptor protein) mice, as well as in murine HSCs expressing D299G and/or T399I SNPs; increased apoptosis in these lines was accompanied by decreased phospho-ERK and Bcl-2. CONCLUSION: TLR4 D299G and T399I SNPs that are associated with protection from hepatic fibrosis reduce TLR4-mediated inflammatory and fibrogenic signaling and lower the apoptotic threshold of activated HSCs. These findings provide a mechanistic link that explains how specific TLR4 SNPs may regulate the risk of fibrosis progression.