Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics.

PURPOSE: Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies. METHODS: Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations. RESULTS: Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases. CONCLUSION: NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Comparative evaluation for the globin gene depletion methods for mRNA sequencing using the whole blood-derived total RNAs.

BACKGROUND: There are challenges in generating mRNA-Seq data from whole-blood derived RNA as globin gene and rRNA are frequent contaminants. Given the abundance of erythrocytes in whole blood, globin genes comprise some 80% or more of the total RNA. Therefore, depletion of globin gene RNA and rRNA are critical steps required to have adequate coverage of reads mapping to the reference transcripts and thus reduce the total cost of sequencing. In this study, we directly compared the performance of probe hybridization (GLOBINClear Kit and Globin-Zero Gold rRNA Removal Kit) and RNAse-H enzymatic depletion (NEBNext® Globin & rRNA Depletion Kit and Ribo-Zero Plus rRNA Depletion Kit) methods from 1 µg of whole blood-derived RNA on mRNA-Seq profiling. All RNA samples were treated with DNaseI for additional cleanup before the depletion step and were processed for poly-A selection for library generation. RESULTS: Probe hybridization revealed a better overall performance than the RNAse-H enzymatic depletion method, detecting a higher number of genes and transcripts without 3' region bias. After depletion, samples treated with probe hybridization showed globin genes at 0.5% (±0.6%) of the total mapped reads; the RNAse-H enzymatic depletion had 3.2% (±3.8%). Probe hybridization showed more junction reads and transcripts compared with RNAse-H enzymatic depletion and also had a higher correlation (R > 0.9) than RNAse-H enzymatic depletion (R > 0.85). CONCLUSION: In this study, our results showed that 1 µg of high-quality RNA from whole blood could be routinely used for transcriptional profiling analysis studies with globin gene and rRNA depletion pre-processing. We also demonstrated that the probe hybridization depletion method is better suited to mRNA sequencing analysis with minimal effect on RNA quality during depletion procedures.

Bacterial Single Cell Whole Transcriptome Amplification in Microfluidic Platform Shows Putative Gene Expression Heterogeneity.

Single cell RNA sequencing is a technology that provides the capability of analyzing the transcriptome of a single cell from a population. So far, single cell RNA sequencing has been focused mostly on human cells due to the larger starting amount of RNA template for subsequent amplification. One of the major challenges of applying single cell RNA sequencing to microbial cells is to amplify the femtograms of the RNA template to obtain sufficient material for downstream sequencing with minimal contamination. To achieve this goal, efforts have been focused on multiround RNA amplification, but would introduce additional contamination and bias. In this work, we for the first time coupled a microfluidic platform with multiple displacement amplification technology to perform single cell whole transcriptome amplification and sequencing of Porphyromonas somerae, a microbe of interest in endometrial cancer, as a proof-of-concept demonstration of using single cell RNA sequencing tool to unveil gene expression heterogeneity in single microbial cells. Our results show that the bacterial single-cell gene expression regulation is distinct across different cells, supporting widespread heterogeneity.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Homozygous/Compound Heterozygous Triadin Mutations Associated With Autosomal-Recessive Long-QT Syndrome and Pediatric Sudden Cardiac Arrest: Elucidation of the Triadin Knockout Syndrome.

BACKGROUND: Long-QT syndrome (LQTS) may result in syncope, seizures, or sudden cardiac arrest. Although 16 LQTS-susceptibility genes have been discovered, 20% to 25% of LQTS remains genetically elusive. METHODS AND RESULTS: We performed whole-exome sequencing child-parent trio analysis followed by recessive and sporadic inheritance modeling and disease-network candidate analysis gene ranking to identify a novel underlying genetic mechanism for LQTS. Subsequent mutational analysis of the candidate gene was performed with polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing on a cohort of 33 additional unrelated patients with genetically elusive LQTS. After whole-exome sequencing and variant filtration, a homozygous p.D18fs*13 TRDN-encoded triadin frameshift mutation was discovered in a 10-year-old female patient with LQTS with a QTc of 500 milliseconds who experienced recurrent exertion-induced syncope/cardiac arrest beginning at 1 year of age. Subsequent mutational analysis of TRDN revealed either homozygous or compound heterozygous frameshift mutations in 4 of 33 unrelated cases of LQTS (12%). All 5 TRDN-null patients displayed extensive T-wave inversions in precordial leads V1 through V4, with either persistent or transient QT prolongation and severe disease expression of exercise-induced cardiac arrest in early childhood (≤3 years of age) and required aggressive therapy. The overall yield of TRDN mutations was significantly greater in patients ≤10 years of age (5 of 10, 50%) compared with older patients (0 of 24, 0%; P=0.0009). CONCLUSIONS: We identified TRDN as a novel underlying genetic basis for recessively inherited LQTS. All TRDN-null patients had strikingly similar phenotypes. Given the recurrent nature of potential lethal arrhythmias, patients fitting this phenotypic profile should undergo cardiac TRDN genetic testing.

Promoter-proximal transcription factor binding is transcriptionally active when coupled with nucleosome repositioning in immediate vicinity.

Previous studies have analyzed patterns of transcription, transcription factor (TF) binding or mapped nucleosome occupancy across the genome. These suggest that the three aspects are genetically connected but the cause and effect relationships are still unknown. For example, physiologic TF binding studies involve many TFs, consequently, it is difficult to assign nucleosome reorganization to the binding site occupancy of any particular TF. Therefore, several aspects remain unclear: does TF binding influence nucleosome (re)organizations locally or impact the chromatin landscape at a more global level; are all or only a fraction of TF binding a result of reorganization in nucleosome occupancy and do all TF binding and associated changes in nucleosome occupancy result in altered gene expression? With these in mind, following characterization of two states (before and after induction of a single TF of choice) we determined: (i) genomic binding sites of the TF, (ii) promoter nucleosome occupancy and (iii) transcriptome profiles. Results demonstrated that promoter-proximal TF binding influenced expression of the target gene when it was coupled to nucleosome repositioning at or close to its binding site in most cases. In contrast, only in few cases change in target gene expression was found when TF binding occurred without local nucleosome reorganization.

Genetic alterations associated with progression from pancreatic intraepithelial neoplasia to invasive pancreatic tumor.

BACKGROUND & AIMS: Increasing grade of pancreatic intraepithelial neoplasia (PanIN) has been associated with progression to pancreatic ductal adenocarcinoma (PDAC). However, the mechanisms that control progression from PanINs to PDAC are not well understood. We investigated the genetic alterations involved in this process. METHODS: Genomic DNA samples from laser-capture microdissected PDACs and adjacent PanIN2 and PanIN3 lesions from 10 patients with pancreatic cancer were analyzed by exome sequencing. RESULTS: Similar numbers of somatic mutations were identified in PanINs and tumors, but the mutational load varied greatly among cases. Ten of the 15 isolated PanINs shared more than 50% of somatic mutations with associated tumors. Mutations common to tumors and clonally related PanIN2 and PanIN3 lesions were identified as genes that could promote carcinogenesis. KRAS and TP53 frequently were altered in PanINs and tumors, but few other recurrently modified genes were detected. Mutations in DNA damage response genes were prevalent in all samples. Genes that encode proteins involved in gap junctions, the actin cytoskeleton, the mitogen-activated protein kinase signaling pathway, axon guidance, and cell-cycle regulation were among the earliest targets of mutagenesis in PanINs that progressed to PDAC. CONCLUSIONS: Early stage PanIN2 lesions appear to contain many of the somatic gene alterations required for PDAC development.

Genome-wide analysis reveals recurrent structural abnormalities of TP63 and other p53-related genes in peripheral T-cell lymphomas.

Peripheral T-cell lymphomas (PTCLs) are aggressive malignancies of mature T lymphocytes with 5-year overall survival rates of only ∼ 35%. Improvement in outcomes has been stymied by poor understanding of the genetics and molecular pathogenesis of PTCL, with a resulting paucity of molecular targets for therapy. We developed bioinformatic tools to identify chromosomal rearrangements using genome-wide, next-generation sequencing analysis of mate-pair DNA libraries and applied these tools to 16 PTCL patient tissue samples and 6 PTCL cell lines. Thirteen recurrent abnormalities were identified, of which 5 involved p53-related genes (TP53, TP63, CDKN2A, WWOX, and ANKRD11). Among these abnormalities were novel TP63 rearrangements encoding fusion proteins homologous to ΔNp63, a dominant-negative p63 isoform that inhibits the p53 pathway. TP63 rearrangements were seen in 11 (5.8%) of 190 PTCLs and were associated with inferior overall survival; they also were detected in 2 (1.2%) of 164 diffuse large B-cell lymphomas. As TP53 mutations are rare in PTCL compared with other malignancies, our findings suggest that a constellation of alternate genetic abnormalities may contribute to disruption of p53-associated tumor suppressor function in PTCL.

Novel TRAF1-ALK fusion identified by deep RNA sequencing of anaplastic large cell lymphoma.

Chromosomal translocations leading to expression of abnormal fusion proteins play a major role in the pathogenesis of various hematologic malignancies. The recent development of high-throughput, "deep" sequencing has allowed discovery of novel translocations leading to a rapid increase in understanding these diseases. Translocations involving the anaplastic lymphoma kinase (ALK) gene leading to ALK fusion proteins originally were discovered in anaplastic large cell lymphomas (ALCLs). Among ALCLs, NPM1-ALK fusions are most common and lead to nuclear localization of the fusion protein. Here, we present a 50-year-old male with ALCL demonstrating cytoplasmic ALK immunoreactivity only, suggesting the presence of a non-NPM1 fusion partner. We performed deep RNA sequencing of tumor tissue from this patient and identified a novel transcript fusing Exon 6 of TRAF1 to Exon 20 of ALK. The TRAF1-ALK fusion transcript was confirmed at the mRNA level by Sanger sequencing and the fusion protein was visualized by Western blot. The discovery of this TRAF1-ALK fusion expands the diversity of known ALK fusion partners and highlights the power of deep sequencing for fusion transcript discovery. © 2013 Wiley Periodicals, Inc.

FKBP5 genetic variation: association with selective serotonin reuptake inhibitor treatment outcomes in major depressive disorder.

OBJECTIVES: FKBP51 (51 kDa immunophilin) acts as a modulator of the glucocorticoid receptor and a negative regulator of the Akt pathway. Genetic variation in FKBP5 plays a role in antidepressant response. The aim of this study was to comprehensively assess the role of genetic variation in FKBP5, identified by both Sanger and Next Generation DNA resequencing, as well as genome-wide single nucleotide polymorphisms (SNPs) associated with FKBP5 expression in the response to the selective serotonin reuptake inhibitor (SSRI) treatment of major depressive disorder. METHODS: We identified 657 SNPs in FKBP5 by Next Generation sequencing of 96 DNA samples from white patients, and 149 SNPs were selected for the genotyping together with 235 SNPs that were trans-associated with variation in FKBP5 expression in lymphoblastoid cells. A total of 529 DNA samples from the Mayo Clinic PGRN-SSRI Pharmacogenomic trial for which genome-wide SNPs had already been obtained were genotyped for these 384 SNPs, and associations with treatment outcomes were determined. The most significant SNPs were genotyped using 96 DNA samples from white non-Hispanic patients of the NIMH-supported Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study to attempt replication, followed by functional genomic studies. RESULTS: Genotype-phenotype association analysis indicated that rs352428 was associated with both 8-week treatment response in the Mayo study (odds ratio=0.49; P=0.003) and 6-week response in the STAR*D replication study (odds ratio=0.74; P=0.05). The electrophoresis mobility shift assay and the reporter gene assay confirmed the possible role of this SNP in transcription regulation. CONCLUSION: This comprehensive FKBP5 sequence study provides insight into the role of common genetic polymorphisms that might influence SSRI treatment outcomes in major depressive disorder patients.

Discovery of recurrent t(6;7)(p25.3;q32.3) translocations in ALK-negative anaplastic large cell lymphomas by massively parallel genomic sequencing.

The genetics of peripheral T-cell lymphomas are poorly understood. The most well-characterized abnormalities are translocations involving ALK, occurring in approximately half of anaplastic large cell lymphomas (ALCLs). To gain insight into the genetics of ALCLs lacking ALK translocations, we combined mate-pair DNA library construction, massively parallel ("Next Generation") sequencing, and a novel bioinformatic algorithm. We identified a balanced translocation disrupting the DUSP22 phosphatase gene on 6p25.3 and adjoining the FRA7H fragile site on 7q32.3 in a systemic ALK-negative ALCL. Using fluorescence in situ hybridization, we demonstrated that the t(6;7)(p25.3;q32.3) was recurrent in ALK-negative ALCLs. Furthermore, t(6;7)(p25.3;q32.3) was associated with down-regulation of DUSP22 and up-regulation of MIR29 microRNAs on 7q32.3. These findings represent the first recurrent translocation reported in ALK-negative ALCL and highlight the utility of massively parallel genomic sequencing to discover novel translocations in lymphoma and other cancers.

Identification of gene mutations in autosomal dominant polycystic kidney disease through targeted resequencing.

Mutations in two large multi-exon genes, PKD1 and PKD2, cause autosomal dominant polycystic kidney disease (ADPKD). The duplication of PKD1 exons 1-32 as six pseudogenes on chromosome 16, the high level of allelic heterogeneity, and the cost of Sanger sequencing complicate mutation analysis, which can aid diagnostics of ADPKD. We developed and validated a strategy to analyze both the PKD1 and PKD2 genes using next-generation sequencing by pooling long-range PCR amplicons and multiplexing bar-coded libraries. We used this approach to characterize a cohort of 230 patients with ADPKD. This process detected definitely and likely pathogenic variants in 115 (63%) of 183 patients with typical ADPKD. In addition, we identified atypical mutations, a gene conversion, and one missed mutation resulting from allele dropout, and we characterized the pattern of deep intronic variation for both genes. In summary, this strategy involving next-generation sequencing is a model for future genetic characterization of large ADPKD populations.

Relationship between BCL2 mutations and follicular lymphoma outcome in the chemoimmunotherapy era.

How to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression is unclear. Building on a prior study demonstrating the early transformation of FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS, and MYC, in 199 newly diagnosed grade 1 and 2 FLs. BCL2 mutations with VAF ≥20% occurred in 52% of cases. Among 97 FL patients who did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF ≥20% were associated with increased transformation risk (HR 3.01, 95% CI 1.04-8.78, p = 0.043) and a trend toward shorter event-free survival (EFS, median 20 months with mutations versus 54 months without, p = 0.052). Other sequenced genes were less frequently mutated and did not increase the prognostic value of the panel. Across the entire population, nonsynonymous BCL2 mutations at VAF ≥20% were associated with decreased EFS (HR 1.55, 95% CI 1.02-2.35, p = 0.043 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.82, 95% CI 1.05-3.17, p = 0.034). Thus, high VAF nonsynonymous BCL2 mutations remain prognostic even in the chemoimmunotherapy era.

Serine hydroxymethyltransferase 1 and 2: gene sequence variation and functional genomic characterization.

Serine hydroxymethyltransferase (SHMT) catalyzes the transfer of a ß-carbon from serine to tetrahydrofolate to form glycine and 5,10-methylene-tetrahydrofolate. This reaction plays an important role in neurotransmitter synthesis and metabolism. We set out to resequence SHMT1 and SHMT2, followed by functional genomic studies. We identified 87 and 60 polymorphisms in SHMT1 and SHMT2, respectively. We observed no significant functional effect of the 13 non-synonymous single-nucleotide polymorphism (SNPs) in these genes, either on catalytic activity or protein quantity. We imputed additional variants across the two genes using '1000 Genomes' data, and identified 14 variants that were significantly associated (p<1.0E-10) with SHMT1 messenger RNA expression in lymphoblastoid cell lines. Many of these SNPs were also significantly correlated with basal SHMT1 protein expression in 268 human liver biopsy samples. Reporter gene assays suggested that the SHMT1 promoter SNP, rs669340, contributed to this variation. Finally, SHMT1 and SHMT2 expression were significantly correlated with those of other Folate and Methionine Cycle genes at both the messenger RNA and protein levels. These experiments represent a comprehensive study of SHMT1 and SHMT2 gene sequence variation and its functional implications. In addition, we obtained preliminary indications that these genes may be co-regulated with other Folate and Methionine Cycle genes.

Gemcitabine metabolic pathway genetic polymorphisms and response in patients with non-small cell lung cancer.

BACKGROUND AND OBJECTIVE: Gemcitabine is widely used to treat non-small cell lung cancer (NSCLC). The aim of this study was to assess the pharmacogenomic effects of the entire gemcitabine metabolic pathway, we genotyped single nucleotide polymorphisms (SNPs) within the 17 pathway genes using DNA samples from patients with NSCLC treated with gemcitabine to determine the effect of genetic variants within gemcitabine pathway genes on overall survival (OS) of patients with NSCLC after treatment of gemcitabine. METHODS: Eight of the 17 pathway genes were resequenced with DNA samples from Coriell lymphoblastoid cell lines (LCLs) using Sanger sequencing for all exons, exon-intron junctions, and 5'-, 3'-UTRs. A total of 107 tagging SNPs were selected on the basis of the resequencing data for the eight genes and on HapMap data for the remaining nine genes, followed by successful genotyping of 394 NSCLC patient DNA samples. Association of SNPs/haplotypes with OS was performed using the Cox regression model, followed by functional studies performed with LCLs and NSCLC cell lines. RESULTS: Five SNPs in four genes (CDA, NT5C2, RRM1, and SLC29A1) showed associations with OS of those patients with NSCLC, as well as nine haplotypes in four genes (RRM1, RRM2, SLC28A3, and SLC29A1) with a P value of less than 0.05. Genotype imputation using the LCLs was performed for a region of 200 kb surrounding those SNPs, followed by association studies with gemcitabine cytotoxicity. Functional studies demonstrated that downregulation of SLC29A1, NT5C2, and RRM1 in NSCLC cell lines altered cell susceptibility to gemcitabine. CONCLUSION: These studies help in identifying biomarkers to predict gemcitabine response in NSCLC, a step toward the individualized chemotherapy of lung cancer.

Human liver methionine cycle: MAT1A and GNMT gene resequencing, functional genomics, and hepatic genotype-phenotype correlation.

The "methionine cycle" plays a critical role in the regulation of concentrations of (S)-adenosylmethionine (AdoMet), the major biological methyl donor. We set out to study sequence variation in genes encoding the enzyme that synthesizes AdoMet in liver, methionine adenosyltransferase 1A (MAT1A) and the major hepatic AdoMet using enzyme, glycine N-methyltransferase (GNMT), as well as functional implications of that variation. We resequenced MAT1A and GNMT using DNA from 288 subjects of three ethnicities, followed by functional genomic and genotype-phenotype correlation studies performed with 268 hepatic biopsy samples. We identified 44 and 42 polymorphisms in MAT1A and GNMT, respectively. Quantitative Western blot analyses for the human liver samples showed large individual variation in MAT1A and GNMT protein expression. Genotype-phenotype correlation identified two genotyped single-nucleotide polymorphisms (SNPs), reference SNP (rs) 9471976 (corrected p = 3.9 × 10(-10)) and rs11752813 (corrected p = 1.8 × 10(-5)), and 42 imputed SNPs surrounding GNMT that were significantly associated with hepatic GNMT protein levels (corrected p values < 0.01). Reporter gene studies showed that variant alleles for both genotyped SNPs resulted in decreased transcriptional activity. Correlation analyses among hepatic protein levels for methionine cycle enzymes showed significant correlations between GNMT and MAT1A (p = 1.5 × 10(-3)) and between GNMT and betaine homocysteine methyltransferase (p = 1.6 × 10(-7)). Our discovery of SNPs that are highly associated with hepatic GNMT protein expression as well as the "coordinate regulation" of methionine cycle enzyme protein levels provide novel insight into the regulation of this important human liver biochemical pathway.

Non-random genetic alterations in the cyanobacterium Nostoc sp. exposed to space conditions.

Understanding the impact of long-term exposure of microorganisms to space is critical in understanding how these exposures impact the evolution and adaptation of microbial life under space conditions. In this work we subjected Nostoc sp. CCCryo 231-06, a cyanobacterium capable of living under many different ecological conditions, and also surviving in extreme ones, to a 23-month stay at the International Space Station (the Biology and Mars Experiment, BIOMEX, on the EXPOSE-R2 platform) and returned it to Earth for single-cell genome analysis. We used microfluidic technology and single cell sequencing to identify the changes that occurred in the whole genome of single Nostoc cells. The variant profile showed that biofilm and photosystem associated loci were the most altered, with an increased variant rate of synonymous base pair substitutions. The cause(s) of these non-random alterations and their implications to the evolutionary potential of single bacterial cells under long-term cosmic exposure warrants further investigation.

Whole genome sequencing of cyanobacterium Nostoc sp. CCCryo 231-06 using microfluidic single cell technology.

The Nostoc sp. strain CCCryo 231-06 is a cyanobacterial strain capable of surviving under extreme conditions and thus is of great interest for the astrobiology community. The knowledge of its complete genome sequence would serve as a guide for further studies. However, a major concern has been placed on the effects of contamination on the quality of sequencing data without a reference genome. Here, we report the use of microfluidic technology combined with single cell sequencing and de novo assembly to minimize the contamination and recover the complete genome of the Nostoc strain CCCryo 231-06 with high quality. 100% of the whole genome was recovered with all contaminants removed and a strongly supported phylogenetic tree. The data reported can be useful for comparative genomics for phylogenetic and taxonomic studies. The method used in this work can be applied to studies that require high-quality assemblies of genomes of unknown microorganisms.