The purpose of a Results section is to present, without interpretation, the key results of your research. Your paper does not need to include every result you obtained during your experiments. Results are “key” when they are relevant to addressing the research questions or hypotheses presented at the beginning of your paper.
We use the Results subsections to show the reader what types of outcomes they can expect when using the methodology that we present. In our papers, we write a “Methods Overview” as the first subsection of the Results section. (We discuss writing the “Methods Overview” subsection in a previous writing tips post.) Remaining subsections in your paper’s Results section present your findings in the form of text, figures, and tables.
Each Results subsection should make a specific point, and the subsection heading should be a succinct description of this message. Effective subsection headings declare a statement that communicates to the reader what the method is capable of doing or what types of data the method can be applied to. For example, in a recent paper published by our group, the heading of a subsection that demonstrates how a new GWAS approach controls for false positive results is: “Phenotype Imputation Controls Type 1 Error.”
Here, a two-paragraph Results subsection has a heading that tells the reader which specific type of analysis is discussed, since the paper presents a method that can be applied toward numerous different analytical tasks.
Cell type composition and diversity
We hypothesized that differences in microbial diversity may be linked to whole blood cell type composition. Since the actual cell counts were not available for these individuals, we used cell-proportion estimates derived from available DNA methylation data to test this hypothesis (Houseman et al. 2012; Aryee et al. 2014; Horvath and Levine 2015).
We assessed methylation data from 65 controls from our replication sample, and compared methylation-derived blood cell proportions to alpha diversity after adjusting for age, gender, RIN, and all technical parameters. We tested whether alpha diversity levels are associated to cell type abundance estimates. Our analysis shows one cell type, CD8+ CD28- CD45RA- cells, to be significantly negatively correlated with alpha diversity after correction for all other cell-count estimates (correlation = -0.41, P=7.3e-4, Figure S6, Table S6). These cells are T cells that lack CD8+ naïve cell markers CD28 and CD45RA and are thought to represent a subpopulation of differentiated CD8+ T cells (Koch et al. 2008; Horvath and Levine 2015). We observed that low alpha diversity correlates with high levels of this population of T cells cell abundance.
Total RNA Sequencing reveals microbial communities in human blood and disease specific effects
For each subsection, we include one figure that illustrates the heading’s message. The figure’s legend (also referred to as a “caption”) can simply be the subsection heading with additional information explaining the methods and data involved in the visual output. It may be helpful to select a figure and write a legend before composing text for the subsection.
At this point, you could probably write an entire paper on each figure! In general, we limit the text in each Results subsection to one to two paragraphs. Here, we use the minimum amount of text that is necessary to walk our reader through the figure. Think about what the reader needs to know in order to start using the method for their own analysis. Relevant information includes the type of data used, analytical steps and parameters, and a summary of conclusions. In many cases, the subsection text and figure legend will be repetitive.
This one-paragraph section provides relevant results in terms of statistical parameters, numerical output, and a supplemental figure. This subsection gives the reader a good idea of what to expect if they want to incorporate this new approach in their own project.
Phenotype Imputation Controls Type I Error
We simulated datasets for multiple phenotypes under the null model where the variant we are testing has no effect (effect size of zero) toward the target phenotype. We computed the type I error under five different significance thresholds: 0.05, 0.01, 0.005, 5 3 10-6, and 5 3 10-8. We generated 100,000,000 simulated datasets that consist of 1,000 individuals. The type I error rates for our imputation method were 0.049, 0.0099, 0.00489, 4.90 3 10-6, and 4.89 3 10-8 for the significance thresholds of 0.05, 0.01, 0.005, 5 3 10-6, and 5 3 10-8, respectively. This indicates that the type I error is correctly controlled in our imputation method. The Northern Finland Birth Cohort dataset 13 was used to show that the type I error is controlled (see Figure S1). We plot the Q-Q plot of the Z-score for the imputed triglyceride (TG) phenotype from the Finland dataset. There is no inflation in the Q-Qplot as shown in Figure S1.
Imputing Phenotypes for Genome-wide Association Studies
Bonus challenge: After you finish writing your paper, try to remove the sentence highlighting the result’s importance from the Figure caption.
The order in which you present your results can be organized in many different ways. Typically, ordering of subsections is not important for initial manuscripts. One simple approach is to order Results subsections sequentially to support the argument that you are building in your paper.
Here, we present another example of a Results subsection, including the description of a relevant figure. The subsection heading is making it clear to the reader that this part of the paper discusses applying ForestPMPlot, a visualization tool for analyzing meta-analysis studies, to eQTL data.
Application to multi-tissue eQTL analysis
One powerful application of our proposed framework is in multi-tissue eQTL analysis in the Genotype-Tissue Expression (GTEx) project. The GTEx project studies human gene expression and genetic regulation in multiple tissues, providing valuable insights into the mechanisms of gene regulation, which can lead to the new discovery of disease-related perturbations. In this project, genetic variation between individuals will be examined for correlation with differences in gene expression level to identify regions of the genome that influence whether, and by how much, a gene is expressed. In particular, examining multiple tissues can give us valuable insights into the genetic architecture of the regulatory mechanism, because many regulatory regions are known to act in a tissue specific manner (Ernst et al. 2011; Encode Project Consortium 2012). Hence, understanding the role of regulatory variants, and the tissues in which they act, is essential for the functional interpretation of GWAS loci and insights into disease etiology.
Figure 2 is an example of the output of ForestPMPlot for a multitissue eQTL study for SEMA3B gene (GTEx Consortium 2015). Examining both the forest plot and the PM-Plot allows us to obtain an insight into the tissue-specific genetics effects in eQTL analysis, which leads to the identification of three significant eQTL tissues (heart left ventricle, stomach, and thyroid). This example clearly shows that examining both the forest plot and the PM-Plot allows us to easily hypothesize that there is a specific group of studies showing tissue differences in eQTL analysis.
ForestPMPlot: A Flexible Tool for Visualizing Heterogeneity Between Studies in Meta-analysis
Below, we provide examples of several different types of figures that can illustrate the point of a Results subsection.