SHEsisPlus, a toolset for genetic studies on polyploid species

Polymorphism of the PPARD Gene and Dynamic Balance Performance in Han Chinese Children

BACKGROUND

Athletic performances are complex traits with heritability of ~66%. Dynamic balance is one of the most important athletic performances, and there has been little studies for it in sports genomics. The candidate PPARD gene was reported to be able to affect muscle development for balance predisposition and influence the athletic performance including skiing triumph in the Caucasian population. This study aims to investigate whether the PPARD gene is a susceptibility gene for dynamic balance performance in Han Chinese children.

RESULTS

A total 2244 children were recruited and their balance beam performances were measured. Five polymorphisms in the PPARD gene were genotyped through the MassARRAY Sequenom platform. Rs2016520 exerted significant association with dynamic balance performance (minor allele C, P = 0.015, Pcorrected < 0.05) and was affirmed in a meta-analysis by combining previously reported Caucasian cohorts (OR = 1.57, 95% CI = [1.30, 1.91], P < 10 -5) . Another polymorphism, rs2267668, was also significantly associated with dynamic balance performance (minor allele G, P = 0.015, Pcorrected < 0.05). In the dichotomous study, 321 cases (61% boys and 39% girls) and 370 controls (49% boys and 51% girls) in our samples were selected as representatives, and the thresholds were the mean velocity (0.737 m/s) ± standard deviation (0.264 m/s), in which rs2016520-C and rs2267668-G still remained significant (CI =1.41 [1.11~1.79], P = 0.004, Pcorrected < 0.016; CI =1.45 [1.14~1.86], P = 0.002, Pcorrected < 0.016). In different genders, consistent OR direction was observed for each variant.

CONCLUSIONS

Our results suggested that the PPARD gene is associated with dynamic balance performance of human being, and further studies to reveal its etiology is strongly suggested.

pSBVB: A Versatile Simulation Tool To Evaluate Genomic Selection in Polyploid Species

Genomic Selection (GS) is the procedure whereby molecular information is used to predict complex phenotypes and it is standard in many animal and plant breeding schemes. However, only a small number of studies have been reported in horticultural crops, and in polyploid species in particular. In this paper, we have developed a versatile forward simulation tool, called polyploid Sequence Based Virtual Breeding (pSBVB), to evaluate GS strategies in polyploids; pSBVB is an efficient gene dropping software that can simulate any number of complex phenotypes, allowing a very flexible modeling of phenotypes suited to polyploids. As input, it takes genotype data from the founder population, which can vary from single nucleotide polymorphisms (SNP) chips up to sequence, a list of causal variants for every trait and their heritabilities, and the pedigree. Recombination rates between homeologous chromosomes can be specified, so that both allo- and autopolyploid species can be considered. The program outputs phenotype and genotype data for all individuals in the pedigree. Optionally, it can produce several genomic relationship matrices that consider exact or approximate genotype values. pSBVB can therefore be used to evaluate GS strategies in polyploid species (say varying SNP density, genetic architecture or population size, among other factors), or to optimize experimental designs for association studies. We illustrate pSBVB with SNP data from tetraploid potato and partial sequence data from octoploid strawberry, and we show that GS is a promising breeding strategy for polyploid species but that the actual advantage critically depends on the underlying genetic architecture. Source code, examples and a complete manual are freely available in GitHub https://github.com/lauzingaretti/pSBVB.

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a 'sweet' crop with 'bitter' genomes

Sugarcane (Saccharum spp.) is a highly energy-efficient crop primarily for sugar and bio-ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high-quality genome-wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in-depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.

Common variants in the SLC28A2 gene are associated with serum uric acid level and hyperuricemia and gout in Han Chinese

Background

Serum uric acid (SUA), hyperuricemia (HUA) and gout are complex traits with relatively high heritability. This study aims to identify whether a candidate gene, SLC28A2, exerts susceptibility for SUA fluctuation and incidence of HUA and gout in the Han Chinese population.

Results

Three sample sets of 1376 gout patients, 1290 long-term HUA subjects (no gout attack) and 1349 normouricemic controls were recruited for this study. Eight polymorphisms in the SLC28A2 gene were genotyped using the ligase detection reaction-polymerase chain reaction (LDR-PCR) technology. Rs16941238 showed the most significant associations with SUA level (minor allele “A”, BETA = − 13.84 μmol/L, P = 0.0041, P_perm = 0.0042) and HUA (OR = 0.7734, P = 0.0033, P_perm = 0.0020), but not with gout (OR = 0.8801, P = 0.1315, P_perm = 0.1491). Rs2271437 was significantly associated with gout (minor allele “G”, OR = 1.387, P = 0.0277, P_perm = 0.0288), and was further confirmed in the meta-analysis with the previously published gout GWAS dataset (OR = 1.3221, P = 0.0089). Each variant basically conferred consistent OR direction on gout and HUA, compared with the normouricemic control.

Conclusions

Our findings support the associations of the SLC28A2 gene with the SUA level, the HUA phenotype and gout in Han Chinese.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0078-0) contains supplementary material, which is available to authorized users.

An association study of dopaminergic (DRD2) and serotoninergic (5-HT2) gene polymorphism and schizophrenia in a North Indian population

The aim of this study was to investigate whether single nucleotide polymorphisms (SNPs) in DRD2 and 5-HT2 receptor genes are associated with schizophrenia in North Indian population. Four hundred forty-three patients who met ICD10-DCR criteria for schizophrenia were enrolled from six participating centers along with 443 genetically related healthy subjects and 150 genetically unrelated healthy participants. A total of 7 gene polymorphisms from DRD2 (rs1800497, rs1079597, rs1800498, rs1801028) and 5-HT2 A (rs6313, rs6311, rs6305) were genotyped for their association with schizophrenia. No significant difference was found in frequency of various genotypes and alleles of the studied markers for DRD2 and 5-HT2 A genes between the cases and their genetic controls. However, significant differences were noted for rs1079597 genotype (Taq1B; p = 0.039) and its allele frequencies (p = 0.029) in persons with schizophrenia and the unrelated healthy controls. The DRD2 (Taq1 A-B-D) and 5-HT2 A (rs6311-rs6313-rs6305) haplotype frequencies differed significantly for A2B1D2 [p = 0.038; OR = 0.685 (95%CI = 0.479-0.981)] and ACC [p = 0.001; OR = 0.621 (95%CI = 0.461-0.838)] for the cases vs genetically related healthy controls. Similarly, significant difference was observed for the frequencies of GCC [p = 0.006; OR = 0.692 (95%CI = 0.532-0.900)] and ACC [p < 0.001; OR = 3.622 (95%CI = 1.73-7.585)] in the cases and unrelated healthy controls. Unlike previous research from India as well as abroad, the predominance of B1 allele of rs1079597 in patients with schizophrenia and absence of Cys311 in all study participants is a salient difference. Concluding, the B2 allele of rs1079597 may increase the risk of schizophrenia while the A2B1D2 haplotype may be protective in North Indian population.

ACTN3 is associated with children's physical fitness in Han Chinese

The ACTN3 gene locates on 11q13-q14 and encodes the α-actinin-3 protein, which is only expressed in human skeletal muscle and influenced muscle function and metabolism. The previous studies reported that SNP rs1815739 is associated with elite power athletes' performance. In this study, we investigated the association between five SNPs within the ACTN3 gene and Chinese children physical fitness. We recruited 2244 Han Chinese children participants, and measured their 25-m run, stand broad jump, 10-m shuttle run, handgrip, BMI (calculated by weight and height) data. SNPs rs1671064, rs2275998, rs2290463, rs10791881, and rs1815739 of ACTN3 gene were genotyped and analyzed in five physical fitness data. QTL analysis on genotype and physical fitness data was carried out in all samples. Furthermore, a dichotomous division of samples into an overweight group (543) and a normal group (1701) was used for an association study of overweight. In the QTL analysis, we found rs2290463 was significantly associated with stand broad jump (corrected P value = 0.009, beta = 2.692). After added age and gender as covariates in the regression test, the association became more significant (P value = 5.80 × 10^- 5, corrected P value = 4.06 × 10^- 4); when we used BMI as a covariate, the association still existed (P value = 4.65 × 10^- 4, corrected P value = 0.001). In the association study of overweight, rs2275998 was found to be significant (OR, 95% CI = 0.733 [0.6-0.895]; P_allele = 0.011, P_genotype = 0.024) after the Bonferroni correction, and the association did not change much after a further correction for gender, age, and stand broad jump performance. Our results showed that common variants in ACTN3 are significantly associated with both stand broad jump performance and overweight in Han Chinese children.

Interactions Among Polymorphisms of Susceptibility Loci for Alzheimer's Disease or Depressive Disorder

Background

Several genetic susceptibility loci for major depressive disorder (MDD) or Alzheimer’s disease (AD) have been described. Interactions among polymorphisms are thought to explain the differences between low- and high-risk groups. We tested for the contribution of interactions between multiple functional polymorphisms in the risk of MDD or AD.

Material/Methods

A genetic association case-control study was performed in 68 MDD cases, 84 AD cases (35 of them with comorbid depression), and 90 controls. The contribution of 7 polymorphisms from 5 genes (APOE, HSPA1A, SLC6A4, HTR2A, and BDNF) related to risk of MDD or AD development was analyzed.

Results

Significant associations were found between MDD and interactions among polymorphisms in HSPA1A, SLC6A4, and BDNF or HSPA1A, BDNF, and APOE genes. For polymorphisms in the APOE gene in AD, significant differences were confirmed on the distributions of alleles and genotype rates compared to the control or MDD. Increased probability of comorbid depression was found in patients with AD who do not carry the ɛ4 allele of APOE.

Conclusions

Assessment of the interactions among polymorphisms of susceptibility loci in both MDD and AD confirmed a synergistic effect of genetic factors influencing inflammatory, serotonergic, and neurotrophic pathways at these heterogenous complex diseases. The effect of interactions was greater in MDD than in AD. A presence of the ɛ4 allele was confirmed as a genetic susceptibility factor in AD. Our findings indicate a role of APOE genotype in onset of comorbid depression in a subgroup of patients with AD who are not carriers of the APOE ɛ4 allele.

Tools for Genetic Studies in Experimental Populations of Polyploids

Polyploid organisms carry more than two copies of each chromosome, a condition rarely tolerated in animals but which occurs relatively frequently in the plant kingdom. One of the principal challenges faced by polyploid organisms is to evolve stable meiotic mechanisms to faithfully transmit genetic information to the next generation upon which the study of inheritance is based. In this review we look at the tools available to the research community to better understand polyploid inheritance, many of which have only recently been developed. Most of these tools are intended for experimental populations (rather than natural populations), facilitating genomics-assisted crop improvement and plant breeding. This is hardly surprising given that a large proportion of domesticated plant species are polyploid. We focus on three main areas: (1) polyploid genotyping; (2) genetic and physical mapping; and (3) quantitative trait analysis and genomic selection. We also briefly review some miscellaneous topics such as the mode of inheritance and the availability of polyploid simulation software. The current polyploid analytic toolbox includes software for assigning marker genotypes (and in particular, estimating the dosage of marker alleles in the heterozygous condition), establishing chromosome-scale linkage phase among marker alleles, constructing (short-range) haplotypes, generating linkage maps, performing genome-wide association studies (GWAS) and quantitative trait locus (QTL) analyses, and simulating polyploid populations. These tools can also help elucidate the mode of inheritance (disomic, polysomic or a mixture of both as in segmental allopolyploids) or reveal whether double reduction and multivalent chromosomal pairing occur. An increasing number of polyploids (or associated diploids) are being sequenced, leading to publicly available reference genome assemblies. Much work remains in order to keep pace with developments in genomic technologies. However, such technologies also offer the promise of understanding polyploid genomes at a level which hitherto has remained elusive.

IL-1RN and IL-1β Polymorphism and ARV-Associated Hepatotoxicity

The severity of hepatic injury depends upon cytokines. Previous studies associated IL-1RN allele 2 with IL-1β production. Hence, we examined the association of IL-1 RN and IL-1β polymorphisms with ARV-associated hepatotoxicity. Genotyping of IL-1RN (VNTR), IL-1β (-511C/T) polymorphisms was done in 162 HIV-infected patients, 34 with ARV hepatotoxicity, 128 without hepatotoxicity, and 152 healthy controls using PCR and PCR-RFLP method. The haplotypes 1T and 2C enhanced the risk for severe hepatotoxicity (OR = 1.41, P = 0.25; OR = 1.67, P = 0.31). IL-1β-511TT genotype significantly represented among tobacco using HIV-infected individuals compared to nonusers (OR = 3.74, P = 0.05). IL-1β-511TT genotype among alcohol users increased the risk for hepatotoxicity (OR = 1.80, P = 0.90). IL-1β-511CT and -511TT genotypes overrepresented in alcohol using HIV-infected individuals (OR = 2.29, P = 0.27; OR = 2.64, P = 0.19). IL-RN 2/2 and 1/3 genotypes represented higher in nevirapine using hepatotoxicity patients (OR = 1.42, P = 0.64, OR = 8.79, P = 0.09). IL-1β-511CT and -511 TT genotypes among nevirapine users enhanced the risk for severe hepatotoxicity (OR = 4.29, P = 0.20; OR = 1.95, P = 0.56). IL-1β-511CT and -511TT genotypes were overrepresented in combined nevirapine and alcohol using HIV-infected individuals as compared to nevirapine users and alcohol nonusers (OR = 2.56, P = 0.26; OR = 2.84, P = 0.24). IL-1β-511TT genotype with tobacco, alcohol, and nevirapine usage revealed a trend of risk for the development of ARV-associated hepatotoxicity and its severity.

Association study of NDST3 gene for schizophrenia, bipolar disorder, major depressive disorder in the Han Chinese population

The NDST3 gene at 4q26 was a functional candidate gene for mental disorders. Recently, a novel genome-wide significant risk locus at chromosome 4q26 was identified and the top single nucleotide polymorphism rs11098403 in the vicinity of NDST3 gene was reported to confer risk of schizophrenia in Caucasian. Nevertheless, association between NDST3 gene polymorphisms and schizophrenia, bipolar disorder, or major depressive disorders has not been well studied in the Han Chinese population. To further investigate whether NDST3 is a risk gene for these mental disorders, we genotyped and analyzed eight tag SNPs (rs11098403, rs10857057, rs2389521, rs4833564, rs6837896, rs7689157, rs3817274, rs609512) covering NDST3 gene in 1,248 schizophrenia cases, 1,056 major depression cases, 1,344 bipolar disorder cases, and 1,248 controls of Chinese origin. However, there was no significant difference in allelic or genotypic frequency observed between each case group and healthy controls. Accordingly, our study does not support that the NDST3 gene plays a major role in schizophrenia, bipolar disorder, and major depressive disorder in the Han Chinese population.

Association of candidate genes with heading date in a diverse Dactylis glomerata population

Flowering occurs in response to cues from both temperature and photoperiod elicitors in cool-season, long-day forage grasses, and genes involved in sensing the elicitors and inducing downstream flowering responses have been associated with heading date and flowering time in perennial forage grasses as well as cereal grasses. In this study we test for association between orchardgrass (Dactylis glomerata L.) heading date and polymorphisms in the CONSTANS (DgCO1), FLOWERING TIME (DgFT1), a VRN1 like MADS-box (DgMADS), and PHOTOPERIOD (DgPPD1-like) containing genes. A diverse population of 150 genotypes was measured for heading date across three years, genotyped, and candidate genes sequenced. Although pairwise population kinship values were generally low, the genotypes fit into a two-group structure model. Linkage disequilibrium decayed rapidly, reaching r² levels below 0.2 within the 500bp of each gene. SNPs significantly associated with heading date were detected in equal-dose and tetraploid dosage models. The DgCO1 gene had the most significant polymorphisms and those with the largest effects, while DgMADS had several significant polymorphisms in its first intron with smaller effects. These polymorphisms can be used for further validation, selection, and development of breeding lines of orchardgrass.

Genetic risk between the CACNA1I gene and schizophrenia in Chinese Uygur population

Background

Schizophrenia (SCZ) is a common mental disorder with high heritability, and genetic factors play a major role in the pathogenesis. Recent researches indicated that the CACNA1I involved in calcium channels probably affect the potential pathogenesis of SCZ.

Results

In this study, we attempted to investigate whether the CACNA1I gene contributes the risk to SCZ in the Uighur Chinese population, and performed a case-control study involving 985 patient samples and 1218 normal controls to analyze nine SNPs within the CACNA1I gene. Among these sites, six SNPs were significantly associated with SCZ in the allele distribution: rs132575 (adjusted P_allele = 0.039, OR = 1.159), rs713860 (adjusted P_allele = 0.039, OR = 0.792), rs738168 (adjusted P_allele = 0.039, OR = 0.785), rs136805 (adjusted P_allele = 0.014, OR = 1.212), rs5757760 (adjusted P_allele = 0.042, OR = 0.873) and rs5750871 (adjusted P_allele = 0.039, OR = 0.859). In addition, two SNPs turned to be risk factors for SCZ not only in the allele distribution, but also in the genotype distribution: rs132575 (adjusted P_genotype = 0.037) and rs136805 (adjusted P_genotype = 0.037).

Conclusions

Overall, the present study provided evidence that significant association exists between the CACNA1I gene and SCZ in the Uighur Chinese population, subsequent validation of functional analysis and genetic association studies are needed to further extend this study.

Common variants in ZMIZ1 and near NGF confer risk for primary dysmenorrhoea

Primary dysmenorrhoea, defined as painful menstrual cramps in the absence of pelvic pathology, is a common problem in women of reproductive age. Its aetiology and pathophysiology remain largely unknown. Here we performed a two-stage genome-wide association study and subsequent replication study to identify genetic factors associated with primary dysmenorrhoea in a total of 6,770 Chinese individuals. Our analysis provided evidence of a significant (P<5 × 10⁻⁸) association at rs76518691 in the gene ZMIZ1 and at rs7523831 near NGF. ZMIZ1 has previously been associated with several autoimmune diseases, and NGF plays a key role in the generation of pain and hyperalgesia and has been associated with migraine. These findings provide future directions for research on susceptibility mechanisms for primary dysmenorrhoea. Furthermore, our genetic architecture analysis provides molecular support for the heritability and polygenic nature of this condition.

Primary dysmenorrhoea, the most common gynaecologic complaint, remains genetically and pathophysiologically elusive. Here, Li and colleagues identify common variants in ZMIZ1 and near NGF conferring risk for primary dysmenorrhoea using genome-wide association study in a Chinese population.