Janette Edson

A comparative study of techniques for differential expression analysis on RNA-Seq data.

Recent advances in next-generation sequencing technology allow high-throughput cDNA sequencing (RNA-Seq) to be widely applied in transcriptomic studies, in particular for detecting differentially expressed genes between groups. Many software packages have been developed for the identification of differentially expressed genes (DEGs) between treatment groups based on RNA-Seq data. However, there is a lack of consensus on how to approach an optimal study design and choice of suitable software for the analysis. In this comparative study we evaluate the performance of three of the most frequently used software tools: Cufflinks-Cuffdiff2, DESeq and edgeR. A number of important parameters of RNA-Seq technology were taken into consideration, including the number of replicates, sequencing depth, and balanced vs. unbalanced sequencing depth within and between groups. We benchmarked results relative to sets of DEGs identified through either quantitative RT-PCR or microarray. We observed that edgeR performs slightly better than DESeq and Cuffdiff2 in terms of the ability to uncover true positives. Overall, DESeq or taking the intersection of DEGs from two or more tools is recommended if the number of false positives is a major concern in the study. In other circumstances, edgeR is slightly preferable for differential expression analysis at the expense of potentially introducing more false positives.

Experience-Dependent Accumulation of N6-Methyladenosine in the Prefrontal Cortex Is Associated with Memory Processes in Mice.

The RNA modification N6-methyladenosine (m6A) influences mRNA stability and cell-type-specific developmental programming, and is highly abundant in the adult brain. However, it has not been determined whether m6A is dynamically regulated by experience. Based on transcriptome-wide profiling of m6A, we report that the level of m6A increases in the medial prefrontal cortex (mPFC) of mice in response to behavioral experience. The modulation was enriched near the stop codon of mRNAs, including genes related to neuronal plasticity. In primary cortical neurons, in vitro, modulation of m6A by the RNA demethylase FTO influenced the degradation profiles of a subset of transcripts with modulated sites. In vivo, the expression of Fto and the m6A methyltransferase, Mettl3 correlated with the observed increase in m6A levels post-training. Furthermore, targeted knockdown of FTO in the mPFC led to enhanced consolidation of cued fear memory. Thus, together with its role in early development, the dynamic regulation of m6A in the adult brain serves as an important epitranscriptomic mechanism associated with behavioral adaptation. SIGNIFICANCE STATEMENT N6-methyladenosine (m6A) is the most prevalent internal modification on RNA, however, its cellular dynamics in vivo remains elusive. Here we provide the first demonstration of m6A upregulation in the mouse medial prefrontal cortex (mPFC) following behavioral training. Knocking down the m6A demethylase FTO in the mPFC, which increases total m6A level, results in enhanced consolidation of fear memory. Our findings suggest that m6A is regulated in an activity-dependent manner in the adult brain, and may function to fine-tune mRNA turnover during memory-related processes.

Towards a molecular characterization of autism spectrum disorders: an exome sequencing and systems approach

The hypothetical ‘AXAS’ gene network model that profiles functional patterns of heterogeneous DNA variants overrepresented in autism spectrum disorder (ASD), X-linked intellectual disability, attention deficit and hyperactivity disorder and schizophrenia was used in this current study to analyze whole exome sequencing data from an Australian ASD cohort. An optimized DNA variant filtering pipeline was used to identify loss-of-function DNA variations. Inherited variants from parents with a broader autism phenotype and de novo variants were found to be significantly associated with ASD. Gene ontology analysis revealed that putative rare causal variants cluster in key neurobiological processes and are overrepresented in functions involving neuronal development, signal transduction and synapse development including the neurexin trans-synaptic complex. We also show how a complex gene network model can be used to fine map combinations of inherited and de novo variations in families with ASD that converge in the L1CAM pathway. Our results provide an important step forward in the molecular characterization of ASD with potential for developing a tool to analyze the pathogenesis of individual affected families.

Reduced reaction volumes and increased Taq DNA polymerase concentration improve STR profiling outcomes from a real-world low template DNA source: telogen hairs

PurposeThe primary method for analysis of low template DNA (LTDNA) is known as the low copy number (LCN) method involving an increased number of PCR cycles (typically 34). In common with other LTDNA methods, LCN profiles are characterized by allelic imbalance, drop in, and drop out that require complicated interpretation rules. They often require replicate PCR reactions to generate a “consensus” profile in a specialized facility. An ideal method for analysis of LTDNA should enhance profiling outcomes without elevated error rates and be performed using standard facilities, with minimum additional cost.MethodsIn this study, we present a comparison of four method variations for the amplification of STRs from LTDNA with a widely used, commercially available kit (AmpFℓSTR® Profiler Plus®): the standard method, the standard method with a post-PCR clean up, the LCN method, and a reduced reaction volume with increased Taq DNA polymerase concentration.ResultsUsing telogen hairs—a common source of LTDNA—and matched reference DNA, the LCN method produced the highest number of concordant and non-concordant (i.e., dropped-in) alleles. In comparison, the reduced reaction volume with increased Taq polymerase yielded more full and concordant DNA profiles (all alleles combined) and less off-ladder alleles from a broad range of input DNA. In addition, this method resulted in less non-concordant alleles than LCN and no more than for standard PCR, which suggests that it may be preferred over increased PCR cycles for LTDNA analysis, either with or without consensus profiling and statistical modelling.ConclusionsOverall, this study highlights the importance and benefit of optimizing PCR conditions and developing improved laboratory methods to amplify and analyze LTDNA.

Extracting DNA from ‘jaws’: high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so‐called bio‐swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high‐yield sources of DNA for genomic‐scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.

Whole exome sequencing and DNA methylation analysis in a clinical amyotrophic lateral sclerosis cohort

Gene discovery has provided remarkable biological insights into amyotrophic lateral sclerosis (ALS). One challenge for clinical application of genetic testing is critical evaluation of the significance of reported variants.

Whole-exome sequencing in amyotrophic lateral sclerosis suggests NEK1 is a risk gene in Chinese

BackgroundAmyotrophic lateral sclerosis (ALS) is a progressive neurological disease characterised by the degeneration of motor neurons, which are responsible for voluntary movement. There remains limited understanding of disease aetiology, with median survival of ALS of three years and no effective treatment. Identifying genes that contribute to ALS susceptibility is an important step towards understanding aetiology. The vast majority of published human genetic studies, including for ALS, have used samples of European ancestry. The importance of trans-ethnic studies in human genetic studies is widely recognised, yet a dearth of studies of non-European ancestries remains. Here, we report analyses of novel whole-exome sequencing (WES) data from Chinese ALS and control individuals.MethodsWES data were generated for 610 ALS cases and 460 controls drawn from Chinese populations. We assessed evidence for an excess of rare damaging mutations at the gene level and the gene set level, considering only singleton variants filtered to have allele frequency less than 5 × 10–5 in reference databases. To meta-analyse our results with a published study of European ancestry, we used a Cochran–Mantel–Haenszel test to compare gene-level variant counts in cases vs controls.ResultsNo gene passed the genome-wide significance threshold with ALS in Chinese samples alone. Combining rare variant counts in Chinese with those from the largest WES study of European ancestry resulted in three genes surpassing genome-wide significance: TBK1 (p = 8.3 × 10–12), SOD1 (p = 8.9 × 10–9) and NEK1 (p = 1.1 × 10–9). In the Chinese data alone, SOD1 and NEK1 were nominally significantly associated with ALS (p = 0.04 and p = 7 × 10–3, respectively) and the case/control frequencies of rare coding variants in these genes were similar in Chinese and Europeans (SOD1: 1.5%/0.2% vs 0.9%/0.1%, NEK1 1.8%/0.4% vs 1.9%/0.8%). This was also true for TBK1 (1.2%/0.2% vs 1.4%/0.4%), but the association with ALS in Chinese was not significant (p = 0.14).ConclusionsWhile SOD1 is already recognised as an ALS-associated gene in Chinese, we provide novel evidence for association of NEK1 with ALS in Chinese, reporting variants in these genes not previously found in Europeans.

An integrative analysis of non-coding regulatory DNA variations associated with autism spectrum disorder

A number of genetic studies have identified rare protein-coding DNA variations associated with autism spectrum disorder (ASD), a neurodevelopmental disorder with significant genetic etiology and heterogeneity. In contrast, the contributions of functional, regulatory genetic variations that occur in the extensive non-protein-coding regions of the genome remain poorly understood. Here we developed a genome-wide analysis to identify the rare single nucleotide variants (SNVs) that occur in non-coding regions and determined the regulatory function and evolutionary conservation of these variants. Using publicly available datasets and computational predictions, we identified SNVs within putative regulatory regions in promoters, transcription factor binding sites, and microRNA genes and their target sites. Overall, we found that the regulatory variants in ASD cases were enriched in ASD-risk genes and genes involved in fetal neurodevelopment. As with previously reported coding mutations, we found an enrichment of the regulatory variants associated with dysregulation of neurodevelopmental and synaptic signaling pathways. Among these were several rare inherited SNVs found in the mature sequence of microRNAs predicted to affect the regulation of ASD-risk genes. We show a paternally inherited miR-873-5p variant with altered binding affinity for several risk-genes including NRXN2 and CNTNAP2 putatively overlay maternally inherited loss-of-function coding variations in NRXN1 and CNTNAP2 to likely increase the genetic liability in an idiopathic ASD case. Our analysis pipeline provides a new resource for identifying loss-of-function regulatory DNA variations that may contribute to the genetic etiology of complex disorders.

Small non-coding RNA expression from anterior cingulate cortex in schizophrenia shows sex specific regulation

MicroRNAs (miRNAs) are known to regulate the expression of genes that are important for brain development and function, but the roles of other classes of small non-coding RNAs (sncRNAs) are less well understood. Additionally, although miRNA expression studies have been conducted in post-mortem brain samples from schizophrenia (SCZ) patients, other classes of sncRNAs are yet to be investigated in SCZ. We profiled the expression of miRNAs, piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs) in SCZ by applying small RNA sequencing (RNA-Seq) to sncRNA isolated from post-mortem anterior cingulate cortex (ACC) of SCZ-affected individuals (n=22) and matched controls (n=22). We identified about one-third of annotated miRNAs, one-quarter of snoRNAs and a small proportion of piRNAs and snRNAs. No sncRNAs were significantly differentially expressed between SCZ and controls, but there was evidence for an interaction between disease status and sex on the expression level of a number of miRNAs and snoRNAs. Many of these transcripts exhibited differential expression between male and female cases, and/or between female cases and controls, suggesting sex based dysregulation in ACC of SCZ. These findings require replication in an independent sample, but our study provides further insights into the potential involvement of sncRNAs in brain function and SCZ.