Replication Study in a Japanese Population of Six Susceptibility Loci for Type 2 Diabetes Originally Identified by a Transethnic Meta-Analysis of Genome-Wide Association Studies

A transcriptome-wide association study of uterine fibroids to identify potential genetic markers and toxic chemicals

Uterine fibroid is one of the most prevalent benign tumors in women, with high socioeconomic costs. Although genome-wide association studies (GWAS) have identified several loci associated with uterine fibroid risks, they could not successfully interpret the biological effects of genomic variants at the gene expression levels. To prioritize uterine fibroid susceptibility genes that are biologically interpretable, we conducted a transcriptome-wide association study (TWAS) by integrating GWAS data of uterine fibroid and expression quantitative loci data. We identified nine significant TWAS genes including two novel genes, RP11-282O18.3 and KBTBD7, which may be causal genes for uterine fibroid. We conducted functional enrichment network analyses using the TWAS results to investigate the biological pathways in which the overall TWAS genes were involved. The results demonstrated the immune system process to be a key pathway in uterine fibroid pathogenesis. Finally, we carried out chemical–gene interaction analyses using the TWAS results and the comparative toxicogenomics database to determine the potential risk chemicals for uterine fibroid. We identified five toxic chemicals that were significantly associated with uterine fibroid TWAS genes, suggesting that they may be implicated in the pathogenesis of uterine fibroid. In this study, we performed an integrative analysis covering the broad application of bioinformatics approaches. Our study may provide a deeper understanding of uterine fibroid etiologies and informative notifications about potential risk chemicals for uterine fibroid.

Evolutionary forces in diabetes and hypertension pathogenesis in Africans

Rates of type 2 diabetes (T2D) and hypertension are increasing rapidly in urbanizing sub-Saharan Africa (SSA). While lifestyle factors drive the increases in T2D and hypertension prevalence, evidence across populations shows that genetic variation, which is driven by evolutionary forces including a natural selection that shaped the human genome, also plays a role. Here we report the evidence for the effect of selection in African genomes on mechanisms underlying T2D and hypertension, including energy metabolism, adipose tissue biology, insulin action and salt retention. Selection effects found for variants in genes PPARA and TCF7L2 may have enabled Africans to respond to nutritional challenges by altering carbohydrate and lipid metabolism. Likewise, African-ancestry-specific characteristics of adipose tissue biology (low visceral adipose tissue [VAT], high intermuscular adipose tissue and a strong association between VAT and adiponectin) may have been selected for in response to nutritional and infectious disease challenges in the African environment. Evidence for selection effects on insulin action, including insulin resistance and secretion, has been found for several genes including MPHOSPH9, TMEM127, ZRANB3 and MC3R. These effects may have been historically adaptive in critical conditions, such as famine and inflammation. A strong correlation between hypertension susceptibility variants and latitude supports the hypothesis of selection for salt retention mechanisms in warm, humid climates. Nevertheless, adaptive genomics studies in African populations are scarce. More work is needed, particularly genomics studies covering the wide diversity of African populations in SSA and Africans in diaspora, as well as further functional assessment of established risk loci.

Accurate Diagnosis of Small Ruminant Lentivirus Infection Is Needed for Selection of Resistant Sheep through TMEM154 E35K Genotyping

Small ruminant lentiviruses (SRLV) cause an incurable multiorganic disease widely spread in sheep and goats that disturbs animal welfare and production. In the absence of a vaccine, control measures have been traditionally based on early diagnosis and breeding with virus-inactivated colostrum with segregation of seropositive animals. However, antigenic heterogeneity, poor antibody production due to low viral load, and single strain design of most available ELISA, pose a threat to SRLV diagnosis. Genome-wide association studies have described TMEM154 E35K polymorphism as a good genetic marker for selection of resistant animals in some American and European breeds. In this study, a multitargeted serological and virological screening of more than 500 animals from four different breeds (latxa, raza Navarra, assaf, and churra) attending to SRLV infection status was performed. Then, animals were genotyped to characterize TMEM154 E35K polymorphism. ELISA procedures, individually considered, only identified a proportion of the seropositive animals, and PCR detected a fraction of seronegative animals, globally offering different animal classifications according to SRLV infection status. TMEM154 allele frequency differed substantially among breeds and a positive association between seroprevalence and TMEM154 genotype was found only in one breed. Selection based on TMEM154 may be suitable for specific ovine breeds or SRLV strains, however generalization to the whole SRLV genetic spectrum, ovine breeds, or epidemiological situation may need further validation.

Accurate Diagnosis of Small Ruminant Lentivirus Infection Is Needed for Selection of Resistant Sheep through TMEM154 E35K Genotyping

Small ruminant lentiviruses (SRLV) cause an incurable multiorganic disease widely spread in sheep and goats that disturbs animal welfare and production. In the absence of a vaccine, control measures have been traditionally based on early diagnosis and breeding with virus-inactivated colostrum with segregation of seropositive animals. However, antigenic heterogeneity, poor antibody production due to low viral load, and single strain design of most available ELISA, pose a threat to SRLV diagnosis. Genome-wide association studies have described TMEM154 E35K polymorphism as a good genetic marker for selection of resistant animals in some American and European breeds. In this study, a multitargeted serological and virological screening of more than 500 animals from four different breeds (latxa, raza Navarra, assaf, and churra) attending to SRLV infection status was performed. Then, animals were genotyped to characterize TMEM154 E35K polymorphism. ELISA procedures, individually considered, only identified a proportion of the seropositive animals, and PCR detected a fraction of seronegative animals, globally offering different animal classifications according to SRLV infection status. TMEM154 allele frequency differed substantially among breeds and a positive association between seroprevalence and TMEM154 genotype was found only in one breed. Selection based on TMEM154 may be suitable for specific ovine breeds or SRLV strains, however generalization to the whole SRLV genetic spectrum, ovine breeds, or epidemiological situation may need further validation.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies

Background

Our longitudinal exome‐wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single‐nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow‐up period, 5 years) in 4884–6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease‐susceptibility loci may help us uncover the pathogenesis of the related complex disorders.

Methods

In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease‐susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project.

Results

Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out‐of‐Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non‐Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease‐associated loci may be affected by demographic events.

Conclusion

Our findings indicate that derived allele frequencies of nine disease‐associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome‐wide association studies largely increased in non‐Africans after Out‐of‐Africa.

Genetic and Environmental Predictors of Adolescent PTSD Symptom Trajectories Following a Natural Disaster

: Genes, environmental factors, and their interplay affect posttrauma symptoms. Although environmental predictors of the longitudinal course of posttraumatic stress disorder (PTSD) symptoms are documented, there remains a need to incorporate genetic risk into these models, especially in youth who are underrepresented in genetic studies. In an epidemiologic sample tornado-exposed adolescents (n = 707, 51% female, Mage = 14.54 years), trajectories of PTSD symptoms were examined at baseline and at 4-months and 12-months following baseline. This study aimed to determine if rare genetic variation in genes previously found in the sample to be related to PTSD diagnosis at baseline (MPHOSPH9, LGALS13, SLC2A2), environmental factors (disaster severity, social support), or their interplay were associated with symptom trajectories. A series of mixed effects models were conducted. Symptoms decreased over the three time points. Elevated tornado severity was associated with elevated baseline symptoms. Elevated recreational support was associated with lower baseline symptoms and attenuated improvement over time. Greater LGLAS13 variants attenuated symptom improvement over time. An interaction between MPHOSPH9 variants and tornado severity was associated with elevated baseline symptoms, but not change over time. Findings suggest the importance of rare genetic variation and environmental factors on the longitudinal course of PTSD symptoms following natural disaster trauma exposure.

A preliminary investigation of rare variants associated with genetic risk for PTSD in a natural disaster-exposed adolescent sample

Background: Posttraumatic stress disorder (PTSD) involves a complex interaction of biological, psychological, and social factors. Numerous studies have demonstrated genetic variation associated with the development of PTSD, primarily in adults. However, the contribution of low frequency and rare genetic variants to PTSD is unknown to date. Moreover, there is limited work on genetic risk for PTSD in child and adolescent populations. Objective: This preliminary study aimed to identify the low frequency and rare genetic variation that contributes to PTSD using an exome array. Method: This post-disaster, adolescent sample (n = 707, 51% females, M _age = 14.54) was assessed for PTSD diagnosis and symptom count following tornado exposure. Results: Gene-based models, covarying for ancestry principal components, age, sex, tornado severity, and previous trauma identified variants in four genes associated with diagnosis and 276 genes associated with symptom count (at p _adj < .001). Functional class analyses suggested an association with variants in the nonsense class (nonsynonymous variant that results in truncation of, and usually non-functional, protein) with both outcomes. An exploratory gene network pathway analysis showed a great number of significant genes involved in brain and immune function, illustrating the usefulness of downstream examination of gene-based findings that may point to relevant biological processes. Conclusions: While further investigation in larger samples is warranted, findings align with extant PTSD literature that has identified variants associated with biological conditions such as immune function.

Exposure to maternal obesity alters gene expression in the preimplantation ovine conceptus

Background

Embryonic and fetal exposure to maternal obesity causes several maladaptive morphological and epigenetic changes in exposed offspring. The timing of these events is unclear, but changes can be observed even after a short exposure to maternal obesity around the time of conception. The hypothesis of this work is that maternal obesity influences the ovine preimplantation conceptus early in pregnancy, and this exposure will affect gene expression in embryonic and extraembryonic tissues.

Results

Obese and lean ewe groups were established by overfeeding or normal feeding, respectively. Ewes were then bred to genetically similar rams. Conceptuses were collected at day 14 of gestation. Morphological assessments were made, conceptuses were sexed by genomic PCR analysis, and samples underwent RNA-sequencing analysis. While no obvious morphological differences existed between conceptuses, differentially expressed genes (≥ 2-fold; ≥ 0.2 RPKM; ≤ 0.05 FDR) were detected based on maternal obesity exposure (n = 21). Also, differential effects of maternal obesity were noted on each conceptus sex (n = 347). A large portion of differentially expressed genes were associated with embryogenesis and placental development.

Conclusions

Findings reveal that the preimplantation ovine conceptus genome responds to maternal obesity in a sex-dependent manner. The sexual dimorphism in response to the maternal environment coupled with changes in placental gene expression may explain aberrations in phenotype observed in offspring derived from obese females.

Electronic supplementary material

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

Two novel susceptibility loci for type 2 diabetes mellitus identified by longitudinal exome-wide association studies in a Japanese population

Recent genome-wide association studies identified genetic variants that confer susceptibility to type 2 diabetes mellitus (T2DM). However, few longitudinal genome-wide association studies of this metabolic disorder have been reported to date. Therefore, we performed a longitudinal exome-wide association study of T2DM, using 24,579 single nucleotide polymorphisms (SNPs) and repeated measurements from 6022 Japanese individuals. The generalized estimating equation model was applied to test relations of SNPs to three T2DM-related parameters: prevalence of T2DM, fasting plasma glucose level, and blood glycosylated hemoglobin content. Three SNPs that passed quality control were significantly (P<2.26×10^-7) associated with two of the three T2DM-related parameters in additive and recessive models. Of the three SNPs, rs6414624 in EVC and rs78338345 in GGA3 were novel susceptibility loci for T2DM. In the present study, the SNP of GGA3 was predicted to be a genetic variant whose minor allele frequency has recently increased in East Asia.