Henrik Bengtsson

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ∼0.1× genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost.

Breast cancers are comprised of molecularly distinct subtypes that may respond differently to pathway-targeted therapies now under development. Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways, and drug response. In a test of 77 therapeutic compounds, nearly all drugs showed differential responses across these cell lines, and approximately one third showed subtype-, pathway-, and/or genomic aberration-specific responses. These observations suggest mechanisms of response and resistance and may inform efforts to develop molecular assays that predict clinical response.

Sex in mammals is genetically determined and is defined at the cellular level by sex chromosome complement (XY males and XX females). The Y chromosome-linked gene sex-determining region Y (SRY) is believed to be the master initiator of male sex determination in almost all eutherian and metatherian mammals, functioning to upregulate expression of its direct target gene Sry-related HMG box-containing gene 9 (SOX9). Data suggest that SRY evolved from SOX3, although there is no direct functional evidence to support this hypothesis. Indeed, loss-of-function mutations in SOX3 do not affect sex determination in mice or humans. To further investigate Sox3 function in vivo, we generated transgenic mice overexpressing Sox3. Here, we report that in one of these transgenic lines, Sox3 was ectopically expressed in the bipotential gonad and that this led to frequent complete XX male sex reversal. Further analysis indicated that Sox3 induced testis differentiation in this particular line of mice by upregulating expression of Sox9 via a similar mechanism to Sry. Importantly, we also identified genomic rearrangements within the SOX3 regulatory region in three patients with XX male sex reversal. Together, these data suggest that SOX3 and SRY are functionally interchangeable in sex determination and support the notion that SRY evolved from SOX3 via a regulatory mutation that led to its de novo expression in the early gonad.