Linkage Disequilibrium Test Research Articles

Towards pharmacogenomics-guided tuberculosis (TB) therapy: N-acetyltransferase-2 genotypes among TB-infected Kenyans of mixed ethnicity

BackgroundThough persons of African descent have one of the widest genetic variability, genetic polymorphisms of drug-metabolising enzymes such as N-Acetyltransferase-2 (NAT2) are understudied. This study aimed to identify prevalent NAT2 single nucleotide polymorphisms (SNPs) and infer their potential effects on enzyme function among Kenyan volunteers with tuberculosis (TB) infection. Genotypic distribution at each SNP and non-random association of alleles were evaluated by testing for Hardy-Weinberg Equilibrium (HWE) and Linkage Disequilibrium (LD).MethodsWe isolated genomic DNA from cryopreserved Peripheral Blood Mononuclear Cells of 79 volunteers. We amplified the protein-coding region of the NAT2 gene by polymerase chain reaction (PCR) and sequenced PCR products using the Sanger sequencing method. Sequencing reads were mapped and aligned to the NAT2 reference using the Geneious software (Auckland, New Zealand). Statistical analyses were performed using RStudio version 4.3.2 (2023.09.1 + 494).ResultsThe most frequent haplotype was the wild type NAT2*4 (37%). Five genetic variants: 282C > T (NAT2*13), 341 T > C (NAT2*5), 803A > G (NAT2*12), 590G > A (NAT2*6) and 481C > T (NAT2*11) were observed with allele frequencies of 29%, 18%, 6%, 6%, and 4% respectively. According to the bimodal distribution of acetylation activity, the predicted phenotype was 76% rapid (mainly consisting of the wildtype NAT2*4 and the NAT2*13A variant). A higher proportion of rapid acetylators were female, 72% vs 28% male (p = 0.022, odds ratio [OR] 3.48, 95% confidence interval [CI] 1.21 to 10.48). All variants were in HWE. NAT2 341 T > C was in strong complete LD with the 590G > A variant (D′ = 1.0, r2 = − 0.39) but not complete LD with the 282C > T variant (D′ = 0.94, r2 = − 0.54).ConclusionThe rapid acetylation haplotypes predominated. Despite the LD observed, none of the SNPs could be termed tag SNP. This study adds to the genetic characterisation data of African populations at NAT2, which may be useful for developing relevant pharmacogenomic tools for TB therapy. To support optimised, pharmacogenomics-guided TB therapy, we recommend genotype-phenotype studies, including studies designed to explore gender-associated differences.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species.

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Drug resistance markers in Plasmodium vivax isolates from a Kanchanaburi province, Thailand between January to May 2023.

Plasmodium vivax has become the predominant species in the border regions of Thailand. The emergence and spread of antimalarial drug resistance in P. vivax is one of the significant challenges for malaria control. Continuous surveillance of drug resistance is therefore necessary for monitoring the development of drug resistance in the region. This study aims to investigate the prevalence of the mutation in the P. vivax multidrug resistant 1 (Pvmdr1), dihydrofolate reductase (Pvdhfr), and dihydropteroate synthetase (Pvdhps) genes conferred resistance to chloroquine (CQ), pyrimethamine (P) and sulfadoxine (S), respectively. 100 P. vivax isolates were obtained between January to May 2023 from a Kanchanaburi province, western Thailand. Nucleotide sequences of Pvmdr1, Pvdhfr, and Pvdhps genes were amplified and sequenced. The frequency of single nucleotide polymorphisms (SNPs)-haplotypes of drug-resistant alleles was assessed. The linkage disequilibrium (LD) tests were also analyzed. In Pvmdr1, T958M, Y976F, and F1076L, mutations were detected in 100%, 21%, and 23% of the isolates, respectively. In Pvdhfr, the quadruple mutant allele (I57R58M61T117) prevailed in 84% of the samples, followed by (L57R58M61T117) in 11%. For Pvdhps, the double mutant allele (G383G553) was detected (48%), followed by the triple mutant allele (G383M512G553) (47%) of the isolates. The most prevalent combination of Pvdhfr (I57R58M61T117) and Pvdhps (G383G553) alleles was sextuple mutated haplotypes (48%). For LD analysis, the association in the SNPs pairs was found between the intragenic and intergenic regions of the Pvdhfr and Pvdhps genes. The study has recently updated the high prevalence of three gene mutations associated with CQ and SP resistance. Genetic monitoring is therefore important to intensify in the regions to further assess the spread of drug resistant. Our data also provide evidence on the distribution of drug resistance for the early warning system, thereby threatening P. vivax malaria treatment policy decisions at the national level.

Effects of Recurrent Selection on Population Structure and Allele Frequencies in the M3S Maize Population

The effects of four cycles of recurrent selection on the allele frequencies of simple sequence repeat (SSR) markers and population structure were examined in the Maksimir 3 Synthetic (M3S) maize population (Zea mays L.). Genotyping of 32 plants from each selection cycle at 38 SSR loci revealed that the mean number of alleles per locus and the mean expected heterozygosity were preserved across selection cycles, indicating the maintenance of sufficient genetic variability in the population required for future genetic gain. The Waples test of selective neutrality revealed that genetic drift was the main force in changing allele frequencies in the population. The proportion of selectively non-neutral loci in single cycles of selection varied between 16% and 37%. Some non-neutral loci shared the same genomic locations with previously published QTLs controlling important agronomic traits. An analysis of molecular variance revealed that 5.6% of the genetic variation occurred among and 94.4% within cycle populations. Between 5% and 29% of loci were found to be in a significant Hardy–Weinberg (HW) disequilibrium, with the majority showing an excess of homozygosity. The excess of homozygosity at several loci was highly consistent across cycle populations, suggesting positive assortative mating as a possible cause of the observed HW disequilibrium. Linkage disequilibrium (LD) tests revealed that the M3S population was essentially in linkage equilibrium. The proportion of pairs of loci in significant LD varied from 0.1% to 1.8% across selection cycles, probably due to the effects of genetic drift and epistatic selection.

Interplay Between Genes from Peroxisome Biogenesis and Integrated Stress Response Pathways is Significantly Associated with Alzheimer’s Disease

AbstractBackgroundPeroxisomes are single‐membrane organelles involved among other things in ether‐phospholipids production that plays an important role in axon myelination in the central nervous system. Disruptions in cytosolic receptor PEX5L (involved in protein import into peroxisome matrix) may lead to peroxisome malfunctioning and disturbance of myelin homeostasis, which, in turn, may contribute to developing Alzheimer’s disease (AD). Experimental studies showed that peroxisome dysfunction is associated with activation of the integrative stress response (ISR). Recent studies provide evidence that the EIF2AK4 gene, a sensor of amino acid deprivation in ISR, is involved in the AD regulation. It is, however, unclear whether interplay between PEX5L and EIF2A4 is associated with AD.MethodTo test the association of the interplay between PEX5L and EIF2AK4 genes with the late‐onset‐AD, we analyzed the Health and Retirement Study (HRS) data using logistic regression model with the interaction term. Education, sex, smoking status, and first five principal components were included as observed covariates. The binary AD outcome was defined as “0” if AD onset was at any age, and as “1” if a person survived to age 85+ and did not get AD during the life. Linkage Disequilibrium test and SNP‐clumping procedure were used to reduce number of tests in SNPxSNP interaction analysis.ResultInteractions between SNPs in fifteen SNP pairs (one SNP in EIF2AK4 gene and another in PEX5L gene) were significantly associated with decreased odds of AD, with effect sizes (β coefficients before the interaction term) ranging between 0.62‐0.69, and p‐values between 1.64E‐04 ‐ 2.12E‐05, with Bonferroni correction threshold = 1.67E‐04. Six SNPs from PEX5L (rs1404271, rs11708154, rs59018122, rs9836757, rs76133869, rs6443674) and five SNPs from EIF2AK4 (rs3207297, rs2279579, rs2291625, rs12591842, rs117584784) participated in such interactions.ConclusionIn this study, interactions between SNPs in EIF2AK4 and PEX5L genes were significantly associated with lower odds of AD in HRS participants. Our findings support the idea that the interplay between genes involved in cellular stress response pathway and peroxisome biogenesis can play a protective role in AD. These results also demonstrate how the knowledge from experimental research can be translated to human applications population‐based studies of AD.

HLA-G 3'UTR haplotype analyses in HCV infection and HCV-derived cirrhosis, hepatocellular carcinoma and fibrosis.

Hepatitis C virus (HCV) infection is a major cause of chronic liver disease. Chronic HCV infection is also an important cause of hepatic fibrosis, cirrhosis and hepatocellular carcinoma (HCC). HCV has the capacity to evade immune surveillance by altering the host immune response. Moreover, variations in immune-related genes can lead to differential susceptibility to HCV infection as well as interfere on the susceptibility to the development of hepatic fibrosis, cirrhosis and HCC. The human leucocyte antigen G (HLA-G) gene codes for an immunomodulatory protein known to be expressed in the maternal-foetal interface and in immune-privileged tissues. The HLA-G 3' untranslated region (3'UTR) is important for mRNA stability, and variants in this region are known to impact gene expression. Studies, mainly focusing in a 14bp insertion/deletion polymorphism, have correlated HLA-G 3'UTR with susceptibility to viral infections, but other polymorphic variants in the HLA-G 3'UTR might also affect HCV infection as they are inherited as haplotypes. The present study evaluated HLA-G 3'UTR polymorphisms and performed linkage disequilibrium test and haplotype assembly in 286 HCV infected patients who have developed fibrosis, cirrhosis or HCC, as well as in 129 healthy control subjects. Haplotypes UTR-1, UTR-2 and UTR-3 were the most observed in HCV+ patients, in the frequencies of 0.276, 0.255 and 0.121, respectively. No statistically significant difference was observed between HCV+ and control subjects, even when patients were grouped according to outcome (HCC, cirrhosis or fibrosis). Despite that, some trends in the results were observed, and therefore, we cannot rule out the possibility that variants associated to high HLA-G expression can be involved in HCV infection susceptibility.

MicroRNA in Alzheimer’s Disease: Results of interaction analysis

AbstractBackgroundMiRNA are multifunctional molecules involved in regulating posttranscriptional expression of genes. Several studies showed that miR‐107 influences expression of the BACE1 gene in the process of Aβ production. Extracellular accumulation of Aβ is well‐known hallmark of Alzheimer’s disease (AD). The product NLRP3 gene is involved in inflammasome that induces neuroinflammation a pathological process accompanying AD. The role of miR‐107 in regulation of NRLP3 is unknown.MethodWe used logistic regression model with the interaction term to test the hypothesis that interaction of SNPs from miR‐107 with SNPs from the NLRP3 gene is associated with AD. For this we used data on genotyped individuals collected in Framingham Health Study (FHS). The Cardiovascular Health Study (CHS) was used to replicate the research findings. Linkage Disequilibrium test and clumping procedure were used to reduce the number of tests and make Bonferroni correction less conservative.ResultIn the analysis of FHS data we found that the association of the rs2926616 SNP from the miR‐107 miRNA with the rs10754557 SNP from the NLRP3 gene is positive and statistically significant. The effect size β = 0.356, p‐value 6.71E‐03. This value is smaller than the p‐value threshold calculated in accordance with Bonferroni correction 7.14E‐03. This result is replicated in the analysis of CHS data with β = 0.289, and p‐value 4.76E‐02. Using FHS data we also tested the role of interaction of SNPs from the miR‐107 with SNPs from the TLR2 gene in AD and found a pair of SNPs whose interaction has statistically significant association with AD. This finding was replicated in the analysis of the CHS data where an interaction of a pair of SNPs also showed significant association with AD. Detected SNPs from this pair were in linkage disequilibrium with the pair SNPs from the CHS data.ConclusionThe results of these analyses indicate that one miRNA may influence several components of AD pathological development. This property may have important implications for developing efficient medications against AD. More studies are needed to evaluate miRNA involvement in different pathways of AD regulation, as well as in identifying factors capable in affecting such miRNA.

SSR Genotypes of the Puccinia triticina in 15 Provinces of China Indicate Regional Migration in One Season from East to West and South to North

Leaf rust of wheat caused by Puccinia triticina (Pt) is one of the most common fungal diseases in the southwest and northwest of China, the middle and lower reaches of the Yangtze River, and the southern part of the Huang-Huai-Hai river basin. Using 13 simple sequence repeat (SSR) markers, we systematically revealed the genotypic diversities, population differentiation and reproduction of Pt isolates in 15 wheat-producing areas in China. A total of 622 isolates were divided into 3 predominant populations, including the eastern Pt populations, consisting of Pt samples from 8 eastern provinces of Beijing, Hebei, Shanxi, Shaanxi, Anhui, Shandong, Henan, and Heilongjiang; the 4 western Pt populations from Gansu, Qinghai, Sichuan, and Inner Mongolia provinces; and the bridge Pt populations including Jiangsu, Hubei, and Yunnan, which communicated the other 2 populations as a “bridge”. The pathogen transmission of eastern Pt populations was more frequent than western Pt populations. The linkage disequilibrium test indicated that the whole Pt population was in a state of linkage disequilibrium. However, populations of Beijing, Hebei, Shaanxi, Jiangsu, Henan, and Heilongjiang provinces showed obvious linkage equilibrium, while the five provinces of Qinghai, Hubei, Anhui, Shandong, and Inner Mongolia supported clonal modes of reproduction.

SSR genotypes of Puccinia triticina in 10 provinces of China reveal a possible source of the wheat leaf rust fungus

AbstractWheat leaf rust is a fungal disease caused by Puccinia triticina and spread by wind. In recent years, with changes in the climate and farming systems, damage caused by wheat leaf rust has become increasingly severe. To study the genetic structure of P. triticina populations in China and provide a basis for control of wheat leaf rust, collections of P. triticina were obtained from Xinjiang, Inner Mongolia, Hebei, Henan, Jiangsu, Anhui, Hubei, Sichuan, Guizhou and Yunnan provinces in 2019. A total of 177 single uredinial isolates were tested with 13 simple‐sequence repeat (SSR) primers. The results showed that there was frequent gene exchange within P. triticina populations in China, but a distinct population divergence was also revealed between those in Xinjiang and Yunnan provinces. Results also suggested that Inner Mongolia might be a source of P. triticina for other regions. In addition, the linkage disequilibrium test indicated that the P. triticina populations reproduced by clonal production of urediniospores.

Genetic diversity and population structure of Tongcheng pigs in China using whole-genome SNP chip.

Tongcheng (TC) pigs, distinguished by their superior meat quality, are a Chinese indigenous pig breed. Recently, the genetic resources of TC pigs are under tremendous threat due to the introduction of cosmopolitan pig breeds and African swine fever disease. To promote their management and conservation, the present study assessed genetic diversity and population structure of TC pigs using single nucleotide polymorphism (SNP) markers. A total of 26, 999 SNPs were screened from 51, 315 SNPs in 68 TC pigs. The multi-dimensional scaling (MDS) analysis and neighbor-joining tree revealed that all 68 pigs were from a purebred population. The effective population size decreased over time, and it was 96 prior to generation 20. Both linkage disequilibrium (LD) and neutrality test indicated a low selection of TC pigs with average LD value of 0.15 ± 0.23. Genetic diversity results exhibited a minor allele frequency (MAF) of 0.23, observed heterozygosity (HO) of 0.32, expected heterozygosity (He) of 0.31, and nucleotide diversity (Pi) of 0.31. All these parameters indicated a remarkably high genetic diversity of TC pigs. Additionally, 184 runs of homozygosity (ROH) segments were detected from the whole genome of TC pigs with an average ROH length of 23.71Mb, ranging from 11.26Mb to 69.02 Mb. The highest ROH coverage was found on chromosome 1 (10.12%), while the lowest was on chromosome 18 (1.49%). The average inbreeding coefficients based on ROH (FROH) was 0.04%. Fourteen ROH islands containing 240 genes were detected on 9 different autosomes. Some of these 240 genes were overlapped with the genes related to biological processes such as immune function, reproduction, muscular development, and fat deposition, including FFAR2, FFAR4, MAPK8, NPY5R, KISS1, and these genes might be associated with such traits as meat quality and disease resistance in TC pigs. Taken together, population structure and genetic diversity results suggested that the TC pig represented a valuable genetic resource. However, TC pig breed conservation program remains to be further optimized to ensure adequate genetic diversity and avoid inbreeding depression. Our findings provide theoretical basis for formulating management and conservation strategies for TC pigs.

High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations.

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals. Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA). Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies. Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

Phytophthora palmivora from Sulawesi and Java Islands, Indonesia, reveals high genotypic diversity and lack of population structure

Phytophthora palmivora is the causal agent of cocoa black pod disease, one of the primary diseases of cocoa in Indonesia. A better understanding of P. palmivora population genetics is needed to aid the development of relevant disease management strategies. This study is the first population genetic study of P. palmivora in Indonesia using microsatellite markers based on the alleles genotyping method. The microsatellite markers were used to determine the genotype of 44 P. palmivora isolates from Sulawesi (24) and Java (20) islands. The total number of observed multilocus genotypes (MLG) from both populations was 34. The genotypic diversity of P. palmivora from Sulawesi (2.90; 16.0; 0.938) and Java (2.76; 14.3; 0.930) islands was high as seen from Shannon's diversity index (H), Stoddart and Taylor's Index (G), and Simpson's Index (λ) respectively. Evenness and Nei's unbiased gene diversity exhibited similarly high levels from both populations. The linkage disequilibrium test indicated that sexual recombination occurred in the Java population (P = 0.312). Analysis of molecular variance (AMOVA) and Bayesian clustering revealed five genetic clusters, and isolates from both islands were evenly distributed across the five gene clusters. All genetic diversity was from within individuals. P. palmivora from Sulawesi and Java showed a high genotypic diversity but a lack of genetic differentiation among populations (Fst = 0.006). Both populations formed one highly diverse group. Minimum spanning network analysis showed no particular grouping of MLGs, and shared MLGs from both populations indicated long-distance migration of P. palmivora facilitated by human activities.

Genetic variations in ABCA1/G1 associated with plasma lipid levels and risk of ischemic stroke

BackgroundATP binding cassette transporters ABCA1 and ABCG1 play a crucial role in cholesterol efflux and reverse cholesterol transport (RCT), thereby rendering ischemic stroke (IS) susceptibility. Variants of ABCA1/G1 have been implicated in etiology of IS. This study aimed to investigate the association between single-nucleotide polymorphisms (SNPs) of ABCA1/G1 with plasma lipid variability and the risk of IS in Chinese Han Population. MethodsTotally 249 IS patients and 226 healthy controls were enrolled and 10 SNPs of ABCA1/G1 were screened for genotyping by kompetitive allele-specific polymerase chain reaction (KASP) and validated by sanger sequencing. The logistic regression analysis was performed to identify risk alleles of IS and appropriate genetic model. The genetic risk scores (GRS) and predicted risks for all individuals was computed. Based on different plasma lipid levels, we applied stratified analyses for subgroups. Linkage disequilibrium (LD) test was used to explore different functional haplotype combinations. Association between specific allele or genotype of the SNPs of ABCA1/G1 and plasma lipid or lipoproteins levels were also investigated. ResultsBesides total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C), significant differences of clinical data were observed between IS and control group. The rare GG genotype frequencies of rs4149338 on ABCA1 was higher in IS patients than those in controls (11.4%, 4.6%, respectively, P = 0.037). Frequencies of rs57137919 on ABCG1 for rare AA genotype was lower in IS group than those in control group (4.6%, 13.3%, respectively, P = 0.030). GRS showed ability to discriminate IS patients and controls (AUC = 0.633, P < 0.001). Haplotype A-A (rs4149339-rs4149338) was correlated with reduced risk of IS (P = 0.023). Association analysis showed that subjects with rare AA genotype of rs57137919 had the lowest LDL-C levels while rare GG genotype of rs4149338 had lower TC level than those with AA genotype. The mRNA expression of ABCG1 was higher in IS patients, especially in the patients with frequent GG genotype of rs57137919, and was positively correlated with higher ABCG1 expression level and plasma LDL-C level. ConclusionsPolymorphisms of ABCA1/G1 associated with varieties of plasma lipid levels and risk of IS.

Human leukocyte antigen allele and haplotype frequencies in Singapore bone marrow donors and cord blood units.

We describe the allele and haplotype frequencies seen in a volunteer unrelated bone marrow donor registry, a public cord blood bank, and donor/recipient samples processed by the Health Sciences Authority (HSA) in Singapore. Historical human leukocyte antigen (HLA) typing reports were anonymized and combined. They were checked for HLA typing nomenclature discrepancies or ambiguities using the HLA-net UNIFORMATE tool, and for analysis, the validated data were subsequently separated into Chinese, Malay, Indian, and "Others," according to the race classification system used in Singapore. Individual ethnic allele and haplotype frequencies were calculated with the HLA-net GENE[RATE] pipeline using basic statistics. The Basic Statistics Tool of HLA-net was used to estimate haplotype frequency using an expectation maximization algorithm, given a set of multi-allelic data pairs for a given HLA locus. The outputs downloaded from the site comprised plain text files with haplotype frequency estimates, results of a global linkage disequilibrium test, and standardized residuals (stdres) corresponding to deviations from expected frequencies. HLA typing results from 59,186 individuals met the inclusion criteria, yielding 118,372 analyzable alleles. In our study population, the haplotype A*33:03-B*58:01-C*03:02-DRB1*03:01～DQB1*02:01:01G with a frequency of 4.91% was the most common. This haplotype was also the most common among Singaporean Chinese donors. Consistent with the predominant Chinese population, haplotypes with a frequency greater than 1% were also the most frequently observed haplotypes in the Singaporean population. In the Malay donor population, the most common haplotype was A*33:03～B*44:03～C*07:01:01G～ DRB1*07:01-DQB1*02:01:01G, with a frequency of 3.41%, whereas within the Indian donor population, the most common haplotype was A*01:01-B*57:01-C*06:02～DRB1*07:01-DQB1*03:03, with a frequency of 3.42%. Haplotype diversity and composition statistics within donor pools provide HLA background data required for the targeted recruitment of donors to support the hematopoietic stem cell donor requirements of the country. These data may be used in the future to devise donor recruitment strategies for optimizing the donor pool through targeted publicity and accruals.

Identification and Characterization of Nine Novel X-Chromosomal Short Tandem Repeats on Xp21.1, Xq21.31, and Xq23 Regions.

The application of X-chromosomal short tandem repeats (X-STRs) has been recognized as a powerful tool in complex kinship testing. To support further development of X-STR analysis in forensic use, we identified nine novel X-STRs, which could be clustered into three linkage groups on Xp21.1, Xq21.31, and Xq23. A multiplex PCR system was built based on the electrophoresis. A total of 198 unrelated Shanghai Han samples along with 168 samples from 43 families was collected to investigate the genetic polymorphism and forensic parameters of the nine loci. Allele numbers ranged from 5 to 12, and amplicon sizes ranged from 146 to 477 bp. The multiplex showed high values for the combined power of discrimination (0.99997977 in males and 0.99999999 in females) and combined mean exclusion chances (0.99997918 and 0.99997821 in trios, 0.99984939 in duos, and 0.99984200 in deficiency cases). The linkage between all pairs of loci was estimated via Kosambi mapping function and linkage disequilibrium test, and further investigated through the family study. The data from 43 families strongly demonstrated an independent transmission between LGs and a tight linkage among loci within the same LG. All these results support that the newly described X-STRs and the multiplex system are highly promising for further forensic use.

SNPInt-GPU: Tool for Epistasis Testing with Multiple Methods and GPU Acceleration.

We present SNPInt-GPU, a software providing several methods for statistical epistasis testing. SNPInt-GPU supports GPU acceleration using the Nvidia CUDA framework, but can also be used without GPU hardware. The software implements logistic regression (as in PLINK epistasis testing), BOOST, log-linear regression, mutual information (MI), and information gain (IG) for pairwise testing as well as mutual information and information gain for third-order tests. Optionally, r2 scores for testing for linkage disequilibrium (LD) can be calculated on-the-fly. SNPInt-GPU is publicly available at GitHub. The software requires a Linux-based operating system and CUDA libraries. This chapter describes detailed installation and usage instructions as well as examples for basic preliminary quality control and analysis of results.

Plasmodium vivax Cysteine-Rich Protective Antigen Polymorphism at Exon-1 Shows Recombination and Signatures of Balancing Selection.

Plasmodium vivax Cysteine-Rich Protective Antigen (CyRPA) is a merozoite protein participating in the parasite invasion of human reticulocytes. During natural P. vivax infection, antibody responses against PvCyRPA have been detected. In children, low anti-CyRPA antibody titers correlated with clinical protection, which suggests this protein as a potential vaccine candidate. This work analyzed the genetic and amino acid diversity of pvcyrpa in Mexican and global parasites. Consensus coding sequences of pvcyrpa were obtained from seven isolates. Other sequences were extracted from a repository. Maximum likelihood phylogenetic trees, genetic diversity parameters, linkage disequilibrium (LD), and neutrality tests were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. In 22 sequences from Southern Mexico, two synonymous and 21 nonsynonymous mutations defined nine private haplotypes. These parasites had the highest LD-R2 index and the lowest nucleotide diversity compared to isolates from South America or Asia. The nucleotide diversity and Tajima’s D values varied across the coding gene. The exon-1 sequence had greater diversity and Rm values than those of exon-2. Exon-1 had significant positive values for Tajima’s D, β-α values, and for the Z (HA: dN > dS) and MK tests. These patterns were similar for parasites of different origin. The polymorphic amino acid residues at PvCyRPA resembled the conformational B-cell peptides reported in PfCyRPA. Diversity at pvcyrpa exon-1 is caused by mutation and recombination. This seems to be maintained by balancing selection, likely due to selective immune pressure, all of which merit further study.

Genetic Characterization of the Yugoslavian Shepherd Dog – Sharplanina, a Livestock Guard Dog from the Western Balkans

Abstract Yugoslavian Shepherd Dog - Sharplanina (YSD) is a livestock guard dog from the Western Balkans present in this region over a long time, but recognized by the Fédération Cynologique Internationale as a distinct breed as late as 1957. However, the information regarding the origin and the size of the breed’s foundation stock is still lacking. In order to contribute towards better understanding of the genetic make-up of the YSD and its foundation stock, we re-analyzed previously generated genetic profiles of 94 registered YSD dogs assessed with nine nuclear microsatellites. Studied individuals comprised 90 unrelated dogs and two pairs of full-sibs, sampled at four sampling sites: three dog shows and at a military training centre for dogs in Serbia. We supported earlier findings on high levels of genetic diversity in YSD (HE =0.728±0.027) and lack of inbreeding, and revealed substructure of the breed because we found two distinct gene pools in the Bayesian clustering analysis, indicated also by the excess of homozygotes (i.e., Wahlund effect) and outcomes of other analyses: linkage disequilibrium tests, Neighbour-Joining tree, principal coordinates and two-dimensional scaling analyses. The two gene pools were almost equally represented at each sampling site. One gene pool was composed of individuals with high genetic integrity, while the other gene pool was characterized with admixed ancestry, developed possibly via hybridization with native breeding stock outside the registry system, other breeds, such as the Caucasian Shepherd, and/or individuals admixed with wolves. Thus, we demonstrate rather complex and diverse ancestry implying a genetically heterogeneous foundation stock of the YSD.

Association between platelet-derived growth factor-C single nucleotide polymorphisms and nonsyndromic cleft lip with or without cleft palate in Western Chinese population

To explore the association between two single nucleotide polymorphisms (SNPs), namely, rs4691383 and rs7667857, in the platelet-derived growth factor-C (PDGF-C) gene, the genotypes, environmental exposure factors, and nonsyndromic cleft lip with or without cleft palate (NSCL/P) in Western Chinese population. A total of 268 case-parent trios were selected, and two SNPs (rs4691383 andrs7667857) were genotyped by using polymerase chain reaction and restriction enzyme fragment length polymorphic method and direct sequencing method. Hardy-Weinberg equilibrium, linkage disequilibrium test, transmission disequilibrium test, and haplotype analysis were conducted to analyze the data. Meanwhile, the questionnaires on the epidemiology of cleft lip and palate filled by the included samples were collected, and the interaction between the genotypes of the two SNPs and environmental exposure factors was assessed by conditional logistic regression. The A allele at rs4691383 and the G allele at rs7667857 of PDGF-C gene were over-transmitted for NSCL/P (P<0.05). No interaction effect was observed between the three environmental exposure factors (history of smoking/passive smoking, folic acid supplementation, and long-term inhalation of harmful environmental gases) and the PDGF-C genotypes among NSCL/P (P>0.05). The rs4691383 and rs7667857 at PDGF-C gene are closely related to the occurrence of NSCL/P in Western Chinese population. However, the interaction between environmental exposure factors and PDGF-C genotypes is not obvious in the occurrence of NSCL/P.

Development and characterization of 16 novel microsatellite markers by Transcriptome sequencing for Angelica dahurica and test for cross-species amplification

BackgroundAngelica dahurica (Apiaceae) is an important herb in traditional Chinese medicine. Because of its important medicinal and economic values, its wild resources were over-exploited and increasingly reduced. Meanwhile, the diversity of cultivars of A. dahurica has decreased as a result of long-term artificial cultivation. However, there are no population genetics studies of natural A. dahurica reported yet, especially for using microsatellite markers (SSRs) to investigate population genetics of the species.ResultsSixteen polymorphic EST-SSR loci were isolated from A. dahurica with transcriptome sequencing technology (RNA-Seq). The number of alleles varied from 2 to 15 per polymorphic locus over populations with the observed and expected heterozygosities ranging from 0.000 to 1.000 and from 0.000 to 0.829, respectively. Significant deviations from Hardy–Weinberg equilibrium were observed at 8 loci. Tests of linkage disequilibrium showed 11 informative locus pairs were significant across all populations. Cross-species amplification showed that 14 out of 16 SSR loci have the transferability in cultivar-A. dahurica cv. ‘Hangbaizhi’ and A. decursiva.ConclusionsThe 16 newly developed loci microsatellite primers with RNA-Seq will be useful for further investigating population genetics of A. dahurica, cultivars and other members of this genus.