A New Baltic Population-Specific Human Genetic Marker in the PMCA4 Gene.

OBJECTIVES: The PMCA gene family consists of 4 genes and at least 21 splice variants; among these, the Ca2+ ATPase 4 (PMCA4) gene encodes a plasma membrane protein abundantly expressed in several tissues, including the kidney, heart, and sperm. Knockout of PMCA4 causes infertility due to immotile sperm in mouse models. We therefore investigated variants in this gene for potential association with infertility in groups of Estonian (n = 191) and Latvian (n = 92) men with reduced sperm motility. METHODS: All exons, exon-intron boundaries, 5' and 3' untranslated regions, and the promoter region of the PMCA4 gene were analysed by direct sequencing for a group of Estonian infertile men. Genotyping of guanine and adenine alleles of rs147729934 was performed, using a custom-designed TaqMan® probe for a group of Latvian infertile men as well as additional groups from Latvia and several groups of people with proven ethnicity from the Baltic region. RESULTS: Although we did not identify any significant associations between variants in the gene and infertility, our results indicated that in all studied Latvian and Estonian groups the adenine allele of the variant rs147729934 was present at a higher frequency than expected. Analysis of additional samples indicated that the adenine allele of rs147729934 likely originated once in the modern-day Baltic or western Russia area, as the frequency of the minor adenine allele observed in this region is remarkably higher than that in the general European population. CONCLUSIONS: Our results revealed no significant difference in frequencies of genetic variants in PMCA4 gene between men with normal and those with reduced sperm motility. The adenine allele of the variant rs147729934 is potentially an informative tool for future population studies concerning ancient Baltic and Finno-Ugric history.

Robust genotyping tool for autosomal recessive type of limb-girdle muscular dystrophies.

BACKGROUND: Limb-girdle muscular dystrophies are characterized by predominant involvement of the shoulder and pelvic girdle and trunk muscle groups. Currently, there are 31 genes implicated in the different forms of limb-girdle muscular dystrophies, which exhibit similar phenotypes and clinical overlap; therefore, advanced molecular techniques are required to achieve differential diagnosis. METHODS: We investigated 26 patients from Latvia and 34 patients from Lithuania with clinical symptoms of limb-girdle muscular dystrophies, along with 565 healthy unrelated controls from general and ethnic populations using our developed test kit based on the Illumina VeraCode GoldenGate genotyping platform, Ion AmpliSeq Inherited Disease Panel and direct sequencing of mutations in calpain 3 (CAPN3), anoctamin 5 (ANO5) and fukutin related protein (FKRP) genes. RESULTS: Analysis revealed a homozygous CAPN3 c.550delA mutation in eight patients and three heterozygous variants in controls: dysferlin (DYSF) c.5028delG, CAPN3 c.2288A > G, and FKRP c.135C > T. Additionally, three mutations within FKRP gene were found: homozygous c.826C > A, and two compound - c.826C > A/c.404_405insT and c.826C > A/c.204_206delCTC mutations, and one mutation within CLCN1 gene - c.2680C > T p.Arg894Ter. ANO5 c.191dupA was not present. CONCLUSIONS: Genetic diagnosis was possible in 12 of 60 patients (20%). The allele frequency of CAPN3 gene mutation c.550delA in Latvia is 0.0016 and in Lithuania - 0.0029. The allele frequencies of CAPN3 gene mutation c.2288A > G and DYSF gene mutation c.4872delG are 0.003.

Y-Chromosomal Lineages of Latvians in the Context of the Genetic Variation of the Eastern-Baltic Region.

Variations of the nonrecombining Y-chromosomal region were investigated in 159 unrelated Baltic-speaking ethnic Latvians from four different geographic regions, using 28 biallelic markers and 12 short tandem repeats. Eleven different haplogroups (hgs) were detected in a regionally homogeneous Latvian population, among which N1c, R1a, and I1 cover more than 85% of its paternal lineages. When compared its closest geographic neighbors, the composition of the Latvian Y-chromosomal gene pool was found to be very similar to those of Lithuanians and Estonians. Despite the comparable frequency distribution of hg N1c in Latvians and Lithuanians with the Finno-Ugric-speaking populations from the Eastern coast of the Baltic Sea, the observed differences in allelic variances of N1c haplotypes between these two groups are in concordance with the previously stated hypothesis of different dispersal ways of this lineage in the region. More than a third of Latvian paternal lineages belong specifically to a recently defined R1a-M558 hg, indicating an influence from a common source within Eastern Slavic populations on the formation of the present-day Latvian Y-chromosome gene pool.

BCL3 gene role in facial morphology.

BACKGROUND: Cleft lip (CL) with or without palate (CLP) and isolated cleft palate (CP) are etiologically complex diseases with interactions among various environmental and genetic factors. The aim of the current study was to identify association with genetic markers and phenotypic craniofacial data in patients with CL/CLP/CP parents. METHODS: Posteroanterior and lateral digital radiographs of the cranium were obtained from 74 parents of patients with CL/CLP/CP. One hundred seventy-three patients with CL/CLP/CP and 190 controls were enrolled in the study for the association test. Five genetic markers of the IRF6 gene and 14 markers of the 19q13 locus were genotyped. Linear regression analysis was performed for the relationship of cephalometric measurements with genotype data adjusted for age, gender, and cleft type. Chi-square and transmission disequilibrium tests were performed to evaluate differences in alleles of the BCL3 gene. Positive findings were replicated in an independent sample (n = 95) of patients with CL/CLP/CP parents. RESULTS: Genetic markers of the BCL3 gene at 19q13, rs7257231, and rs1979377 in the familial association test and rs10401176 in the case-control association test, were associated with craniofacial phenotype. Carriers of BCL3 allele rs7257231T had longer posterior cranial bases than noncarriers (p(adjusted) = 0.0028), and in the familial-based association test showed the statistically strongest relationship (p(adjusted) = 0.05) to phenotype. Relation of rs7257231 to facial formation was confirmed in the replication group (p = 0.0024). CONCLUSIONS: The results indicate that BCL3, which has functions related to cell adhesion and whose downregulation can cause disruption of ectodermal development, is likely to be important in facial formation. Birth Defects Research (Part A), 2012.

Variation in FGF1, FOXE1, and TIMP2 genes is associated with nonsyndromic cleft lip with or without cleft palate.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (CL/P) is a common complex birth defect caused by the interaction between multiple genes and environmental factors. METHODS: Five hundred and eighty-seven single nucleotide polymorphisms in 40 candidate genes related to orofacial clefting were tested for association with CL/P in a clefting sample composed of 300 patients and 606 controls from Estonian, Latvian, and Lithuanian populations. RESULTS: In case-control comparisons, the minor alleles of FGF1 rs34010 (p = 4.56 × 10(-4) ), WNT9B rs4968282 (p = 0.0013), and FOXE1 rs7860144 (p = 0.0021) were associated with a decreased risk of CL/P. Multiple haplotypes in FGF1, FOXE1, and TIMP2 and haplotypes in WNT9B, PVRL2, and LHX8 were associated with CL/P. The strongest association was found for protective haplotype rs250092/rs34010 GT in the FGF1 gene (p = 5.01 × 10(-4) ). The strongest epistatic interaction was observed between the COL2A1 and WNT3 genes. CONCLUSIONS: Our results provide for the first time evidence implicating FGF1 in the occurrence of CL/P, and support TIMP2 and WNT9B as novel loci predisposing to CL/P. We have also replicated recently reported significant associations between variants in or near FOXE1 and CL/P. It is likely that variation in FOXE1, TIMP2, and the FGF and Wnt signaling pathway genes confers susceptibility to nonsyndromic CL/P in Northeastern European populations.

Association studies of candidate genes and cleft lip and palate taking into consideration geographical origin.

Isolated cleft lip and/or palate (CL/CLP) is a complex congenital anomaly with many contributing factors. There are several genes involved in the aetiology of CL/CLP, they are different in selected populations. In a previous study, the mitochondrial haplotypes of Latvian subjects with CL/CLP were characterized. Latvian subjects with CL/CLP have mostly mitochondrial haplogroups U4/U5 compared with the ethnic population of Latvia. The aim of this study was to stratify the results of genotyping based on European mitochondrial DNA (mtDNA) haplotypes. DNA samples from 108 patients with CL/CLP and from 182 unrelated and unaffected individuals selected randomly in Latvia (used as controls) were obtained for investigation. In this study, we analysed the data taking into consideration mitochondrial haplogroups and found that gene associations depended on the genetic origin of the population. The phenotype of patients with non-U haplotypes was associated with markers in wingless-type MMTV integration site family, member 3 (WNT3), collagen, type XI, alpha 2 (COL11A2), and fibroblast growth factor receptor 1 (FGFR1), whereas patients with U4 and U5 haplotypes showed significant association with WNT3 and COL11A2. It is unlikely that mtDNA variants play a direct role in the development of CL/CLP; rather, they may be a surrogate for population substructure and provide a tool to increase homogeneity and statistical power.

Genetic variants in COL2A1, COL11A2, and IRF6 contribute risk to nonsyndromic cleft palate.

BACKGROUND: Orofacial clefts are among the most common birth defects with a strong genetic component. Nonsyndromic cleft palate (NSCP) is a complex malformation determined by the interaction between multiple genes and environmental risk factors. METHODS: We conducted a case-control association study to investigate the role of 40 candidate genes in predisposition to orofacial clefting. Five hundred ninety-one haplotype tagging single nucleotide polymorphism (tagSNPs) were genotyped in a clefting sample from the Baltic region, composed of 104 patients with nonsyndromic cleft palate and 606 controls from an Estonian, Latvian, and Lithuanian population. RESULTS: In case-control comparisons, the minor alleles of IRF6 rs17389541 (p = 5.45 × 10(-4)) and COL2A1 rs1793949 (p = 7.26 × 10(-4)) were associated with increased risk of NSCP. Multiple haplotypes in COL2A1 and COL11A2 and haplotypes in WNT3, FGFR1, and CLPTM1were associated with NSCP. The strongest associations were found for IRF6 haplotype rs17389541/rs9430018 GT (p = 2.23 × 10(-4)) and COL2A1 haplotype rs12822608/rs6823 GC (p = 3.68 × 10(-4)). The strongest epistatic interactions were observed between MSX1 and BMP2, FGF1 and PVRL2, and COL2A1 and FGF2 genes. CONCLUSIONS: This study provides for the first time evidence of the implication of IRF6, COL2A1, and WNT3 in the occurrence of NSCP. It is likely that variation in cartilage collagen II and XI genes, IRF6, and the Wnt and FGF signaling pathway genes contributes susceptibility to nonsyndromic cleft palate in Northeastern European populations.

Genetic variation spectrum in ATP7B gene identified in Latvian patients with Wilson disease.

BACKGROUND: Wilson disease (WD) is an autosomal recessive disorder of copper metabolism caused by allelic variants in ATP7B gene. More than 500 distinct variants have been reported, the most common WD causing allelic variant in the patients from Central, Eastern, and Northern Europe is H1069Q. METHODS: All Latvian patients with clinically confirmed WD were screened for the most common mutation p.H1069Q by PCR Bi-PASA method. Direct DNA sequencing of gene ATP7B (all 21 exons) was performed for the patients with WD symptoms, being either heterozygous for H1069Q or without it on any allele. RESULTS: We identified 15 different allelic variants along with eight non-disease-causing allelic variants. Based on the gene molecular analysis and patients' clinical data variant p.His1069Gln was found in 66.9% of WD alleles. Wide clinical variability was observed among individuals with the same ATP7B genotype. The results of our study confirm that neurological manifestations of WD are typically present later than the liver disease but no significant association between the presence/absence of the most common genetic variant and mode of initial WD presentation or age at presentation was identified. CONCLUSIONS: (1) The most prevalent mutation in Latvian patients with Wilson disease was c.3207C>A (p.His1069Gln); (2) No significant phenotype-genotype correlation was found in Latvian patients with Wilson disease; (3) The estimated prevalence of Wilson disease in Latvia is 1 of 24,000 cases which is higher than frequently quoted prevalence of 1: 30,000.

Mutation analysis of the MSX1 gene exons and intron in patients with nonsyndromic cleft lip and palate.

UNLABELLED: Cleft lip with or without cleft palate and cleft palate (CL/CLP/CP) is one of the most common malformations among newborns. The estimated prevalence in Latvia is 1/700. Nonsyndromic CL/CLP/CP is a complex trait determined by multiple, interacting genetic and environmental factors. MSX1 gene is one of the most important candidate-genes, which had been analyzed in relation with nonsyndromic CL/CLP/CP. The objective of our study was to examine the etiologic role of MSX1 gene mutations in the development of nonsyndromic CL/CLP/CP in Latvian population. MATERIALS AND METHODS: DNA was extracted from venous blood of 53 patients with cleft lip with or without palate. Polymerase chain reaction (PCR) was performed of selected segments of MSX1 gene. These were sequenced and analysed by comparison with reference sequence, accession Nr. AF426432 (NCBI). RESULTS: 16 DNA sequence variations were identified in 53 patient samples; 6 of them have not been previously described. Identified sequence variations localized in coding regions do not cause amino acid substitutions, therefore they are not considered as mutations with an etiological role in CL/CLP/CP development. Baltic-Taiwan joint research project "Identification of genes involved in craniofacial morphogenesis and susceptibility to orofacial clefting in a human genome scan 2004-2006".

Comparison of telomere length between population-specific mitochondrial haplogroups among different age groups in a Latvian population.

Population studies have demonstrated that telomere length (TL) displays great diversity among different populations. Previously described controversial findings associated longevity with specific mitochondrial DNA haplogroups (hgs) (e.g., J and U). These observations may be influenced by population diversity, geographic location, and/or specific historic background. The aims of this study were to identify a specific hg which correlates with aging in a Latvian populating and to evaluate the possible association of TL variability with specific mitochondrial hgs. The results show no significant correlation between TL, mitochondrial DNA hgs and longevity. A slight increase in frequency was observed among centenarians of hg H; however, these findings were not statistically significant. TL did not show any statically significant difference, only hg W had slightly longer telomeres among others. An insignificant increase in TL was observed in the 55-89 age group of hg W but in the <90 age group for hg J which also had the longest TL in the 20-45 age group. In conclusion this study indicates that specific mitochondrial DNA hgs do not have a significant, if any, influence on the variation of TL in Latvians.

Association between inherited monogenic liver disorders and chronic hepatitis C.

AIM: To determine the frequencies of mutations that cause inherited monogenic liver disorders in patients with chronic hepatitis C. METHODS: This study included 86 patients with chronic hepatitis C (55 men, 31 women; mean age at diagnosis, 38.36 ± 14.52 years) who had undergone antiviral therapy comprising pegylated interferon and ribavirin. Viral load, biochemical parameter changes, and liver biopsy morphological data were evaluated in all patients. The control group comprised 271 unrelated individuals representing the general population of Latvia for mutation frequency calculations. The most frequent mutations that cause inherited liver disorders [gene (mutation): ATP7B (H1069Q), HFE (C282Y, H63D), UGT1A1 (TA)7, and SERPINA1 (PiZ)] were detected by polymerase chain reaction (PCR), bidirectional PCR allele-specific amplification, restriction fragment length polymorphism analysis, and sequencing. RESULTS: The viral genotype was detected in 80 of the 86 patients. Viral genotypes 1, 2, and 3 were present in 61 (76%), 7 (9%), and 12 (15%) patients, respectively. Among all 86 patients, 50 (58%) reached an early viral response and 70 (81%) reached a sustained viral response. All 16 patients who did not reach a sustained viral response had viral genotype 1. Case-control analysis revealed a statistically significant difference in only the H1069Q mutation between patients and controls (patients, 0.057; controls, 0.012; odds ratio, 5.514; 95%CI: 1.119-29.827, P = 0.022). However, the H1069Q mutation was not associated with antiviral treatment outcomes or biochemical indices. The (TA) 7 mutation of the UGT1A1 gene was associated with decreased ferritin levels (beta regression coefficient = -295.7, P = 0.0087). CONCLUSION: Genetic mutations that cause inherited liver diseases in patients with hepatitis C should be studied in detail.

Association of BMP4 polymorphisms with non-syndromic cleft lip with or without cleft palate and isolated cleft palate in Latvian and Lithuanian populations.

Cleft lip with or without cleft palate (CLP and CL, respectively) and isolated cleft palate (CP) represent one of the most common human birth defects, with a prevalence of approximately 1 in 300-2500 depending on the population. Formation of non-syndromic CL/CLP and CP arises from the interaction of environmental and genetic factors. The objective of this study was to investigate the association between the BMP4 gene (encoding bone morphogenetic protein 4) and non-syndromic CL/CLP and CP in order to clarify the role of this gene in the aetiology of the malformation in Latvian and Lithuanian populations. We genotyped three markers of the BMP4 gene (rs17563, rs2071047 and rs1957860) in order to perform single marker and haplotype association analyses for Latvian and Lithuanian non-syndromic CL/CLP and CP patients and controls. Transmission disequilibrium test was also conducted for Latvian and Lithuanian proband-parent trios. The case-control analysis revealed that SNP rs2071047 allele A was associated with a decreased risk of CL/CLP in the Latvian population, which was confirmed by the haplotype analysis. A modest association was detected between SNP rs1957860 and CP in the Lithuanian population, where allele C was associated with a decreased risk of this cleft phenotype, corroborating haplotype analysis data. Our findings support a role of the BMP4 gene in the aetiology of non-syndromic CL/CLP and CP in the studied populations.

Dynamics of telomere length in different age groups in a Latvian population.

The shortening of telomeres with ageing is a well-documented observation; however, the reported number of nucleotides in telomeres varies between different laboratories and studies. Such variability is likely caused by ethnic differences between the populations studied. Until now, there were no studies that investigated the variability of telomere length in a senescent Latvian population of the most common mitochondrial haplogroups, defined as H (45%), U (25%), Y chromosomal N1c (40%) and R1a1 (40%). Telomere length was determined in 121 individuals in different age groups, including a control group containing individuals of 20-40 years old and groups of individuals between 60-70 years old, 71-80 years old, 81-90 years old, and above 90 years old. Telomere length was determined using the Southern blot telomeric restriction fragment assay (TRF). Decreased telomere length with ageing was confirmed, but a comparison of centenarians and individuals between 60-90 years of age did not demonstrate a significant difference in telomere length. However, significant variability in telomere length was observed in the control group, indicating probable rapid telomere shortening in some individuals that could lead up to development of health status decline appearing with ageing. Telomere length measured in mononuclear blood cells (MNC) was compared with the telomere length measured in whole peripheral white blood cells (WBC) using TRF. Telomere length in MNC was longer than in WBC for the control group with individuals 20 to 40 years old; in contrast, for the group of individuals aged 65 to 85 years old, measured telomere length was shorter in MNC when compared to WBC.

Mitochondrial DNA origins of the Latvian clefting population.

Latvia has one of the highest prevalence of isolated cleft lip with or without cleft palate (CL/P) in Europe. To clarify the genetic origins of the Latvian cleft population and establish a method for genetic mapping, mitochondrial DNA variation was studied in a population affected with clefting. One-hundred and seven subjects and 351 samples from unrelated healthy volunteers representing four anthropologically, archaeologically and ethno-linguistically different regions of Latvia were selected. The case group showed a higher frequency of haplogroups U4 (p=0.02) and U5 (p=0.0003) than in non-U haplogroups. We hypothesize that U4 and U5 mtDNA haplotype carriers may also carry susceptibility genes for clefts. Future studies will take into consideration these definitions based on mtDNA haplotypes when analyzing genetic variations and their possible contribution to CL/P susceptibility.

The link between mitochondrial DNA hypervariable segment I heteroplasmy and ageing among genetically unrelated Latvians.

Various studies have demonstrated that mitochondrial DNA (mtDNA) heteroplasmy tends to increase with age and that the observed frequency of heteroplasmy among populations mostly depends on the way it is measured. Therefore, we investigated age-related association on the presence of mtDNA heteroplasmy within the hypervariable segment 1 (HVS-I) in a selected study group. The study group consisted of 300 maternally unrelated Latvians ranging in age from 18 to over 90 years. To determine the optimal method for mtDNA heteroplasmy detection, three approaches were used: (i) SURVEYOR Mutation Detection Kit, (ii) sequencing and (iii) denaturing gradient-gel electrophoresis (DGGE). Among the studied individuals, 30.3% were found to be heteroplasmic. The distribution of heteroplasmy statistically significantly increased with individuals' age (17%; 95% confidence interval [CI] 0.095-0.244 in the 18-40 year age group vs. 39%; [CI] 0.294-0.487 in the >90 year age group). Heteroplasmy occurred in a total of 21 different positions within HVS-I, and was the most frequent at fast-mutated positions 16189, 16304 and 16311. The results indicate that heteroplasmy in HVS-I is relatively common and occurs in a broad spectrum of sites. The above is supported by evidence to eventual increase of the probability of heteroplasmy with age due to specific mitochondrial haplogroup background.

The western and eastern roots of the Saami--the story of genetic "outliers" told by mitochondrial DNA and Y chromosomes.

The Saami are regarded as extreme genetic outliers among European populations. In this study, a high-resolution phylogenetic analysis of Saami genetic heritage was undertaken in a comprehensive context, through use of maternally inherited mitochondrial DNA (mtDNA) and paternally inherited Y-chromosomal variation. DNA variants present in the Saami were compared with those found in Europe and Siberia, through use of both new and previously published data from 445 Saami and 17,096 western Eurasian and Siberian mtDNA samples, as well as 127 Saami and 2,840 western Eurasian and Siberian Y-chromosome samples. It was shown that the "Saami motif" variant of mtDNA haplogroup U5b is present in a large area outside Scandinavia. A detailed phylogeographic analysis of one of the predominant Saami mtDNA haplogroups, U5b1b, which also includes the lineages of the "Saami motif," was undertaken in 31 populations. The results indicate that the origin of U5b1b, as for the other predominant Saami haplogroup, V, is most likely in western, rather than eastern, Europe. Furthermore, an additional haplogroup (H1) spread among the Saami was virtually absent in 781 Samoyed and Ob-Ugric Siberians but was present in western and central European populations. The Y-chromosomal variety in the Saami is also consistent with their European ancestry. It suggests that the large genetic separation of the Saami from other Europeans is best explained by assuming that the Saami are descendants of a narrow, distinctive subset of Europeans. In particular, no evidence of a significant directional gene flow from extant aboriginal Siberian populations into the haploid gene pools of the Saami was found.

Phylogeography of Y-chromosome haplogroup I reveals distinct domains of prehistoric gene flow in europe.

To investigate which aspects of contemporary human Y-chromosome variation in Europe are characteristic of primary colonization, late-glacial expansions from refuge areas, Neolithic dispersals, or more recent events of gene flow, we have analyzed, in detail, haplogroup I (Hg I), the only major clade of the Y phylogeny that is widespread over Europe but virtually absent elsewhere. The analysis of 1,104 Hg I Y chromosomes, which were identified in the survey of 7,574 males from 60 population samples, revealed several subclades with distinct geographic distributions. Subclade I1a accounts for most of Hg I in Scandinavia, with a rapidly decreasing frequency toward both the East European Plain and the Atlantic fringe, but microsatellite diversity reveals that France could be the source region of the early spread of both I1a and the less common I1c. Also, I1b*, which extends from the eastern Adriatic to eastern Europe and declines noticeably toward the southern Balkans and abruptly toward the periphery of northern Italy, probably diffused after the Last Glacial Maximum from a homeland in eastern Europe or the Balkans. In contrast, I1b2 most likely arose in southern France/Iberia. Similarly to the other subclades, it underwent a postglacial expansion and marked the human colonization of Sardinia approximately 9,000 years ago.