Australian Genomics: Outcomes of a 5-year national program to accelerate the integration of genomics in healthcare.

Australian Genomics is a national collaborative partnership of more than 100 organizations piloting a whole-of-system approach to integrating genomics into healthcare, based on federation principles. In the first five years of operation, Australian Genomics has evaluated the outcomes of genomic testing in more than 5,200 individuals across 19 rare disease and cancer flagship studies. Comprehensive analyses of the health economic, policy, ethical, legal, implementation and workforce implications of incorporating genomics in the Australian context have informed evidence-based change in policy and practice, resulting in national government funding and equity of access for a range of genomic tests. Simultaneously, Australian Genomics has built national skills, infrastructure, policy, and data resources to enable effective data sharing to drive discovery research and support improvements in clinical genomic delivery.

MTR3D: identifying regions within protein tertiary structures under purifying selection.

The identification of disease-causal variants is non-trivial. By mapping population variation from over 448,000 exome and genome sequences to over 81,000 experimental structures and homology models of the human proteome, we have calculated both regional intolerance to missense variation (Missense Tolerance Ratio, MTR), using a sliding window of 21-41 codons, and introduce a new 3D spatial intolerance to missense variation score (3D Missense Tolerance Ratio, MTR3D), using spheres of 5-8 Å. We show that the MTR3D is less biased by regions with limited data and more accurately identifies regions under purifying selection than estimates relying on the sequence alone. Intolerant regions were highly enriched for both ClinVar pathogenic and COSMIC somatic missense variants (Mann-Whitney U test P < 2.2 × 10-16). Further, we combine sequence- and spatial-based scores to generate a consensus score, MTRX, which distinguishes pathogenic from benign variants more accurately than either score separately (AUC = 0.85). The MTR3D server enables easy visualisation of population variation, MTR, MTR3D and MTRX scores across the entire gene and protein structure for >17,000 human genes and >42,000 alternative alternate transcripts, including both Ensembl and RefSeq transcripts. MTR3D is freely available by user-friendly web-interface and API at http://biosig.unimelb.edu.au/mtr3d/.

aCLImatise: automated generation of tool definitions for bioinformatics workflows.

SUMMARY: aCLImatise is a utility for automatically generating tool definitions compatible with bioinformatics workflow languages, by parsing command-line help output. aCLImatise also has an associated database called the aCLImatise Base Camp, which provides thousands of pre-computed tool definitions. AVAILABILITY AND IMPLEMENTATION: The latest aCLImatise source code is available within a GitHub organisation, under the GPL-3.0 license: https://github.com/aCLImatise. In particular, documentation for the aCLImatise Python package is available at https://aclimatise.github.io/CliHelpParser/, and the aCLImatise Base Camp is available at https://aclimatise.github.io/BaseCamp/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A head-to-head evaluation of the diagnostic efficacy and costs of trio versus singleton exome sequencing analysis.

Diagnostic exome sequencing (ES) can be performed on the proband only (singleton; sES) or with additional samples, often including both biological parents with the proband (trio; tES). In this study we sought to compare the efficiencies of exome sequencing (ES) by trio (tES) versus singleton (sES) approach, determine costs, and identify factors to consider when deciding on optimal implementation strategies for the diagnosis of monogenic disorders. We undertook ES in 30 trios and analysed each proband's sES and tES data in parallel. Two teams were randomly allocated to either sES or tES analysis for each case and blinded to each other's work. Each task was timed and cost analyses were based on time taken and diagnostic yield. We modelled three scenarios to determine the factors to consider in the implementation of tES. sES diagnosed 11/30 (36.7%) cases and tES identified one additional diagnosis (12/30 (40.0%)). tES obviated the need for Sanger segregation, reduced the number of variants for curation, and had lower cost-per-diagnosis when considering analysis alone. When sequencing costs were included, tES nearly doubled the cost of sES. Reflexing to tES in those who remain undiagnosed after sES was cost-saving over tES in all as first-line. This approach requires a large differential in diagnostic yield between sES and tES for maximal benefit given current sequencing costs. tES may be preferable when scaling up laboratory throughput due to efficiency gains and opportunity cost considerations. Our findings are relevant to clinicians, laboratories and health services considering tES over sES.

Exome sequencing has higher diagnostic yield compared to simulated disease-specific panels in children with suspected monogenic disorders.

As test costs decline, whole-exome sequencing (WES) has become increasingly used for clinical diagnosis, and now represents the primary alternative to gene panel testing for patients with a suspected genetic disorder. We sought to compare the diagnostic yield of singleton-WES with simulated application of commercial gene panels in children suspected of having a genetically heterogeneous condition. Recruitment, singleton-WES and phenotype-driven variant analysis was completed for 145 paediatric patients. At recruitment, clinicians were required to propose commercial gene panel tests as an alternative to WES and nominate a phenotype-driven candidate gene list. In WES-diagnosed children, three commercial options for each proposed panel were identified and evaluated for hypothetical diagnostic yield assuming 100% analytical sensitivity and specificity. We compared the price of WES with the least costly panel in WES-diagnosed children. In WES-undiagnosed children, we evaluated the exonic coverage of their phenotype-driven gene list using aggregate data. WES diagnoses were made in genes not included in at least one-of-three commercial panels in 42% of cases. Had a panel been selected instead, 23% of WES-diagnosed children would not have been diagnosed. In 26% of cases, the least costly panel option would have been more expensive than WES. Evaluation of WES coverage found that at the most stringent level of 20× coverage, the likelihood of missing a clinically relevant variant in a candidate gene list was maximally 8%. The broader coverage of WES makes it a superior alternative to gene panel testing at similar financial cost for children with suspected complex monogenic phenotypes.

Evaluation of computational programs to predict HLA genotypes from genomic sequencing data.

Motivation: Despite being essential for numerous clinical and research applications, high-resolution human leukocyte antigen (HLA) typing remains challenging and laboratory tests are also time-consuming and labour intensive. With next-generation sequencing data becoming widely accessible, on-demand in silico HLA typing offers an economical and efficient alternative. Results: In this study we evaluate the HLA typing accuracy and efficiency of five computational HLA typing methods by comparing their predictions against a curated set of > 1000 published polymerase chain reaction-derived HLA genotypes on three different data sets (whole genome sequencing, whole exome sequencing and transcriptomic sequencing data). The highest accuracy at clinically relevant resolution (four digits) we observe is 81% on RNAseq data by PHLAT and 99% accuracy by OptiType when limited to Class I genes only. We also observed variability between the tools for resource consumption, with runtime ranging from an average of 5 h (HLAminer) to 7 min (seq2HLA) and memory from 12.8 GB (HLA-VBSeq) to 0.46 GB (HLAminer) per sample. While a minimal coverage is required, other factors also determine prediction accuracy and the results between tools do not correlate well. Therefore, by combining tools, there is the potential to develop a highly accurate ensemble method that is able to deliver fast, economical HLA typing from existing sequencing data.

A clinically driven variant prioritization framework outperforms purely computational approaches for the diagnostic analysis of singleton WES data.

Rapid identification of clinically significant variants is key to the successful application of next generation sequencing technologies in clinical practice. The Melbourne Genomics Health Alliance (MGHA) variant prioritization framework employs a gene prioritization index based on clinician-generated a priori gene lists, and a variant prioritization index (VPI) based on rarity, conservation and protein effect. We used data from 80 patients who underwent singleton whole exome sequencing (WES) to test the ability of the framework to rank causative variants highly, and compared it against the performance of other gene and variant prioritization tools. Causative variants were identified in 59 of the patients. Using the MGHA prioritization framework the average rank of the causative variant was 2.24, with 76% ranked as the top priority variant, and 90% ranked within the top five. Using clinician-generated gene lists resulted in ranking causative variants an average of 8.2 positions higher than prioritization based on variant properties alone. This clinically driven prioritization approach significantly outperformed purely computational tools, placing a greater proportion of causative variants top or in the top 5 (permutation P-value=0.001). Clinicians included 40 of the 49 WES diagnoses in their a priori list of differential diagnoses (81%). The lists generated by PhenoTips and Phenomizer contained 14 (29%) and 18 (37%) of these diagnoses respectively. These results highlight the benefits of clinically led variant prioritization in increasing the efficiency of singleton WES data analysis and have important implications for developing models for the funding and delivery of genomic services.

Diagnostic Impact and Cost-effectiveness of Whole-Exome Sequencing for Ambulant Children With Suspected Monogenic Conditions.

Importance: Optimal use of whole-exome sequencing (WES) in the pediatric setting requires an understanding of who should be considered for testing and when it should be performed to maximize clinical utility and cost-effectiveness. Objectives: To investigate the impact of WES in sequencing-naive children suspected of having a monogenic disorder and evaluate its cost-effectiveness if WES had been available at different time points in their diagnostic trajectory. Design, Setting, and Participants: This prospective study was part of the Melbourne Genomics Health Alliance demonstration project. At the ambulatory outpatient clinics of the Victorian Clinical Genetics Services at the Royal Children's Hospital, Melbourne, Australia, children older than 2 years suspected of having a monogenic disorder were prospectively recruited from May 1 through November 30, 2015, by clinical geneticists after referral from general and subspecialist pediatricians. All children had nondiagnostic microarrays and no prior single-gene or panel sequencing. Exposures: All children underwent singleton WES with targeted phenotype-driven analysis. Main Outcomes and Measures: The study examined the clinical utility of a molecular diagnosis and the cost-effectiveness of alternative diagnostic trajectories, depending on timing of WES. Results: Of 61 children originally assessed, 44 (21 [48%] male and 23 [52%] female) aged 2 to 18 years (mean age at initial presentation, 28 months; range, 0-121 months) were recruited, and a diagnosis was achieved in 23 (52%) by singleton WES. The diagnoses were unexpected in 8 of 23 (35%), and clinical management was altered in 6 of 23 (26%). The mean duration of the diagnostic odyssey was 6 years, with each child having a mean of 19 tests and 4 clinical genetics and 4 nongenetics specialist consultations, and 26 (59%) underwent a procedure while under general anesthetic for diagnostic purposes. Economic analyses of the diagnostic trajectory identified that WES performed at initial tertiary presentation resulted in an incremental cost savings of A$9020 (US$6838) per additional diagnosis (95% CI, A$4304-A$15 404 [US$3263-US$11 678]) compared with the standard diagnostic pathway. Even if WES were performed at the first genetics appointment, there would be an incremental cost savings of A$5461 (US$4140) (95% CI, A$1433-A$10 557 [US$1086- US$8004]) per additional diagnosis compared with the standard diagnostic pathway. Conclusions and Relevance: Singleton WES in children with suspected monogenic conditions has high diagnostic yield, and cost-effectiveness is maximized by early application in the diagnostic pathway. Pediatricians should consider early referral of children with undiagnosed syndromes to clinical geneticists.

Diagnostic and cost utility of whole exome sequencing in peripheral neuropathy.

OBJECTIVE: To explore the diagnostic utility and cost effectiveness of whole exome sequencing (WES) in a cohort of individuals with peripheral neuropathy. METHODS: Singleton WES was performed in individuals recruited though one pediatric and one adult tertiary center between February 2014 and December 2015. Initial analysis was restricted to a virtual panel of 55 genes associated with peripheral neuropathies. Patients with uninformative results underwent expanded analysis of the WES data. Data on the cost of prior investigations and assessments performed for diagnostic purposes in each patient was collected. RESULTS: Fifty patients with a peripheral neuropathy were recruited (median age 18 years; range 2-68 years). The median time from initial presentation to study enrollment was 6 years 9 months (range 2 months-62 years), and the average cost of prior investigations and assessments for diagnostic purposes AU$4013 per patient. Eleven individuals received a diagnosis from the virtual panel. Eight individuals received a diagnosis following expanded analysis of the WES data, increasing the overall diagnostic yield to 38%. Two additional individuals were diagnosed with pathogenic copy number variants through SNP microarray. CONCLUSIONS: This study provides evidence that WES has a high diagnostic utility and is cost effective in patients with a peripheral neuropathy. Expanded analysis of WES data significantly improves the diagnostic yield in patients in whom a diagnosis is not found on the initial targeted analysis. This is primarily due to diagnosis of conditions caused by newly discovered genes and the resolution of complex and atypical phenotypes.

An Emerging Female Phenotype with Loss-of-Function Mutations in the Aristaless-Related Homeodomain Transcription Factor ARX.

The devastating clinical presentation of X-linked lissencephaly with abnormal genitalia (XLAG) is invariably caused by loss-of-function mutations in the Aristaless-related homeobox (ARX) gene. Mutations in this X-chromosome gene contribute to intellectual disability (ID) with co-morbidities including seizures and movement disorders such as dystonia in affected males. The detection of affected females with mutations in ARX is increasing. We present a family with multiple affected individuals, including two females. Two male siblings presenting with XLAG were deceased prior to full-term gestation or within the first few weeks of life. Of the two female siblings, one presented with behavioral disturbances, mild ID, a seizure disorder, and complete agenesis of the corpus callosum (ACC), similar to the mother's phenotype. A novel insertion mutation in Exon 2 of ARX was identified, c.982delCinsTTT predicted to cause a frameshift at p.(Q328Ffs* 37). Our finding is consistent with loss-of-function mutations in ARX causing XLAG in hemizygous males and extends the findings of ID and seizures in heterozygous females. We review the reported phenotypes of females with mutations in ARX and highlight the importance of screening ARX in male and female patients with ID, seizures, and in particular with complete ACC.

Real-world utility of whole exome sequencing with targeted gene analysis for focal epilepsy.

OBJECTIVE: Driven by advances in genomic technology and reduction in costs, next-generation sequencing (NGS) is venturing into routine clinical care. The 'real-world' clinical utility of NGS remains to be determined in focal epilepsies, which account for 60% of all epilepsies and for which the importance of genetic factors is just beginning to emerge. We investigated the diagnostic yield and management implications of whole exome sequencing (WES)-based screening of selected genes in the routine care of common focal epilepsies suspected to have a genetic basis. METHODS: We performed WES, followed by targeted analysis of 64 epilepsy genes, on 40 consecutive children and adults enrolled prospectively from routine clinical practice who had MRI-negative focal epilepsy and a family history of febrile seizures or any type of epilepsy in at least one first- or second-degree relative. Exclusion criteria were previous genetic testing, severe intellectual disability and benign focal epilepsies of childhood. RESULTS: 5/40 (12.5%) patients had a pathogenic or likely pathogenic variant, detected in SCN1A, DEPDC5, PCDH19, GABRG2 or NPRL2. Identifying a pathogenic SCN1A variant in a patient with drug-resistant epilepsy prompted to halt presurgical investigations due to concern of unfavorable post-surgical outcome. It also led in the same patient to discontinue long-standing carbamazepine therapy (a potentially aggravating drug in epilepsies due to SCN1A mutations), resulting in complete seizure control. Patients with pathogenic or likely pathogenic variants had a younger median age of seizure onset (range) compared to those without [18 months (8 months-18 years) vs 18 years (18 months-70 years), p=0.02]. SIGNIFICANCE: Our data demonstrate that WES with targeted gene analysis is an effective diagnostic tool for patients with common focal epilepsies in whom a genetic etiology is suspected. It can also influence clinical decision-making, including antiepileptic drug selection and consideration of epilepsy surgery, hence supporting its incorporation in the routine clinical care of this patient group.

A prospective evaluation of whole-exome sequencing as a first-tier molecular test in infants with suspected monogenic disorders.

PURPOSE: To prospectively evaluate the diagnostic and clinical utility of singleton whole-exome sequencing (WES) as a first-tier test in infants with suspected monogenic disease. METHODS: Singleton WES was performed as a first-tier sequencing test in infants recruited from a single pediatric tertiary center. This occurred in parallel with standard investigations, including single- or multigene panel sequencing when clinically indicated. The diagnosis rate, clinical utility, and impact on management of singleton WES were evaluated. RESULTS: Of 80 enrolled infants, 46 received a molecular genetic diagnosis through singleton WES (57.5%) compared with 11 (13.75%) who underwent standard investigations in the same patient group. Clinical management changed following exome diagnosis in 15 of 46 diagnosed participants (32.6%). Twelve relatives received a genetic diagnosis following cascade testing, and 28 couples were identified as being at high risk of recurrence in future pregnancies. CONCLUSIONS: This prospective study provides strong evidence for increased diagnostic and clinical utility of singleton WES as a first-tier sequencing test for infants with a suspected monogenic disorder. Singleton WES outperformed standard care in terms of diagnosis rate and the benefits of a diagnosis, namely, impact on management of the child and clarification of reproductive risks for the extended family in a timely manner.Genet Med 18 11, 1090-1096.

Personalised Medicine Possible With Real-Time Integration of Genomic and Clinical Data To Inform Clinical Decision-Making.

Despite widespread use of genomic sequencing in research, there are gaps in our understanding of the performance and provision of genomic sequencing in clinical practice. The Melbourne Genomics Health Alliance (the Alliance), has been established to determine the feasibility, performance and impact of using genomic sequencing as a diagnostic tool. The Alliance has partnered with BioGrid Australia to enable the linkage of genomic sequencing, clinical treatment and outcome data for this project. This integrated dataset of genetic, clinical and patient sourced information will be used by the Alliance to evaluate the potential diagnostic value of genomic sequencing in routine clinical practice. This project will allow the Alliance to provide recommendations to facilitate the integration of genomic sequencing into clinical practice to enable personalised disease treatment.

Cpipe: a shared variant detection pipeline designed for diagnostic settings.

The benefits of implementing high throughput sequencing in the clinic are quickly becoming apparent. However, few freely available bioinformatics pipelines have been built from the ground up with clinical genomics in mind. Here we present Cpipe, a pipeline designed specifically for clinical genetic disease diagnostics. Cpipe was developed by the Melbourne Genomics Health Alliance, an Australian initiative to promote common approaches to genomics across healthcare institutions. As such, Cpipe has been designed to provide fast, effective and reproducible analysis, while also being highly flexible and customisable to meet the individual needs of diverse clinical settings. Cpipe is being shared with the clinical sequencing community as an open source project and is available at http://cpipeline.org.

High-resolution characterization of sequence signatures due to non-random cleavage of cell-free DNA.

BACKGROUND: High-throughput sequencing of cell-free DNA fragments found in human plasma has been used to non-invasively detect fetal aneuploidy, monitor organ transplants and investigate tumor DNA. However, many biological properties of this extracellular genetic material remain unknown. Research that further characterizes circulating DNA could substantially increase its diagnostic value by allowing the application of more sophisticated bioinformatics tools that lead to an improved signal to noise ratio in the sequencing data. METHODS: In this study, we investigate various features of cell-free DNA in plasma using deep-sequencing data from two pregnant women (>70X, >50X) and compare them with matched cellular DNA. We utilize a descriptive approach to examine how the biological cleavage of cell-free DNA affects different sequence signatures such as fragment lengths, sequence motifs at fragment ends and the distribution of cleavage sites along the genome. RESULTS: We show that the size distributions of these cell-free DNA molecules are dependent on their autosomal and mitochondrial origin as well as the genomic location within chromosomes. DNA mapping to particular microsatellites and alpha repeat elements display unique size signatures. We show how cell-free fragments occur in clusters along the genome, localizing to nucleosomal arrays and are preferentially cleaved at linker regions by correlating the mapping locations of these fragments with ENCODE annotation of chromatin organization. Our work further demonstrates that cell-free autosomal DNA cleavage is sequence dependent. The region spanning up to 10 positions on either side of the DNA cleavage site show a consistent pattern of preference for specific nucleotides. This sequence motif is present in cleavage sites localized to nucleosomal cores and linker regions but is absent in nucleosome-free mitochondrial DNA. CONCLUSIONS: These background signals in cell-free DNA sequencing data stem from the non-random biological cleavage of these fragments. This sequence structure can be harnessed to improve bioinformatics algorithms, in particular for CNV and structural variant detection. Descriptive measures for cell-free DNA features developed here could also be used in biomarker analysis to monitor the changes that occur during different pathological conditions.

Harnessing gene expression networks to prioritize candidate epileptic encephalopathy genes.

We apply a novel gene expression network analysis to a cohort of 182 recently reported candidate Epileptic Encephalopathy genes to identify those most likely to be true Epileptic Encephalopathy genes. These candidate genes were identified as having single variants of likely pathogenic significance discovered in a large-scale massively parallel sequencing study. Candidate Epileptic Encephalopathy genes were prioritized according to their co-expression with 29 known Epileptic Encephalopathy genes. We utilized developing brain and adult brain gene expression data from the Allen Human Brain Atlas (AHBA) and compared this to data from Celsius: a large, heterogeneous gene expression data warehouse. We show replicable prioritization results using these three independent gene expression resources, two of which are brain-specific, with small sample size, and the third derived from a heterogeneous collection of tissues with large sample size. Of the nineteen genes that we predicted with the highest likelihood to be true Epileptic Encephalopathy genes, two (GNAO1 and GRIN2B) have recently been independently reported and confirmed. We compare our results to those produced by an established in silico prioritization approach called Endeavour, and finally present gene expression networks for the known and candidate Epileptic Encephalopathy genes. This highlights sub-networks of gene expression, particularly in the network derived from the adult AHBA gene expression dataset. These networks give clues to the likely biological interactions between Epileptic Encephalopathy genes, potentially highlighting underlying mechanisms and avenues for therapeutic targets.

Use of copy number deletion polymorphisms to assess DNA chimerism.

BACKGROUND: We describe a novel approach that harnesses the ubiquity of copy number deletion polymorphisms in human genomes to definitively detect and quantify chimeric DNA in clinical samples. Unlike other molecular approaches to chimerism analysis, the copy number deletion (CND) method targets genomic loci (>50 base pairs in length) that are wholly absent from wild-type (i.e., self) background DNA sequences in a sex-independent manner. METHODS: Bespoke quantitative PCR (qPCR) CND assays were developed and validated using a series of DNA standards and chimeric plasma DNA samples collected from 2 allogeneic kidney transplant recipients and 12 pregnant women. Assay performance and informativeness were assessed using appropriate statistical methods. RESULTS: The CND qPCR assays showed high sensitivity, precision, and reliability for linear quantification of DNA chimerism down to 16 genomic equivalents (i.e., 106 pg). Fetal fraction (%) in 12 singleton male pregnancies was calculated using the CND qPCR approach, which showed closer agreement with single-nucleotide polymorphism-based massively parallel sequencing than the SRY (sex determining region Y) (Y chromosome) qPCR assay. The latter consistently underestimated the fetal fraction relative to the other methods. We also were able to measure biological changes in plasma nonself DNA concentrations in 2 renal transplant recipients. CONCLUSIONS: The CND qPCR technique is suitable for measurement of chimerism for monitoring of rejection in allogeneic organ transplantation and quantification of the cell-free fetal DNA fraction in maternal plasma samples used for noninvasive prenatal genetic testing.

Investigating and correcting plasma DNA sequencing coverage bias to enhance aneuploidy discovery.

Pregnant women carry a mixture of cell-free DNA fragments from self and fetus (non-self) in their circulation. In recent years multiple independent studies have demonstrated the ability to detect fetal trisomies such as trisomy 21, the cause of Down syndrome, by Next-Generation Sequencing of maternal plasma. The current clinical tests based on this approach show very high sensitivity and specificity, although as yet they have not become the standard diagnostic test. Here we describe improvements to the analysis of the sequencing data by reducing GC bias and better handling of the genomic repeats. We show substantial improvements in the sensitivity of the standard trisomy 21 statistical tests, which we measure by artificially reducing read coverage. We also explore the bias stemming from the natural cleavage of plasma DNA by examining DNA motifs and position specific base distributions. We propose a model to correct this fragmentation bias and observe that incorporating this bias does not lead to any further improvements in the detection of fetal trisomy. The improved bias corrections that we demonstrate in this work can be readily adopted into existing fetal trisomy detection protocols and should also lead to improvements in sub-chromosomal copy number variation detection.

The cost of reducing starting RNA quantity for Illumina BeadArrays: a bead-level dilution experiment.

BACKGROUND: The demands of microarray expression technologies for quantities of RNA place a limit on the questions they can address. As a consequence, the RNA requirements have reduced over time as technologies have improved. In this paper we investigate the costs of reducing the starting quantity of RNA for the Illumina BeadArray platform. This we do via a dilution data set generated from two reference RNA sources that have become the standard for investigations into microarray and sequencing technologies. RESULTS: We find that the starting quantity of RNA has an effect on observed intensities despite the fact that the quantity of cRNA being hybridized remains constant. We see a loss of sensitivity when using lower quantities of RNA, but no great rise in the false positive rate. Even with 10 ng of starting RNA, the positive results are reliable although many differentially expressed genes are missed. We see that there is some scope for combining data from samples that have contributed differing quantities of RNA, but note also that sample sizes should increase to compensate for the loss of signal-to-noise when using low quantities of starting RNA. CONCLUSIONS: The BeadArray platform maintains a low false discovery rate even when small amounts of starting RNA are used. In contrast, the sensitivity of the platform drops off noticeably over the same range. Thus, those conducting experiments should not opt for low quantities of starting RNA without consideration of the costs of doing so. The implications for experimental design, and the integration of data from different starting quantities, are complex.