The need to select mediators from a high dimensional data source, such as neuroimaging data and genetic data, arises in much scientific research. In this work, we formulate a multiple-hypothesis testing framework for mediator selection from a high-dimensional candidate set, and propose a method, which extends the recent development in false discovery rate (FDR)-controlled variable selection with knockoff to select mediators with FDR control. We show that the proposed method and algorithm achieved finite sample FDR control. We present extensive simulation results to demonstrate the power and finite sample performance compared with the existing method. Lastly, we demonstrate the method for analyzing the Adolescent Brain Cognitive Development (ABCD) study, in which the proposed method selects several resting-state functional magnetic resonance imaging connectivity markers as mediators for the relationship between adverse childhood events and the crystallized composite score in the NIH toolbox.

BACKGROUND: Fetal sex affects fetal and maternal health outcomes in pregnancy, but this connection remains poorly understood. As the placenta is the route of fetomaternal communication and derives from the fetal genome, placental gene expression sex differences may explain these outcomes.
OBJECTIVE: We utilized next generation sequencing to study the normal human placenta in both sexes in first and third trimester to generate a normative transcriptome based on sex and gestation.
METHODS: We analyzed 124 first trimester (T1, 59 female and 65 male) and 43 third trimester (T3, 18 female and 25 male) samples for sex differences within each trimester and sex-specific gestational differences.
RESULTS: Placenta shows more significant sexual dimorphism in T1, with 94 T1 and 26 T3 differentially expressed genes (DEGs). The sex chromosomes contributed 60.6% of DEGs in T1 and 80.8% of DEGs in T3, excluding X/Y pseudoautosomal regions. There were 6 DEGs from the pseudoautosomal regions, only significant in T1 and all upregulated in males. The distribution of DEGs on the X chromosome suggests genes on Xp (the short arm) may be particularly important in placental sex differences. Dosage compensation analysis of X/Y homolog genes shows expression is primarily contributed by the X chromosome. In sex-specific analyses of first versus third trimester, there were 2815 DEGs common to both sexes upregulated in T1, and 3263 common DEGs upregulated in T3. There were 7 female-exclusive DEGs upregulated in T1, 15 female-exclusive DEGs upregulated in T3, 10 male-exclusive DEGs upregulated in T1, and 20 male-exclusive DEGs upregulated in T3.
CONCLUSIONS: This is the largest cohort of placentas across gestation from healthy pregnancies defining the normative sex dimorphic gene expression and sex common, sex specific and sex exclusive gene expression across gestation. The first trimester has the most sexually dimorphic transcripts, and the majority were upregulated in females compared to males in both trimesters. The short arm of the X chromosome and the pseudoautosomal region is particularly critical in defining sex differences in the first trimester placenta. As pregnancy is a dynamic state, sex specific DEGs across gestation may contribute to sex dimorphic changes in overall outcomes.

Neuropeptides represent a unique class of signaling molecules that have garnered much attention but require special consideration when identifications are gleaned from mass spectra. With highly variable sequence lengths, neuropeptides must be analyzed in their endogenous state. Further, neuropeptides share great homology within families, differing by as little as a single amino acid residue, complicating even routine analyses and necessitating optimized computational strategies for confident and accurate identifications. We present EndoGenius, a database searching strategy designed specifically for elucidating neuropeptide identifications from mass spectra by leveraging optimized peptide-spectrum matching approaches, an expansive motif database, and a novel scoring algorithm to achieve broader representation of the neuropeptidome and minimize reidentification. This work describes an algorithm capable of reporting more neuropeptide identifications at 1% false-discovery rate than alternative software in five Callinectes sapidus neuronal tissue types.

Top-down proteomics (TDP) directly analyzes intact proteins and thus provides more comprehensive qualitative and quantitative proteoform-level information than conventional bottom-up proteomics (BUP) that relies on digested peptides and protein inference. While significant advancements have been made in TDP in sample preparation, separation, instrumentation, and data analysis, reliable and reproducible data analysis still remains one of the major bottlenecks in TDP. A key step for robust data analysis is the establishment of an objective estimation of proteoform-level false discovery rate (FDR) in proteoform identification. The most widely used FDR estimation scheme is based on the target-decoy approach (TDA), which has primarily been established for BUP. We present evidence that the TDA-based FDR estimation may not work at the proteoform-level due to an overlooked factor, namely the erroneous deconvolution of precursor masses, which leads to incorrect FDR estimation. We argue that the conventional TDA-based FDR in proteoform identification is in fact protein-level FDR rather than proteoform-level FDR unless precursor deconvolution error rate is taken into account. To address this issue, we propose a formula to correct for proteoform-level FDR bias by combining TDA-based FDR and precursor deconvolution error rate.

Phosphorylation is a post-translational modification of great interest to researchers due to its relevance in many biological processes. LC-MS/MS techniques have enabled high-throughput data acquisition, with studies claiming identification and localization of thousands of phosphosites. The identification and localization of phosphosites emerge from different analytical pipelines and scoring algorithms, with uncertainty embedded throughout the pipeline. For many pipelines and algorithms, arbitrary thresholding is used, but little is known about the actual global false localization rate in these studies. Recently, it has been suggested to use decoy amino acids to estimate global false localization rates of phosphosites, among the peptide-spectrum matches reported. Here, we describe a simple pipeline aiming to maximize the information extracted from these studies by objectively collapsing from peptide-spectrum match to the peptidoform-site level, as well as combining findings from multiple studies while maintaining track of false localization rates. We show that the approach is more effective than current processes that use a simpler mechanism for handling phosphosite identification redundancy within and across studies. In our case study using eight rice phosphoproteomics data sets, 6368 unique sites were confidently identified using our decoy approach compared to 4687 using traditional thresholding in which false localization rates are unknown.

Neuropeptides are a class of endogenous peptides that have key regulatory roles in biochemical, physiological, and behavioral processes. Mass spectrometry analyses of neuropeptides often rely on protein informatics tools for database searching and peptide identification. As neuropeptide databases are typically experimentally built and comprised of short sequences with high sequence similarity to each other, we developed a novel database searching tool, HyPep, which utilizes sequence homology searching for peptide identification. HyPep aligns de novo sequenced peptides, generated through PEAKS software, with neuropeptide database sequences and identifies neuropeptides based on the alignment score. HyPep performance was optimized using LC-MS/MS measurements of peptide extracts from various Callinectes sapidus neuronal tissue types and compared with a commercial database searching software, PEAKS DB. HyPep identified more neuropeptides from each tissue type than PEAKS DB at 1% false discovery rate, and the false match rate from both programs was 2%. In addition to identification, this report describes how HyPep can aid in the discovery of novel neuropeptides.

Protein methylation is a widespread post-translational modification (PTM) involved in several important biological processes including, but not limited to, RNA splicing, signal transduction, translation, and DNA repair. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is considered today the most versatile and accurate technique to profile PTMs with high precision and proteome-wide depth; however, the identification of protein methylations by MS is still prone to high false discovery rates. In this chapter, we describe the heavy methyl SILAC metabolic labeling strategy that allows high-confidence identification of in vivo methyl-peptides by MS-based proteomics. We provide a general protocol that covers the steps of heavy methyl labeling of cultured cells, protein sample preparation, LC-MS/MS analysis, and downstream computational analysis of the acquired MS data.

BACKGROUND: Genome-wide association studies (GWAS) are limited in power to detect associations that exceed the stringent genome-wide significance threshold. This limitation can be alleviated by leveraging relevant auxiliary data, such as functional genomic data. Frameworks utilising the conditional false discovery rate have been developed for this purpose, and have been shown to increase power for GWAS discovery whilst controlling the false discovery rate. However, the methods are currently only applicable for continuous auxiliary data and cannot be used to leverage auxiliary data with a binary representation, such as whether SNPs are synonymous or non-synonymous, or whether they reside in regions of the genome with specific activity states.
RESULTS: We describe an extension to the cFDR framework for binary auxiliary data, called \"Binary cFDR\". We demonstrate FDR control of our method using detailed simulations, and show that Binary cFDR performs better than a comparator method in terms of sensitivity and FDR control. We introduce an all-encompassing user-oriented CRAN R package ( https://annahutch.github.io/fcfdr/ ; https://cran.r-project.org/web/packages/fcfdr/index.html ) and demonstrate its utility in an application to type 1 diabetes, where we identify additional genetic associations.
CONCLUSIONS: Our all-encompassing R package, fcfdr, serves as a comprehensive toolkit to unite GWAS and functional genomic data in order to increase statistical power to detect genetic associations.

Irritable bowel syndrome with diarrhoea (IBS-D) and functional diarrhoea (FDr) are the two major functional bowel disorders characterized by diarrhoea. In spite of their high prevalence, IBS-D and FDr are associated with major uncertainties, especially regarding their optimal diagnostic work-up and management. A Delphi consensus was performed with experts from 10 European countries who conducted a literature summary and voting process on 31 statements. Quality of evidence was evaluated using the grading of recommendations, assessment, development, and evaluation criteria. Consensus (defined as >80% agreement) was reached for all the statements. The panel agreed with the potential overlapping of IBS-D and FDr. In terms of diagnosis, the consensus supports a symptom-based approach also with the exclusion of alarm symptoms, recommending the evaluation of full blood count, C-reactive protein, serology for coeliac disease, and faecal calprotectin, and consideration of diagnosing bile acid diarrhoea. Colonoscopy with random biopsies in both the right and left colon is recommended in patients older than 50 years and in presence of alarm features. Regarding treatment, a strong consensus was achieved for the use of a diet low fermentable oligo-, di-, monosaccharides and polyols, gut-directed psychological therapies, rifaximin, loperamide, and eluxadoline. A weak or conditional recommendation was achieved for antispasmodics, probiotics, tryciclic antidepressants, bile acid sequestrants, 5-hydroxytryptamine-3 antagonists (i.e. alosetron, ondansetron, or ramosetron). A multinational group of European experts summarized the current state of consensus on the definition, diagnosis, and management of IBS-D and FDr.

Proteogenomics refers to the integrated analysis of the genome and proteome that leverages mass-spectrometry (MS)-based proteomics data to improve genome annotations, understand gene expression control through proteoforms and find sequence variants to develop novel insights for disease classification and therapeutic strategies. However, proteogenomic studies often suffer from reduced sensitivity and specificity due to inflated database size. To control the error rates, proteogenomics depends on the target-decoy search strategy, the de-facto method for false discovery rate (FDR) estimation in proteomics. The proteogenomic databases constructed from three- or six-frame nucleotide database translation not only increase the search space and compute-time but also violate the equivalence of target and decoy databases. These searches result in poorer separation between target and decoy scores, leading to stringent FDR thresholds. Understanding these factors and applying modified strategies such as two-pass database search or peptide-class-specific FDR can result in a better interpretation of MS data without introducing additional statistical biases. Based on these considerations, a user can interpret the proteogenomics results appropriately and control false positives and negatives in a more informed manner. In this review, first, we briefly discuss the proteogenomic workflows and limitations in database construction, followed by various considerations that can influence potential novel discoveries in a proteogenomic study. We conclude with suggestions to counter these challenges for better proteogenomic data interpretation.