Leveraging Prior Information to Detect Causal Variants via Multi-Variant Regression

Although many methods are available to test sequence variants for association with complex diseases and traits, methods that specifically seek to identify causal variants are less developed. Here we develop and evaluate a Bayesian hierarchical regression method that incorporates prior information on the likelihood of variant causality through weighting of variant effects. By simulation studies using both simulated and real sequence variants, we compared a standard single variant test for analyzing variant-disease association with the proposed method using different weighting schemes. We found that by leveraging linkage disequilibrium of variants with known GWAS signals and sequence conservation (phastCons), the proposed method provides a powerful approach for detecting causal variants while controlling false positives.     

Genetic Variability for Pathogenicity, Isozyme, Ribosomal DNA and Colony Color Variants in Populations of Rhynchosporium Secalis

Samples of Rhynchosporium secalis were collected from two experimental barley populations known to carry a diverse array of alleles for resistance to this fungal pathogen. Classification of 163 isolates for four putative isozyme systems, a colony color dimorphism and 20 ribosomal DNA restriction fragment length variants revealed 49 different multilocus phenotypes (haplotypes). The six most common haplotypes differed significantly in pathogenicity. Genetic analyses of the data indicated that effective population sizes of the fungus were very large, that the effects of genetic drift were small, and that negligible recombination occurred in the populations studied. Frequency dependent selection was suggested as an explanation for the maintenance of variation in pathogenicity in the fungus.     

Identifying Genetic Variants for Heart Rate Variability in the Acetylcholine Pathway

Heart rate variability is an important risk factor for cardiovascular disease and all-cause mortality. The acetylcholine pathway plays a key role in explaining heart rate variability in humans. We assessed whether 443 genotyped and imputed common genetic variants in eight key genes (CHAT, SLC18A3, SLC5A7, CHRNB4, CHRNA3, CHRNA, CHRM2 and ACHE) of the acetylcholine pathway were associated with variation in an established measure of heart rate variability reflecting parasympathetic control of the heart rhythm, the root mean square of successive differences (RMSSD) of normal RR intervals. The association was studied in a two stage design in individuals of European descent. First, analyses were performed in a discovery sample of four cohorts (n = 3429, discovery stage). Second, findings were replicated in three independent cohorts (n = 3311, replication stage), and finally the two stages were combined in a meta-analysis (n = 6740). RMSSD data were obtained under resting conditions. After correction for multiple testing, none of the SNPs showed an association with RMSSD. In conclusion, no common genetic variants for heart rate variability were identified in the largest and most comprehensive candidate gene study on the acetylcholine pathway to date. Future gene finding efforts for RMSSD may want to focus on hypothesis free approaches such as the genome-wide association study.     

Hidden Genetic Variability within Electromorphs in Finite Populations

The amount of hidden genetic variability within electromorphs in finite populations is studied by using the infinite site model and stepwise mutation model simultaneously. A formula is developed for the bivariate probability generating function for the number of codon differences and the number of electromorph state differences between two randomly chosen cistrons. Using this formula, the distribution as well as the mean and variance of the number of codon differences between two identical or nonidentical electromorphs are studied. The distribution of the number of codon differences between two randomly chosen identical electromorphs is similar to the geometric distribution but more leptokurtic. Studies are also made on the number of codon differences between two electromorphs chosen at random one from each of two populations which have been separated for an arbitrary number of generations. It is shown that the amount of hidden genetic variability is very large if the product of effective population size and mutation rate is large.     

Genetic Variants Contribute to Gene Expression Variability in Humans

Expression quantitative trait loci (eQTL) studies have established convincing relationships between genetic variants and gene expression. Most of these studies focused on the mean of gene expression level, but not the variance of gene expression level (i.e., gene expression variability). In the present study, we systematically explore genome-wide association between genetic variants and gene expression variability in humans. We adapt the double generalized linear model (dglm) to simultaneously fit the means and the variances of gene expression among the three possible genotypes of a biallelic SNP. The genomic loci showing significant association between the variances of gene expression and the genotypes are termed expression variability QTL (evQTL). Using a data set of gene expression in lymphoblastoid cell lines (LCLs) derived from 210 HapMap individuals, we identify cis-acting evQTL involving 218 distinct genes, among which 8 genes, ADCY1, CTNNA2, DAAM2, FERMT2, IL6, PLOD2, SNX7, and TNFRSF11B, are cross-validated using an extra expression data set of the same LCLs. We also identify ∼300 trans-acting evQTL between >13,000 common SNPs and 500 randomly selected representative genes. We employ two distinct scenarios, emphasizing single-SNP and multiple-SNP effects on expression variability, to explain the formation of evQTL. We argue that detecting evQTL may represent a novel method for effectively screening for genetic interactions, especially when the multiple-SNP influence on expression variability is implied. The implication of our results for revealing genetic mechanisms of gene expression variability is discussed.     

Identification of New Genetic Risk Variants for Type 2 Diabetes

Although more than 20 genetic susceptibility loci have been reported for type 2 diabetes (T2D), most reported variants have small to moderate effects and account for only a small proportion of the heritability of T2D, suggesting that the majority of inter-person genetic variation in this disease remains to be determined. We conducted a multistage, genome-wide association study (GWAS) within the Asian Consortium of Diabetes to search for T2D susceptibility markers. From 590,887 SNPs genotyped in 1,019 T2D cases and 1,710 controls selected from Chinese women in Shanghai, we selected the top 2,100 SNPs that were not in linkage disequilibrium (r²<0.2) with known T2D loci for in silico replication in three T2D GWAS conducted among European Americans, Koreans, and Singapore Chinese. The 5 most promising SNPs were genotyped in an independent set of 1,645 cases and 1,649 controls from Shanghai, and 4 of them were further genotyped in 1,487 cases and 3,316 controls from 2 additional Chinese studies. Consistent associations across all studies were found for rs1359790 (13q31.1), rs10906115 (10p13), and rs1436955 (15q22.2) with P-values (per allele OR, 95%CI) of 6.49×10⁻⁹ (1.15, 1.10–1.20), 1.45×10⁻⁸ (1.13, 1.08–1.18), and 7.14×10⁻⁷ (1.13, 1.08–1.19), respectively, in combined analyses of 9,794 cases and 14,615 controls. Our study provides strong evidence for a novel T2D susceptibility locus at 13q31.1 and the presence of new independent risk variants near regions (10p13 and 15q22.2) reported by previous GWAS.     

Genetic Variability in Natural Populations of Eurytopic Ostracod Candona neglecta Sars

An electrophoretic survey of allozyme variation was conducted in four, highly polymorphic loci on nine populations of ostracod Candona neglecta Sars from three different environments: the profundal of post-glacial lakes, deep muddy bottom of the Baltic Sea and small astatic water bodies. The results suggest lack of genetic isolation between populations from lake profundal and the Baltic Sea. On the other hand a very distinct founder effect can be noted in the case of young, isolated populations from small astatic basins. It is suggested that a population inhabiting a large lake may be genetically subdivided due to differentiated eutrophication.     

Search for A-factor-Dependent Variants in Actinomycete Populations

A study of 28 nocardia-like, asporogenous, and oligosporous spontaneous morphological variants belonging to 23 species of streptomycetes revealed five strains producing regulators of the A-factor group. Streptomyces griseus 1439, which forms aerial mycelium and spores only in the presence of exogenous A-factor, was used as the test strain. Among the 28 spontaneous variants, three new A-factor-dependent strains were revealed, which represented the species Streptomyces griseus, S. citreofluorescens, and S. viridovulgaris subsp. albomarinus. These weakly differentiated variants did not produce A-factor and behaved as its recipients, responding by changes in their morphological characteristics at a concentration of this regulator in the medium of 0.01 g/ml or higher. The original collection strains in whose populations the variants were selected produced substances of the A-factor group. The A-factor-dependent variants differed in the level of the regulator required for maximal expression of the morphological characteristics: it was necessary to introduce the A-factor at a concentration of 1 g/ml for S. citreofluorescens and S. viridovulgaris subsp. albomarinus and at 10 g/ml for S. griseus.     

Use of the IRAP Marker to Study Genetic Variability in Pseudocercospora fijiensis Populations

Pseudocercospora fijiensis is the etiological agent of black Sigatoka, which is currently considered as one of the most destructive banana diseases in all locations where it occurs. It is estimated that a large portion of the P. fijiensis genome consists of transposable elements, which allows researchers to use transposon-based molecular markers in the analysis of genetic variability in populations of this pathogen. In this context, the inter-retrotransposon-amplified polymorphism (IRAP) was used to study the genetic variability in P. fijiensis populations from different hosts and different geographical origins in Brazil. A total of 22 loci were amplified and 77.3 % showed a polymorphism. Cluster analysis revealed two major groups in Brazil. The observed genetic diversity (H _E) was 0.22, and through molecular analysis of variance, it was determined that the greatest genetic variability occurs within populations. The discriminant analysis of principal components revealed no structuring related to the geographical origin of culture of the host. The IRAP-based marker system is a suitable tool for the study of genetic variability in P. fijiensis.     

Association of Genetic Variants with Primary Angle Closure Glaucoma in Two Different Populations

Purpose

A recent large genome-wide association study (GWAS) identified multiple variants associated with primary angle-closure glaucoma (PACG). The present study investigated the role of these variants in two cohorts with PACG recruited from Australia and Nepal.

Method

Patients with PACG and appropriate controls were recruited from eye clinics in Australia (n = 232 cases and n = 288 controls) and Nepal (n = 106 cases and 204 controls). Single nucleotide polymorphisms (SNPs) rs3753841 (COL11A1), rs1015213 (located between PCMTD1 and ST18), rs11024102 (PLEKHA7), and rs3788317 (TXNRD2) were selected and genotyped on the Sequenom. Analyses were conducted using PLINK and METAL.

Results

After adjustment for age and sex, SNP rs3753841 was found to be significantly associated with PACG in the Australian cohort (p = 0.017; OR = 1.34). SNPs rs1015213 (p = 0.014; OR 2.35) and rs11024102 (p = 0.039; OR 1.43) were significantly associated with the disease development in the Nepalese cohort. None of these SNPs survived Bonferroni correction (p = 0.05/4 = 0.013). However, in the combined analysis, of both cohorts, rs3753841 and rs1015213 showed significant association with p-values of 0.009 and 0.004, respectively both surviving Bonferroni correction. SNP rs11024102 showed suggestive association with PACG (p-value 0.035) and no association was found with rs3788317.

Conclusion

The present results support the initial GWAS findings, and confirm the SNP’s contribution to PACG. This is the first study to investigate these loci in both Australian Caucasian and Nepalese populations.     

Comparison of Genetic Variants in Cancer-Related Genes between Chinese Hui and Han Populations

Background

The Chinese Hui population, as the second largest minority ethnic group in China, may have a different genetic background from Han people because of its unique demographic history. In this study, we aimed to identify genetic differences between Han and Hui Chinese from the Ningxia region of China by comparing eighteen single nucleotide polymorphisms in cancer-related genes.

Methods

DNA samples were collected from 99 Hui and 145 Han people from the Ningxia Hui Autonomous Region in China, and SNPs were detected using an improved multiplex ligase detection reaction method. Genotyping data from six 1000 Genomes Project population samples (99 Utah residents with northern and western European ancestry (CEU), 107 Toscani in Italy (TSI), 108 Yoruba in Ibadan (YRI), 61 of African ancestry in the southwestern US (ASW), 103 Han Chinese in Beijing (CHB), and 104 Japanese in Tokyo (JPT)) were also included in this study. Differences in the distribution of alleles among the populations were assessed using χ² tests, and F_ST was used to measure the degree of population differentiation.

Results

We found that the genetic diversity of many SNPs in cancer-related genes in the Hui Chinese in Ningxia was different from that in the Han Chinese in Ningxia. For example, the allele frequencies of four SNPs (rs13361707, rs2274223, rs465498, and rs753955) showed different genetic distributions (p<0.05) between Chinese Ningxia Han and Chinese Ningxia Hui. Five SNPs (rs730506, rs13361707, rs2274223, rs465498 and rs753955) had different F_ST values (F_ST >0.000) between the Hui and Han populations.

Conclusions

These results suggest that some SNPs associated with cancer-related genes vary among different Chinese ethnic groups. We suggest that population differences should be carefully considered in evaluating cancer risk and prognosis as well as the efficacy of cancer therapy.     

Identification of msp1 Gene Variants in Populations of Meloidogyne incognita Using PCR-DGGE

Effectors of root-knot nematodes are essential for parasitism and prone to recognition by adapted variants of the host plants. This selective pressure initiates hypervariability of effector genes. Diversity of the gene variants within nematode populations might correlate with host preferences. In this study we developed a method to compare the distribution of variants of the effector gene msp1 among populations of Meloidogyne incognita. Primers were designed to amplify a 234-bp fragment of msp1. Sequencing of cloned PCR products revealed five msp1 variants from seven populations that were distinguishable in their reproduction on five host plants. A protocol for denaturing gradient gel electrophoresis (DGGE) was developed to separate these msp1 variants. DGGE for replicated pools of juveniles from the seven populations revealed ten variants of msp1. A correlation between the presence of a particular gene variant and the reproductive potential on particular hosts was not evident. Especially race 3 showed substantial variation within the population. DGGE fingerprints of msp1 tended to cluster the populations according to their reproduction rate on pepper. The developed method could be useful for analyzing population heterogeneity and epidemiology of M. incognita.     

Causal Inference for Spatial Constancy across Saccades

Our ability to interact with the environment hinges on creating a stable visual world despite the continuous changes in retinal input. To achieve visual stability, the brain must distinguish the retinal image shifts caused by eye movements and shifts due to movements of the visual scene. This process appears not to be flawless: during saccades, we often fail to detect whether visual objects remain stable or move, which is called saccadic suppression of displacement (SSD). How does the brain evaluate the memorized information of the presaccadic scene and the actual visual feedback of the postsaccadic visual scene in the computations for visual stability? Using a SSD task, we test how participants localize the presaccadic position of the fixation target, the saccade target or a peripheral non-foveated target that was displaced parallel or orthogonal during a horizontal saccade, and subsequently viewed for three different durations. Results showed different localization errors of the three targets, depending on the viewing time of the postsaccadic stimulus and its spatial separation from the presaccadic location. We modeled the data through a Bayesian causal inference mechanism, in which at the trial level an optimal mixing of two possible strategies, integration vs. separation of the presaccadic memory and the postsaccadic sensory signals, is applied. Fits of this model generally outperformed other plausible decision strategies for producing SSD. Our findings suggest that humans exploit a Bayesian inference process with two causal structures to mediate visual stability.     

Genetic Variability in Markers of HLA-C Expression in Two Diverse South African Populations

An insertion-deletion (indel) polymorphism within the 3′ untranslated region (UTR) of HLA-C has been shown to be involved in the regulation of HLA-C expression. Individuals who carry a deletion at this position exhibit increased HLA-C expression, which associates with lower viral set point in HIV-1 infected individuals. This 263 indel (rs67384697) is reported to be in strong linkage disequilibrium (LD) with a single nucleotide polymorphism (SNP) 35 kilobases upstream of HLA-C (-35T/C; rs9264942) in Caucasian individuals, making this SNP a potential marker for both HLA-C expression and HIV-1 disease progression. We therefore examined genetic variation within the HLA-C 3′ UTR of 265 Black and Caucasian South Africans by direct sequencing and identified haplotypes encompassing the 263 indel and another indel at position 230 in both populations. Concomitant evaluation of variability at the −35 SNP revealed this polymorphism to be an inappropriate marker for the 263 indel in these populations. These findings provide important insights into genetic variability within the regulatory regions of HLA-C that have potential implications for our understanding of the regulation of HLA-C expression and its impact on HIV-1 disease progression.     

Genetic Variability in European Populations of an Invasive American Crayfish: Preliminary Results

Since 1973, the red swamp crayfish, Procambarus clarkii, native to south-central United States and northeastern Mexico, has spread throughout Europe. Here, we surveyed the genetic variability of five European populations of the species using RAPD markers. Genetic variation was found to be so high as to uniquely fingerprint most of the surveyed individuals. Analysis of molecular variance (AMOVA) of the RAPD markers showed that 1) a large part of the genetic variation can be attributed to the differentiation among localities, and 2) the differentiation was mainly due to the separation of the samples from Louisiana with respect to the European set. A hypothesis emerged in which subsequent introductions of crayfish from different sources were performed. This hypothesis might explain the high genetic diversity observed within each population and the genetic differentiation among populations, as the result, respectively, of the introduction of different sets of crayfish and the casual bias of introductions. Although preliminary, our results suggest that RAPDs could be helpful in providing information about human-mediated introduced populations.