Rare Variants Association Analysis in Large-Scale Sequencing Studies at the Single Locus Level

GWAS for identification of genomic regions and candidate genes in vegetable crops

Genome-wide association Studies (GWAS), initially developed for human genetics, have been highly effective in plant research, particularly for vegetable crops. GWAS is a robust tool for identifying genes associated with key traits such as yield, nutritional value, disease resistance, adaptability, and bioactive compound biosynthesis. Unlike traditional methods, GWAS does not require prior biological knowledge and can accurately pinpoint loci, minimizing false positives. The process involves developing a diverse panel, rigorous phenotyping and genotyping, and sophisticated statistical analysis using various models and software tools. By scanning the entire genome, GWAS identifies specific loci or single nucleotide polymorphisms (SNPs) linked to target traits. When a causal SNP variant is not directly genotyped, GWAS identifies SNPs in linkage disequilibrium (LD) with the causal variant, mapping the genetic interval. The method begins with careful panel selection, phenotyping, and genotyping, controlling for environmental effects and utilizing Best Linear Unbiased Prediction (BLUP). High-correlation, high-heritability traits are prioritized. Various genotyping methods address confounders like population structure and kinship. Bonferroni correction (BC) prevents false positives, and significant associations are shown in Manhattan plots. Candidate genes are identified through LD analysis and fine mapping, followed by functional validation. GWAS offers critical insights for enhancing vegetable crop breeding efficiency and precision, driving breakthroughs through advanced methods.

Spatially adaptive variable screening in presurgical functional magnetic resonance imaging data analysis

Accurate delineation of functional brain regions adjacent to tumors is imperative for planning neurosurgery that preserves critical functions. Functional magnetic resonance imaging (fMRI) plays an increasingly pivotal role in presurgical counseling and planning. In the analysis of presurgical fMRI data, the impact of false negatives on patients surpasses that of false positives because failure to identify functional regions and unintentionally resecting critical tissues can result in severe harm to patients. This paper introduces a novel metric, the Bayesian missed discovery rate (BMDR), designed for controlling false negatives within the voxel-specific mixture model. Building on the BMDR metric, we propose a new variable screening procedure that not only ensures effective control of false negatives but also capitalizes on the spatial structure of fMRI data. In comparison to existing statistical methods in fMRI data analysis, our new procedure directly regulates false negatives at a desirable level and is entirely data-driven. Moreover, it significantly differs from current false-negative control procedures by incorporating spatial information. Numerical examples demonstrate that the new method outperforms several state-of-the-art methods in retaining signal voxels, particularly the subtle ones at the boundaries of functional regions, while achieving a cleaner separation of functional regions from background noise. These findings hold promising implications for planning function-preserving neurosurgery.

Effective SNP ranking improves the performance of eQTL mapping

Genome-wide expression quantitative trait loci (eQTLs) mapping explores the relationship between gene expression and DNA variants, such as single-nucleotide polymorphism (SNPs), to understand genetic basis of human diseases. Due to the large number of genes and SNPs that need to be assessed, current methods for eQTL mapping often suffer from low detection power, especially for identifying trans-eQTLs. In this paper, we propose the idea of performing SNP ranking based on the higher criticism statistic, a summary statistic developed in large-scale signal detection. We illustrate how the HC-based SNP ranking can effectively prioritize eQTL signals over noise, greatly reduce the burden of joint modeling, and improve the power for eQTL mapping. Numerical results in simulation studies demonstrate the superior performance of our method compared to existing methods. The proposed method is also evaluated in HapMap eQTL data analysis and the results are compared to a database of known eQTLs.

Identifying individual risk rare variants using protein structure guided local tests (POINT)

Rare variants are of increasing interest to genetic association studies because of their etiological contributions to human complex diseases. Due to the rarity of the mutant events, rare variants are routinely analyzed on an aggregate level. While aggregation analyses improve the detection of global-level signal, they are not able to pinpoint causal variants within a variant set. To perform inference on a localized level, additional information, e.g., biological annotation, is often needed to boost the information content of a rare variant. Following the observation that important variants are likely to cluster together on functional domains, we propose a protein structure guided local test (POINT) to provide variant-specific association information using structure-guided aggregation of signal. Constructed under a kernel machine framework, POINT performs local association testing by borrowing information from neighboring variants in the 3-dimensional protein space in a data-adaptive fashion. Besides merely providing a list of promising variants, POINT assigns each variant a p-value to permit variant ranking and prioritization. We assess the selection performance of POINT using simulations and illustrate how it can be used to prioritize individual rare variants in PCSK9, ANGPTL4 and CETP in the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial data.

A modified association test for rare and common variants based on affected sib-pair design

Current genome-wide association analysis has identified a great number of rare and common variants associated with common complex traits, however, more effective approaches for detecting associations between rare and common variants with common diseases are still demanded. Approaches for detecting rare variant association analysis will compromise the power when detecting the effects of rare and common variants simultaneously. In this paper, we extend an existing method of testing for rare variant association based on affected sib pairs (TOW-sib) and propose a variable weight test for rare and common variants association based on affected sib pairs (abbreviated as VW-TOWsib). The VW-TOWsib can be used to achieve the purpose of detecting the association of rare and common variants with complex diseases. Simulation results in various scenarios show that our proposed method is more powerful than existing methods for detecting effects of rare and common variants. At the same time, the VW-TOWsib also performs well as a method for rare variant association analysis.

Large-Scale Analysis of Genetic and Clinical Patient Data

Biomedical data science has experienced an explosion of new data over the past decade. Abundant genetic and genomic data are increasingly available in large, diverse data sets due to the maturation of modern molecular technologies. Along with these molecular data, dense, rich phenotypic data are also available on comprehensive clinical data sets from health care provider organizations, clinical trials, population health registries, and epidemiologic studies. The methods and approaches for interrogating these large genetic/genomic and clinical data sets continue to evolve rapidly, as our understanding of the questions and challenges continue to emerge. In this review, the state-of-the-art methodologies for genetic/genomic analysis along with complex phenomics will be discussed. This field is changing and adapting to the novel data types made available, as well as technological advances in computation and machine learning. Thus, I will also discuss the future challenges in this exciting and innovative space. The promises of precision medicine rely heavily on the ability to marry complex genetic/genomic data with clinical phenotypes in meaningful ways.

Whole-Exome Sequencing to Identify Novel Biological Pathways Associated With Infertility After Pelvic Inflammatory Disease

BACKGROUND

Ideal management of sexually transmitted infections (STI) may require risk markers for pathology or vaccine development. Previously, we identified common genetic variants associated with chlamydial pelvic inflammatory disease (PID) and reduced fecundity. As this explains only a proportion of the long-term morbidity risk, we used whole-exome sequencing to identify biological pathways that may be associated with STI-related infertility.

METHODS

We obtained stored DNA from 43 non-Hispanic black women with PID from the PID Evaluation and Clinical Health Study. Infertility was assessed at a mean of 84 months. Principal component analysis revealed no population stratification. Potential covariates did not significantly differ between groups. Sequencing kernel association test was used to examine associations between aggregates of variants on a single gene and infertility. The results from the sequencing kernel association test were used to choose "focus genes" (P < 0.01; n = 150) for subsequent Ingenuity Pathway Analysis to identify "gene sets" that are enriched in biologically relevant pathways.

RESULTS

Pathway analysis revealed that focus genes were enriched in canonical pathways including, IL-1 signaling, P2Y purinergic receptor signaling, and bone morphogenic protein signaling.

CONCLUSIONS

Focus genes were enriched in pathways that impact innate and adaptive immunity, protein kinase A activity, cellular growth, and DNA repair. These may alter host resistance or immunopathology after infection. Targeted sequencing of biological pathways identified in this study may provide insight into STI-related infertility.