MIDAS: software for analysis and visualisation of interallelic disequilibrium between multiallelic markers

Methylenetetrahydrofolate reductase polymorphic variants C677T and A1298C in rectal cancer in Slavic population: significance for cancer risk and response to chemoradiotherapy

Background: Methylenetetrahydrofolate reductase (MTHFR) single nucleotide polymorphisms (SNPs) have been suggested as risk, prognostic, and predictive factors for colorectal cancer in various populations, but have not been validated so far. The aim of this study was to examine the association of MTHFR C677T (rs1801133) and A1298C (rs1801131) single nucleotide polymorphisms with the risk of rectal cancer as well as the response to neoadjuvant chemoradiotherapy (nCRT) based on 5-Fluorouracil (5-FU)/leucovorin (LV) in the locally advanced setting. Patients and methods: This case-control study included 119 healthy controls and 97 patients with locally advanced rectal cancer (LARC). For MTHFR genotyping, restriction fragment length polymorphism analysis (PCR-RFLP) was employed. Results: In silico analysis highlighted that SNPs C677T and A1298T correlate with MTHFR gene expression, and that gene expression profile correlates with cancer risk and stage. Using dominant and recessive models, it was found that the MTHFR 677CC vs. 677CT+677TT have increased risk of cancer development (odds ratio (OR): 2.27; 95% confidence interval (CI): 1.30-3.95, p = 0.002) as well as 677CC+677CT compared to 677TT (OR: 4.18, 95% CI: 1.16-14.99, p = 0.014). MTHFR 1298AA also shown increased risk for cancer development compared to 1298AC+1298CC (OR:2.0, 95% CI: 1.20-3.59, p = 0.035) Statistical analysis of combined genotypes highlighted the protective role of CT/AC combined genotype (OR: 3.15 95% CI: 1.576-6.279, p = 0.002) while the CC/AA genotype showed an increased risk for rectal cancer development (OR: 2.499, 95% CI: 1.246-5.081, p = 0.016) The carriers of the 677C/1298A haplotype had the highest risk for developing rectal cancer (OR: 1.74; 95% CI: 1.198-2.530, p = 0.002) while the 677T/1298C haplotype seems to provide a protective effect. (OR: 0.44; 95%CI 0.248-0.795, p = 0.003). No significant association with response to chemoradiotherapy was found. Conclusion: Our data point to MTHFR 667C allele and 1298A alleles as low-penetrance risk factors for rectal cancer in our population. To the best of our knowledge, this is the first study of this type performed on the Slavic population in the Western Balkan, as various population-based factors might also be significant our findings can be used for future meta-analyses and the construction of genetic cancer risk prediction panels.

Human NCR3 gene variants rs2736191 and rs11575837 alter longitudinal risk for development of pediatric malaria episodes and severe malarial anemia

BACKGROUND

Plasmodium falciparum malaria is a leading cause of pediatric morbidity and mortality in holoendemic transmission areas. Severe malarial anemia [SMA, hemoglobin (Hb) < 5.0 g/dL in children] is the most common clinical manifestation of severe malaria in such regions. Although innate immune response genes are known to influence the development of SMA, the role of natural killer (NK) cells in malaria pathogenesis remains largely undefined. As such, we examined the impact of genetic variation in the gene encoding a primary NK cell receptor, natural cytotoxicity-triggering receptor 3 (NCR3), on the occurrence of malaria and SMA episodes over time.

METHODS

Susceptibility to malaria, SMA, and all-cause mortality was determined in carriers of NCR3 genetic variants (i.e., rs2736191:C > G and rs11575837:C > T) and their haplotypes. The prospective observational study was conducted over a 36 mos. follow-up period in a cohort of children (n = 1,515, aged 1.9-40 mos.) residing in a holoendemic P. falciparum transmission region, Siaya, Kenya.

RESULTS

Poisson regression modeling, controlling for anemia-promoting covariates, revealed a significantly increased risk of malaria in carriers of the homozygous mutant allele genotype (TT) for rs11575837 after multiple test correction [Incidence rate ratio (IRR) = 1.540, 95% CI = 1.114-2.129, P = 0.009]. Increased risk of SMA was observed for rs2736191 in children who inherited the CG genotype (IRR = 1.269, 95% CI = 1.009-1.597, P = 0.041) and in the additive model (presence of 1 or 2 copies) (IRR = 1.198, 95% CI = 1.030-1.393, P = 0.019), but was not significant after multiple test correction. Modeling of the haplotypes revealed that the CC haplotype had a significant additive effect for protection against SMA (i.e., reduced risk for development of SMA) after multiple test correction (IRR = 0.823, 95% CI = 0.711-0.952, P = 0.009). Although increased susceptibility to SMA was present in carriers of the GC haplotype (IRR = 1.276, 95% CI = 1.030-1.581, P = 0.026) with an additive effect (IRR = 1.182, 95% CI = 1.018-1.372, P = 0.029), the results did not remain significant after multiple test correction. None of the NCR3 genotypes or haplotypes were associated with all-cause mortality.

CONCLUSIONS

Variation in NCR3 alters susceptibility to malaria and SMA during the acquisition of naturally-acquired malarial immunity. These results highlight the importance of NK cells in the innate immune response to malaria.

Polymorphisms Related to Iron Homeostasis Associate with Liver Disease in Chronic Hepatitis C

Dysregulation of iron metabolism in chronic hepatitis C (CHC) is a significant risk factor for hepatic cirrhosis and cancer. We studied if known genetic variants related to iron homeostasis associate with liver disease progression in CHC. Retrospective analysis included 249 CHC patients qualified for antiviral therapy between 2004 and 2014. For all patients, nine SNPs within HFE, TFR2, HDAC2, HDAC3, HDAC5, TMPRSS6, and CYBRD1 genes were genotyped. Expression of selected iron-related genes, was determined with qRT-PCR in 124 liver biopsies, and mRNA expression of co-inhibitory receptors (PD-1, Tim3, CTLA4) was measured in 79 liver samples. CYBRD1 rs884409, HDAC5 rs368328, TFR2 rs7385804, and TMPRSS6 rs855791 associated with histopathological changes in liver tissue at baseline. The combination of minor allele in HDAC3 rs976552 and CYBRD1 rs884409 linked with higher prevalence of hepatocellular carcinoma (HCC) during follow up (OR 8.1 CI 2.2-29.2; p = 0.001). Minor allele in HDAC3 rs976552 associated with lower hepatic expression of CTLA4. Tested polymorphisms related to iron homeostasis associate with histopathological changes in the liver. The presence of both HDAC3 rs976552 G and CYBRD1 rs884409 G alleles correlates with HCC occurrence, especially in the group of patients with elevated AST (>129 IU/L). rs976552 in HDAC3 could impact immunological processes associated with carcinogenesis in CHC.

Some New Aspects of Genetic Variability in Patients with Cutaneous T-Cell Lymphoma

AIM

Cutaneous T-cell lymphoma (CTCL) is a group of T-cell malignancies that develop in the skin. Though studied intensively, the etiology and pathogenesis of CTCL remain elusive. This study evaluated the survival of CTCL patients in the 1st Department of Dermatovenereology of St. Anne's University Hospital Brno. It included analysis of 19 polymorphic gene variants based on their expected involvement in CTCL severity.

MATERIAL AND METHODS

75 patients with CTCL, evaluated and treated at the 1st Department of Dermatovenereology of St. Anne´s University Hospital Brno, Faculty of Medicine, Masaryk University, were recruited for the study over the last 28 years (44 men and 31 women, average age 58 years, range 20-82 years). All patients were genotyped for 19 chosen gene polymorphisms by the conventional PCR method with restriction analysis. A multivariate Cox regression model was calculated to reveal genetic polymorphisms and other risk factors for survival.

RESULTS

The model identified MDR Ex21 2677 (rs2032582) as a significant genetic factor influencing the survival of the patients, with the T-allele playing a protective role. A multivariate stepwise Cox regression model confirmed the following as significant independent risk factors for overall survival: increased age at admission, clinical staging of the tumor, and male sex.

CONCLUSION

We showed that the TT genotype at position 2677 of the MDR1 gene exhibited statistically significant longer survival in CTCL patients. As such, the TT genotype of MDR1 confers a significant advantage for the CTCL patients who respond to treatment.

Nonsynonymous amino acid changes in the α-chain of complement component 5 influence longitudinal susceptibility to Plasmodium falciparum infections and severe malarial anemia in kenyan children

Background: Severe malarial anemia (SMA; Hb < 5.0 g/dl) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission regions such as western Kenya. Methods: We investigated the relationship between two novel complement component 5 (C5) missense mutations [rs17216529:C>T, p(Val145Ile) and rs17610:C>T, p(Ser1310Asn)] and longitudinal outcomes of malaria in a cohort of Kenyan children (under 60 mos, n = 1,546). Molecular modeling was used to investigate the impact of the amino acid transitions on the C5 protein structure. Results: Prediction of the wild-type and mutant C5 protein structures did not reveal major changes to the overall structure. However, based on the position of the variants, subtle differences could impact on the stability of C5b. The influence of the C5 genotypes/haplotypes on the number of malaria and SMA episodes over 36 months was determined by Poisson regression modeling. Genotypic analyses revealed that inheritance of the homozygous mutant (TT) for rs17216529:C>T enhanced the risk for both malaria (incidence rate ratio, IRR = 1.144, 95%CI: 1.059-1.236, p = 0.001) and SMA (IRR = 1.627, 95%CI: 1.201-2.204, p = 0.002). In the haplotypic model, carriers of TC had increased risk of malaria (IRR = 1.068, 95%CI: 1.017-1.122, p = 0.009), while carriers of both wild-type alleles (CC) were protected against SMA (IRR = 0.679, 95%CI: 0.542-0.850, p = 0.001). Conclusion: Collectively, these findings show that the selected C5 missense mutations influence the longitudinal risk of malaria and SMA in immune-naïve children exposed to holoendemic P. falciparum transmission through a mechanism that remains to be defined.

Complement component 3 mutations alter the longitudinal risk of pediatric malaria and severe malarial anemia

Severe malarial anemia (SMA) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission regions. To gain enhanced understanding of predisposing factors for SMA, we explored the relationship between complement component 3 (C3) missense mutations [rs2230199 (2307C>G, Arg>Gly¹⁰²) and rs11569534 (34420G>A, Gly>Asp¹²²⁴)], malaria, and SMA in a cohort of children (n = 1617 children) over 36 months of follow-up. Variants were selected based on their ability to impart amino acid substitutions that can alter the structure and function of C3. The 2307C>G mutation results in a basic to a polar residue change (Arg to Gly) at position 102 (β-chain) in the macroglobulin-1 (MG1) domain, while 34420G>A elicits a polar to acidic residue change (Gly to Asp) at position 1224 (α-chain) in the thioester-containing domain. After adjusting for multiple comparisons, longitudinal analyses revealed that inheritance of the homozygous mutant (GG) at 2307 enhanced the risk of SMA (RR = 2.142, 95%CI: 1.229-3.735, P = 0.007). The haplotype containing both wild-type alleles (CG) decreased the incident risk ratio of both malaria (RR = 0.897, 95%CI: 0.828-0.972, P = 0.008) and SMA (RR = 0.617, 95%CI: 0.448-0.848, P = 0.003). Malaria incident risk ratio was also reduced in carriers of the GG (Gly¹⁰²Gly¹²²⁴) haplotype (RR = 0.941, 95%CI: 0.888-0.997, P = 0.040). Collectively, inheritance of the missense mutations in MG1 and thioester-containing domain influence the longitudinal risk of malaria and SMA in children exposed to intense Plasmodium falciparum transmission.

Type 2 diabetes is associated with the MTNR1B gene, a genetic bridge between circadian rhythm and glucose metabolism, in a Turkish population

Type 2 diabetes (T2D) is a complicated public health problem in Turkey as well as worldwide. Genome-wide approaches have been guiding in very challenging situations, such as the elucidation of genetic variations underlying complex diseases such as T2D. Despite intensive studies worldwide, few studies have determined the genetic susceptibility to T2D in Turkish populations. In this study, we investigated the effect of genes that are strongly associated with T2D in genome-wide association (GWA) studies, including MTNR1B, CDKAL1, THADA, ADAMTS9 and ENPP1, on T2D and its characteristic traits in a Turkish population. In 824 nonobese individuals (454 T2D patients and 370 healthy individuals), prominent variants of these GWA genes were genotyped by real-time PCR using the LightSNiP Genotyping Assay System. The SNP rs1387153 C/T, which is located 28 kb upstream of the MTNR1B gene, was significantly associated with T2D and fasting blood glucose levels (P < 0.05). The intronic SNP rs10830963 C/G in the MTNR1B gene was not associated with T2D, but it was associated with fasting blood glucose, HbA1C and LDL levels (P < 0.05). The other important GWA loci investigated in our study were not found to be associated with T2D or its traits. Only the SNP rs1044498 (A/C variation) in the ENPP1 gene was determined to be related to fasting blood glucose (P < 0.05). Our study suggests, consistent with the literature, that the MTNR1B locus, which has a prominent role in glucose regulation, is associated with T2D development by affecting blood glucose levels in our population.

Effects of thymidylate synthase polymorphisms on toxicities associated with high-dose methotrexate in childhood acute lymphoblastic leukemia

BACKGROUND

High-dose methotrexate (HD- MTX) is the cornerstone of chemotherapy for acute lymphoblastic leukemia (ALL), and one of its target enzymes is Thymidylate Synthase (TYMS). We hypothesized that genetic polymorphisms of TYMS gene would be associated with MTX toxicity in ALL children.

METHODS

64 children with ALL were included in this study. Genotyping analysis was conducted on three common polymorphisms: tandem repeats in the promoter-enhancer region (VNTR), 6 bp ins/del (1494del6) in the 5'UTR, and rs2790 A > G in the 3'-untranslated region (3'-UTR). The association between genetic polymorphisms and MTX toxicity was studied.

RESULTS

Genetic polymorphism of TYMS was associated with hematological toxicities but not with non-hematological adverse events. A significant association between TYMS 1494del6 genotypes and incidence of neutropenia (ANC < 1700 mm³), infection and leukopenia was observed. Carriers of the dominant allele (Del) were 6 times more likely to develop neutropenia compared to minor genotype carriers (OR (95% CI) 6 (1.2-31.1); p = 0.04), and 4.2 times less likely to have infection, as compared to Ins/Ins carriers (OR  4.2, 95% CI (1.1-16); p = 0.04). Carriers of Del allele were 9.2 times more likely to develop grade 3 and 4 leukopenia, p = 0.02, 95% CI (1.1-75.6). Significant association was found between 28 bp VNTR and thrombocytopenia; (OR  3.3, 95% CI (1.1-10), p = 0.04). No significant association was found between TYMS rs2790 A > G genetic polymorphisms and MTX hematologic toxicities.

CONCLUSION

Genetic polymorphism of TYMS1494del6 may modulate susceptibility to MTX toxicity.

Association of Fetuin-B with Subclinical Atherosclerosis in Obese Chinese Adults

Aim: We aimed to explore the independent associations of serum Fetuin-B and common genetic variants in FETUB locus with subclinical atherosclerosis.

Methods: A cross-sectional study of 1,140 obese adults, who underwent serum Fetuin-B testing, hepatic ultrasonography scanning, genotyping on four tagging single nucleotide polymorphisms (SNPs) in FETUB locus and atherosclerosis detection, was conducted in Xiamen, China.

Results: Increasing tertiles of brachial ankle pulse wave velocity (ba-PWV) were significantly associated with higher prevalence of nonalcoholic fatty liver disease (NAFLD) (48.8%, 61.5%, and 70.5% for tertiles of 1–3, respectively, p < 0.001) and serum Fetuin-B (3.85 ± 1.39, 4.09 ± 1.40, and 4.27 ± 1.46 µg/ml, p = 0.047). Multivariable linear regression analyses with adjustment for potential confounding factors, even NAFLD per se, showed that serum Fetuin-B were significantly and positively associated with ba-PWV, with standardized regression coefficients (β) ranging from 0.055 to 0.075 (all p-values < 0.05) in different models. However, the significant relationship between serum Fetuin-B and ba-PWV disappeared with further adjustment for insulin resistance. Serum Fetuin-B was not significantly associated with ankle-brachial index (ABI). All genotypes of the four tested FETUB tagging SNPs were not significantly associated with either ba-PWV or ABI with adjustment for potential confounding factors

Conclusion: Serum Fetuin-B was positively associated with ba-PWV and may link liver fat accumulation to subclinical atherosclerosis via insulin resistance.

Estimating linkage disequilibrium from genotypes under Hardy-Weinberg equilibrium

Background

Measures of linkage disequilibrium (LD) play a key role in a wide range of applications from disease association to demographic history estimation. The true population LD cannot be measured directly and instead can only be inferred from genetic samples, which are unavoidably subject to measurement error. Previous studies of r² (a measure of LD), such as the bias due to finite sample size and its variance, were based on the special case that the true population-wise LD is zero. These results generally do not hold for non-zero \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {r}_{true}^2 $$\end{document}rtrue2 values, which are more common in real genetic data.

Results

This work generalises the estimation of r² to all levels of LD, and for both phased and unphased data. First, we provide new formulae for the effect of finite sample size on the observed r² values. Second, we find a new empirical formula for the variance of the observed r², equals to 2E[r²](1 − E[r²])/n, where n is the diploid sample size. Third, we propose a new routine, Constrained ML, a likelihood-based method to directly estimate haplotype frequencies and r² from diploid genotypes under Hardy-Weinberg Equilibrium. While serving the same purpose as the pre-existing Expectation-Maximisation algorithm, the new routine can have better convergence and is simpler to use. A new likelihood-ratio test is also introduced to test for the absence of a particular haplotype. Extensive simulations are run to support these findings.

Conclusion

Most inferences on LD will benefit from our new findings, from point and interval estimation to hypothesis testing. Genetic analyses utilising r² information will become more accurate as a result.

Genetic Differentiation in Hatchery and Stocked Populations of Sea Trout in the Southern Baltic: Selection Evidence at SNP Loci

The impacts and interactions between hatchery-bred fish and wild fish populations has been a topic of active investigation in recent decades. In some instances, the benefits of stocking can be overshadowed by negative effects such as genetic introgression with natural populations, loss of genetic diversity, and dilution of local adaptations. Methods that facilitate the identification of stocked fish enable us to estimate not only the effectiveness of stocking but also the level of natural reproduction and the degree of hybridization. The longest Baltic river, the Vistula, also has the second highest discharge. Historically, it hosted numerous populations of the anadromous form of brown trout (sea trout); however, dam construction has since interfered with and reduced spawning migration to a rate that is much lower than before. Reduced spawning has resulted in a population collapse and a negative flow-on effect on commercial catches. In response, Poland (along with many other Baltic countries) initiated an intensive stocking program which continues today and which sees the average annual release of 700,000 smolts. As a consequence, today’s main-river and inshore catches come from stock-enhanced populations. High-throughput single-nucleotide polymorphism (SNP) genotyping was performed on samples of sea trout from southern Baltic populations; results suggest that a significant portion of the sea trout catches in the Vistula mouth region have direct hatchery origin and indicate the presence of Pomeranian specimens. SNP loci identified as outliers indicate a potential selection pressure that may be related with effects of hatchery breeding and mixing with natural populations. The brown trout SNP array applied in this study showed high effectiveness not only for population differentiation, but more importantly, it emerged as a sensitive tool to provide evidence of detection selection.

Cyclooxygenase-2 haplotypes influence the longitudinal risk of malaria and severe malarial anemia in Kenyan children from a holoendemic transmission region

Cyclooxygenase-2 [(COX-2) or prostaglandin endoperoxide H2 synthase-2 (PTGS-2)] induces the production of prostaglandins as part of the host-immune response to infections. Although a number of studies have demonstrated the effects of COX-2 promoter variants on autoimmune and inflammatory diseases, their role in malaria remains undefined. As such, we investigated the relationship between four COX-2 promoter variants (COX-2 -512 C > T, -608 T > C, -765 G > C, and -1195 A > G) and susceptibility to malaria and severe malarial anemia (SMA) upon enrollment and longitudinally over a 36-month follow-up period. All-cause mortality was also explored. The investigation was carried out in children (n = 1081, age; 2-70 months) residing in a holoendemic Plasmodium falciparum transmission region of western Kenya. At enrollment, genotypes/haplotypes (controlling for anemia-promoting covariates) did not reveal any strong effects on susceptibility to either malaria or SMA. Longitudinal analyses showed decreased malaria episodes in children who inherited the -608 CC mutant allele (RR = 0.746, P = 1.811 × 10-4) and -512C/-608T/-765G/-1195G (CTGG) haplotype (RR = 0.856, P = 0.011), and increased risk in TTCA haplotype carriers (RR = 1.115, P = 0.026). Over the follow-up period, inheritance of the rare TTCG haplotype was associated with enhanced susceptibility to both malaria (RR = 1.608, P = 0.016) and SMA (RR = 5.714, P = 0.004), while carriage of the rare TTGG haplotype increased the risk of malaria (RR = 1.755, P = 0.007), SMA (RR = 8.706, P = 3.97 × 10-4), and all-cause mortality (HR = 110.000, P = 0.001). Collectively, these results show that SNP variations in the COX-2 promoter, and their inherited combinations, are associated with the longitudinal risk of malaria, SMA, and all-cause mortality among children living in a high transmission area for P. falciparum.

Host genetic background affects the course of infection and treatment response in patients with chronic hepatitis B

BACKGROUND

Hepatitis B virus (HBV) utilizes proteins encoded by the host to infect hepatocytes and replicate. Recently, several novel host factors have been identified and described as important to the HBV lifecycle. The influence of host genetic background on chronic hepatitis B (CHB) pathogenesis is still poorly understood.

OBJECTIVES

Here, we aimed to investigate the association of NTCP, FXRα, HNF1α, HNF4α, and TDP2 genetic polymorphisms with the natural course of CHB and antiviral treatment response.

STUDY DESIGN

We genotyped 18 single-nucleotide polymorphisms using MALDI-TOF mass spectrometry in 136 patients with CHB and 100 healthy individuals. We investigated associations of the selected polymorphisms with biochemical, serological and hepatic markers of disease progression and treatment response.

RESULTS

No significant differences in genotypic or allelic distribution between CHB and control groups were observed. Within TDP2, rs3087943 variations were associated with treatment response, and rs1047782 modified the risk of advanced liver inflammation. Rs7154439 within NTCP was associated with HBeAg seroconversion after 48 weeks of nucleos(t)ide analogue treatment. HNF1α genotypes were associated with treatment response, liver damage and baseline HBeAg presence. HNF4α rs1800961 predicted PEG-IFNα treatment-induced HBsAg clearance in long-term follow up.

CONCLUSIONS

This study indicates host genetic background relevance in the course of CHB and confirms the role of recently described genes for HBV infection. The obtained results might serve as a starting point for validation studies on the clinical application of selected genetic variants to predict individual risks of CHB-induced liver failure and treatment response.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Metagenomic Approaches for Public Health Surveillance of Foodborne Infections: Opportunities and Challenges

Foodborne disease surveillance in the United States is at a critical point. Clinical and diagnostic laboratories are using culture-independent diagnostic tests (CIDTs) to identify the pathogen causing foodborne illness from patient specimens. CIDTs are molecular tests that allow doctors to rapidly identify the bacteria causing illness within hours. CIDTs, unlike previous gold standard methods such as bacterial culture, do not produce an isolate that can be subtyped as part of the national molecular subtyping network for foodborne disease surveillance, PulseNet. Without subtype information, cases can no longer be linked using molecular data to identify potentially related cases that are part of an outbreak. In this review, we discuss the public health needs for a molecular subtyping approach directly from patient specimen and highlight different approaches, including amplicon and shotgun metagenomic sequencing.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

Folate pathway genetic polymorphisms modulate methotrexate-induced toxicity in childhood acute lymphoblastic leukemia

BACKGROUND

Acute lymphoblastic leukemia (ALL) is one of the major malignancies affecting children in Jordan. Methotrexate (MTX) is the cornerstone of chemotherapy for ALL, and works by targeting enzymes involved in the folate pathway. We hypothesize that genetic polymorphisms of the folate pathway are associated with MTX toxicity in children with ALL.

METHODS

A total of 64 children with ALL were included in this study; 31 (48.4%) boys and 33 (51.6%) girls aged 2-16 years. The folate pathway genes were genotyped using polymerase chain reaction followed by sequencing and studying the association between genetic polymorphisms and MTX toxicity.

RESULTS

The immunophenotype was B-lineage in 55 patients (85.9%) and T-lineage in nine patients (14.1%). All genetic polymorphisms, except for dihydropyrimidine dehydrogenase polymorphisms, were associated with hematological toxicities and did not appear to precipitate any non-hematological adverse events. Patients with ALL carrying dominant alleles of methylene tetrahydrofolate (MTHFR) C677T and dihydrofolate reductase 19 bp deletion were at a higher risk of developing severe leucopenia [OR (95% CI) = 4.5 (1.2-17), p = 0.03; 5.4 (1.6-17.8); p = 0.006] while minor allele carriers of MTHFR A1298C were more likely to develop neutropenia [OR (95% CI) = 6.1 (1.3-29.5); 0.04]. Furthermore, dominant allele carriers of thymidylate synthase 1494 del6 were at a higher risk of developing neutropenia [OR (95% CI) = 6 (1.2-31.1); p = 0.04].

CONCLUSION

Genetic polymorphisms of the folate pathway may modulate MTX-induced toxicity in childhood ALL.

The rs4686434 variant in the fetuin B (FETUB) locus is associated with intrahepatic triglyceride content in obese Chinese adults

BACKGROUND

This study explored associations of genetic variants in the fetuin B (FETUB) locus with intrahepatic triglyceride (IHTG) content.

METHODS

Four tagging single-nucleotide polymorphisms (SNPs) of the FETUB locus and patatin-like phospholipase domain containing 3 (PNPLA3) rs738409 and transmembrane 6 super family member 2 (TM6SF2) rs58542926 were genotyped in 418 obese Chinese adults in whom serum FETUB and IHTG were measured.

RESULTS

Subjects carrying the minor G allele for FETUB rs4686434 (AG/GG) had lower serum FETUB levels (mean [±SD] 3.89 ± 1.36 vs 4.22 ± 1.46 μg/mL; P = 0.021) and IHTG content (12.7 ± 9.4% vs 14.6 ± 9.8%; P = 0.045) than their controls (AA), whereas IHTG content was higher in those carrying the minor G allele for PNPLA3 rs738409 (CG/GG) than in their controls (CC; 14.5 ± 10.1% vs 12.0 ± 8.6%; P = 0.012). After adjusting for potential confounders, IHTG content was lower in carriers of the minor G allele for FETUB rs4686434 (AG/GG vs AA, β -2.27 ± 0.91, P = 0.012), but was higher in carriers of the minor G allele for PNPLA3 rs738409 (CG/GG vs CC, β 2.65 ± 0.97, P = 0.006). There was a significant joint effect between FETUB rs4686434 and PNPLA3 rs738409 on IHTG content, with increasing genetic risk score (counting the risk allele of A in rs4686434 and G in rs738409) being associated with higher IHTG content (β 1.85 ± 0.48, P <0.001).

CONCLUSIONS

Carrying the minor G allele for FETUB rs4686434 was significantly associated with decreased IHTG content and may affect hepatic triglyceride accumulation in individuals at high risk of non-alcoholic fatty liver disease.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Catechol-O-Methyltransferase moderates effect of stress mindset on affect and cognition

There is evidence that altering stress mindset-the belief that stress is enhancing vs. debilitating-can change cognitive, affective and physiological responses to stress. However individual differences in responsiveness to stress mindset manipulations have not been explored. Given the previously established role of catecholamines in both placebo effects and stress, we hypothesized that genetic variation in catechol-O-methyltransferase (COMT), an enzyme that metabolizes catecholamines, would moderate responses to an intervention intended to alter participants' mindsets about stress. Participants (N = 107) were exposed to a stress mindset manipulation (videos highlighting either the enhancing or debilitating effects of stress) prior to engaging in a Trier Social Stress task and subsequent cognitive tasks. The associations of the COMT rs4680 polymorphism with the effect of stress mindset video manipulations on cognitive and affective responses were examined. Genetic variation at rs4680 modified the effects of stress mindset on affective and cognitive responses to stress. Individuals homozygous for rs4680 low-activity allele (met/met) were responsive to the stress-is-enhancing mindset manipulation as indicated by greater increases in positive affect, improved cognitive functioning, and happiness bias in response to stress. Conversely, individuals homozygous for the high-activity allele (val/val) were not as responsive to the stress mindset manipulation. These results suggest that responses to stress mindset intervention may vary with COMT genotype. These findings contribute to the understanding of gene by environment interactions for mindset interventions and stress reactivity and therefore warrant further investigations.

The Regulation of Monoamine Oxidase A Gene Expression by Distinct Variable Number Tandem Repeats

The monoamine oxidase A (MAOA) uVNTR (upstream variable number tandem repeat) is one of the most often cited examples of a gene by environment interaction (GxE) in relation to behavioral traits. However, MAOA possesses a second VNTR, 500 bp upstream of the uVNTR, which is termed d- or distal VNTR. Furthermore, genomic analysis indicates that there are a minimum of two transcriptional start sites (TSSs) for MAOA, one of which encompasses the uVNTR within the 5' untranslated region of one of the isoforms. Through expression analysis in semi-haploid HAP1 cell lines genetically engineered in order to knockout (KO) either the uVNTR, dVNTR, or both VNTRs, we assessed the effect of the two MAOA VNTRs, either alone or in combination, on gene expression directed from the different TSSs. Complementing our functional analysis, we determined the haplotype variation of these VNTRs in the general population. The expression of the two MAOA isoforms was differentially modulated by the two VNTRs located in the promoter region. The most extensively studied uVNTR, previously considered a positive regulator of the MAOA gene, did not modulate the expression of what it is considered the canonical isoform, while we found that the dVNTR positively regulated this isoform in our model. In contrast, both the uVNTR and the dVNTR were found to act as negative regulators of the second less abundant MAOA isoform. The haplotype analysis for these two VNTRs demonstrated a bias against the presence of one of the potential variants. The uVNTR and dVNTR differentially affect expression of distinct MAOA isoforms, and thus, their combined profiling offers new insights into gene-regulation, GxE interaction, and ultimately MAOA-driven behavior.

Common genetic variants in the FETUB locus, genetically predicted fetuin-B levels, and risk of insulin resistance in obese Chinese adults

Elevated serum fetuin-B is suggested to be associated with insulin resistance, but it is unknown if this association is causal. The aim of this study was to explore the potential causal relationship between fetuin-B and insulin resistance.We used Mendelian randomization analysis by incorporating information of genetic variants in FETUB and serum fetuin-B concentrations with insulin resistance in 1148 obese Chinese adults.Common genetic variants (FETUB rs4686434, rs6785067, and rs3733159) were significantly associated with serum fetuin-B concentrations but not with insulin resistance. Higher serum fetuin-B levels were significantly associated with increased homeostasis model assessment of insulin resistance (HOMA-IR) (0.17 [95%CI: 0.01 to 0.32, P = .037] 10 mol IU L higher per SD). However, Mendelian randomization analysis using 3 single-nucleotide polymorphisms as instrumental variables did not support a significant association between genetically predicted fetuin-B levels and HOMA-IR (-0.09 [95%CI: -0.62 to 0.44, P = .738] 10 mol IU L lower per SD). The regression coefficients for measured and genetically predicted fetuin-B concentrations on HOMA-IR were significantly different (P <.001).This study suggests the association between fetuin-B and insulin resistance may not be causal. Future studies on the nongenetic determinants of serum fetuin-B concentration to assess if such unmeasured factors may confound the association between fetuin-B and insulin resistance as well as more pathway analysis for this association are warranted.

Single nucleotide polymorphisms and microsatellites in the canine glutathione S-transferase pi 1 (GSTP1) gene promoter

BACKGROUND

Genetic polymorphisms within the glutathione S-transferase P1 (GSTP1) gene affect the elimination of toxic xenobiotics by the GSTP1 enzyme. In dogs, exposure to environmental chemicals that may be GSTP1 substrates is associated with cancer. The objectives of this study were to investigate the genetic variability in the GSTP1 promoter in a diverse population of 278 purebred dogs, compare the incidence of any variants found between breeds, and predict their effects on gene expression. To provide information on ancestral alleles, a number of wolves, coyotes, and foxes were also sequenced.

RESULTS

Fifteen single nucleotide polymorphisms (SNPs) and two microsatellites were discovered. Three of these loci were only polymorphic in dogs while three other SNPs were unique to wolves and coyotes. The major allele at c.-46 is T in dogs but is C in the wild canids. The c.-185 delT variant was unique to dogs. The microsatellite located in the 5' untranslated region (5'UTR) was a highly polymorphic GCC tandem repeat, consisting of simple and compound alleles that varied in size from 10 to 22-repeat units. The most common alleles consisted of 11, 16, and 17-repeats. The 11-repeat allele was found in 10% of dogs but not in the other canids. Unequal recombination and replication slippage between similar and distinct alleles may be the mechanism for the multiple microsatellites observed. Twenty-eight haplotypes were constructed in the dog, and an additional 8 were observed in wolves and coyotes. While the most common haplotype acrossbreeds was the wild-type *1A(17), other prevalent haplotypes included *3A(11) in Greyhounds, *6A(16) in Labrador Retrievers, *9A(16) in Golden Retrievers, and *8A(19) in Standard Poodles. Boxers and Siberian Huskies exhibited minimal haplotypic diversity. Compared to the simple 16*1 allele, the compound 16*2 allele (found in 12% of dogs) may interfere with transcription factor binding and/or the stability of the GSTP1 transcript.

CONCLUSIONS

Dogs and other canids exhibit extensive variation in the GSTP1 promoter. Genetic polymorphisms within distinct haplotypes prevalent in certain breeds can affect GSTP1 expression and carcinogen detoxification, and thus may be useful as genetic markers for cancer in dogs.

Population genetic structure of Guiana dolphin (Sotalia guianensis) from the southwestern Atlantic coast of Brazil

Sotalia guianensis is a small dolphin that is vulnerable to anthropogenic impacts. Along the Brazilian Atlantic coast, this species is threatened with extinction. A prioritized action plan for conservation strategies relies on increased knowledge of the population. The scarcity of studies about genetic diversity and assessments of population structure for this animal have precluded effective action in the region. Here, we assessed, for the first time, the genetic differentiation at 14 microsatellite loci in 90 S. guianensis specimens stranded on the southeastern Atlantic coast of the State of Espírito Santo, Brazil. We estimated population parameters and structure, measured the significance of global gametic disequilibrium and the intensity of non-random multiallelic interallelic associations and constructed a provisional synteny map using Bos taurus, the closest terrestrial mammal with a reference genome available. All microsatellite loci were polymorphic, with at least three and a maximum of ten alleles each. Allele frequencies ranged from 0.01 to 0.97. Observed heterozygosity ranged from 0.061 to 0.701. The mean inbreeding coefficient was 0.103. Three loci were in Hardy-Weinberg disequilibrium even when missing genotypes were inferred. Although 77 of the 91 possible two-locus associations were in global gametic equilibrium, we unveiled 13 statistically significant, sign-based, non-random multiallelic interallelic associations in 10 two-locus combinations with either coupling (D' values ranging from 0.782 to 0.353) or repulsion (D' values -0.517 to -1.000) forces. Most of the interallelic associations did not involve the major alleles. Thus, for either physically or non-physically linked loci, measuring the intensity of non-random interallelic associations is important for defining the evolutionary forces at equilibrium. We uncovered a small degree of genetic differentiation (F_ST = 0.010; P-value = 0.463) with a hierarchical clustering into one segment containing members from the southern and northern coastal regions. The data thus support the scenario of little genetic structure in the population of S. guianensis in this geographic area.

Interferon lambda polymorphisms associate with body iron indices and hepatic expression of interferon-responsive long non-coding RNA in chronic hepatitis C

Single nucleotide polymorphisms (SNPs) within DNA region containing interferon lambda 3 (IFNL3) and IFNL4 genes are prognostic factors of treatment response in chronic hepatitis C (CHC). Iron overload, frequently diagnosed in CHC, is associated with unfavorable disease course and a risk of carcinogenesis. Its etiology and relationship with the immune response in CHC are not fully explained. Our aim was to determine whether IFNL polymorphisms in CHC patients associate with body iron indices, and whether they are linked with hepatic expression of genes involved in iron homeostasis and IFN signaling. For 192 CHC patients, four SNPs within IFNL3-IFNL4 region (rs12979860, rs368234815, rs8099917, rs12980275) were genotyped. In 185 liver biopsies, histopathological analyses were performed. Expression of five mRNAs and three long non-coding RNAs (lncRNAs) was determined with qRT-PCR in 105 liver samples. Rs12979860 TT or rs8099917 GG genotypes as well as markers of serum and hepatocyte iron overload associated with higher activity of gamma-glutamyl transpeptidase and liver steatosis. The presence of two minor alleles in any of the tested SNPs predisposed to abnormally high serum iron concentration and correlated with higher hepatic expression of lncRNA NRIR. On the other hand, homozygosity in any major allele associated with higher viral load. Patients bearing rs12979860 CC genotype had lower hepatic expression of hepcidin (HAMP; P = 0.03). HAMP mRNA level positively correlated with serum iron indices and degree of hepatocyte iron deposits. IFNL polymorphisms influence regulatory pathways of cellular response to IFN and affect body iron balance in chronic hepatitis C virus infection.

Influence of the CYP3A4/5 genetic score and ABCB1 polymorphisms on tacrolimus exposure and renal function in Brazilian kidney transplant patients

BACKGROUND

Polymorphisms in genes encoding transport proteins and metabolizing enzymes involved in tacrolimus (TAC) disposition may be important sources of individual variability during treatment.

OBJECTIVE

The aim of this study was to investigate the effect of combined CYP3A4 and CYP3A5 variants, using a CYP3A4/5 genetic score, and ABCB1 polymorphisms on therapeutic TAC monitoring and their relationship with clinical outcomes.

MATERIAL AND METHODS

Brazilian kidney transplant recipients (n=151), who received TAC over 3 months after transplantation, were genotyped for CYP3A4 rs2242480 (g.20230G>A), CYP3A5 rs15524 (g.31611C>T) and rs776746 (g.6986A>G), ABCB1 rs1128503 (c.1236C>T), rs1045642 (c.3435C>T), and rs2032582 (c.2677G>T/A) polymorphisms.

RESULTS

Frequencies of CYP3A4 g.20230A, CYP3A5 g.31611C, and g.6986A were 0.37, 0.26, and 0.28, respectively. These alleles were associated with TAC rapid metabolization and were used for CYP3A4/5 genetic score construction. A higher CYP3A4/5 genetic score was associated with higher TAC dose and lower concentrations for dose administered (Co/D, P<0.05). Ninety days after transplantation, the presence of two or more rapid metabolization alleles contributed toward 27.7% of Co/D variability and was associated with a lower estimated glomerular filtration rate values (P<0.05). For ABCB1, the frequencies of c.1236T, c.3435T, and c.2677T/A alleles were 0.42, 0.42, and 0.33/0.04. At 30 days after transplantation, patients carrying ABCB1 c.1236TT+c.3435TT+(c.2677TT+TA) genotypes had higher TAC Co/D than those with common or heterozygous genotypes (P<0.05).

CONCLUSION

The results show the impact of the CYP3A4/5 genetic score on TAC exposure and renal function in Brazilian patients. Furthermore, ABCB1 polymorphisms, in a combined analysis, influenced TAC Co/D at 30 days after transplantation.

Single Marker and Haplotype-Based Association Analysis of Semolina and Pasta Colour in Elite Durum Wheat Breeding Lines Using a High-Density Consensus Map

Association mapping is usually performed by testing the correlation between a single marker and phenotypes. However, because patterns of variation within genomes are inherited as blocks, clustering markers into haplotypes for genome-wide scans could be a worthwhile approach to improve statistical power to detect associations. The availability of high-density molecular data allows the possibility to assess the potential of both approaches to identify marker-trait associations in durum wheat. In the present study, we used single marker- and haplotype-based approaches to identify loci associated with semolina and pasta colour in durum wheat, the main objective being to evaluate the potential benefits of haplotype-based analysis for identifying quantitative trait loci. One hundred sixty-nine durum lines were genotyped using the Illumina 90K Infinium iSelect assay, and 12,234 polymorphic single nucleotide polymorphism (SNP) markers were generated and used to assess the population structure and the linkage disequilibrium (LD) patterns. A total of 8,581 SNPs previously localized to a high-density consensus map were clustered into 406 haplotype blocks based on the average LD distance of 5.3 cM. Combining multiple SNPs into haplotype blocks increased the average polymorphism information content (PIC) from 0.27 per SNP to 0.50 per haplotype. The haplotype-based analysis identified 12 loci associated with grain pigment colour traits, including the five loci identified by the single marker-based analysis. Furthermore, the haplotype-based analysis resulted in an increase of the phenotypic variance explained (50.4% on average) and the allelic effect (33.7% on average) when compared to single marker analysis. The presence of multiple allelic combinations within each haplotype locus offers potential for screening the most favorable haplotype series and may facilitate marker-assisted selection of grain pigment colour in durum wheat. These results suggest a benefit of haplotype-based analysis over single marker analysis to detect loci associated with colour traits in durum wheat.

Thymidylate Synthase Polymorphisms and Risk of Lung Cancer among the Jordanian Population: a Case Control Study

BACKGROUND

Thymidylate synthase (TS) catalyzes the methylation of deoxyuridylate to deoxythymidylate and is involved in DNA methylation, synthesis and repair. Two common polymorphisms have been reported, tandem repeats in the promoter-enhancer region (TSER), and 6bp ins/del in the 5'UTR, that are implicated in a number of human diseases, including cancer. The association between the two polymorphisms in risk for lung cancer (LC) was here investigated in the Jordanian population.

MATERIALS AND METHODS

An age, gender, and smoking-matched case-control study involving 84 lung cancer cases and 71 controls was conducted. The polymerase chain reaction/restriction fragment length polymorphism (PCR-RFLP) technique was used to detect the polymorphism of interest.

RESULTS

Individuals bearing the ins/ins genotype were 2.5 times more likely to have lung cancer [(95%CI: 0.98-6.37), p=0.051]. Individuals who were less than or equal to 57 years and carrying ins/ins genotype were 4.6 times more susceptible to lung cancer [OR<57 vs >57years: 4.6 (95%CI: 0.93-22.5), p=0.059)]. Genotypes and alleles of TSER were distributed similarly between cases and controls. Weak linkage disequilibrium existed between the two loci of interest (Lewontin's coefficient [D']) (LC: D' =0.03, r2: 0. 001, p= 0.8;

CONTROLS

D' =0.29, r2: 0.08, p=0.02). Carriers of the "3 tandem repeats_insertion" haplotype (3R_ins) were 2 times more likely to have lung cancer [2 (95%CI: 1.13-3.48), p=0.061].

CONCLUSIONS

Genetic polymorphism of TS at 3` UTR and its haplotype analysis may modulate the risk of lung cancer in Jordanians. The 6bp ins/del polymorphism of TS at 3 `UTR is more informative than TSER polymorphism in predicting increased risk.

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life-history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness-related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O3+4 arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.

The Relation between eNOS -786 C/T, 4 a/b, MMP-13 rs640198 G/T, Eotaxin 426 C/T, -384 A/G, and 67 G/A Polymorphisms and Long-Term Outcome in Patients with Coronary Artery Disease

Aim. The purpose of this study is to determine the association between eotaxin 426 C/T, −384 A/G, 67 G/A, eNOS −786 T/C, 4 a/b, and MMP-13 rs640198 G/T and prognosis of patients with known CAD. Methods. From total of 1161 patients referred to coronary angiography, 532 patients with angiographically confirmed CAD were selected. Their long-term outcome was followed up using hospital database. Subsequent events were assessed in this study: death or combined endpoint-myocardial infarction, unstable angina pectoris, revascularization, heart failure hospitalization, and cardioverter-defibrillator implantation. Results. The multivariate Cox regression model identified age, smoking, and 3-vessel disease as significant predictors of all-cause death. Further analysis showed that eotaxin 67 G/A (GA + AA versus GG) and eotaxin −384 A/G (GG versus GA + AA) were significant independent prognostic factors when added into the model: HR (95% CI) 2.81 (1.35–5.85), p = 0.006; HR (95% CI) 2.63 (1.19–5.83), p = 0.017; eotaxin −384 A/G was significantly associated with the event-free survival, but it did not provide the prognostic information above the effect of two- or three-vessel disease. Conclusion. The A allele in eotaxin 67 G/A polymorphism is associated with worse survival in CAD patients.

Lack of association of rs3798220 with small apolipoprotein(a) isoforms and high lipoprotein(a) levels in East and Southeast Asians

OBJECTIVE

The variant allele of rs3798220 in the apolipoprotein(a) gene (LPA) is used to assess the risk for coronary artery disease (CAD) in Europeans, where it is associated with short alleles of the Kringle IV-2 (KIV-2) copy number variation (CNV) and high lipoprotein(a) (Lp(a)) concentrations. No association of rs3798220 with CAD was detected in a GWAS of East Asians. Our study investigated the association of rs3798220 with Lp(a) concentrations and KIV-2 CNV size in non-European populations to explain the missing association of the variant with CAD in Asians.

METHODS

We screened three populations from Africa and seven from Asia by TaqMan Assay for rs3798220 and determined KIV-2 CNV sizes of LPA alleles by pulsed-field gel electrophoresis (PFGE). Additionally, CAD cases from India were analysed. To investigate the phylogenetic origin of rs3798220, 40 LPA alleles from Chinese individuals were separated by PFGE and haplotyped for further SNPs.

RESULTS

The variant was not found in Africans. Allele frequencies in East and Southeast Asians ranged from 2.9% to 11.6%, and were very low (0.15%) in CAD cases and controls from India. The variant was neither associated with short KIV-2 CNV alleles nor elevated Lp(a) concentrations in Asians.

CONCLUSION

Our study shows that rs3798220 is no marker for short KIV-2 CNV alleles and high Lp(a) in East and Southeast Asians, although the haplotype background is shared with Europeans. It appears unlikely that this SNP confers atherogenic potential on its own. Furthermore, this SNP does not explain Lp(a) attributed risk for CAD in Asian Indians.

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the 'species trait' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory.

Patients with chronic three-vessel disease in a 15-year follow-up study: genetic and non-genetic predictors of survival

Genetic and non-genetic predictors of 15-year survival in patients with chronic three-vessel disease (3VD) were investigated. Coronary angiography was performed on 810 subjects with symptoms of stable ischemic heart disease in 1998. The patients with 3VD were genotyped for 23 candidate polymorphisms covering the PPAR-RXR pathway, matrix metalloproteinase-2, renin-angiotensin-aldosterone system, endothelin-1, cytokine genes, MTHFR and APO E variants. Fifteen-year survival data were obtained from the national insurance registry. All data were available in the case of 150 patients with 3VD. Statistical analysis used stepwise Cox regression with dominant, recessive, or additive mode of genetic expression. Involved variables included age, sex, BMI, blood pressure, diabetes, ejection fraction, left main stenosis, previously diagnosed coronary stenosis, myocardial infarction in personal history, and coronary bypass along with polymorphisms pre-selected by log-rank tests. Out of the 23 polymorphisms, four were included in the model construction. SNP in the IL-6 gene rs1800795 (-174 G/C) has been found to be a significant predictor of survival. This SNP was in a linkage disequilibrium with rs1800797 (-597 G/A) in the same gene (D'=1.0), which was also found to constitute a significant predictor of survival when rs1800795 was not included in the model construction. Age, increased BMI, diabetes, low EF, and left main stenosis were also significant predictors in all models. Age, increased BMI, diabetes, low ejection fraction, left main stenosis, and genetic variation in the IL-6 promoter were established as significant independent risk factors for the survival of patients with three-vessel disease.

Maternal prenatal anxiety and child brain-derived neurotrophic factor (BDNF) genotype: Effects on internalizing symptoms from 4 to 15 years of age

Multiple behavioral and health outcomes, including internalizing symptoms, may be predicted from prenatal maternal anxiety, depression, or stress. However, not all children are affected, and those that are can be affected in different ways. Here we test the hypothesis that the effects of prenatal anxiety are moderated by genetic variation in the child's brain-derived neurotrophic factor (BDNF) gene, using the Avon Longitudinal Study of Parents and Children population cohort. Internalizing symptoms were assessed from 4 to 13 years of age using the Strengths and Difficulties Questionnaire (n = 8,584); a clinical interview with the adolescents was conducted at age 15 years (n = 4,704). Obstetric and psychosocial risk and postnatal maternal symptoms were included as covariates. Results show that prenatal maternal anxiety predicted internalizing symptoms, including with the diagnostic assessment at 15 years. There was a main effect of two BDNF polymorphisms (rs6265 [val66met] and rs11030104) on internalizing symptoms up to age 13. There was also genetic moderation of the prenatal anxiety effect by different BDNF polymorphisms (rs11030121 and rs7124442), although significant effects were limited to preadolescence. The findings suggest a role for BDNF gene–environment interactions in individual vulnerability to the effects of prenatal anxiety on child internalizing symptoms.

Variation and linkage disequilibrium between a structurally polymorphic Alu located near the OR12D2 gene of the extended major histocompatibility complex class I region and HLA-A alleles

We investigated the genetic structure and population frequency of an Alu repeat dimorphism (absence or presence) located near the OR12D2 gene within the olfactory receptor gene region telomeric of the alpha HLA class I region (HLA-J, -A, -G, -F). The structurally polymorphic Alu insertion (POALIN) locus rs33972478 that we designated as AluOR and its allele and haplotype frequencies and association with HLA-A and six other structurally polymorphic retroelements (3 Alu, 2 SVA and an HERVK9) were determined in 100 Japanese, 174 Caucasians and 100 African American DNA samples. The AluOR insertion varied in population frequency between 14.4% and 31.5% with significant differences between the Japanese and Caucasians, but not between the Caucasian and African Americans. Although AluOR is located 600 kb from the HLA-A gene, there was a significant linkage disequilibrium between the two loci and a high percentage association of the AluOR insertion with HLA-A29 (79%) in Caucasians and HLA-A31 (69.4%) in Japanese. Inferred haplotypes among three-locus to eight-locus haplotype structures showed maximum differences between the populations with the eight-locus haplotypes. The most frequent multilocus haplotype shared between the populations was the HLA-A2 allele in combination with the AluHG insertion. The AluOR whether investigated alone or together with the HLA class I alleles and other dimorphic retroelements is an informative ancestral marker for the identification of lineages and variations within the same and/or different populations and for examining the linkage and crossing-over between the HLA and OR genomic regions in the extended MHC.

Microsatellite polymorphism and its association with body weight and selected morphometrics of farm red fox (Vulpes vulpes L.)

Polymorphism of 30 canine-derived microsatellites was studied in a group of 200 red foxes kept on 2 Polish farms. 22 out of 30 microsatellites were selected to study association between marker genotypes and body weight (BW), body length (BL), body circumference (BC), tail length (TL), ear height (EH), length of the right front limb (FRLL), length of the right rear limb (RRLL), length of the right front foot (FRFL) and length of the right rear foot (RRFL). A total of 112 alleles and 243 genotypes were found at 22 autosomal microsatellite loci. Three monomorphic loci deemed as uninformative were excluded from the study. The association between marker genotypes and the studied traits was analysed using general linear model (GLM) procedure and least squares means (LSM). Linkage disequilibrium (LD) was estimated to assess non-random association between microsatellite loci. Out of 19 microsatellites studied four markers showed no association with the studied traits, three markers had a significant effect on one trait, and another three markers had significant effect on two traits. Among ten microsatellites with significant effect on four economically important traits (BW, BL, BC, TL) four were associated with two characters: marker FH2613 with BW and BC, marker FH2097withBL and BC, marker ZUBECA6 with BW and BC, whereas marker REN75M10 was associated with BL and TL. The strongest LD (r² ranged from 0.15 to 0.33) was estimated between nine loci with significant effect on economically important traits (BW, BL, BC, TL).

Emerging paradigms in genomics-based crop improvement

Next generation sequencing platforms and high-throughput genotyping assays have remarkably expedited the pace of development of genomic tools and resources for several crops. Complementing the technological developments, conceptual shifts have also been witnessed in designing experimental populations. Availability of second generation mapping populations encompassing multiple alleles, multiple traits, and extensive recombination events is radically changing the phenomenon of classical QTL mapping. Additionally, the rising molecular breeding approaches like marker assisted recurrent selection (MARS) that are able to harness several QTLs are of particular importance in obtaining a "designed" genotype carrying the most desirable combinations of favourable alleles. Furthermore, rapid generation of genome-wide marker data coupled with easy access to precise and accurate phenotypic screens enable large-scale exploitation of LD not only to discover novel QTLs via whole genome association scans but also to practise genomic estimated breeding value (GEBV)-based selection of genotypes. Given refinements being experienced in analytical methods and software tools, the multiparent populations will be the resource of choice to undertake genome wide association studies (GWAS), multiparent MARS, and genomic selection (GS). With this, it is envisioned that these high-throughput and high-power molecular breeding methods would greatly assist in exploiting the enormous potential underlying breeding by design approach to facilitate accelerated crop improvement.

Patterns of cyto-nuclear linkage disequilibrium in Silene latifolia: genomic heterogeneity and temporal stability

Non-random association of alleles in the nucleus and cytoplasmic organelles, or cyto-nuclear linkage disequilibrium (LD), is both an important component of a number of evolutionary processes and a statistical indicator of others. The evolutionary significance of cyto-nuclear LD will depend on both its magnitude and how stable those associations are through time. Here, we use a longitudinal population genetic data set to explore the magnitude and temporal dynamics of cyto-nuclear disequilibria through time. We genotyped 135 and 170 individuals from 16 and 17 patches of the plant species Silene latifolia in Southwestern VA, sampled in 1993 and 2008, respectively. Individuals were genotyped at 14 highly polymorphic microsatellite markers and a single-nucleotide polymorphism (SNP) in the mitochondrial gene, atp1. Normalized LD (D') between nuclear and cytoplasmic loci varied considerably depending on which nuclear locus was considered (ranging from 0.005-0.632). Four of the 14 cyto-nuclear associations showed a statistically significant shift over approximately seven generations. However, the overall magnitude of this disequilibrium was largely stable over time. The observed origin and stability of cyto-nuclear LD is most likely caused by the slow admixture between anciently diverged lineages within the species' newly invaded range, and the local spatial structure and metapopulation dynamics that are known to structure genetic variation in this system.

Effect of genetic polymorphisms in CA6 gene on the expression and catalytic activity of human salivary carbonic anhydrase VI

Carbonic anhydrase isoenzyme VI (CA VI) plays an important role in the homeostasis of oral tissues participating in the processes of taste, protection of dental tissues against the loss of minerals, caries, and possibly in the formation of dental calculus in periodontal disease. This study aimed to verify the correlation between changes in the expression and activity of human salivary carbonic anhydrase VI and genetic polymorphisms in its gene (CA6). The study population consisted of 182 healthy volunteers (female and male, aged 18-22). Samples of total saliva were assayed for CA VI concentrations using a specific time-resolved immunofluorometric assay. CA VI catalytic activity was detected by a modified protocol of Kotwica et al. [J Physiol Pharmacol 2006;57(suppl 8):107-123], adapted to CA VI in saliva. Samples of genomic DNA were genotyped for polymorphisms rs2274327 (C/T), rs2274328 (A/C) and rs2274333 (A/G) by TaqMan® SNP Genotyping Assays. The concentration and catalytic activity of the salivary CA VI obtained for the different genotypes were analyzed using the Kruskal-Wallis nonparametric test and the Dunn test. The results showed that individuals with TT genotype (rs2274327) had significantly lower CA VI concentrations than the individuals with genotypes CT or CC (p < 0.05). There was also an association between polymorphism rs2274333 and salivary CA VI concentrations. There were no associations between the three polymorphisms analyzed and variations in CA VI activity. Our results suggest that polymorphisms in the CA6 gene are associated with the concentrations of secreted CA VI.

Common genetic variation in the human FNDC5 locus, encoding the novel muscle-derived 'browning' factor irisin, determines insulin sensitivity

Aims/hypothesis

Recently, the novel myokine irisin was described to drive adipose tissue ‘browning’, to increase energy expenditure, and to improve obesity and insulin resistance in high fat-fed mice. Here, we assessed whether common single nucleotide polymorphisms (SNPs) in the FNDC5 locus, encoding the irisin precursor, contribute to human prediabetic phenotypes (overweight, glucose intolerance, insulin resistance, impaired insulin release).

Methods

A population of 1,976 individuals was characterized by oral glucose tolerance tests and genotyped for FNDC5 tagging SNPs. Subgroups underwent hyperinsulinaemic-euglycaemic clamps, magnetic resonance imaging/spectroscopy, and intravenous glucose tolerance tests. From 37 young and 14 elderly participants recruited in two different centres, muscle biopsies were obtained for the preparation of human myotube cultures.

Results

After appropriate adjustment and Bonferroni correction for the number of tested variants, SNPs rs16835198 and rs726344 were associated with in vivo measures of insulin sensitivity. Via interrogation of publicly available data from the Meta-Analyses of Glucose and Insulin-related traits Consortium, rs726344’s effect on insulin sensitivity was replicated. Moreover, novel data from human myotubes revealed a negative association between FNDC5 expression and appropriately adjusted in vivo measures of insulin sensitivity in young donors. This finding was replicated in myotubes from elderly men.

Conclusions/interpretation

This study provides evidence that the FNDC5 gene, encoding the novel myokine irisin, determines insulin sensitivity in humans. Our gene expression data point to an unexpected insulin-desensitizing effect of irisin.

Do pathogens reduce genetic diversity of their hosts? Variable effects of sylvatic plague in black-tailed prairie dogs

Introduced diseases can cause dramatic declines in-and even the loss of-natural populations. Extirpations may be followed by low recolonization rates, leading to inbreeding and a loss of genetic variation, with consequences on population viability. Conversely, extirpations may create vacant habitat patches that individuals from multiple source populations can colonize, potentially leading to an influx of variation. We tested these alternative hypotheses by sampling 15 colonies in a prairie dog metapopulation during 7 years that encompassed an outbreak of sylvatic plague, providing the opportunity to monitor genetic diversity before, during and after the outbreak. Analysis of nine microsatellite loci revealed that within the metapopulation, there was no change in diversity. However, within extirpated colonies, patterns varied: In half of the colonies, allelic richness after recovery was less than the preplague conditions, and in the other half, richness was greater than the preplague conditions. Finally, analysis of variation within individuals revealed that prairie dogs present in recolonized colonies had higher heterozygosity than those present before plague. We confirmed plague survivorship in six founders; these individuals had significantly higher heterozygosity than expected by chance. Collectively, our results suggest that high immigration rates can maintain genetic variation at a regional scale despite simultaneous extirpations in spatially proximate populations. Thus, virulent diseases may increase genetic diversity of host populations by creating vacant habitats that allow an influx of genetic diversity. Furthermore, even highly virulent diseases may not eliminate individuals randomly; rather, they may selectively remove the most inbred individuals.

The genetic content of chromosomal inversions across a wide latitudinal gradient

There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.

IL-1RN +2018T>C polymorphism is correlated with colorectal cancer

Epidemiological and experimental evidence indicates chronic inflammation as a risk factor for colorectal cancer. We investigated whether IL-1B -511C>T (rs16944), IL-1B +3954C>T (rs1143634) and IL1-RN +2018T>C (rs419598) cytokine polymorphisms are correlated with colorectal cancer. Blood samples were obtained from 377 Romanian subjects: 144 patients with sporadic colorectal cancer and 233 healthy controls. Polymorphisms were analyzed by allelic discrimination TaqMan PCR assays with specific probes. The results of our study showed that IL-1RN +2018T>C polymorphism is associated with colorectal cancer. We found that there was a significant difference in the frequency of CC genotype between patients with colorectal cancer and the control group (OR 2.42, 95 % CI: 1.06-5.53, p = 0,034) when TT genotype was used as reference. Furthermore, in a stratified analysis, a positive association was found only for IL-1RN +2018CC genotype, that was limited to early I and II stages (OR 2.72, 95 % CI: 1.05-7.03, p = 0,033). We did not find any association between any of the IL-1B polymorphisms and colorectal cancer. In conclusion this study found that IL-1RN +2018T>C polymorphism is associated with colorectal cancer, mainly for localized disease.

Catechol-O-methyltransferase val158met polymorphism predicts placebo effect in irritable bowel syndrome

Identifying patients who are potential placebo responders has major implications for clinical practice and trial design. Catechol-O-methyltransferase (COMT), an important enzyme in dopamine catabolism plays a key role in processes associated with the placebo effect such as reward, pain, memory and learning. We hypothesized that the COMT functional val158met polymorphism, was a predictor of placebo effects and tested our hypothesis in a subset of 104 patients from a previously reported randomized controlled trial in irritable bowel syndrome (IBS). The three treatment arms from this study were: no-treatment (“waitlist”), placebo treatment alone (“limited”) and, placebo treatment “augmented” with a supportive patient-health care provider interaction. The primary outcome measure was change from baseline in IBS-Symptom Severity Scale (IBS-SSS) after three weeks of treatment. In a regression model, the number of methionine alleles in COMT val158met was linearly related to placebo response as measured by changes in IBS-SSS (p = .035). The strongest placebo response occurred in met/met homozygotes treated in the augmented placebo arm. A smaller met/met associated effect was observed with limited placebo treatment and there was no effect in the waitlist control. These data support our hypothesis that the COMT val158met polymorphism is a potential biomarker of placebo response.

Microsatellite variation in the equine MHC

Genes within the major histocompatibility complex (MHC) encode proteins involved in innate and adaptive immune responses. Genetic variation in this region can influence the immune response of an individual animal to challenges from a variety of pathogens; however, a complete documentation of genetic variation in the MHC is lacking for most domestic animals, including horses. To provide additional genetic markers for study of the horse MHC, or ELA (equine lymphocyte antigen), we identified 37 polymorphic microsatellite repeats in ELA and used these variations separately and together with published SNPs to investigate linkage disequilibrium (LD) and haplotype structure in a sample of Thoroughbred horses. ELA SNPs alone detected little LD, but microsatellites, either separately or combined with SNPs, revealed substantially more LD. A subset of markers in very high LD across the breadth of ELA may be predictive of structural polymorphisms or linked epistases that are important drivers of haplotype structure in Thoroughbreds.