Large-scale next generation sequencing based analysis of SLCO1B1 pharmacogenetics variants in the Saudi population.

BACKGROUND: SLCO1B1 plays an important role in mediating hepatic clearance of many different drugs including statins, angiotensin-converting enzyme inhibitors, chemotherapeutic agents and antibiotics. Several variants in SLCO1B1 have been shown to have a clinically significant impact, in relation to efficacy of these medications. This study provides a comprehensive overview of SLCO1B1 variation in Saudi individuals, one of the largest Arab populations in the Middle East. METHODS: The dataset of 11,889 (9,961 exomes and 1,928 pharmacogenetic gene panel) Saudi nationals, was used to determine the presence and frequencies of SLCO1B1 variants, as described by the Clinical Pharmacogenetic Implementation Consortium (CPIC). RESULTS: We identified 141 previously described SNPs, of which rs2306283 (50%) and rs4149056 (28%), were the most common. In addition, we observed six alleles [*15 (24.7%) followed by *20 (8.04%), *14 (5.86%), *5 (3.84%), *31 (0.21%) and *9 (0.03%)] predicted to be clinically actionable. Allele diplotype to phenotype conversion revealed 41 OATP1B1 diplotypes. We estimated the burden of rare, and novel predicted deleterious variants, resulting from 17 such alterations. CONCLUSIONS: The data we present, from one of the largest Arab cohorts studied to date, provides the most comprehensive overview of SLCO1B1 variants, and the subsequent OATP1B1 activity of this ethnic group, which thus far remains relatively underrepresented in available international genomic databases. We believe that the presented data provides a basis for further clinical investigations and the application of personalized statin drug therapy guidance in Arabs.

Impact of Four Ovine TMEM154 Haplotypes on Ewes during Multiyear Lentivirus Exposure.

Polypeptide variation encoded by the ovine transmembrane protein 154 gene (TMEM154) is associated with susceptibility to ovine lentivirus, the causative agent of Ovine Progressive Pneumonia (OPP) and Visna/Maedi. Our aim was to compare the four most prevalent TMEM154 haplotypes on the incidence of infection and ewe productivity during natural multiyear virus exposure. Prospective cohort studies were designed to test gene action and estimate effects of TMEM154 haplotypes encoding distinctive variant residues: K35 ("1"), I70 ("2"), ancestral ("3"), and A4del/M44 ("4"). Exposure consisted of co-mingling infected ewes at a rate greater than 30% with serological status evaluated every four months. For ewes with one or two copies of the highly susceptible haplotypes "2" and "3", the infection prevalence steadily increased to nearly 100% at 55 months. Haplotypes "2" and "3" were equally susceptible and dominant to haplotype "1". A difference was not detected (p < 0.53) in the magnitude of effect with haplotype combinations of "1" and "4". The ewe infection prevalence with "1,1"; "1,4"; and "4,4" was 10% to 40% at 55 months. The latter suggested that two copies of the K35 amino acid substitution ("1") were as effective as a homozygous TMEM154 "knockout" with the frame-shift deletion mutation ("4") in reducing infection susceptibility. When considering ewe reproductive performance, a difference was not detected when comparing haplotypes "2", and "3" to each other, or "1" and "4" to each other. Our study indicated that ewes with two copies of the severely truncated versions of TMEM154 ("4,4") had normal lamb productivity. Without complete understanding of the natural function of TMEM154 our recommendations to producers interested in using TMEM154 selection to reduce their flock's genetic predisposition to OPP are encouraged to increase the frequency of TMEM154 haplotype K35 ("1") since it encodes a full-length protein with minimal difference to the ancestral polypeptide.

Association of three SNPs in adiponectin gene with lipid traits of Tianzhu Black Muscovy (Cairina moschata).

Adiponectin plays a critically biological role in atherosclerosis, glucose utilization, lipid and carbohydrate metabolism, and triglyceride synthesis in animals and humans. However, little is known about the effect of adiponectin on lipid metabolism of the avian species. The aim of the preset study was to investigate the potential associations between adiponectin gene single nucleotide polymorphisms (SNPs) and the lipid traits in 348 females of Tianzhu Black Muscovy. Three novel SNPs (167G>A, 290T>C and 711G>A) were detected in adiponectin gene. 167G>A and 290T>C has linked very closely, and then 711G>A with 167G>A and 290T>C has no strong linkage disequilibrium, respectively. The Chi square test showed that allelic frequency and genotype frequency of two SNPs (167G>A and 711G>A) didn't agree with the Hardy-Weinberg equilibrium (P>0.05). Four haplotypes and nine diplotypes were formed on the three SNPs of adiponectin gene. Association analysis indicated that the 167G>A genotypes were strongly associated with intramuscular fat (IMF) of chest muscle and serum total cholesterol (TC) (P < 0.01); the 290T>C genotypes were strongly associated with IMF, TC, and serum triglyceride (TG) (P < 0.01); furthermore, the 711G>A genotypes were significantly associated with TG and TC (P < 0.05); the diplotypes were strongly associated with IMF, TC, and TG (P < 0.01). Therefore, three SNPs in adiponectin were potential markers for improving IMF in Muscovy ducks.

Polymorphism identification in GDF9 gene and its association analysis with reproduction traits in Jinghai Yellow chicken.

GDF9 (growth differentiation factor 9) belongs to the transforming growth factor-ß (TGF-ß) superfamily and plays an irreplaceable role in female fertility. To reveal its genetic effects on productivity performance in chickens, 373 Jinghai Yellow chickens were chosen randomly to detect SNPs in GDF9 by PCR-SSCP and DNA sequencing methods. Eventually, four SNPs (g.2053G > A, g.2275T > C, g.2338C > T, g.2420T > C) in total had been detected. Amongst them, g.2420T > C was first found significantly associated with reproduction trait in chickens and heterozygous type C2T2 had higher average egg weight at 300 days of age (AEWD300) than T2T2 (p < 0.01). Least squares analysis showed that age at first laying (AFE) of H1 and H1H1 chickens were significantly earlier than that of H7 and H7H7 ones, respectively (p < 0.05). H1H5 hens showed higher AEWD300 than H4H7 ones (p < 0.05). For total egg number at 300 days of age (END300), mean of H5H5 was significantly higher than that of H4H4 (p < 0.05). Hence, the study suggested that hybrid vigor at g.2420T > C could be utilized in practice. H1H1, H1H5 and H5H5 could be the dominant diplotypes for chicken breeding. The study may contribute to the breeding progress of productive chickens and supply reference for oviparous animal production practice.

Characterization of CYP2C19 pharmacogenetic variation in African populations and comparison with other global populations.

Background: CYP2C19 is important in the metabolism of clopidogrel and several antidepressants. This study aimed to characterize the distribution of CYP2C19 star alleles (haplotypes) across diverse African populations compared with global populations. Methods: CYP2C19 star alleles and diplotypes were called from high coverage genomes using the StellarPGx pipeline. Results: CYP2C19*1 (51%), *2 (17%) and *17 (22%) were the most common star alleles across African populations in this study. It was observed that 3% of African participants had potentially novel CYP2C19 haplotypes. Conclusion: This study supports the necessity for CYP2C19 pharmacogenetic testing in African and global clinical settings, as well as the importance of comprehensive star allele characterization in the African context.

Analysis of 1276 Haplotype-Resolved Genomes Allows Characterization of Cis- and Trans-Abundant Genes.

Many methods for haplotyping have materialized, but their application on a significant scale has been rare to date. Here we summarize analyses that were carried out in 1092 genomes from the 1000 Genomes Consortium and validated in an unprecedented number of 184 PGP genomes that have been experimentally haplotype-resolved by application of the Long-Fragment Read (LFR) technology. These analyses provided first insights into the diplotypic nature of human genomes and its potential functional implications. Thus, protein-changing variants were not randomly distributed between the two homologues of 18,121 autosomal protein-coding genes but occurred significantly more frequently in cis than in trans configurations in virtually each of the 1276 phased genomes. This resulted in global cis/trans ratios of ~60:40, establishing "cis abundance" as a universal characteristic of diploid human genomes. This phenomenon was based on two different classes of genes, a larger one exhibiting cis configurations of protein-changing variants in excess, so-called "cis-abundant" genes, and a smaller one of "trans-abundant" genes. These two gene classes, which together constitute a common diplotypic exome, were further functionally distinguished by means of gene ontology (GO) and pathway enrichment analysis. Moreover, they were distinguishable in terms of their effects on the human interactome, where they constitute distinct cis and trans modules, as shown with network propagation on a large integrated protein-protein interaction network. These analyses, recently performed with updated database and analysis tools, further consolidated the characterization of cis- and trans-abundant genes while expanding previous results. In this chapter, we present the key results along with the materials and methods to motivate readers to investigate these findings independently and gain further insights into the diplotypic nature of genes and genomes.

Characterization of POR haplotype distribution in African populations and comparison with other global populations.

Background & aim: POR is an enzyme that mediates electron transfer to enable the drug-metabolizing activity of CYP450 proteins. However, POR has been understudied in pharmacogenomics despite this vital role. This study aimed to characterize the genetic variation in POR across African populations and to compare the star allele (haplotype) distribution with that in other global populations. Materials & methods: POR star alleles were called from whole-genome sequencing data using the StellarPGx pipeline. Results: In addition to the common POR*1 and *28 (defined by rs1057868), five novel rare haplotypes were computationally inferred. No significant frequency differences were observed among the majority of African populations. However, POR*28 was observed at a higher frequency in individuals of non-African ancestry. Conclusion: This study highlights the distribution of POR alleles in Africa and across global populations with a view toward informing future precision medicine implementation.

HLA-G 3'UTR polymorphism diplotypes and soluble HLA-G plasma levels impact cervical cancer susceptibility and prognosis.

Background: Human leukocyte antigen G (HLA-G) is an immune checkpoint molecule with relevance in several cancers. The aim of this study was to evaluate the potential role of soluble HLA-G (sHLA-G), its genetic polymorphisms and its haplotype structure in the susceptibility and prognosis of primary cervical cancer in a Chinese Han population. Methods: We investigated sHLA-G plasma levels and 3' untranslated region (3'UTR) polymorphisms through ELISA and direct DNA sequencing, respectively, in cervical cancer patients (120 cases) and healthy control women (96 cases). The data were analyzed for associations using PowerMarker, Haploview, and GraphPad Prism. Results: In this study, 8 polymorphic sites, 16 haplotypes and 23 diplotypes in the HLA-G 3'UTR were identified in our study population. We observed that each pair of 8 polymorphic sites exhibited linkage disequilibrium. The heterozygote CT genotype at position +3422 (rs17875408) was more common in cervical cancer patients than in healthy women (OR=5.285, P<0.05). Haplotypes UTR-1, UTR-3, and UTR-7 accounted for more than 85% of both groups, but no significant difference was found. The frequency of the UTR-1/UTR-3 diplotype in patients was significantly higher than that in controls (P<0.05). In addition, we further observed that HLA-G 3'UTR polymorphisms may influence the sHLA-G plasma level in patients' peripheral blood, especially 14 bp Ins/Del (rs371194629) and +3142 C/G (rs1063320). A receiver operating characteristic (ROC) curve analysis showed that the sHLA-G level had good diagnostic performance in differentiating patients with cervical cancer from healthy women (AUC>0.7). Among patients, mean sHLA-G levels increased with increasing FIGO stages but were not related to the overall survival time. Conclusions: The results of the present study enhance our understanding of how HLA-G 3'UTR polymorphisms can influence the peripheral sHLA-G plasma level and play a key role in cervical carcinogenesis. This study further confirmed that sHLA-G may represent a novel plasma biomarker for the prognosis and potential therapeutic target of cervical cancer.

Genetic variation in CSF2 (5q31.1) is associated with longitudinal susceptibility to pediatric malaria, severe malarial anemia, and all-cause mortality in a high-burden malaria and HIV region of Kenya.

Plasmodium falciparum infections remain among the leading causes of morbidity and mortality in holoendemic transmission areas. Located within region 5q31.1, the colony-stimulating factor 2 gene (CSF2) encodes granulocyte-macrophage colony-stimulating factor (GM-CSF), a hematopoietic growth factor that mediates host immune responses. Since the effect of CSF2 variation on malaria pathogenesis remains unreported, we investigated the impact of two genetic variants in the 5q31.1 gene region flanking CSF2:g-7032 G > A (rs168681:G > A) and CSF2:g.64544T > C (rs246835:T > C) on the rate and timing of malaria and severe malarial anemia (SMA, Hb < 5.0 g/dL) episodes over 36 months of follow-up. Children (n = 1654, aged 2-70 months) were recruited from a holoendemic P. falciparum transmission area of western Kenya. Decreased incidence rate ratio (IRR) for malaria was conferred by inheritance of the CSF2:g.64544 TC genotype (P = 0.0277) and CSF2 AC/GC diplotype (P = 0.0015). Increased IRR for malaria was observed in carriers of the CSF2 AT/GC diplotype (P = 0.0237), while the inheritance of the CSF2 AT haplotype increased the IRR for SMA (P = 0.0166). A model estimating the longitudinal risk of malaria showed decreased hazard rates among CSF2 AC haplotype carriers (P = 0.0045). Investigation of all-cause mortality revealed that inheritance of the GA genotype at CSF2:g-7032 increased the risk of mortality (P = 0.0315). Higher risk of SMA and all-cause mortality were observed in younger children (P < 0.0001 and P = 0.0015), HIV-1(+) individuals (P < 0.0001 and P < 0.0001), and carriers of HbSS (P = 0.0342 and P = 0.0019). Results from this holoendemic P. falciparum area show that variation in gene region 5q31.1 influences susceptibility to malaria, SMA, and mortality, as does age, HIV-1 status, and inheritance of HbSS.

Haplotype-aware diplotyping from noisy long reads.

Current genotyping approaches for single-nucleotide variations rely on short, accurate reads from second-generation sequencing devices. Presently, third-generation sequencing platforms are rapidly becoming more widespread, yet approaches for leveraging their long but error-prone reads for genotyping are lacking. Here, we introduce a novel statistical framework for the joint inference of haplotypes and genotypes from noisy long reads, which we term diplotyping. Our technique takes full advantage of linkage information provided by long reads. We validate hundreds of thousands of candidate variants that have not yet been included in the high-confidence reference set of the Genome-in-a-Bottle effort.

Estrogen receptor α is not a candidate gene for metabolic syndrome in Caucasian elderly subjects.

OBJECTIVE: Variants of estrogen receptor α (ERα) have been associated with obesity, dyslipidemia, diabetes and blood pressure. The Middle East registers some of the highest rate of metabolic syndrome worldwide. The aim of this study is to investigate the relationship between metabolic syndrome, a clustered combination of these metabolic factors, and polymorphisms PvuII and XbaI of ERα in Lebanese Caucasian elderly overweight subjects. MATERIAL/METHODS: 250 Caucasian Lebanese unrelated elderly men and women, median age 71 years, were studied. ERα intronic polymorphisms variants, PvuII and XbaI diplotypes and genotypes, were examined. Associations with metabolic syndrome, defined by the American Heart Association/National Heart, Lung, and Blood Institute (AHA/NHLBI), and its components, namely high density lipoprotein (HDL), fasting glucose levels, blood pressure, and waist circumference were evaluated in regression models. RESULTS: ER α diplotypes and genotypes distributions were similar between participants with and without metabolic syndrome, in the overall group of subjects, and by gender. No consistent associations between the diplotypes and genotypes tested and metabolic syndrome, or its components, could be detected. CONCLUSIONS: Genetic variants in ERα were not associated with metabolic syndrome or its components, in a group of 250 Lebanese Caucasian elderly participants, a group with a high prevalence of metabolic syndrome.