Genomic analyses inform on migration events during the peopling of Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Functional divergence of Helicobacter pylori related to early gastric cancer.

Helicobacter pylori is an extra macro- and microdiverse bacterial species, but where and when diversity arises is not well-understood. To test whether a new environment accelerates H. pylori genetic changes for quick adaptation, we have examined the genetic and phenotypic changes in H. pylori obtained from different locations of the stomach from patients with early gastric cancer (ECG) or chronic gastritis (CG). Macroarray analysis did not detect differences in genetic content among all of the isolates obtained from different locations within the same stomach of patients with EGC or CG. The extent and types of functional diversity of H. pylori isolates were characterized by 2-D difference gel electrophoresis (2D DIGE). Our analysis revealed 32 differentially expressed proteins in H. pylori related to EGC and 14 differentially expressed proteins in H. pylori related to CG. Most of the differentially expressed proteins belong to the antioxidant protein group (SodB, KatA, AphC/TsaA, TrxA, Pfr), tricarbon acid cycle proteins (Idh, FrdA, FrdB, FldA, AcnB) and heat shock proteins (GroEL and ClpB). We conclude that H. pylori protein expression variability is mostly associated with microorganism adaptation to morphologically different parts of the stomach, which has histological features and morphological changes due to pathological processes; gene loss or acquisition is not involved in the adaptation process.

Application of matrix-assisted laser desorption/ionization time-of-flight mass spectrometry for the study of Helicobacter pylori.

The characteristics of matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry based investigation of extremely variable bacteria such as Helicobacter pylori were studied. H. pylori possesses a very high natural variability. Accurate tools for species identification and epidemiological characterization could help the scientific community to better understand the transmission pathways and virulence mechanisms of these bacteria. Seventeen clinical as well as two laboratory strains of H. pylori were analyzed by the MALDI Biotyper method for rapid species identification. Mass spectra collected were found containing 7-13 significant peaks per sample, and only six protein signals were identical for more than half of the strains. Four of them could be assigned to ribosomal proteins RL32, RL33, RL34, and RL36. The reproducible peak with m/z 6948 was identified as a histidine-rich metal-binding polypeptide by tandem mass spectrometry (MS/MS). In spite of the evident protein heterogeneity of H. pylori the mass spectra collected for a particular strain under several cultivations were highly reproducible. Moreover, all clinical strains were perfectly identified as H. pylori species through comparative analysis using the MALDI Biotyper software (Bruker Daltonics, Germany) by pattern matching against a database containing mass spectra from different microbial strains (n = 3287) including H. pylori 26695 and J99. The results of this study allow the conclusion that the MALDI-TOF direct bacterial profiling is suited for H. pylori identification and could be supported by mass spectra fragmentation of the observed polypeptide if necessary.

Comparative genome analysis of Ureaplasma parvum clinical isolates.

Ureaplasma parvum colonizes human mucosal surfaces, primarily in the respiratory and urogenital tracts, causing a wide spectrum of diseases, from non-gonococcal urethritis to pneumonitis in immunocompromised hosts. Although the basis for these diverse clinical outcomes is not yet understood, more severe disease may be associated with strains harboring a certain set of strain-specific genes. To investigate this, whole genome DNA macroarrays were constructed and used to assess genomic diversity in 10 U. parvum clinical strains. We found that 7.6% of U. parvum genes were dispersed into one or more strains, thus defining a minimal functional core of 538 U. parvum genes. Most of the strain-specific genes (79%) were of unknown function and were unique to U. parvum. Four hypervariable plasticity regions were identified in the genome containing 93% of the variability in the gene pool (UU32-UU33, UU145-UU170, UU440-UU447 and UU527-UU529). We hypothesized that one of them (UU145-UU170) was a pathogenicity island in U. parvum and we characterized it. Thus, we propose that the clinical outcome of U. parvum infection is probably associated with this newly identified pathogenicity island.

Coexpressionfinder: a new algorithm for finding groups of coexpressed genes.

RESULTS: A new algorithm is developed which is intended to find groups of genes whose expression values change in a concordant manner in a series of experiments with DNA arrays. This algorithm is named as CoexpressionFinder. It can find more complete and internally coordinated groups of gene expression vectors than hierarchical clustering. Also, it finds more genes having coordinated expression. The algorithm's design allows parallel execution. AVAILABILITY: The algorithm is implemented as a Java application which is freely available at: http://www.bioinformatics.ru/cf/index.jsp and http://bioinformatics.ru/cf/index.jsp.

Phenylethynylpyrene-labeled oligonucleotide probes for excimer fluorescence SNP analysis of 23S rRNA gene in clarithromycin-resistant Helicobacter pylori strains.

The use of phenylethynylpyrene excimer forming pair in the design of specific fluorescent probes for determination of A2144G (A2143G and/or A2143C) mutations in 23S rRNA gene of Helicobacter pylori is described. Analysis of fluorescence spectra of model duplexes revealed optimal positions of fluorophore residues in the probe sequences for maximum efficiency of SNP detection. Application of excimer forming probes for analysis of DNA samples isolated from natural bacterial strains of H. pylori was demonstrated.

A Novel Mutation in a Kazakh Family with X-Linked Alport Syndrome.

Alport syndrome is a genetic condition that results in hematuria, progressive renal impairment, hearing loss, and occasionally lenticonus and retinopathy. Approximately 80% of Alport syndrome cases are caused by X-linked mutations in the COL4A5 gene encoding type IV collagen. The objective of this study was to define the SNP profiles for COL4A5 in patients with hereditary nephritis and hematuria. For this, we examined four subjects from one Kazakh family clinically affected with X-linked Alport syndrome due to COL4A5 gene mutations. All 51 exons of the COL4A5 gene were screened by linkage analysis and direct DNA sequencing, resulting in the identification of a novel mutation (G641E) in exon 25. The mutation was found only in two affected family individuals but was not present in healthy family members or 200 unrelated healthy controls. This result demonstrates that this novel mutation is pathogenic and has meaningful implications for the diagnosis of patients with Alport syndrome.

Behçet's disease in HLA-B*51 negative Germans and Turks shows association with HLA-Bw4-80I.

INTRODUCTION: Behçet's disease (BD) as systemic vasculitis of unknown etiology is associated with HLA-B*51 in European and Asian populations. HLA-A*26 was claimed as an additional BD susceptibility marker in Japanese and Greek patients. This study was performed to test for HLA associations in HLA-B*51 negative German and Turkish BD populations. METHODS: In total, 65 German and 46 Turkish patients lacking HLA-B*51 were analyzed in comparison to healthy HLA-B*51 negative Germans (n = 1500) and Turks (n = 130). HLA-A/B genotypes were determined by SSOP. P-values with correction for multiple testing (p(c)), χ2-test and odds ratio (OR) were used for statistical evaluation. RESULTS: HLA-A*26 was significantly more frequent in HLA-B*51- German patients [p(c) = 0.0076, OR = 3.23, 95% CI 1.63 to 6.39] than in respective controls. HLA-A*26 was also elevated in a smaller group of Turkish patients versus the controls. Significant association of HLA-Bw4 with isoleucine at amino-acid position 80 (HLA-Bw4-80I) was found in the HLA-B*51(-) German cohort of BD patients [p(c) = 0.0042, OR = 2.35, 95% CI 1.41 to 3.93) and in the Turkish patients in comparison to the respective controls [p = 0.025, OR = 2.17, 95% CI 1.09 to 4.31]. On the contrary, HLA-Bw4-80 T was reduced in both HLA-B*51(-) BD patient cohorts. CONCLUSIONS: The study shows a significant association of HLA-Bw4-80I present on HLA-B*51 as well as on other B-locus molecules with BD. This indicates that distinctive Bw4 epitopes on HLA-B locus molecules could play a role in BD pathogenesis. The study also indicates an association with HLA-A*26 in German and Turkish BD patients as a genetic risk factor independent of HLA-B*51.

Polymorphisms of HLA-DRB1, -DQA1 and -DQB1 in inhabitants of Astana, the capital city of Kazakhstan.

BACKGROUND: Kazakhstan has been inhabited by different populations, such as the Kazakh, Kyrgyz, Uzbek and others. Here we investigate allelic and haplotypic polymorphisms of human leukocyte antigen (HLA) genes at DRB1, DQA1 and DQB1 loci in the Kazakh ethnic group, and their genetic relationship between world populations. METHODOLOGY/PRINCIPAL FINDINGS: A total of 157 unrelated Kazakh ethnic individuals from Astana were genotyped using sequence based typing (SBT-Method) for HLA-DRB1, -DQA1 and -DQB1 loci. Allele frequencies, neighbor-joining method, and multidimensional scaling analysis have been obtained for comparison with other world populations. Statistical analyses were performed using Arlequin v3.11. Applying the software PAST v. 2.17 the resulting genetic distance matrix was used for a multidimensional scaling analysis (MDS). Respectively 37, 17 and 19 alleles were observed at HLA-DRB1, -DQA1 and -DQB1 loci. The most frequent alleles were HLA-DRB1*07:01 (13.1%), HLA-DQA1*03:01 (13.1%) and HLA-DQB1*03:01 (17.6%). In the observed group of Kazakhs DRB1*07:01-DQA1*02:01-DQB1*02:01 (8.0%) was the most common three loci haplotype. DRB1*10:01-DQB1*05:01 showed the strongest linkage disequilibrium. The Kazakh population shows genetic kinship with the Kazakhs from China, Uyghurs, Mongolians, Todzhinians, Tuvinians and as well as with other Siberians and Asians. CONCLUSIONS/SIGNIFICANCE: The HLA-DRB1, -DQA1 and -DQB1 loci are highly polymorphic in the Kazakh population, and this population has the closest relationship with other Asian and Siberian populations.

The need for standardized biobanks in Kazakhstan.

Biobanks are an important tool for clinical and research studies conducted on biomarkers of genetic therapy, diagnostic tests and new drugs; however, most biobanks remain incomplete and are often used without uniform standards and criteria. There is also a a lack of high-quality biological samples and many bioethical problems are often overlooked. Currently, Kazakhstan has no standard requirements and protocols for biomedical organizations. However, .an analysis of published data shows that possibly hundreds of samples are analyzed. Therefore, an establishment of biobank with standardized requirements could create better quality research. The National Center for Biotechnology has already started a biobank with more than 1,500 blood samples, with the ultimate goal of creating a biobank including around 10,000 blood samples of healthy volunteers, the same number of samples obtained from individuals with cardiovascular and endocrine diseases with samples stored under special conditions. The database contains demographic characteristics of donor's medical history. Informed consent for research received from all donors. This biobank can be considered as a national resource for scientific research.

Frequency of NAT2 and GSTP1 polymorphisms in the Kazakh population.

INTRODUCTION: Phase II xenobiotic biotransformation enzymes perform detoxification of hydrophilic and often toxic Phase I products through glutathionetransferase (GST), UDP-glucuronosyltransferase (UDF), N-acetyltransferase (NAT) families and other enzymes. GST protein families metabolize a large number of electrophilic xenobiotics, by conjugating fusing them with glutathione. Arylamine-N-acetyltransferase (NAT) catalyzes the acetylation of the aromatic and heterocyclic amines. MATERIALS AND METHODS: This study assesses the frequency of NAT2 and GSTP1 gene polymorphisms in 326 healthy individuals from different regions of Kazakhstan by using Real-Time PCR and direct sequencing methods. RESULTS: The allele frequencies were calculated for NAT2*5 (0.54) and GSTP1 (0.27). GSTP1 alleles were in the Hardy-Weinberg equilibrium (p > 0.05), while NAT2*5 (p = 0.00) were not. The population differences between North, Northeast and South Kazakhstan regions were also analyzed. No statistically significant differences in the frequency of genotypes were found. CONCLUSION: Allelic polymorphisms of NAT2*5 and GSTP1 genes greatly varied indifferent populations. The Kazakh population was significantly different from the Asian, Caucasoid, African-American and Hispanic populations by NAT2*5 and GSTP1 genes. Allelic variants of the NAT2*5 had a low frequency in Asian populations. Allelic frequency in other world populations varied from 30 to 50%. The differences between Kazakh (0.54) and the world population were statistically significant (p < 0.05). The frequency of GSTP1 (rs1695) in the African American population was 42%. The frequency of GSTP1 in Asian populations varied from 11% to 23%. The frequency in Caucasoid populations was around 30%. The differences between Kazakh population (0.27) and other populations selected were statistically significant (p < 0.05).The study of mutations in GSTP1 and NAT2 genes is necessary in assessing the risk of the development of various diseases, such as cancer. Information on allelic polymorphisms might also be useful for personal perscriptions such as cyclophosphamide, cisplatin, methotrexate, isoniazid, pyrazinamide, and rifampin.

Pharmacogenetic research in Kazakhstan.

INTRODUCTION: Pharmacogenomics is an emerging field of medicine that combines genetics and pharmacology. Pharmacogenomic research is relatively new in Kazahkstan, but, in recent years, significant progress has been made in this field. The National Scientific Laboratory for Biotechnology has launched several government-funded research projects focused on finding genetic markers that determine susceptibility to various drugs. Another goal of pharmacogenetic research in the laboratory is to find the pharmacogenomic markers that target cardiovascular diseases, accounting for allelic frequencies in selected genes in the Kazakh population. In addition, pharmacogenomic testing kits allow patients to choose the drug dosage. For example, the drug Warfarin has been developed within the framework of the "Technology Commercialization Project," funded jointly by the Ministry of Education and Science of the Republic of Kazakhstan and the World Bank. MATERIAL AND METHODS: The pharmacogenomic studies were conducted using the real-time PCR and direct DNA sequencing. DNA was isolated from venous blood or buccal cells, collected from patients. RESULTS: To date, we have identified the most promising areas of research in the field of pharmacogenomics in Kazakhstan. The allelic frequencies of a number of polymorphisms in the Kazakh population have been calculated (CYP2C9, CYP2C19, CYP3A4, VKORC1, CYP4F2, GGCX, CYP2D6, CYP1A2, NAT2, GSTP1, SLC47A1). A unique repository of DNA samples was established and is being replenished during the implementation of aforementioned projects. Development of the testing kit for individual selection of Warfarin dosage is nearing completion. A patent, named "Method of Selection Based Dose Warfarin Genotyping for the Kazakh Population" has been recently obtained. An application for another patent, titled "Express Method of Correction of Warfarin Dosing, Based on Real-time PCR" has received positive evaluation. The results of domestic pharmacogenomic studies will allow a more rational selection of drugs and their dosage regimens specific to the Kazakh population.

Draft Genome Sequence of Rhodococcus erythropolis DN1, a Crude Oil Biodegrader.

We report the 6,548-Mb genome sequence of Rhodococcus erythropolis strain DN1, isolated from the oil-contaminated soil in the Karagandy region of Kazakhstan. The draft genome sequence of strain DN1 might provide new insights into the genetic mechanisms of crude oil biodegradation.

Complete Genome Sequences of Helicobacter pylori Rifampin-Resistant Strains.

Here we present the complete genome sequences of two Helicobacter pylori rifampin-resistant (Rif(r)) strains (Rif1 and Rif2). Rif(r) strains were obtained by in vitro selection of H. pylori 26695 on agar plates with 20 µg/ml rifampin. The genome data provide insights on the genomic diversity of H. pylori under selection by rifampin.

Estrogen Receptor Gene (ESR1) PVUII and XBAI Polymorphisms and Bone Mineral Density in Kazakh Women.

INTRODUCTION: Osteoporosis is a common age-related disease that is strongly influenced by genetics. Polymorphisms of the estrogen receptor gene alpha (ESR1) are consistently been associated with bone mineral density (BMD) and fracture.The purpose of this investigation was to evaluate potential association of single nucleotide polymorphism (SNP) variants of the ESR1 gene and bone mineral density (BMD) of the lumbar spine in Kazakh women. METHODS: 140 female participants in Pavlodar clinics with varying measures of BMD. We are examined the potential association of BMD with 2 SNPs from the ESR1 gene (rs2234693 [PvuII] and rs9340799 [XbaI]). Genotyping of the PvuII and XbaI polymorphisms was performed by direct sequencing of the gene fragments containing restriction sites with the identification of genotypes PP, Pp, pp and XX, Xx, xx respectively. RESULTS: Unadjusted mean BMD values ranged from 1.14±0.14 g/cm2 in Caucasian women and 1.03±0.11 g/cm2 in Asian women. The association between PvuII polymorphism and BMD at the lumbar spine (p= 0.04 for PP=Pp=pp) was statistically significant in all women. The XbaI polymorphism was not associated with BMD at lumbar spine. The relative risk for low BMD was higher for the marker PvuII (RR=1.51) than for the marker XbaI (RR=1.35). CONCLUSION: The PvuII polymorphism had a weak association with lumbar spine BMD. XbaI polymorphism was unlikely to be a predictor of lumbar spine BMD in Kazakh women. These conclusions could help to determine the genetic risk factors for osteoporosis; however, further studies on the association between gene polymorphisms and BMD are needed including larger numbers of participants and genes to clarify genetic risks.

Genetic profile and determinants of homocysteine levels in Kazakhstan patients with breast cancer.

AIM: To analyze associations between homocysteine level, MTHFR and FTO rs1477196 polymorphisms and folate status in patients with breast cancer (BC) in order to clarify determinants of hyperhomocysteinemia. PATIENTS AND METHODS: The study included 315 BC cases and 604 controls. RESULTS: The MTHFRC677T genotype was associated with an increased incidence of BC [Odds ratio (OR)=1.71; 95% Confidential interval (CI)=1.21-2.43]. The MTHFR A1298C genotype was associated with a decreased risk of BC [OR=0.68; 95% CI: 0.49-0.95]. The homocysteine level was not associated with either MTHFR C677T or A1298C, nor with FTO rs1477196, but was inversely correlated with folate status in cancer cases (p<0.0001) and tended to be higher in patients with the MTHFR 677TT genotype. The folate level (p<0.0005) was an independent predictor of hyper-homocysteinemia in patients with BC. CONCLUSION: These results suggest an important role of homocysteine in breast tumorigenesis. Further studies are warranted to investigate how combined MTHFR genotypes exert their effects on cancer susceptibility.

Vitamin D Receptor Gene Polymorphisms and Breast Cancer Risk in Kazakhstan.

INTRODUCTION: The steroid hormone 1,25-dihydroxyvitamin D3 is thought to protect against breast cancer. The activity of 1,25-dihydroxyvitamin D3 is mediated via the vitamin D receptor (VDR), and a number of polymorphisms in the VDR gene have been identified. These result in distinct genotypes, some of which may alter susceptibility to breast cancer. Two common single nucleotide polymorphisms (SNP) in the VDR gene (VDR), rs1544410 (BsmI) and rs2228570 (FokI), have been inconsistently associated with breast cancer risk. Increased risk has been reported for the FokI ff genotype, which encodes a less transcriptionally active isoform of VDR. A reduced risk has been reported for the BsmI BB genotype which may influence VDR mRNA stability. AIM: We have investigated whether specific VDR gene polymorphisms are associated with breast cancer risk in Kazakhstan women. MATERIAL AND METHODS: In a case-control study, female breast cancer patients (315) and a female control group (n=604) were tested for two VDR polymorphisms. Statistical analysis was conducted using SPSS19.0. RESULTS: : The VDR rs2228570 (FokI) polymorphism was associated with an increased occurence of BC [rs2228570 (folk) ff vs. FF genotype: OR=1.71; 95% CI=1.21-2.43]. No association was noted between rs1544410 (BsmI) BB and breast cancer risk [OR=0.68; 95% CI=0.49-0.95]. CONCLUSION: : Although the factors that increase breast cancer susceptibility remain uncertain, future large studies should integrate genetic variation in VDR with biomarkers of vitamin D status. Additional testing on the effect of varying genotypes on the functional mechanisms of the VDR could help to improve future testing and treatment of woman at risk for breast cancer.

Draft Genome Sequence of the Live Vaccine Strain Brucella abortus 82.

Vaccination is a crucial part of the brucellosis eradication programs worldwide. A live vaccine strain of Brucella abortus 82 has been successfully used for the vaccination of cattle against brucellosis in the former Soviet republics for the last 39 years. Here, we report the genome sequence of Brucella abortus 82.

BRCA1 and BRCA2 Gene Mutations Screening In Sporadic Breast Cancer Patients In Kazakhstan.

BACKGROUND: A large number of distinct mutations in the BRCA1 and BRCA2 genes have been reported worldwide, but little is known regarding the role of these inherited susceptibility genes in breast cancer risk among Kazakhstan women. AIM: To evaluate the role of BRCA1/2 mutations in Kazakhstan women presenting with sporadic breast cancer. METHODS: We investigated the distribution and nature of polymorphisms in BRCA1 and BRCA2 entire coding regions in 156 Kazakhstan sporadic breast cancer cases and 112 age-matched controls using automatic direct sequencing. RESULTS: We identified 22 distinct variants, including 16 missense mutations and 6 polymorphisms in BRCA1/2 genes. In BRCA1, 9 missense mutations and 3 synonymous polymorphisms were observed. In BRCA2, 7 missense mutations and 3 polymorphisms were detected. There was a higher prevalence of observed mutations in Caucasian breast cancer cases compared to Asian cases (p<0.05); higher frequencies of sequence variants were observed in Asian controls. No recurrent or founder mutations were observed in BRCA1/2 genes. There were no statistically significant differences in age at diagnosis, tumor histology, size of tumor, and lymph node involvement between women with breast cancer with or without the BRCA sequence alterations. CONCLUSIONS: Considering the majority of breast cancer cases are sporadic, the present study will be helpful in the evaluation of the need for the genetic screening of BRCA1/2 mutations and reliable genetic counseling for Kazakhstan sporadic breast cancer patients. Evaluation of common polymorphisms and mutations and breast cancer risk in families with genetic predisposition to breast cancer is ongoing in another current investigation.