On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools.

BACKGROUND: In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation. RESULTS: We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed [Formula: see text] profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds. CONCLUSIONS: Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

Genetic structure of Iranian indigenous sheep breeds: insights for conservation.

Genetic structure and genetic diversity levels of indigenous Iranian sheep breeds are not clear, despite the interest this region has in itself as an important center for domestication of livestock. Early population genetic studies have reported high levels of diversity among Iranian sheep breeds until recently, when high admixture levels and genetic homogeneity have been detected. The rapid reduction of diversity observed in Iranian breeds might be due to an increasing trend of intensive crossbreeding practices or even total replacement of native breeds by highly specialized and productive ones. From a conservative perspective, this situation is highly concerning; thus, it might be wise to consider a conservation program in Iran to preserve the original genetic diversity in native sheep breeds. In this study, a total of 1065 animals with the purest morphological features representing 24 Iranian indigenous sheep breeds were sampled, corresponding to ancestral breed diversity. These samples were genotyped for 17 microsatellite loci in order to (1) determine the native ancestral diversity of Iranian breeds, (2) define the degree of genetic relationship among studied breeds, and (3) assess conservation priorities among defined groups. Our results showed no recent loss of diversity, but high genetic diversity levels for indigenous sheep breeds in Iran. Indeed, the analysis of conservation priorities pointed out the importance of 8 breeds for maintaining Iranian sheep breeds' maximum genetic diversity. Thus, under a genetic perspective, these 8 breeds should be the ones included into conservation programs for restocking endangered areas.

Association between combinations of genetic polymorphisms and epidemiopathogenic forms of bovine paratuberculosis.

Control of major mycobacterial diseases affecting livestock is a challenging issue that requires different approaches. The use of genetic markers for improving resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle has been explored as a promising population strategy We performed paratuberculosis epidemiopathogenic phenotypic and genotypic characterization involving 24 SNPs in six candidate genes (NOD2, CD209, SLC11A1, SP110, TLR2 and TLR4) on 502 slaughtered Friesian cows. In the current study, we investigate whether recently proposed paratuberculosis (PTB) epidemiopathogenic (EP) forms (apparently free-AF, latent-LAT and patent-PAT) could be associated with some combination of these 24 SNPs. Best EP form grouping was obtained using a combination of 5 SNPs in four genes (CD209: rs210748127; SLC11A1: rs110090506; SP110: rs136859213 and rs110480812; and TLR2: rs41830058). These groups were defined according to the level of infection progression risk to patent epidemiopathogenic forms and showed the following distributions: LOWIN (low) with 39 (8%) cases (94.9% AF/5.1% LAT/0% PAT); LATIN (low) with 17 (3%) cases (5.9% AF/94.1% LAT/0% PAT); AVERIN (average) with 413 (82%) cases (52.1% AF/38.5% LAT/9.4% PAT) and PATIN (patent) with 33 (7%) cases (36.4% AF/24.2% LAT/39.4% PAT). Age of slaughter was significantly higher for LATIN (88.3 months) compared to AVERIN (65.3 months; p = 0.0007) and PATIN (59.1 months; p = 0.0004), and for LOWIN (73.9 months) compared to PATIN (p = 0.0233), and nearly significant compared to AVERIN (p = 0.0572) These results suggest that some selected genetic polymorphisms have a potential use as markers of PTB EP forms and thus add a new tool for the control of this widespread infection.

Genomic selection signatures in sheep from the Western Pyrenees.

BACKGROUND: The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades. RESULTS: Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index Fst. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program. CONCLUSIONS: Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Genetic association study of systemic lupus erythematosus and disease subphenotypes in European populations.

Epidemiological studies suggest a strong contribution of genetic factors in the pathogenesis of systemic lupus erythematosus (SLE). In the last decades, many risk loci have been identified in several genetic association studies following both candidate gene and genome-wide approaches. The present work was conducted by GAPAID (Genes And Proteins for AutoImmunity Diagnostics) consortium with a dual aim: to replicate the association of several previously reported SLE susceptibility loci in an independent European sample and to explore their relation with some disease subphenotypes. A total of 48 single nucleotide polymorphisms (SNP) from 40 associated loci were typed in a cohort of 208 SLE patients and 152 controls from Rheumatology Units of the University Hospital of Pisa (Italy) and University of Pécs Medical Center (Hungary). Regression analyses were performed to detect disease susceptibility loci and to identify genes affecting specific disease manifestations (renal, neurological, or skin involvement; arthritis; secondary Sjögren syndrome; and secondary antiphospholipid syndrome). Association of previously described risk alleles from HLA locus has been replicated, while IRF5, BLK, ITGAM, and IRF8 loci have been found to be consistent with previous published results. In addition, two new subphenotype-specific associations have been detected: SNP rs5754217 (UBE2L3) with skin involvement and rs3093030 (ICAM1-ICAM4-ICAM5) with hematological disorders. Overall, results from GAPAID project are consistent with previously established associations for HLA, IRF5, BLK, ITGAM, and IRF8 SLE susceptibility loci and report for the first time two subphenotype-specific associations.

Serological and Genetic Evidence for Altered Complement System Functionality in Systemic Lupus Erythematosus: Findings of the GAPAID Consortium.

Systemic lupus erythematosus is a chronic autoimmune disease with multifactorial ethiopathogenesis. The complement system is involved in both the early and late stages of disease development and organ damage. To better understand autoantibody mediated complement consumption we examined ex vivo immune complex formation on autoantigen arrays. We recruited patients with SLE (n = 211), with other systemic autoimmune diseases (n = 65) and non-autoimmune control subjects (n = 149). Standard clinical and laboratory data were collected and serum complement levels were determined. The genotype of SNP rs1143679 in the ITGAM gene was also determined. Ex vivo formation of immune complexes, with respect to IgM, IgG, complement C4 and C3 binding, was examined using a functional immunoassay on autoantigen microarray comprising nucleic acids, proteins and lipids. Complement consumption of nucleic acids increased upon binding of IgM and IgG even when serum complement levels were decreased due to consumption in SLE patients. A negative correlation between serum complement levels and ex vivo complement deposition on nucleic acid autoantigens is demonstrated. On the contrary, complement deposition on tested protein and lipid autoantigens showed positive correlation with C4 levels. Genetic analysis revealed that the non-synonymous variant rs1143679 in complement receptor type 3 is associated with an increased production of anti-dsDNA IgG antibodies. Notwithstanding, homozygous carriers of the previously reported susceptible allele (AA) had lower levels of dsDNA specific IgM among SLE patients. Both the non-synonymous variant rs1143679 and the high ratio of nucleic acid specific IgG/IgM were associated with multiple organ involvement. In summary, secondary complement deficiency in SLE does not impair opsonization of nucleic-acid-containing autoantigens but does affect other antigens and potentially other complement dependent processes. Dysfunction of the receptor recognizing complement opsonized immune complexes promotes the development of class-switched autoantibodies targeting nucleic acids.

Genetic variants associated with rheumatoid arthritis patients and serotypes in European populations.

OBJECTIVES: To replicate the association of rheumatoid arthritis (RA) susceptibility loci in an independent European sample and to assess their specificity with anti-citrullinated protein antibodies (ACPA) status. METHODS: A selection of 64 SNP previously associated with RA have been typed in a cohort of 267 RA patients (169 ACPA-positive and 98 ACPA-negative) and 152 controls from the Rheumatology Units of the University Hospital of Pisa (Italy) and the University of Pécs Medical Center (Hungary). Regression analyses were performed first considering overall RA patients and secondly, taking both serotype subgroups as different disease entities. The results have been adjusted for age, gender and origin of individuals. RESULTS: The well-known CD2, REL, TNFAIP3, IRF5, PTPRC, and CCR6 have been confirmed as RA disease associated loci together with recently discovered BACH2, RASGRP1, and IKZF3 loci, taking all RA patients as a unique phenotype. Results from both serological subgroups separately reflect the specificity of these susceptibility loci and show additional ACPA-positive specific associations for variants at IL6R, IL2RA, BLK, DDX6, IL6, and TLE3 genes. CONCLUSIONS: The results from GAPAID project are consistent with previously established RA disease associations for CD2, PTPRC, REL, CCR6, TNFAIP3, IRF5, BLK, IL2RA, and DDX6 loci. In addition, IL6R, BACH2, RASGRP1, TLE3, and IKZF3 are replicated for the first time in an independent European population and IL6 appears to be a suggestive new RA associated locus. The stratified analysis based on ACPA status provides further support for distinct genetic aetiologies of RA subsets, which might have therapeutic implications.

Influence of MTHFR C677T polymorphism on methotrexate monotherapy discontinuation in rheumatoid arthritis patients: results from the GAPAID European project.

OBJECTIVES: Methotrexate (MTX) is the most widely prescribed drug for rheumatoid arthritis (RA) patients, but 45% of them discontinue therapy within two years, either due to inefficacy or toxicity. Several authors have reported contradictory results related to C677T polymorphism in the MTHFR gene and response to MTX in RA. The purpose of this study was to further explore this genotype-response association in a European RA population. METHODS: This retrospective longitudinal study included a total of 269 RA patients from Italy and Hungary, of whom 73.2% had available data on MTX treatment (197 patients). C677T polymorphism (rs1801133) was genotyped by quantitative PCR using TaqMan assays. Genotype association analysis and Kaplan-Meier method were used for statistical comparisons between patients continuing and patients who abandoned MTX treatment. RESULTS: A total of 85 out of the 197 RA patients (43%) abandoned MTX treatment by the time of analysis. No significant genotype-MTX discontinuation association was found for the overall population, either at the end of the study (p=0.375), or during the follow-up (p=0.324). When the analysis was restricted to the 68 patients on MTX monotherapy, a borderline association (OR 3.15, 95% CI 0.93-10.67, p=0.057) was noted with the recessive genetic model. In agreement with that, a Kaplan-Meier analysis showed a significantly shorter time-to-discontinuation of MTX monotherapy for homozygous carriers of the T-allele (p=0.042). CONCLUSIONS: These results demonstrate that the C677T polymorphism in the MTHFR gene is involved in MTX monotherapy discontinuation in a multicentre European patient cohort, confirming previous results.

Genetic association analysis of paratuberculosis forms in holstein-friesian cattle.

A genetic susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infections in ruminants has been longtime suspected to exist. Recently, natural infections in cattle have been reclassified into latent and patent forms based on histopathological findings and their associations with immunological and microbiological variables. This study aims to explore whether these newly defined phenotypes are associated with twenty-four single-nucleotide polymorphisms (SNPs) in six bovine candidate genes: nucleotide-binding oligomerization domain 2 (NOD2), solute carrier family 11 member A1 (SLC11A1), nuclear body protein SP110 (SP110), toll-like receptors (TLRs) 2 and 4, and CD209 (also known as DC-SIGN, dendritic cell-specific ICAM3-grabbing nonintegrin). SNPs were genotyped for 772 Holstein-Friesian animals (52.6% apparently free; 38.1% latent; 9.3% patent) by TaqMan OpenArray technology. Genotypic-phenotypic associations were assessed by logistic regression analysis adjusted for age at slaughter, under five models (codominant, dominant, recessive, overdominant, and log-additive), and corrected for multiple testing. The rs208222804 C allele (CD209 gene) was found to be associated with latent paratuberculosis (log-additive model: P < 0.0034 after permutation procedure; OR = 0.64, 95% CI = 0.48-0.86). No significant association was detected between any SNP and the patent phenotype. Consequently, CD209 gene may play a key role in the pathogenesis of bovine paratuberculosis.