Background: Copy Number Variation (CNV) of the human CNTN6 gene (encoding the contactin-6 protein), caused by deletions or duplications, is responsible for severe neurodevelopmental impairments, often in combination with facial dysmorphias. Conversely, deleterious point mutations of this gene do not show any clinical phenotypes. The aim of this study is to generate mice carrying large deletions, duplications and inversions involving the Cntn6 gene as a new experimental model to study CNV of the human CNTN6 locus.

Results: To generate large chromosomal rearrangements on mouse chromosome 6, we applied CRISPR/Cas9 technology in zygotes. Two guide RNAs (gRNAs) (flanking a DNA fragment of 1137 Mb) together with Cas9 mRNA and single-stranded DNA oligonucleotides (ssODN) were microinjected into the cytoplasm of 599 zygotes of F1 (C57BL x CBA) mice, and 256 of them were transplanted into oviducts of CD-1 females. As a result, we observed the birth of 41 viable F0 offspring. Genotyping of these mice was performed by PCR analysis and sequencing of PCR products. Among the 41 F0 offspring, we identified seven mice with deletions, two animals carrying duplications of the gene and four carrying inversions. Interestingly, two F0 offspring had both deletions and duplications. It is important to note that while three of seven deletion carriers showed expected sequences at the new joint sites, in another three, we identified an absence of 1-10 nucleotides at the CRISPR/Cas9 cut sites, and in one animal, 103 bp were missing, presumably due to error-prone non-homologous end joining. In addition, we detected the absence of 5 and 13 nucleotides at these sites in two F0 duplication carriers. Similar sequence changes at CRISPR/Cas9 cut sites were observed at the right and left boundaries of inversions. Thus, megabase-scale deletions, duplications and inversions were identified in 11 F0 offspring among 41 analyzed, i.e., approximately 25% efficiency. All genetically modified F0 offspring were viable and able to transmit these large chromosomal rearrangements to the next generation.

Conclusions: Using CRISPR/Cas9 technology, we created mice carrying megabase-scale deletions, duplications, and inversions involving the full-sized Cntn6 gene. These mice became founders of new mouse lines, which may be more appropriate experimental models of CNV in the human 3p26.3 region than Сntn6 knockout mice.

Epoxyeicosatrienoic acids (EETs) are important vasoactive products of arachidonic acid metabolism with a wide range of biological actions in the cardiovascular system. The present study investigated whether single nucleotide polymorphisms (SNP) of genes coding cytochrome P450 2C subfamily, enzymes involved in biosynthesis of EETs, are associated with the risk of coronary heart disease (CHD). A total of 1255 unrelated Russian subjects comprising 561 patients with angiographically diagnosed CHD and 694 age- and sex-matched healthy subjects were included in the study. DNA samples from all study participants were genotyped for six common SNPs rs7909236, rs1934953 of CYP2C8, rs9332242, rs4918758 and rs61886769 of CYP2C9 and rs4244285 of CYP2C19 using by the Mass-ARRAY 4 system. SNP rs4918758 of CYP2C9 was associated with decreased risk of CHD (codominant model) at a borderline significance with odds ratio adjusted for sex and age 0.61 (95% CI: 0.41-0.92, P=0.038, Q=0.20). SNP rs9332242 of CYP2C9 showed a trend towards association with increased CHD risk in cigarette smokers (P=0.049, Q=0.29). Log-likelihood ratio test (LRT) pointed out epistatic interactions between rs9332242 and rs61886769 of CYP2C9 (codominant model, Pinteraction=0.02), however, this P-value did not survive after correction for multiple tests. Bioinformatic analysis revealed a regulatory potential for a majority of the investigated SNPs. Our preliminary results demonstrate that polymorphisms of genes encoding CYP2C subfamily represent potential genetic markers of CHD susceptibility. Further studies are required to substantiate the contribution of these genes to the disease risk.

Human mitochondrial DNA haplogroup U is among the initial maternal founders in Southwest Asia and Europe and one that best indicates matrilineal genetic continuity between late Pleistocene hunter-gatherer groups and present-day populations of Europe. While most haplogroup U subclades are older than 30 thousand years, the comparatively recent coalescence time of the extant variation of haplogroup U7 (~16-19 thousand years ago) suggests that its current distribution is the consequence of more recent dispersal events, despite its wide geographical range across Europe, the Near East and South Asia. Here we report 267 new U7 mitogenomes that - analysed alongside 100 published ones - enable us to discern at least two distinct temporal phases of dispersal, both of which most likely emanated from the Near East. The earlier one began prior to the Holocene (~11.5 thousand years ago) towards South Asia, while the later dispersal took place more recently towards Mediterranean Europe during the Neolithic (~8 thousand years ago). These findings imply that the carriers of haplogroup U7 spread to South Asia and Europe before the suggested Bronze Age expansion of Indo-European languages from the Pontic-Caspian Steppe region.

Numerous studies demonstrated an importance of cytochrome P-450 epoxygenase pathway of arachidonic acids metabolism for the pathogenesis of essential hypertension (EH). The present study was designed to investigate whether common single-nucleotide polymorphisms (SNP) of CYP2C gene subfamily such as CYP2C8 (rs7909236 and rs1934953), CYP2C9 (rs9332242), and CYP2C19 (rs4244285) are associated with susceptibility to EH in Russian population. A total of 816 unrelated Russian individuals comprising 425 EH patients and 391 normotensive controls were included into the study. Genotyping of SNPs was performed using the MassARRAY 4 system. SNP rs7909236 of CYP2C8 was significantly associated with increased risk of EH (OR adjusted for sex and age was 2.99 95% CI 1.39-6.44, P = 0.005). SNPs rs1934953 CYP2C8 and rs4244285 of CYP2C19 showed association with EH risk but at a borderline statistical significance (P ≤ 0.04). Combination of genotypes CYP2C8 rs7909236 TT and CYP2C19 rs4244285 GG was associated with increased EH risk (OR 3.34 95%CI 1.48-7.51, P = 0.004). Genotype-phenotype correlation analysis showed that the levels of CYP2C8 mRNA were significantly correlated with SNP rs7909236 (P = 0.01). in silico functional prediction analysis revealed the functionality of majority of investigated SNPs. Thus, genes of CYP2C subfamily are important genetic determinants of susceptibility to essential hypertension in Russians.

Background: The colonization of Eurasia and Australasia by African modern humans has been explained, nearly unanimously, as the result of a quick southern coastal dispersal route through the Arabian Peninsula, the Indian subcontinent, and the Indochinese Peninsula, to reach Australia around 50 kya. The phylogeny and phylogeography of the major mitochondrial DNA Eurasian haplogroups M and N have played the main role in giving molecular genetics support to that scenario. However, using the same molecular tools, a northern route across central Asia has been invoked as an alternative that is more conciliatory with the fossil record of East Asia. Here, we assess as the Eurasian macrohaplogroup R fits in the northern path.

Results: Haplogroup U, with a founder age around 50 kya, is one of the oldest clades of macrohaplogroup R in western Asia. The main branches of U expanded in successive waves across West, Central and South Asia before the Last Glacial Maximum. All these dispersions had rather overlapping ranges. Some of them, as those of U6 and U3, reached North Africa. At the other end of Asia, in Wallacea, another branch of macrohaplogroup R, haplogroup P, also independently expanded in the area around 52 kya, in this case as isolated bursts geographically well structured, with autochthonous branches in Australia, New Guinea, and the Philippines.

Conclusions: Coeval independently dispersals around 50 kya of the West Asia haplogroup U and the Wallacea haplogroup P, points to a halfway core area in southeast Asia as the most probable centre of expansion of macrohaplogroup R, what fits in the phylogeographic pattern of its ancestor, macrohaplogroup N, for which a northern route and a southeast Asian origin has been already proposed.

Purpose: We aimed to investigate the association between the polymorphism and expression patterns of multiple drug resistance genes (MDR) in breast cancer (BC).

Materials and methods: The MDR gene expression levels were measured in tumor tissues of 106 breast cancer patients using quantitative real-time PCR. Affymetrix CytoScan™ HD Array chips were used to assess genotypes. Pairwise correlation analysis for ABCB1, ABCC1, ABCC2 and ABCG2 gene expression levels was carried out to reveal co-expression clusters. Associations between SNPs of MDR genes and their preoperative expression levels were assessed using analysis of covariance adjusting for covariates.

Results: The SNPs associated with the expression of the ABCB1, ABCC1, ABCC2 and ABCG2 genes before NAC were detected. In addition, 21 SNPs associated with the expression of four ABC-transporter genes and involved in the expression regulation were identified. Validation in an independent sample confirmed the association between the MDR cluster genes and 11 SNPs.

Conclusions: Four MDR genes: ABCB1, ABCC1, ABCC2 and ABCG2 were shown to form the functional expression cluster in breast tumor. Further studies are required to discover precise mechanisms of the cluster regulation, thereby providing new approaches and targets to combat the development of the MDR phenotype during chemotherapy.

L-ficolin encoded by FCN2 gene is a crucial factor of defence against infection in humans. We studied the prevalence of the two common variants (rs17549193 and rs7851696) in aboriginal and alien populations of the Taymyr-Dolgan-Nenets region of Krasnoyarskiy Kray, East Siberia, Russia (Nenets, Dolgans, Nganasans, Russians). We found a decreased prevalence of the rs17549193*T allele in all aboriginal populations as compared to Russians. Also, its frequency was the lowest in the Nenets among the studied populations, while frequency of the rs7851696*T allele was increased in this population. The results suggest that the Arctic populations of East Siberia are characterised by specificity of genetic make-up responsible for the activity of L-ficolin. Clinical and epidemiological studies are required to discover if these genetic features correlate with the infant infectious morbidity in East Siberian populations.

Background: Antipsychotic drugs can cause hyperprolactinemia. However, hyperprolactinemia was also observed in treatment-naive patients with a first schizophrenic episode. This phenomenon might be related to the role of prolactin as a cytokine in autoimmune diseases. Extrapituitary prolactin production is regulated by an alternative promoter, which contains the functional single nucleotide polymorphism -1149 G/T (rs1341239). We examined whether this polymorphism was associated with hyperprolactinemia in patients with schizophrenia.

Method: We recruited 443 patients with schizophrenia and 126 healthy controls. The functional polymorphism -1149 G/T (rs1341239) in the prolactin gene was genotyped with multiplexed primer extension, combined with MALDI-TOF mass spectrometry. Genotype and allele frequencies were compared between groups with the χ2 test and logistic regression models adjusting for covariates.

Results: The frequency of genotypes and alleles in patients with schizophrenia did not differ from those in control subjects. A comparison between patients with schizophrenia with and without hyperprolactinemia revealed significantly higher frequency of the G allele in patients with hyperprolactinemia than in patients without it (χ2=7.25; p=0.007; OR=1.44 [1.10-1.89]). Accordingly, patients with hyperprolactinemia carried the GG genotype more frequently than patients without hyperprolactinemia (χ2=9.49; p=0.009). This association remained significant after adjusting the estimates for such covariates as sex, age, duration of the diseases and the dose of chlorpromazine equivalents.

Conclusion: This study revealed a significant association between the polymorphic variant rs1341239 and the development of hyperprolactinemia in patients with schizophrenia. The serum prolactin concentration in patients with schizophrenia treated with antipsychotics may provide an indication of the activity of the gene that regulates extrapituitary prolactin production which is believed to play a role in the immune system.

Hyperprolactinaemia (HPRL) is a classical side effect of antipsychotic drugs primarily attributed to blockade of dopamine D2 subtype receptors in the pituitary gland. Although dopamine is considered the primary factor inhibiting prolactin release, the activity of prolactin-producing lactotrophs is also regulated by the secretagogues thyrotrophin releasing hormone, vasoactive intestinal polypeptide and serotonin (5-hydroxytryptamine; 5-HT).

Methods: We describe the association between HPRL and a set of 29 SNPs from 5-HT receptor genes HTR1A, HTR1B, HTR2A, HTR2C, HTR3A, HTR3B and HTR6 in a population of 446 Caucasians (221 males/225 females) with a clinical diagnosis of schizophrenia (according to ICD-10: F20) who were treated with classical and/or atypical antipsychotic drugs.

Results: None of the studied autosomal markers were found to be associated with HPRL. However, a significant association was established between various HTR2C polymorphisms and HPRL.

Conclusions: This study revealed an association between HPRL and X-chromosome haplotypes comprised of the rs569959 and rs17326429 polymorphisms.

Objective: The present study was designed to investigate whether genetic polymorphisms of the aryl hydrocarbon receptor (AHR) signaling pathway are involved in the molecular basis of essential hypertension (EH).

Methods: A total of 2160 unrelated Russian individuals comprising 1341 EH patients and 819 healthy controls were recruited into the study. Seven common AHR pathway single-nucleotide polymorphisms (SNPs) such as rs2066853, rs2292596, rs2228099, rs1048943, rs762551, rs1056836, and rs1800566 were genotyped by TaqMan-based allele discrimination assays.

Results: We found that SNP rs2228099 of ARNT is associated with an increased risk of EH (odds ratio=1.20 95% confidence interval: 1.01-1.44, P=0.043) in a dominant genetic model, whereas polymorphism rs762551 of CYP1A2 showed an association with a decreased risk of disease in a recessive genetic model (odds ratio=0.68, 95% confidence interval: 0.52-0.89, P=0.006). A log-likelihood ratio test enabled identification of epistatic interaction effects on EH susceptibility for all SNPs. MB-MDR analysis showed that cigarette smoking, rs1048943, rs762551, rs1056836, and rs2228099 were significant contributing factors in 19, 18, 13, 13, and 11 interaction models, respectively. The best MDR model associated with EH risk included rs1048943, rs762551, rs1056836, and cigarette smoking (cross-validation consistency 100%, prediction error 45.7%, Ppermutation<0.0001). The mRNA expression and in-silico function prediction analyses have confirmed a regulatory potential for a majority of SNPs associated with EH susceptibility.

Conclusion: Our pilot study was the first to show that gene-gene and gene-environment interactions in the AHR signaling pathway represent important determinants for the development of EH, and the pathway may become an attractive target for a pharmacological intervention in hypertensive patients in the future.

Intratumor morphological heterogeneity in breast cancer is represented by different morphological structures (tubular, alveolar, solid, trabecular, and discrete) and contributes to poor prognosis; however, the mechanisms involved remain unclear. In this study, we performed 3D imaging, laser microdissection-assisted array comparative genomic hybridization and gene expression microarray analysis of different morphological structures and examined their association with the standard immunohistochemistry scorings and CD44+CD24- cancer stem cells. We found that the intratumor morphological heterogeneity is not associated with chromosomal aberrations. By contrast, morphological structures were characterized by specific gene expression profiles and signaling pathways and significantly differed in progesterone receptor and Ki-67 expression. Most importantly, we observed significant differences between structures in the number of expressed genes of the epithelial and mesenchymal phenotypes and the association with cancer invasion pathways. Tubular (tube-shaped) and alveolar (spheroid-shaped) structures were transcriptionally similar and demonstrated co-expression of epithelial and mesenchymal markers. Solid (large shapeless) structures retained epithelial features but demonstrated an increase in mesenchymal traits and collective cell migration hallmarks. Mesenchymal genes and cancer invasion pathways, as well as Ki-67 expression, were enriched in trabecular (one/two rows of tumor cells) and discrete groups (single cells and/or arrangements of 2-5 cells). Surprisingly, the number of CD44+CD24- cells was found to be the lowest in discrete groups and the highest in alveolar and solid structures. Overall, our findings indicate the association of intratumor morphological heterogeneity in breast cancer with the epithelial-mesenchymal transition and CD44+CD24- stemness and the appeal of this heterogeneity as a model for the study of cancer invasion.

The objective of this study was to identify genes targeted by both copy number and copy-neutral changes in the right coronary arteries in the area of advanced atherosclerotic plaques and intact internal mammary arteries derived from the same individuals with comorbid coronary artery disease and metabolic syndrome. The artery samples from 10 patients were screened for genomic imbalances using array comparative genomic hybridization. Ninety high-confidence, identical copy number variations (CNVs) were detected. We also identified eight copy-neutral changes (cn-LOHs) > 1.5 Mb in paired arterial samples in 4 of 10 individuals. The frequencies of the two gains located in the 10q24.31 (ERLIN1) and 12q24.11 (UNG, ACACB) genomic regions were evaluated in 33 paired arteries and blood samples. Two patients contained the gain in 10q24.31 (ERLIN1) and one patient contained the gain in 12q24.11 (UNG, ACACB) that affected only the blood DNA. An additional two patients harboured these CNVs in both the arteries and blood. In conclusion, we discovered and confirmed a gain of the 10q24.31 (ERLIN1) and 12q24.11 (UNG, ACACB) genomic regions in patients with coronary artery disease and metabolic comorbidity. Analysis of DNA extracted from blood indicated a possible somatic origin for these CNVs.

Background: Genome-wide association studies have identified novel genetic associations for asthma, but without taking into account the role of active tobacco smoking. This study aimed to identify novel genes that interact with ever active tobacco smoking in adult onset asthma.

Methods: We performed a genome-wide interaction analysis in six studies participating in the GABRIEL consortium following two meta-analyses approaches based on 1) the overall interaction effect and 2) the genetic effect in subjects with and without smoking exposure. We performed a discovery meta-analysis including 4,057 subjects of European descent and replicated our findings in an independent cohort (LifeLines Cohort Study), including 12,475 subjects.

Results: First approach: 50 SNPs were selected based on an overall interaction effect at p<10-4. The most pronounced interaction effect was observed for rs9969775 on chromosome 9 (discovery meta-analysis: ORint = 0.50, p = 7.63*10-5, replication: ORint = 0.65, p = 0.02). Second approach: 35 SNPs were selected based on the overall genetic effect in exposed subjects (p <10-4). The most pronounced genetic effect was observed for rs5011804 on chromosome 12 (discovery meta-analysis ORint = 1.50, p = 1.21*10-4; replication: ORint = 1.40, p = 0.03).

Conclusions: Using two genome-wide interaction approaches, we identified novel polymorphisms in non-annotated intergenic regions on chromosomes 9 and 12, that showed suggestive evidence for interaction with active tobacco smoking in the onset of adult asthma.

Tuberculosis is a global healthcare challenge. Host genetic factors were proven to modify risk of the disease. Genome-wide association studies revealed a number of loci associated with TB in different populations. However, no systematic analysis of genetic bases of susceptibility to different clinical stages, such as primary TB and reactivation, was carried out. We set out to validate the results of GWASs in Russians of West Siberia with a consideration of primary and secondary TB. We chose 45 SNP from five large GWASs and genotyped 445 healthy individuals and 323 TB patients including 74 with primary TB and 249 with reactivation. We found that the rs7821565 and rs40363 SNPs were associated with primary TB in Russians (p = 0.019 and 1.4e − 3, respectively), while rs10515787 and rs2837857 were associated with secondary TB (p = 1.2e − 3 and 0.039, respectively). The results suggest genetic basis of primary and secondary TB differs.

В исследование включены больные хроническим вирусным гепатитом С (ХВГС) (n = 184) и популяционная выборка жителей г.Томска (n = 285). Генотипирование 58 полиморфных вариантов выполнено методом масс-спектрометрии на приборе Sequenom MassARRAY® (США). Статистическая обработка данных проводилась в программной среде R с помощью стандартного пакета «stats» (версия 3.0.3) [http://www.R-project.org/]. Больные ХВГС характеризуются более высокой частотой аллелей «А» rs679620 гена MMP3 (OR = 1,74, 95%CI: 1,25-2,42; p = 0,0001) и «С» rs3739998 гена KIAA1462 (OR = 1,33, 95%CI: 1,01-1,75; р = 0,044) и более низкой частотой аллелей «Т» rs20579 гена LIG1 (OR = 0,55, 95%CI: 0,36-0,84; р = 0,004), «G» rs3765124 гена ADAMDEC1 (OR = 0,68, 95%CI: 0,51-0,84; р = 0,006) и «С» rs1007856 гена ITGB5 (OR = 0,75, 95%CI: 0,57-0,99; p = 0,045), по сравнению с контрольной группой. Предрасполагающими к развитию ХВГС являются генотипы «СС» rs10087305 (OR = 5,83, 95%CI: 1,74-19,02; р = 0,002) и «АА» rs3765124 (OR = 1,78, 95%CI: 1,16-8,46; р = 0,008) гена ADAMDEC1 ; «АА» rs679620 гена MMP3 (OR = 2,40, 95%CI: 1,40-4,12; р = 0,001); «TT» rs1007856 гена ITGB5 (OR = 1,74, 95%CI: 1,10-2,73; р = 0,015); «СС» rs20579 гена LIG1 (OR = 1,92, 95%CI: 1,18-3,12; р = 0,007); «СС» rs3739998 гена KIAA1462 (OR = 1,89, 95%CI: 1,16-3,10; р = 0,009). Таким образом, в предрасположенность ХВГС вносят вклад гены регуляции метаболизма экстрацеллюлярного матрикса ( ADAMDEC1 ; MMP3; ITGB5 ), непосредственно влияющие на процессы фиброгенеза, а также гены, ответственные за функционирование эндотелия ( KIAA1462 ) и репарацию ДНК ( LIG1 ).