Multitest Correction Research Articles

Identifying the Genetic Basis of Mineral Elements in Rice Grain Using Genome-Wide Association Mapping.

Mineral malnutrition is a major problem in many rice-consuming countries. It is essential to know the genetic mechanisms of accumulation of mineral elements in the rice grain to provide future solutions for this issue. This study was conducted to identify the genetic basis of six mineral elements (Cu, Fe, K, Mg, Mn, and Zn) by using three models for single-locus and six models for multi-locus analysis of a genome-wide association study (GWAS) using 174 diverse rice accessions and 6565 SNP markers. To declare a SNP as significant, -log10(P) ≥ 3.0 and 15% FDR significance cut-off values were used for single-locus models, while LOD ≥ 3.0 was used for multi-locus models. Using these criteria, 147 SNPs were detected by one or two GWAS methods at -log10(P) ≥ 3.0, 48 of which met the 15% FDR significance cut-off value. Single-locus models outperformed multi-locus models before applying multi-test correction, but once applied, multi-locus models performed better. While 14 (~29%) of the identified quantitative trait loci (QTLs) after multiple test correction co-located with previously reported genes/QTLs and marker associations, another 34 trait-associated SNPs were novel. After mining genes within 250 kb of the 48 significant SNP loci, in silico and gene enrichment analyses were conducted to predict their potential functions. These shortlisted genes with their functions could guide future experimental validation, helping us to understand the complex molecular mechanisms controlling rice grain mineral elements.

Is it possible ABC transporters genetic variants influence the outcomes of a weight-loss diet in obese women?

ATP-Binding Cassette (ABC) transporters are involved in cholesterol metabolism and their dysfunctions could lead to obesity-associated complications. It was investigated whether SNPs in the ABCA1 (rs1800977 and rs2230806), ABCA7 (rs2279796) and ABCG1 (rs692383 and rs3827225) genes can modulate the responsiveness of 137 obese women to a weight-loss dietary intervention. Thus, anthropometric and lipid profiles were collected at baseline and after nine weeks of a calorie-restricted diet of 600kcal per day and participants were genotyped for the ABC genes SNPs. Regarding the transversal analysis, the ABCA7 rs2279796 GG genotype was associated with higher levels of total cholesterol and LDL-c at baseline (p = 0.044 for both). Association between ABCG1 rs692383 AA genotype and lower BMI were found in the post-diet moment, however, statistical significance was lost after multi-test correction. Regarding the longitudinal analysis, after multi-test correction, the association remained between ABCG1 rs692383 G allele and HDL-c levels: G allele carriers had a lower HDL-c reduction (p = 0.043). Results suggest the standard weight-loss diet applied in this study could attenuate the ABCA7 rs2279796 GG genotype effects found at baseline and non-dyslipidemic obese women with ABCG1 rs692383 G allele are benefitting from the diet with a lower reduction in HDL-c levels.

Tissue-Specific eQTL in Zebrafish.

Copy number variants (CNVs) refer to the loss or gain of copies of a genomic DNA region. While some CNVs may play a role in species evolution by enriching the diversity of an organism, CNVs may also be linked to certain diseases such as neurological disorders, early onset obesity, and cancer. CNVs may affect gene expression by direct overlap of the genic region or by an indirect effect where the CNV is located outside the gene location. These indirect CNV regions may contain regulatory elements such as transcription enhancers or repressors as well as regulators such as miRNAs which may work at the level of transcription or translation. Danio rerio (zebrafish) is an excellent model organism for CNV studies. Zebrafish genomes contain a large amount of variation with 14.6% of the zebrafish reference genome found to be copy number variable. This level of variation is more than four times the percentage of reference genome sequence covered by similarly common CNVs in humans. It is this high level of variation that makes zebrafish interesting to investigate the effects of CNV on gene expression. Additionally, zebrafish share 70% of genetic similarities with humans, and 84% of genes associated with human disease are also found in zebrafish. Expressive quantitative trait loci (eQTL) analysis may be used in zebrafish to explore how CNVs may affect gene expression in both a direct and indirect manner. eQTL analysis may be performed for cis associations with a 1-Mb (megabase) window upstream and downstream from the transcription probe midpoint to CGH midpoint. Trans associations (variants that are located beyond the 1-Mb window of the gene either on the same chromosome as the gene or on a different chromosome) may be investigated as well through eQTL analysis; however, trans associations require more tests to be performed than cis associations, which limits power to detect associations. Pairwise associations between each pair of copy number variant and gene will be investigated separately from the same individual using Spearman rank correlations with significant associations found being followed with a multi-test correction technique to assess significance of those CNV gene expression associations. An association between a CNV to a gene expression phenotype should be considered significant only if the p value from the analysis of the observed data is lower than the 0.001 tail threshold from a distribution of the minimal p values (which are found from all comparisons for a given gene from 10,000 permutations of the expression phenotypes). Associations between CNVs and genes may be found to be direct or indirect as well as positive (increased copy number-increased expression) or negative (increased copy number-decreased expression, decreased copy number-increased expression). Ongoing analyses with these associations will investigate the impact of CNVs on gene functionality including immune function and potential disease susceptibility.

Genetic differences between wild and hatchery-bred brown trout (Salmo trutta L.) in single nucleotide polymorphisms linked to selective traits.

To study effects from natural selection acting on brown trout in a natural stream habitat compared with a hatchery environment, 3,781 single nucleotide polymorphism (SNP) markers were analyzed in three closely related groups of brown trout (Salmo trutta L.). Autumn (W/0+, n = 48) and consecutive spring (W/1+, n = 47) samples of brown trout individuals belonging to the same cohort and stream were retrieved using electrofishing. A third group (H/1+, n = 48) comprised hatchery‐reared individuals, bred from a mixture of wild parents of the strain of the two former groups and from a neighboring stream. Pairwise analysis of FST outliers and analysis under a hierarchical model by means of ARLEQUIN software detected 421 (10.8%) candidates of selection, before multitest correction. BAYESCAN software detected 10 candidate loci, all of which were included among the ARLEQUIN candidate loci. Body length was significantly different across genotypes at 10 candidate loci in the W/0+, at 34 candidate loci in the W/1+ and at 21 candidate loci in the H/1+ group. The W/1+ sample was tested for genotype‐specific body length at all loci, and significant differences were found in 10.6% of all loci, and of these, 14.2% had higher frequency of the largest genotype in the W/1+ sample than in W/0+. The corresponding proportion among the candidate loci of W/1+ was 22.7% with genotype‐specific body length, and 88.2% of these had increased frequency of the largest genotype from W/0+ to W/1+, indicating a linkage between these loci and traits affecting growth and survival under this stream's environmental conditions. Bayesian structuring of all loci, and of the noncandidate loci suggested two (K = 2), alternatively four clusters (K = 4). This differed from the candidate SNPs, which suggested only two clusters. In both cases, the hatchery fish dominated one cluster, and body length of W/1+ fish was positively correlated with membership of one cluster both from the K = 2 and the K = 4 structure. Our analysis demonstrates profound genetic differentiation that can be linked to differential selection on a fitness‐related trait (individual growth) in brown trout living under natural vs. hatchery conditions. Candidate SNP loci linked to genes affecting individual growth were identified and provide important inputs into future mapping of the genetic basis of brown trout body size selection.

Exaggerated CpH methylation in the autism-affected brain

BackgroundThe etiology of autism, a complex, heritable, neurodevelopmental disorder, remains largely unexplained. Given the unexplained risk and recent evidence supporting a role for epigenetic mechanisms in the development of autism, we explored the role of CpG and CpH (H = A, C, or T) methylation within the autism-affected cortical brain tissue.MethodsReduced representation bisulfite sequencing (RRBS) was completed, and analysis was carried out in 63 post-mortem cortical brain samples (Brodmann area 19) from 29 autism-affected and 34 control individuals. Analyses to identify single sites that were differentially methylated and to identify any global methylation alterations at either CpG or CpH sites throughout the genome were carried out.ResultsWe report that while no individual site or region of methylation was significantly associated with autism after multi-test correction, methylated CpH dinucleotides were markedly enriched in autism-affected brains (~2-fold enrichment at p < 0.05 cutoff, p = 0.002).ConclusionsThese results further implicate epigenetic alterations in pathobiological mechanisms that underlie autism.

A pharmacogenomic study revealed an association between SLC6A4 and risperidone-induced weight gain in Chinese Han population.

We carried out a pharmacogenomic study in order to identify susceptible genes for antipsychotics induced weight gain within the Chinese Han population. We enrolled 216 patients with schizophrenia in our study. All of them underwent risperidone monotherapy, and fulfilled 4-week follow-up. Weight gain was measured before treatment and 4 weeks later. Seven hundred and sixty-eight SNPs from 85 genes were calculated for association with weight gain percentage. Fifty-seven SNPs located at 16 genes with a p-value less than 0.05.4 SNPs located on serotonin transporter gene (solute carrier family 6, member 4, SLC6A4) remained significant after multitest correction (rs3813034, p = 0.000357, q = 0.08, rs1042173, rs4325622, rs9303628, p = 0.000451, q = 0.08). SLC6A4 might be susceptible gene for risperidone-induced weight gain within the Chinese Han population.

Hepatitis B viraemia: its heritability and association with common genetic variation in the interferon γ signalling pathway

ObjectiveHigh viraemia of hepatitis B virus (HBV) influences all phases in the development of hepatocellular carcinoma (HCC). This study was designed to estimate the overall contribution of host genetics to...

Comparing three different methods to detect selective loci using dominant markers

We carried out a simulation study to compare the efficiency of three alternative programs (DFDIST, DETSELD and BAYESCAN) to detect loci under directional selection from genome-wide scans using dominant markers. We also evaluated the efficiency of correcting for multiple testing those methods that use a classical probability approach. Under a wide range of scenarios, we conclude that BAYESCAN appears to be more efficient than the other methods, detecting a usually high percentage of true selective loci as well as less than 1% of outliers (false positives) under a fully neutral model. In addition, the percentage of outliers detected by this software is always correlated with the true percentage of selective loci in the genome. Our results show, nevertheless, that false positives are common even with a combination of methods and multitest correction, suggesting that conclusions obtained from this approach should be taken with extreme caution.

A new multitest correction (SGoF) that increases its statistical power when increasing the number of tests

BackgroundThe detection of true significant cases under multiple testing is becoming a fundamental issue when analyzing high-dimensional biological data. Unfortunately, known multitest adjustments reduce their statistical power as the number of tests increase. We propose a new multitest adjustment, based on a sequential goodness of fit metatest (SGoF), which increases its statistical power with the number of tests. The method is compared with Bonferroni and FDR-based alternatives by simulating a multitest context via two different kinds of tests: 1) one-sample t-test, and 2) homogeneity G-test.ResultsIt is shown that SGoF behaves especially well with small sample sizes when 1) the alternative hypothesis is weakly to moderately deviated from the null model, 2) there are widespread effects through the family of tests, and 3) the number of tests is large.ConclusionTherefore, SGoF should become an important tool for multitest adjustment when working with high-dimensional biological data.

Identification of novel susceptibility genes in childhood‐onset systemic lupus erythematosus using a uniquely designed candidate gene pathway platform

Childhood-onset systemic lupus erythematosus (SLE) presents a unique subgroup of patients for genetic study. The present study was undertaken to identify susceptibility genes contributing to SLE, using a novel candidate gene pathway microarray platform to investigate gene expression in patients with childhood-onset SLE and both of their parents. Utilizing bioinformatic tools, a platform of 9,412 single-nucleotide polymorphisms (SNPs) from 1,204 genes was designed and validated. Molecular inversion probes and high-throughput SNP technologies were used for assay development. Seven hundred fifty three subjects, corresponding to 251 full trios of childhood-onset SLE families, were genotyped and analyzed using transmission disequilibrium testing (TDT) and multitest corrections. Family-based TDT showed a significant association of SLE with a N673S polymorphism in the P-selectin gene (SELP) (P = 5.74 x 10(-6)) and a C203S polymorphism in the interleukin-1 receptor-associated kinase 1 gene (IRAK1) (P = 9.58 x 10(-6)). These 2 SNPs had a false discovery rate for multitest correction of <0.05, and therefore a >95% probability of being considered as proven. Furthermore, 7 additional SNPs showed q values of <0.5, suggesting association with SLE and providing a direction for followup studies. These additional genes notably included TNFRSF6 (Fas) and IRF5, supporting previous findings of their association with SLE pathogenesis. SELP and IRAK1 were identified as novel SLE-associated genes with a high degree of significance, suggesting new directions in understanding the pathogenesis of SLE. The overall design and results of this study demonstrate that the candidate gene pathway microarray platform used provides a novel and powerful approach that is generally applicable in identifying genetic foundations of complex diseases.

A critical evaluation of discography in patients with lumbar intervertebral disc disease

Study design The study is a prospective observational study of 48 continuous patients with symptomatic lumbar degenerative disk disease. Each patient underwent discography, MRI, and a biochemical analysis of disk lavage fluid. Objectives The purpose of this study was to correlate concordant pain on discography with MRI grade and biochemical markers of inflammation in a clinical setting. Summary of background data The pathophysiology of degenerative disk disease is complex. Discography is used to differentiate symptomatic from asymptomatic levels. MRI is used to image changes in disk water content. Biochemical assays have identified molecular markers of inflammation. To date, no study has correlated concordant pain on discography with MRI findings and biochemical markers. Methods Forty-eight (48) continuous patients with symptomatic lumbar degenerative disk disease gave informed consent for study entry. Patient sex, age, insurance, work status and visual analog score (VAS) were recorded. MRI was obtained and Pfirrmann grading was performed by a single spine surgeon. Discography with disc lavage was performed by a single anesthesiologist. Lavage samples were tested for inflammatory markers with high resolution multi-plex bead immunoassays and ELISA with >5 pg/ml resolution. Results None of demographic variables was significantly related to concordant pain on discogram by chi-squared tests and Mann-Whitney U-test. The Pfirrmann score was significantly different for patients with and without concordant pain at L3-L4 (p<0.001), but was insignificant at other levels after multitest correction. Pfirrmann scores were significantly different at any level in patients with and without concordant pain. VAS scores were not significantly correlated with opening pressures at any level. Despite the presence of serum proteins in the disk lavage fluid, none of the tested inflammatory mediators was identified by multi-plex bead immunoassays and ELISA. Conclusions There are only weak correlations between demographic, discogram, and radiographic variables. Response to discogram cannot be predicted by non-invasive means. The disk lavage method was unable to identify the presence of specific inflammatory peptides with multi-plex immunoassays and ELISA.