Multiplex accurate sensitive quantitation (MASQ) with application to minimal residual disease in acute myeloid leukemia.

Measuring minimal residual disease in cancer has applications for prognosis, monitoring treatment and detection of recurrence. Simple sequence-based methods to detect nucleotide substitution variants have error rates (about 10-3) that limit sensitive detection. We developed and characterized the performance of MASQ (multiplex accurate sensitive quantitation), a method with an error rate below 10-6. MASQ counts variant templates accurately in the presence of millions of host genomes by using tags to identify each template and demanding consensus over multiple reads. Since the MASQ protocol multiplexes 50 target loci, we can both integrate signal from multiple variants and capture subclonal response to treatment. Compared to existing methods for variant detection, MASQ achieves an excellent combination of sensitivity, specificity and yield. We tested MASQ in a pilot study in acute myeloid leukemia (AML) patients who entered complete remission. We detect leukemic variants in the blood and bone marrow samples of all five patients, after induction therapy, at levels ranging from 10-2 to nearly 10-6. We observe evidence of sub-clonal structure and find higher target variant frequencies in patients who go on to relapse, demonstrating the potential for MASQ to quantify residual disease in AML.

SMASH, a fragmentation and sequencing method for genomic copy number analysis.

Copy number variants (CNVs) underlie a significant amount of genetic diversity and disease. CNVs can be detected by a number of means, including chromosomal microarray analysis (CMA) and whole-genome sequencing (WGS), but these approaches suffer from either limited resolution (CMA) or are highly expensive for routine screening (both CMA and WGS). As an alternative, we have developed a next-generation sequencing-based method for CNV analysis termed SMASH, for short multiply aggregated sequence homologies. SMASH utilizes random fragmentation of input genomic DNA to create chimeric sequence reads, from which multiple mappable tags can be parsed using maximal almost-unique matches (MAMs). The SMASH tags are then binned and segmented, generating a profile of genomic copy number at the desired resolution. Because fewer reads are necessary relative to WGS to give accurate CNV data, SMASH libraries can be highly multiplexed, allowing large numbers of individuals to be analyzed at low cost. Increased genomic resolution can be achieved by sequencing to higher depth.

The contribution of de novo coding mutations to autism spectrum disorder.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

De novo gene disruptions in children on the autistic spectrum.

Exome sequencing of 343 families, each with a single child on the autism spectrum and at least one unaffected sibling, reveal de novo small indels and point substitutions, which come mostly from the paternal line in an age-dependent manner. We do not see significantly greater numbers of de novo missense mutations in affected versus unaffected children, but gene-disrupting mutations (nonsense, splice site, and frame shifts) are twice as frequent, 59 to 28. Based on this differential and the number of recurrent and total targets of gene disruption found in our and similar studies, we estimate between 350 and 400 autism susceptibility genes. Many of the disrupted genes in these studies are associated with the fragile X protein, FMRP, reinforcing links between autism and synaptic plasticity. We find FMRP-associated genes are under greater purifying selection than the remainder of genes and suggest they are especially dosage-sensitive targets of cognitive disorders.

An oncogenomics-based in vivo RNAi screen identifies tumor suppressors in liver cancer.

Cancers are highly heterogeneous and contain many passenger and driver mutations. To functionally identify tumor suppressor genes relevant to human cancer, we compiled pools of short hairpin RNAs (shRNAs) targeting the mouse orthologs of genes recurrently deleted in a series of human hepatocellular carcinomas and tested their ability to promote tumorigenesis in a mosaic mouse model. In contrast to randomly selected shRNA pools, many deletion-specific pools accelerated hepatocarcinogenesis in mice. Through further analysis, we identified and validated 13 tumor suppressor genes, 12 of which had not been linked to cancer before. One gene, XPO4, encodes a nuclear export protein whose substrate, EIF5A2, is amplified in human tumors, is required for proliferation of XPO4-deficient tumor cells, and promotes hepatocellular carcinoma in mice. Our results establish the feasibility of in vivo RNAi screens and illustrate how combining cancer genomics, RNA interference, and mosaic mouse models can facilitate the functional annotation of the cancer genome.

Topoisomerase levels determine chemotherapy response in vitro and in vivo.

Topoisomerase poisons are chemotherapeutic agents that are used extensively for treating human malignancies. These drugs can be highly effective, yet tumors are frequently refractory to treatment or become resistant upon tumor relapse. Using a pool-based RNAi screening approach and a well characterized mouse model of lymphoma, we explored the genetic basis for heterogeneous responses to topoisomerase poisons in vitro and in vivo. These experiments identified Top2A expression levels as major determinants of response to the topoisomerase 2 poison doxorubicin and showed that suppression of Top2A produces resistance to doxorubicin in vitro and in vivo. Analogously, using a targeted RNAi approach, we demonstrated that suppression of Top1 produces resistance to the topoisomerase 1 poison camptothecin yet hypersensitizes cancer cells to doxorubicin. Importantly, lymphomas relapsing after treatment display spontaneous changes in topoisomerase levels as predicted by in vitro gene knockdown studies. These results highlight the utility of pooled shRNA screens for identifying genetic determinants of chemotherapy response and suggest strategies for improving the effectiveness of topoisomerase poisons in the clinic.

A role for Yin Yang-1 (YY1) in the assembly of snRNA transcription complexes.

The RNA polymerase (pol) II and III human small nuclear RNA (snRNA) genes have very similar promoters and recruit a number of common factors. In particular, both types of promoters utilize the small nuclear RNA activating protein complex (SNAP(c)) and the TATA box binding protein (TBP) for basal transcription, and are activated by Oct-1. We find that SNAP(c) purified from cell lines expressing tagged SNAP(c) subunits is associated with Yin Yang-1 (YY1), a factor implicated in both activation and repression of transcription. Recombinant YY1 accelerates the binding of SNAP(c) to the proximal sequence element, its target within snRNA promoters. Moreover, it enhances the formation of a complex on the pol III U6 snRNA promoter containing all the factors (SNAP(c), TBP, TFIIB-related factor 2 (Brf2), and B double prime 1 (Bdp1)) that are sufficient to direct in vitro U6 transcription when complemented with purified pol III, as well as that of a subcomplex containing TBP, Brf2, and Bdp1. YY1 is found on both the RNA polymerase II U1 and the RNA polymerase III U6 promoters as determined by chromatin immunoprecipitations. Thus, YY1 represents a new factor that participates in transcription complexes formed on both pol II and III promoters.

Structure-function analysis of the human TFIIB-related factor II protein reveals an essential role for the C-terminal domain in RNA polymerase III transcription.

The transcription factors TFIIB, Brf1, and Brf2 share related N-terminal zinc ribbon and core domains. TFIIB bridges RNA polymerase II (Pol II) with the promoter-bound preinitiation complex, whereas Brf1 and Brf2 are involved, as part of activities also containing TBP and Bdp1 and referred to here as Brf1-TFIIIB and Brf2-TFIIIB, in the recruitment of Pol III. Brf1-TFIIIB recruits Pol III to type 1 and 2 promoters and Brf2-TFIIIB to type 3 promoters such as the human U6 promoter. Brf1 and Brf2 both have a C-terminal extension absent in TFIIB, but their C-terminal extensions are unrelated. In yeast Brf1, the C-terminal extension interacts with the TBP/TATA box complex and contributes to the recruitment of Bdp1. Here we have tested truncated Brf2, as well as Brf2/TFIIB chimeric proteins for U6 transcription and for assembly of U6 preinitiation complexes. Our results characterize functions of various human Brf2 domains and reveal that the C-terminal domain is required for efficient association of the protein with U6 promoter-bound TBP and SNAP(c), a type 3 promoter-specific transcription factor, and for efficient recruitment of Bdp1. This in turn suggests that the C-terminal extensions in Brf1 and Brf2 are crucial to specific recruitment of Pol III over Pol II.

Redundant cooperative interactions for assembly of a human U6 transcription initiation complex.

The core human U6 promoter consists of a proximal sequence element (PSE) located upstream of a TATA box. The PSE is recognized by the snRNA-activating protein complex (SNAP(c)), which consists of five types of subunits, SNAP190, SNAP50, SNAP45, SNAP43, and SNAP19. The TATA box is recognized by TATA box binding protein (TBP). In addition, basal U6 transcription requires the SANT domain protein Bdp1 and the transcription factor IIB-related factor Brf2. SNAP(c) and mini-SNAP(c), which consists of just SNAP43, SNAP50, and the N-terminal third of SNAP190, bind cooperatively with TBP to the core U6 promoter. By generating complexes smaller than mini-SNAP(c), we have identified a 50-amino-acid region within SNAP190 that is (i) required for cooperative binding with TBP in the context of mini-SNAP(c) and (ii) sufficient for cooperative binding with TBP when fused to a heterologous DNA binding domain. We show that derivatives of mini-SNAP(c) lacking this region are active for transcription and that with such complexes, TBP can still be recruited to the U6 promoter through cooperative interactions with Brf2. Our results identify complexes smaller than mini-SNAP(c) that are transcriptionally active and show that there are at least two redundant mechanisms to stably recruit TBP to the U6 transcription initiation complex.