Helen McDonald

Frequent mutation of histone-modifying genes in non-Hodgkin lymphoma

Follicular lymphoma (FL) and diffuse large B-cell lymphoma (DLBCL) are the two most common non-Hodgkin lymphomas (NHLs). Here we sequenced tumour and matched normal DNA from 13 DLBCL cases and one FL case to identify genes with mutations in B-cell NHL. We analysed RNA-seq data from these and another 113 NHLs to identify genes with candidate mutations, and then re-sequenced tumour and matched normal DNA from these cases to confirm 109 genes with multiple somatic mutations. Genes with roles in histone modification were frequent targets of somatic mutation. For example, 32% of DLBCL and 89% of FL cases had somatic mutations in MLL2, which encodes a histone methyltransferase, and 11.4% and 13.4% of DLBCL and FL cases, respectively, had mutations in MEF2B, a calcium-regulated gene that cooperates with CREBBP and EP300 in acetylating histones. Our analysis suggests a previously unappreciated disruption of chromatin biology in lymphomagenesis.

Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells

MicroRNAs (miRNAs) are emerging as important, albeit poorly characterized, regulators of biological processes. Key to further elucidation of their roles is the generation of more complete lists of their numbers and expression changes in different cell states. Here, we report a new method for surveying the expression of small RNAs, including microRNAs, using Illumina sequencing technology. We also present a set of methods for annotating sequences deriving from known miRNAs, identifying variability in mature miRNA sequences, and identifying sequences belonging to previously unidentified miRNA genes. Application of this approach to RNA from human embryonic stem cells obtained before and after their differentiation into embryoid bodies revealed the sequences and expression levels of 334 known plus 104 novel miRNA genes. One hundred seventy-one known and 23 novel microRNA sequences exhibited significant expression differences between these two developmental states. Owing to the increased number of sequence reads, these libraries represent the deepest miRNA sampling to date, spanning nearly six orders of magnitude of expression. The predicted targets of those miRNAs enriched in either sample shared common features. Included among the high-ranked predicted gene targets are those implicated in differentiation, cell cycle control, programmed cell death, and transcriptional regulation.

Profiling the HeLa S3 transcriptome using randomly primed cDNA and massively parallel short-read sequencing

Sequence-based methods for transcriptome characterization have typically relied on generation of either serial analysis of gene expression tags or expressed sequence tags. Although such approaches have the potential to enumerate transcripts by counting sequence tags derived from them, they typically do not robustly survey the majority of transcripts along their entire length. Here we show that massively parallel sequencing of randomly primed cDNAs, using a next-generation sequencing-by-synthesis technology, offers the potential to generate relative measures of mRNA and individual exon abundance while simultaneously profiling the prevalence of both annotated and novel exons and exon-splicing events. This technique identifies known single nucleotide polymorphisms (SNPs) as well as novel single-base variants. Analysis of these variants, and previously unannotated splicing events in the HeLa S3 cell line, reveals an overrepresentation of gene categories including those previously implicated in cancer.

Next-generation tag sequencing for cancer gene expression profiling

We describe a new method, Tag-seq, which employs ultra high-throughput sequencing of 21 base pair cDNA tags for sensitive and cost-effective gene expression profiling. We compared Tag-seq data to LongSAGE data and observed improved representation of several classes of rare transcripts, including transcription factors, antisense transcripts, and intronic sequences, the latter possibly representing novel exons or genes. We observed increases in the diversity, abundance, and dynamic range of such rare transcripts and took advantage of the greater dynamic range of expression to identify, in cancers and normal libraries, altered expression ratios of alternative transcript isoforms. The strand-specific information of Tag-seq reads further allowed us to detect altered expression ratios of sense and antisense (S-AS) transcripts between cancer and normal libraries. S-AS transcripts were enriched in known cancer genes, while transcript isoforms were enriched in miRNA targeting sites. We found that transcript abundance had a stronger GC-bias in LongSAGE than Tag-seq, such that AT-rich tags were less abundant than GC-rich tags in LongSAGE. Tag-seq also performed better in gene discovery, identifying >98% of genes detected by LongSAGE and profiling a distinct subset of the transcriptome characterized by AT-rich genes, which was expressed at levels below those detectable by LongSAGE. Overall, Tag-seq is sensitive to rare transcripts, has less sequence composition bias relative to LongSAGE, and allows differential expression analysis for a greater range of transcripts, including transcripts encoding important regulatory molecules.

Alternative expression analysis by RNA sequencing

In alternative expression analysis by sequencing (ALEXA-seq), we developed a method to analyze massively parallel RNA sequence data to catalog transcripts and assess differential and alternative expression of known and predicted mRNA isoforms in cells and tissues. As proof of principle, we used the approach to compare fluorouracil-resistant and -nonresistant human colorectal cancer cell lines. We assessed the sensitivity and specificity of the approach by comparison to exon tiling and splicing microarrays and validated the results with reverse transcription–PCR, quantitative PCR and Sanger sequencing. We observed global disruption of splicing in fluorouracil-resistant cells characterized by expression of new mRNA isoforms resulting from exon skipping, alternative splice site usage and intron retention. Alternative expression annotation databases, source code, a data viewer and other resources to facilitate analysis are available at http://www.alexaplatform.org/alexa_seq/.

In-depth characterization of the microRNA transcriptome in a leukemia progression model

MicroRNAs (miRNAs) have been shown to play important roles in physiological as well as multiple malignant processes, including acute myeloid leukemia (AML). In an effort to gain further insight into the role of miRNAs in AML, we have applied the Illumina massively parallel sequencing platform to carry out an in-depth analysis of the miRNA transcriptome in a murine leukemia progression model. This model simulates the stepwise conversion of a myeloid progenitor cell by an engineered overexpression of the nucleoporin 98 (NUP98)-homeobox HOXD13 fusion gene (ND13), to aggressive AML inducing cells upon transduction with the oncogenic collaborator Meis1. From this data set, we identified 307 miRNA/miRNA species in the ND13 cells and 306 miRNA/miRNA species in ND13+Meis1 cells, corresponding to 223 and 219 miRNA genes. Sequence counts varied between two and 136,558, indicating a remarkable expression range between the detected miRNA species. The large number of miRNAs expressed and the nature of differential expression suggest that leukemic progression as modeled here is dictated by the repertoire of shared, but differentially expressed miRNAs. Our finding of extensive sequence variations (isomiRs) for almost all miRNA and miRNA species adds additional complexity to the miRNA transcriptome. A stringent target prediction analysis coupled with in vitro target validation revealed the potential for miRNA-mediated release of oncogenes that facilitates leukemic progression from the preleukemic to leukemia inducing state. Finally, 55 novel miRNAs species were identified in our data set, adding further complexity to the emerging world of small RNAs.

Use of mutation profiles to refine the classification of endometrial carcinomas.

The classification of endometrial carcinomas is based on pathological assessment of tumour cell type; the different cell types (endometrioid, serous, carcinosarcoma, mixed, undifferentiated, and clear cell) are associated with distinct molecular alterations. This current classification system for high‐grade subtypes, in particular the distinction between high‐grade endometrioid (EEC‐3) and serous carcinomas (ESC), is limited in its reproducibility and prognostic abilities. Therefore, a search for specific molecular classifiers to improve endometrial carcinoma subclassification is warranted. We performed target enrichment sequencing on 393 endometrial carcinomas from two large cohorts, sequencing exons from the following nine genes: ARID1A, PPP2R1A, PTEN, PIK3CA, KRAS, CTNNB1, TP53, BRAF, and PPP2R5C. Based on this gene panel, each endometrial carcinoma subtype shows a distinct mutation profile. EEC‐3s have significantly different frequencies of PTEN and TP53 mutations when compared to low‐grade endometrioid carcinomas. ESCs and EEC‐3s are distinct subtypes with significantly different frequencies of mutations in PTEN, ARID1A, PPP2R1A, TP53, and CTNNB1. From the mutation profiles, we were able to identify subtype outliers, ie cases diagnosed morphologically as one subtype but with a mutation profile suggestive of a different subtype. Careful review of these diagnostically challenging cases suggested that the original morphological classification was incorrect in most instances. The molecular profile of carcinosarcomas suggests two distinct mutation profiles for these tumours: endometrioid‐type (PTEN, PIK3CA, ARID1A, KRAS mutations) and serous‐type (TP53 and PPP2R1A mutations). While this nine‐gene panel does not allow for a purely molecularly based classification of endometrial carcinoma, it may prove useful as an adjunct to morphological classification and serve as an aid in the classification of problematic cases. If used in practice, it may lead to improved diagnostic reproducibility and may also serve to stratify patients for targeted therapeutics. Copyright

Involvement of the X chromosome in non-Hodgkin lymphoma.

Gain of an X chromosome is observed as a secondary, acquired karyotypic alteration in a significant proportion of malignant lymphomas. To determine the potential involvement of X‐linked genes in neoplastic development, we have analyzed the inactivation status of the supernumerary X chromosome in lymphomas in both male and female patients. In males, neither methylation of FMR1 nor expression of XIST was detected, demonstrating that the duplicated chromosome was not subject to inactivation. In females, both expressed polymorphisms and polymorphisms associated with methylation differences between the active and inactive X chromosome were analyzed to determine whether the duplicated chromosome was active or inactive. To facilitate this analysis, allele‐specific PCR primers were designed for detection of previously described polymorphisms in the IDSX and G6PD genes. The female lymphomas were shown to be clonal in origin, and duplication of either the active (5 cases) or inactive (4 cases) X chromosome was observed. Correlations between clinical status and the inactivation status of the X chromosome involved in the duplication were not observed in our relatively small sample, although 4/4 informative cases with a t(14;18) showed duplication of the active X chromosome. In the course of these studies, we detected hypermethylation of the androgen receptor (AR) locus in an extremely high proportion of both male (7/9) and female (9/10) samples. These results are discussed with respect to whether sex chromosome aneuploidies in tumors are involved in, or simply the result of, the neoplastic process. Genes Chromosomes Cancer 28:246–257, 2000.

Genomic organization of the human α-adducin gene and its alternately spliced isoforms

Abstract The cDNA for the human α-adducin gene has been cloned, and different alternately spliced forms have been identified. We report the complete genomic organization of the human α-adducin gene and these alternately spliced forms. The human α-adducin gene, spanning approximately 85 kb, consists of 16 exons ranging in size from 34 to 1892 bp. One of the spliced forms of the human α-adducin gene results from alternate use of the 5′ splice donor site for exon 10, while another results in a truncated protein following insertion of 34 bp comprising exon 15, followed by a premature stop codon. This alternate spliced form of α-adducin is predicted to result in an altered carboxyl terminus that would eliminate a protein kinase and calmodulin binding site. Seven nucleotide substitutions and 4 insertion/deletions were also identified. The 5′ region of the human α-adducin gene contains one Sp1 site, two AP2 sites, and two CAAT boxes. No TATA box was apparent, consistent with features of a housekeeping gene. We have mapped another cDNA within the first intron of the human α-adducin gene, suggesting overlapping genes in this 4p16.3 genomic region.

Lack of expression of XIST from a small ring X chromosome containing the XIST locus in a girl with short stature, facial dysmorphism and developmental delay

A 46,X,r(X) karyotype was found in a three and a half year old girl with short stature, facial dysmorphism and developmental delay. The clinical findings were consistent with the phenotype described in a limited number of patients with small ring X chromosomes lacking the XIST locus, a critical player in the process of X chromosome inactivation. Surprisingly, in our patient, fluorescent in situ hybridisation demonstrated that the XIST locus was present on the ring X. However, expression studies showed that there was no XIST transcript in peripheral blood cells, suggesting that the ring X had not been inactivated. This was confirmed by the demonstration that both of the patients alleles for the androgen receptor gene were unmethylated, and that both of the patients ZXDA alleles were expressed. The active nature of the ring X would presumably result in overexpression of genes that may account for the developmental delay observed for the patient. Using polymorphic markers along the X chromosome, the ring X was determined to be of paternal origin with one breakpoint in the long arm between DXS8037 and XIST and one in the short arm in Xp11.2 between DXS1126 and DXS991. To attempt to determine why the XIST gene failed to be expressed, the promoter region was sequenced and found to have a base change at the same location as a variant previously associated with nonrandom X chromosome inactivation. This mutation was not seen in over one hundred normal X chromosomes examined; however, it was observed in the paternal grandmother who did not show substantial skewing of X chromosome inactivation.