Changes in environmental exposures over decades may influence the genetic architecture of severe spermatogenic failure

STUDY QUESTION

Do the genetic determinants of idiopathic severe spermatogenic failure (SPGF) differ between generations?

SUMMARY ANSWER

Our data support that the genetic component of idiopathic SPGF is impacted by dynamic changes in environmental exposures over decades.

WHAT IS KNOWN ALREADY

The idiopathic form of SPGF has a multifactorial etiology wherein an interaction between genetic, epigenetic, and environmental factors leads to the disease onset and progression. At the genetic level, genome-wide association studies (GWASs) allow the analysis of millions of genetic variants across the genome in a hypothesis-free manner, as a valuable tool for identifying susceptibility risk loci. However, little is known about the specific role of non-genetic factors and their influence on the genetic determinants in this type of conditions.

STUDY DESIGN, SIZE, DURATION

Case-control genetic association analyses were performed including a total of 912 SPGF cases and 1360 unaffected controls.

PARTICIPANTS/MATERIALS, SETTING, METHODS

All participants had European ancestry (Iberian and German). SPGF cases were diagnosed during the last decade either with idiopathic non-obstructive azoospermia (n = 547) or with idiopathic non-obstructive oligozoospermia (n = 365). Case-control genetic association analyses were performed by logistic regression models considering the generation as a covariate and by in silico functional characterization of the susceptibility genomic regions.

MAIN RESULTS AND THE ROLE OF CHANCE

This analysis revealed 13 novel genetic association signals with SPGF, with eight of them being independent. The observed associations were mostly explained by the interaction between each lead variant and the age-group. Additionally, we established links between these loci and diverse non-genetic factors, such as toxic or dietary habits, respiratory disorders, and autoimmune diseases, which might potentially influence the genetic architecture of idiopathic SPGF.

LARGE SCALE DATA

GWAS data are available from the authors upon reasonable request.

LIMITATIONS, REASONS FOR CAUTION

Additional independent studies involving large cohorts in ethnically diverse populations are warranted to confirm our findings.

WIDER IMPLICATIONS OF THE FINDINGS

Overall, this study proposes an innovative strategy to achieve a more precise understanding of conditions such as SPGF by considering the interactions between a variable exposome through different generations and genetic predisposition to complex diseases.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the "Plan Andaluz de Investigación, Desarrollo e Innovación (PAIDI 2020)" (ref. PY20_00212, P20_00583), the Spanish Ministry of Economy and Competitiveness through the Spanish National Plan for Scientific and Technical Research and Innovation (ref. PID2020-120157RB-I00 funded by MCIN/ AEI/10.13039/501100011033), and the 'Proyectos I+D+i del Programa Operativo FEDER 2020' (ref. B-CTS-584-UGR20). ToxOmics-Centre for Toxicogenomics and Human Health, Genetics, Oncology and Human Toxicology, is also partially supported by the Portuguese Foundation for Science and Technology (Projects: UIDB/00009/2020; UIDP/00009/2020). The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Contribution of TEX15 genetic variants to the risk of developing severe non-obstructive oligozoospermia

Background: Severe spermatogenic failure (SPGF) represents one of the most relevant causes of male infertility. This pathological condition can lead to extreme abnormalities in the seminal sperm count, such as severe oligozoospermia (SO) or non-obstructive azoospermia (NOA). Most cases of SPGF have an unknown aetiology, and it is known that this idiopathic form of male infertility represents a complex condition. In this study, we aimed to evaluate whether common genetic variation in TEX15, which encodes a key player in spermatogenesis, is involved in the susceptibility to idiopathic SPGF. Materials and Methods: We designed a genetic association study comprising a total of 727 SPGF cases (including 527 NOA and 200 SO) and 1,058 unaffected men from the Iberian Peninsula. Following a tagging strategy, three tag single-nucleotide polymorphisms (SNPs) of TEX15 (rs1362912, rs323342, and rs323346) were selected for genotyping using TaqMan probes. Case-control association tests were then performed by logistic regression models. In silico analyses were also carried out to shed light into the putative functional implications of the studied variants. Results: A significant increase in TEX15-rs1362912 minor allele frequency (MAF) was observed in the group of SO patients (MAF = 0.0842) compared to either the control cohort (MAF = 0.0468, OR = 1.90, p = 7.47E-03) or the NOA group (MAF = 0.0472, OR = 1.83, p = 1.23E-02). The genotype distribution of the SO population was also different from those of both control (p = 1.14E-02) and NOA groups (p = 4.33-02). The analysis of functional annotations of the human genome suggested that the effect of the SO-associated TEX15 variants is likely exerted by alteration of the binding affinity of crucial transcription factors for spermatogenesis. Conclusion: Our results suggest that common variation in TEX15 is involved in the genetic predisposition to SO, thus supporting the notion of idiopathic SPGF as a complex trait.

A GWAS in Idiopathic/Unexplained Infertile Men Detects a Genomic Region Determining Follicle-Stimulating Hormone Levels

CONTEXT

Approximately 70% of infertile men are diagnosed with idiopathic (abnormal semen parameters) or unexplained (normozoospermia) infertility, with the common feature of lacking etiologic factors. Follicle-stimulating hormone (FSH) is essential for initiation and maintenance of spermatogenesis. Certain single-nucleotide variations (SNVs; formerly single-nucleotide polymorphisms [SNPs]) (ie, FSHB c.-211G > T, FSHR c.2039A > G) are associated with FSH, testicular volume, and spermatogenesis. It is unknown to what extent other variants are associated with FSH levels and therewith resemble causative factors for infertility.

OBJECTIVE

We aimed to identify further genetic determinants modulating FSH levels in a cohort of men presenting with idiopathic or unexplained infertility.

METHODS

We retrospectively (2010-2018) selected 1900 men with idiopathic/unexplained infertility. In the discovery study (n = 760), a genome-wide association study (GWAS) was performed (Infinium PsychArrays) in association with FSH values (Illumina GenomeStudio, v2.0). Minor allele frequencies (MAFs) were analyzed for the discovery and an independent normozoospermic cohort. In the validation study (n = 1140), TaqMan SNV polymerase chain reaction was conducted for rs11031005 and rs10835638 in association with andrological parameters.

RESULTS

Imputation revealed 9 SNVs in high linkage disequilibrium, with genome-wide significance (P < 4.28e-07) at the FSHB locus 11p.14.1 being associated with FSH. The 9 SNVs accounted for up to a 4.65% variance in FSH level. In the oligozoospermic subgroup, this was increased up to 6.95% and the MAF was enhanced compared to an independent cohort of normozoospermic men. By validation, a significant association for rs11031005/rs10835638 with FSH (P = 4.71e-06/5.55e-07) and FSH/luteinizing hormone ratio (P = 2.08e-12/6.4e-12) was evident.

CONCLUSIONS

This GWAS delineates the polymorphic FSHB genomic region as the main determinant of FSH levels in men with unexplained or idiopathic infertility. Given the essential role of FSH, molecular detection of one of the identified SNVs that causes lowered FSH and therewith decreases spermatogenesis could resolve the idiopathic/unexplained origin by this etiologic factor.

Common genetic variation in KATNAL1 non-coding regions is involved in the susceptibility to severe phenotypes of male infertility

Background

Previous studies in animal models evidenced that genetic mutations of KATNAL1, resulting in dysfunction of its encoded protein, lead to male infertility through disruption of microtubule remodelling and premature germ cell exfoliation. Subsequent studies in humans also suggested a possible role of KATNAL1 single‐nucleotide polymorphisms in the development of male infertility as a consequence of severe spermatogenic failure.

Objectives

The main objective of the present study is to evaluate the effect of the common genetic variation of KATNAL1 in a large and phenotypically well‐characterised cohort of infertile men because of severe spermatogenic failure.

Materials and methods

A total of 715 infertile men because of severe spermatogenic failure, including 210 severe oligospermia and 505 non‐obstructive azoospermia patients, as well as 1058 unaffected controls were genotyped for three KATNAL1 single‐nucleotide polymorphism taggers (rs2077011, rs7338931 and rs2149971). Case–control association analyses by logistic regression assuming different models and in silico functional characterisation of risk variants were conducted.

Results

Genetic associations were observed between the three analysed taggers and different severe spermatogenic failure groups. However, in all cases, the haplotype model (rs2077011*C | rs7338931*T | rs2149971*A) better explained the observed associations than the three risk alleles independently. This haplotype was associated with non‐obstructive azoospermia (adjusted p = 4.96E‐02, odds ratio = 2.97), Sertoli‐cell only syndrome (adjusted p = 2.83E‐02, odds ratio = 5.16) and testicular sperm extraction unsuccessful outcomes (adjusted p = 8.99E‐04, odds ratio = 6.13). The in silico analyses indicated that the effect on severe spermatogenic failure predisposition could be because of an alteration of the KATNAL1 splicing pattern.

Conclusions

Specific allelic combinations of KATNAL1 genetic polymorphisms may confer a risk of developing severe male infertility phenotypes by favouring the overrepresentation of a short non‐functional transcript isoform in the testis.

O-118 New insight into the genetic contribution of common variants to the development of extreme phenotypes of unexplained male infertility: a multicenter genome-wide association study

Study question

What is the contribution of the common genetic variation to the development of unexplained male infertility due to severe spermatogenic failure (SPGF)?

Summary answer

Genetic polymorphisms of key immune and spermatogenesis loci are involved in the etiology of the most severe SPGF cases, defined by Sertoli cell-only (SCO) phenotype.

What is known already

Male infertility is a rising worldwide concern that affects millions of couples. Non-obstructive azoospermia (NOA) and severe oligospermia (SO) are two extreme manifestations characterized by SPGF. A genetic cause can be established in only around 20% of affected men, with the remaining cases being classified as otherwise unexplained. To date, the genome-wide association study (GWAS) strategy, although already successfully applied in several other complex traits and diseases, was less fruitful in studies that attempted to decipher the genetic component of unexplained SPGF, mainly due to both a lack of well-powered samples in different ancestries and limitations in study design.

Study design, size, duration

We designed a GWAS for unexplained male infertility due to SPGF including a total of 1,274 affected cases and 1,951 fertile controls from the Iberian Peninsula (Spain and Portugal) and Germany. Different biostatistics and bioinformatics approaches were used to evaluate the possible effect of single-nucleotide polymorphisms (SNPs) across the whole genome in the susceptibility to specific subtypes of unexplained SPGF.

Participants/materials, setting, methods

The case cohort comprised 502 SO and 772 NOA patients, who were subdivided according to histological phenotypes (SCO, maturation arrest, and hypospermatogenesis) and the outcome of testicular sperm extraction techniques (TESE) from testis biopsies. Genotyping was performed with the GSA platform (Illumina). After quality-control and genotype imputation, 6,539,982 SNPs remained for the analysis, which was performed by logistic regression models. The datasets went through a meta-analysis by the inverse variance weighted method under fixed effects.

Main results and the role of chance

Genetic associations with SCO at the genome-wide-level of significance were identified in the major histocompatibility (MHC) class II region (rs1136759, OR = 1.80, P = 1.32E-08) and in a regulatory region of chromosome 14 nearby the vaccinia-related kinase 1 (VRK1) gene (rs115054029, OR = 3.14, P = 4.37-08). VRK1 is a relevant proliferative factor for spermatogenesis that causes progressive loss of spermatogonia when disrupted in mouse models. The role of the MHC system in SCO susceptibility was comprehensively evaluated through a validated imputation method that infers classical MHC alleles and polymorphic amino acid positions. A serine at position 13 of the HLA-DRβ1 protein (defined by the risk allele of the lead variant rs1136759) explained most of the SCO association signals within the MHC class II region. This residue is located in the binding pocket of the HLA-DR molecule and interacts directly with the presented antigen. Interestingly, position 13 of HLA-DRβ1 is the most relevant risk amino acid position for a wide spectrum of immune-mediated disorders. The HLA-DRB1*13 haplotype (which includes the serine at position 13 and represents the strongest NOA-associated marker in Asians to date) was the strongest signal amongst the classical MHC alleles in our study cohort (OR = 1.93, P = 9.90E-07).

Limitations, reasons for caution

Although the statistical power for the overall analysis was appropriate, the subphenotype analyses performed had considerably lower counts, which may influence the identification of genetic variants conferring low to moderate risk effects. Independent studies in larger SCO study cohorts should be performed to confirm our findings.

Wider implications of the findings

The molecular mechanisms underlying unexplained SPGF are largely unknown. Our data suggest a relevant role of common genetic variation in the development of SCO, the most extreme histological phenotype of NOA. SCO is characterized by the loss of germ cells and, therefore, implies a considerably higher probability of unsuccessful TESE.

Trial registration number

N/A

Genetic Association of a Gain-of-Function IFNGR1 Polymorphism and the Intergenic Region LNCAROD/DKK1 With Behçet's Disease

OBJECTIVE

Behçet's disease is a complex systemic inflammatory vasculitis of incompletely understood etiology. This study was undertaken to investigate genetic associations with Behçet's disease in a diverse multiethnic population.

METHODS

A total of 9,444 patients and controls from 7 different populations were included in this study. Genotyping was performed using an Infinium ImmunoArray-24 v.1.0 or v.2.0 BeadChip. Analysis of expression data from stimulated monocytes, and epigenetic and chromatin interaction analyses were performed.

RESULTS

We identified 2 novel genetic susceptibility loci for Behçet's disease, including a risk locus in IFNGR1 (rs4896243) (odds ratio [OR] 1.25; P = 2.42 × 10-9 ) and within the intergenic region LNCAROD/DKK1 (rs1660760) (OR 0.78; P = 2.75 × 10-8 ). The risk variants in IFNGR1 significantly increased IFNGR1 messenger RNA expression in lipopolysaccharide-stimulated monocytes. In addition, our results replicated the association (P < 5 × 10-8 ) of 6 previously identified susceptibility loci in Behçet's disease: IL10, IL23R, IL12A-AS1, CCR3, ADO, and LACC1, reinforcing the notion that these loci are strong genetic factors in Behçet's disease shared across ancestries. We also identified >30 genetic susceptibility loci with a suggestive level of association (P < 5 × 10-5 ), which will require replication. Finally, functional annotation of genetic susceptibility loci in Behçet's disease revealed their possible regulatory roles and suggested potential causal genes and molecular mechanisms that could be further investigated.

CONCLUSION

We performed the largest genetic association study in Behçet's disease to date. Our findings reveal novel putative functional variants associated with the disease and replicate and extend the genetic associations in other loci across multiple ancestries.

Evaluation of Male Fertility-Associated Loci in a European Population of Patients with Severe Spermatogenic Impairment

Infertility is a growing concern in developed societies. Two extreme phenotypes of male infertility are non-obstructive azoospermia (NOA) and severe oligospermia (SO), which are characterized by severe spermatogenic failure (SpF). We designed a genetic association study comprising 725 Iberian infertile men as a consequence of SpF and 1058 unaffected controls to evaluate whether five single-nucleotide polymorphisms (SNPs), previously associated with reduced fertility in Hutterites, are also involved in the genetic susceptibility to idiopathic SpF and specific clinical entities. A significant difference in the allele frequencies of USP8-rs7174015 was observed under the recessive model between the NOA group and both the control group (p = 0.0226, OR = 1.33) and the SO group (p = 0.0048, OR = 1.78). Other genetic associations for EPSTI1-rs12870438 and PSAT1-rs7867029 with SO and between TUSC1-rs10966811 and testicular sperm extraction (TESE) success in the context of NOA were observed. In silico analysis of functional annotations demonstrated cis-eQTL effects of such SNPs likely due to the modification of binding motif sites for relevant transcription factors of the spermatogenic process. The findings reported here shed light on the molecular mechanisms leading to severe phenotypes of idiopathic male infertility, and may help to better understand the contribution of the common genetic variation to the development of these conditions.

Intronic variation of the SOHLH2 gene confers risk to male reproductive impairment

OBJECTIVE

To evaluate whether SOHLH2 intronic variation contributes to the genetic predisposition to male infertility traits, including severe oligospermia (SO) and different nonobstructive azoospermia (NOA) clinical phenotypes.

DESIGN

Genetic association study.

SETTING

Not applicable.

PATIENT(S)

Five hundred five cases (455 infertile patients diagnosed with NOA and 50 with SO) and 1,050 healthy controls from Spain and Portugal.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Genomic DNA extraction from peripheral blood mononuclear cells, genotyping of the SOHLH2 polymorphisms rs1328626 and rs6563386 using the TaqMan allelic discrimination technology, case-control association analyses using logistic regression models, and exploration of functional annotations in publicly available databases.

RESULT(S)

Evidence of association was observed for both rs6563386 with SO and rs1328626 with unsuccessful sperm retrieval after testicular sperm extraction (TESE-) in the context of NOA. A dominant effect of the minor alleles was suggested in both associations, either when the subset of patients with the manifestation were compared against the control group (rs6563386/SO: P=.021, odds ratio [OR] = 0.51; rs1328626/TESE-: P=.066, OR = 1.46) or against the group of patients without the manifestation (rs6563386/SO: P=.014, OR = 0.46; rs1328626/TESE-: P=.012, OR = 2.43). The haplotype tests suggested a combined effect of both polymorphisms. In silico analyses evidenced that this effect could be due to alteration of the isoform population.

CONCLUSION(S)

Our data suggest that intronic variation of SOHLH2 is associated with spermatogenic failure. The genetic effect is likely caused by different haplotypes of rs6563386 and rs1328626, which may predispose to SO or TESE- depending on the specific allelic combination.

Deficiency of the onco-miRNA cluster, miR-106b∼25, causes oligozoospermia and the cooperative action of miR-106b∼25 and miR-17∼92 is required to maintain male fertility

The identification of new genes involved in sexual development and gonadal function as potential candidates causing male infertility is important for both diagnostic and therapeutic purposes. Deficiency of the onco-miRNA cluster miR-17∼92 has been shown to disrupt spermatogenesis, whereas mutations in its paralog cluster, miR-106b∼25, that is expressed in the same cells, were reported to have no effect on testis development and function. The aim of this work is to determine the role of these two miRNA clusters in spermatogenesis and male fertility. For this, we analyzed miR-106b∼25 and miR-17∼92 single and double mouse mutants and compared them to control mice. We found that miR-106b∼25 knock out testes show reduced size, oligozoospermia and altered spermatogenesis. Transcriptomic analysis showed that multiple molecular pathways are deregulated in these mutant testes. Nevertheless, mutant males conserved normal fertility even when early spermatogenesis and other functions were disrupted. In contrast, miR-17∼92+/-; miR-106b∼25-/- double mutants showed severely disrupted testicular histology and significantly reduced fertility. Our results indicate that miR-106b∼25 and miR-17∼92 ensure accurate gene expression levels in the adult testis, keeping them within the required thresholds. They play a crucial role in testis homeostasis and are required to maintain male fertility. Hence, we have identified new candidate genetic factors to be screened in the molecular diagnosis of human males with reproductive disorders. Finally, considering the well-known oncogenic nature of these two clusters and the fact that patients with reduced fertility are more prone to testicular cancer, our results might also help to elucidate the molecular mechanisms linking both pathologies.

Genetic Landscape of Nonobstructive Azoospermia and New Perspectives for the Clinic

Nonobstructive azoospermia (NOA) represents the most severe expression of male infertility, involving around 1% of the male population and 10% of infertile men. This condition is characterised by the inability of the testis to produce sperm cells, and it is considered to have an important genetic component. During the last two decades, different genetic anomalies, including microdeletions of the Y chromosome, karyotype defects, and missense mutations in genes involved in the reproductive function, have been described as the primary cause of NOA in many infertile men. However, these alterations only explain around 25% of azoospermic cases, with the remaining patients showing an idiopathic origin. Recent studies clearly suggest that the so-called idiopathic NOA has a complex aetiology with a polygenic inheritance, which may alter the spermatogenic process. Although we are far from a complete understanding of the molecular mechanisms underlying NOA, the use of the new technologies for genetic analysis has enabled a considerable increase in knowledge during the last years. In this review, we will provide a comprehensive and updated overview of the genetic basis of NOA, with a special focus on the possible application of the recent insights in clinical practice.

GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways

Systemic sclerosis (SSc) is an autoimmune disease that shows one of the highest mortality rates among rheumatic diseases. We perform a large genome-wide association study (GWAS), and meta-analysis with previous GWASs, in 26,679 individuals and identify 27 independent genome-wide associated signals, including 13 new risk loci. The novel associations nearly double the number of genome-wide hits reported for SSc thus far. We define 95% credible sets of less than 5 likely causal variants in 12 loci. Additionally, we identify specific SSc subtype-associated signals. Functional analysis of high-priority variants shows the potential function of SSc signals, with the identification of 43 robust target genes through HiChIP. Our results point towards molecular pathways potentially involved in vasculopathy and fibrosis, two main hallmarks in SSc, and highlight the spectrum of critical cell types for the disease. This work supports a better understanding of the genetic basis of SSc and provides directions for future functional experiments.

Identification of a 3'-Untranslated Genetic Variant of RARB Associated With Carotid Intima-Media Thickness in Rheumatoid Arthritis: A Genome-Wide Association Study

Objective

To investigate the genetic background influencing the development of cardiovascular (CV) disease in patients with rheumatoid arthritis (RA).

Methods

We performed a genome‐wide association study (GWAS) in which, after quality control and imputation, a total of 6,308,944 polymorphisms across the whole genome were analyzed in 2,989 RA patients of European origin. Data on subclinical atherosclerosis, obtained through assessment of carotid intima‐media thickness (CIMT) and presence/absence of carotid plaques by carotid ultrasonography, were available for 1,355 individuals.

Results

A genetic variant of the RARB gene (rs116199914) was associated with CIMT values at the genome‐wide level of significance (minor allele [G] β coefficient 0.142, P = 1.86 × 10⁻⁸). Interestingly, rs116199914 overlapped with regulatory elements in tissues related to CV pathophysiology and immune cells. In addition, biologic pathway enrichment and predictive protein–protein relationship analyses, including suggestive GWAS signals of potential relevance, revealed a functional enrichment of the collagen biosynthesis network related to the presence/absence of carotid plaques (Gene Ontology no. 0032964; false discovery rate–adjusted P = 4.01 × 10⁻³). Furthermore, our data suggest potential influences of the previously described candidate CV risk loci NFKB1,MSRA, and ZC3HC1 (P = 8.12 × 10⁻⁴, P = 5.94 × 10⁻⁴, and P = 2.46 × 10⁻⁴, respectively).

Conclusion

The present findings strongly suggest that genetic variation within RARB contributes to the development of subclinical atherosclerosis in patients with RA.

The potential of PTPN22 as a therapeutic target for rheumatoid arthritis

INTRODUCTION

PTPN22 encodes a lymphoid-specific tyrosine phosphatase (LYP) that is a master regulator of the immune response. This gene is a major susceptibility factor for a wide range of autoimmune conditions, including rheumatoid arthritis (RA) for which it represents the strongest non-HLA contributor to disease risk. A missense PTPN22 allele (R620W) affecting the protein-protein interaction of LYP with other relevant players was described as the functional variant of the association. This review will focus on the role of PTPN22 in the pathogenic mechanisms underlying RA predisposition and discuss the possibility of developing LYP-based treatment strategies with a potential application in clinical practice. Areas covered: This review covers the literature showing how PTPN22 is implicated in signalling pathways involved in the autoimmune and autoinflammatory processes underlying RA. Insights obtained from studies aimed at developing novel selective LYP suppressors for treating RA are summarized. Expert opinion: Targeting key risk factors during the early steps of the disease may represent a good strategy to accomplish complete disease remission. As cumulating evidences suggest that PTPN22 R620W is a gain-of-function variant, a growing interest in developing LYP inhibitors has arisen. The potential efficacy and possible application of such compounds are discussed.

Comprehensive analysis of three TYK2 gene variants in the susceptibility to Chagas disease infection and cardiomyopathy

Tyrosine kinase 2 (TYK2) is a member of the Janus kinases family implicated in the signal transduction of type I interferons and several interleukins. It has been described that genetic mutations within TYK2 lead to multiple deleterious effects in the immune response. In this work, we have analyzed three functional independent variants from the frequency spectrum on the TYK2 gene (common and low-frequency variants) suggested to reduce the function of the gene in mediating cytokine signaling and the susceptibility to infections by Trypanosoma cruzi and/or the development of Chagas cardiomyopathy in the Colombian population. A total of 1,323 individuals from a Colombian endemic region for Chagas disease were enrolled in the study. They were classified as seronegative (n = 445), seropositive asymptomatic (n = 336), and chronic Chagas Cardiomyopathy subjects (n = 542). DNA samples were genotyped using TaqMan probes. Our results showed no statistically significant differences between the allelic frequencies of the three analyzed variants when seropositive and seronegative individuals were compared, therefore these variants were not associated with susceptibility to Chagas disease. Moreover, when Chagas cardiomyopathy patients were compared to asymptomatic patients, no significant associations were found. Previous reports highlighted the association of this gene in immune-related disorders under an autoimmunity context, but not predisposing patients to infectious diseases, which is consistent with our findings. Therefore, according to our results, TYK2 gene variants do not seem to play an important role in Chagas disease susceptibility and/or chronic Chagas cardiomyopathy.

A genome-wide association study suggests the HLA Class II region as the major susceptibility locus for IgA vasculitis

The genetic component of Immunoglobulin-A (IgA) vasculitis is still far to be elucidated. To increase the current knowledge on the genetic component of this vasculitis we performed the first genome-wide association study (GWAS) on this condition. 308 IgA vasculitis patients and 1,018 healthy controls from Spain were genotyped by Illumina HumanCore BeadChips. Imputation of GWAS data was performed using the 1000 Genomes Project Phase III dataset as reference panel. After quality control filters and GWAS imputation, 285 patients and 1,006 controls remained in the datasets and were included in further analysis. Additionally, the human leukocyte antigen (HLA) region was comprehensively studied by imputing classical alleles and polymorphic amino acid positions. A linkage disequilibrium block of polymorphisms located in the HLA class II region surpassed the genome-wide level of significance (OR = 0.56, 95% CI = 0.46–0.68). Although no polymorphic amino acid positions were associated at the genome-wide level of significance, P-values of potential relevance were observed for the positions 13 and 11 of HLA-DRB1 (P = 6.67E-05, P = 1.88E-05, respectively). Outside the HLA, potential associations were detected, but none of them were close to the statistical significance. In conclusion, our study suggests that IgA vasculitis is an archetypal HLA class II disease.

A Genome-wide Association Study Identifies Risk Alleles in Plasminogen and P4HA2 Associated with Giant Cell Arteritis

Giant cell arteritis (GCA) is the most common form of vasculitis in individuals older than 50 years in Western countries. To shed light onto the genetic background influencing susceptibility for GCA, we performed a genome-wide association screening in a well-powered study cohort. After imputation, 1,844,133 genetic variants were analyzed in 2,134 case subjects and 9,125 unaffected individuals from ten independent populations of European ancestry. Our data confirmed HLA class II as the strongest associated region (independent signals: rs9268905, p = 1.94 × 10^-54, per-allele OR = 1.79; and rs9275592, p = 1.14 × 10^-40, OR = 2.08). Additionally, PLG and P4HA2 were identified as GCA risk genes at the genome-wide level of significance (rs4252134, p = 1.23 × 10^-10, OR = 1.28; and rs128738, p = 4.60 × 10^-9, OR = 1.32, respectively). Interestingly, we observed that the association peaks overlapped with different regulatory elements related to cell types and tissues involved in the pathophysiology of GCA. PLG and P4HA2 are involved in vascular remodelling and angiogenesis, suggesting a high relevance of these processes for the pathogenic mechanisms underlying this type of vasculitis.

Evaluation of VDR gene polymorphisms in Trypanosoma cruzi infection and chronic Chagasic cardiomyopathy

Vitamin D is an important modulator of the immune response. It acts over several immune cell types where the Vitamin D receptor (VDR) is expressed. Due to the high relevance of this signaling pathway, several studies have investigated the possible influence of genes involved in the metabolism of Vitamin D and its receptor in different human diseases. Here, we analyzed whether four single-nucleotide polymorphisms of the VDR gene (rs731236, rs7975232, rs1544410 and rs2228570) are involved in the susceptibility to infection by Trypanosoma cruzi and/or to chronic Chagas cardiomyopathy (CCC) in a Colombian endemic population for this parasite. Our results showed that the rs2228570*A allele is associated with CCC development (P = 4.46E-03, OR = 1.51). In summary, the data presented in this report suggest that variation within the VDR gene may affect the immune response against T. cruzi, increasing the probability of cardiac complications in infected individuals.

Analysis of Systemic Sclerosis-associated Genes in a Turkish Population

OBJECTIVE

To evaluate the genetic background of systemic sclerosis (SSc) in the Turkish population.

METHODS

There were 354 cases and 718 unaffected controls from Turkey genotyped for the most relevant SSc genetic markers (IRF5-rs10488631, STAT4-rs3821236, CD247-rs2056626, DNASE1L3-rs35677470, IL12A-rs77583790, and ATG5-rs9373839). Association tests were conducted to identify possible associations.

RESULTS

Except for ATG5, all the analyzed genes showed either significant associations (IRF5: p = 1.32E-05, OR 1.76; CD247: p = 2.20E-03, OR 0.75) or trends of association (STAT4: p = 0.066, OR 1.21; IL12A: p = 0.079, OR 4.07; DNASE1L3: p = 0.097, OR 1.41) with the overall disease or with specific phenotypes.

CONCLUSION

The genetic component of SSc seems to be similar between Turks and Europeans.

Protective Role of the Interleukin 33 rs3939286 Gene Polymorphism in the Development of Subclinical Atherosclerosis in Rheumatoid Arthritis Patients

OBJECTIVES

To determine whether the interleukin-33 (IL-33)-interleukin-1 receptor like 1 (IL-1RL1) signaling pathway is implicated in the risk of subclinical atherosclerosis in patients with rheumatoid arthritis (RA).

METHODS

A total of 576 Spanish RA patients from Northern Spain were genotyped for 6 well-known IL33-IL1RL1 polymorphisms (IL33 rs3939286, IL33 rs7025417, IL33 rs7044343, IL1RL1 rs2058660, IL1RL1 rs2310173 and IL1RL1 rs13015714) by TaqMan genotyping assay. The presence of subclinical atherosclerosis was determined by the assessment of carotid intima-media thickness (cIMT) by carotid ultrasound (US).

RESULTS

RA patients carrying the TT genotype of the IL33 rs3939286 polymorphism had lower cIMT values than those homozygous for the CC genotype (mean ± standard deviation (SD): 0.71 ± 0.14 mm versus 0.76 ± 0.16 mm, respectively) while patients carrying the CT genotype had intermediate cIMT values (mean ± SD: 0.73 ± 0.17 mm). Moreover, RA patients carrying the mutant allele T of the IL33 rs3939286 polymorphism exhibited significantly lower cIMT values than those carrying the wild allele C (mean ± SD: 0.72 ± 0.16 mm versus 0.75 ± 0.18 mm respectively; p = 0.04). The association of both genotype and allele frequencies of IL33 rs3939286 and cIMT levels remained statistically significant after adjustment for sex, age at the time of US study, follow-up and center (p = 0.006 and p = 0.0023, respectively), evidencing that the potential effect conferred by IL33 rs3939286 may be independent of confounder factors. No association with other IL33-IL1RL1 genetic variants was observed.

CONCLUSIONS

In conclusion, our results may suggest a potential protective effect of the IL33 rs3939286 allele T in the risk of subclinical atherosclerosis in patients with RA.

PTGER4 gene variant rs76523431 is a candidate risk factor for radiological joint damage in rheumatoid arthritis patients: a genetic study of six cohorts

Introduction

Prostaglandin E receptor 4 (PTGER4) is implicated in immune regulation and bone metabolism. The aim of this study was to analyze its role in radiological joint damage in rheumatoid arthritis (RA).

Methods

Six independent cohorts of patients with RA of European or North American descent were included, comprising 1789 patients with 5083 sets of X-rays. The Hospital Clínico San Carlos Rheumatoid Arthritis, Princesa Early Arthritis Register Longitudinal study, and Hospital Universitario de La Paz early arthritis (Spain) cohorts were used as discovery cohorts, and the Leiden Early Arthritis Clinic (The Netherlands), Wichita (United States), and National Databank for Rheumatic Diseases (United States and Canada) cohorts as replication cohorts. First, the PTGER4 rs6896969 single-nucleotide polymorphism (SNP) was genotyped using TaqMan assays and available Illumina Immunochip data and studied in the discovery and replication cohorts. Second, the PTGER4 gene and adjacent regions were analyzed using Immunochip genotyping data in the discovery cohorts. On the basis of pooled p values, linkage disequilibrium structure of the region, and location in regions with transcriptional properties, SNPs were selected for replication. The results from discovery, replication, and overall cohorts were pooled using inverse-variance–weighted meta-analysis. Influence of the polymorphisms on the overall radiological damage (constant effect) and on damage progression over time (time-varying effect) was analyzed.

Results

The rs6896969 polymorphism showed a significant association with radiological damage in the constant effect pooled analysis of the discovery cohorts, although no significant association was observed in the replication cohorts or the overall pooled analysis. Regarding the analysis of the PTGER4 region, 976 variants were analyzed in the discovery cohorts. From the constant and time-varying effect analyses, 12 and 20 SNPs, respectively, were selected for replication. Only the rs76523431 variant showed a significant association with radiographic progression in the time-varying effect pooled analysis of the discovery, replication, and overall cohorts. The overall pooled effect size was 1.10 (95 % confidence interval 1.05–1.14, p = 2.10 × 10⁻⁵), meaning that radiographic yearly progression was 10 % greater for each copy of the minor allele.

Conclusions

The PTGER4 gene is a candidate risk factor for radiological progression in RA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0830-z) contains supplementary material, which is available to authorized users.

New insights into the pathogenesis of giant cell arteritis and hopes for the clinic

Giant cell arteritis is a complex immune-mediated disease that involves large blood vessels in individuals older than 50 years. Recent studies have confirmed a strong association of this form of vasculitis with the HLA region, particularly with HLA class II genes. However, other non-HLA loci, such as protein tyrosine phosphatase non-receptor type 22, may also account for the susceptibility to giant cell arteritis. In addition, genetic variants located in genes encoding proinflammatory cytokines seem to influence the phenotypic expression of the disease, including the risk of severe ischemic complications, the presence of polymyalgia rheumatica and the higher incidence of relapses observed in some patients. The identification of putative genetic markers of disease severity could have clear therapeutic implications, as it may allow us to identify patients who are potentially responders to specific treatments.

Identification of the PTPN22 functional variant R620W as susceptibility genetic factor for giant cell arteritis

OBJECTIVE

To analyse the role of the PTPN22 and CSK genes, previously associated with autoimmunity, in the predisposition and clinical phenotypes of giant cell arteritis (GCA).

METHODS

Our study population was composed of 911 patients diagnosed with biopsy-proven GCA and 8136 unaffected controls from a Spanish discovery cohort and three additional independent replication cohorts from Germany, Norway and the UK. Two functional PTPN22 polymorphisms (rs2476601/R620W and rs33996649/R263Q) and two variants of the CSK gene (rs1378942 and rs34933034) were genotyped using predesigned TaqMan assays.

RESULTS

The analysis of the discovery cohort provided evidence of association of PTPN22 rs2476601/R620W with GCA (PFDR=1.06E-04, OR=1.62, CI 95% 1.29 to 2.04). The association did not appear to follow a specific GCA subphenotype. No statistically significant differences between allele frequencies for the other PTPN22 and CSK genetic variants were evident either in the case/control or in stratified case analysis. To confirm the detected PTPN22 association, three replication cohorts were genotyped, and a consistent association between the PTPN22 rs2476601/R620W variant and GCA was evident in the overall meta-analysis (PMH=2.00E-06, OR=1.51, CI 95% 1.28 to 1.79).

CONCLUSIONS

Our results suggest that the PTPN22 polymorphism rs2476601/R620W plays an important role in the genetic risk to GCA.

New insight on the Xq28 association with systemic sclerosis

OBJECTIVE

To evaluate whether the systemic sclerosis (SSc)-associated IRAK1 non-synonymous single-nucleotide polymorphism rs1059702 is responsible for the Xq28 association with SSc or whether there are other independent signals in the nearby methyl-CpG-binding protein 2 gene (MECP2).

METHODS

We analysed a total of 3065 women with SSc and 2630 unaffected controls from five independent Caucasian cohorts. Four tag single-nucleotide polymorphisms of MECP2 (rs3027935, rs17435, rs5987201 and rs5945175) and the IRAK1 variant rs1059702 were genotyped using TaqMan predesigned assays. A meta-analysis including all cohorts was performed to test the overall effect of these Xq28 polymorphisms on SSc.

RESULTS

IRAK1 rs1059702 and MECP2 rs17435 were associated specifically with diffuse cutaneous SSc (PFDR=4.12×10(-3), OR=1.27, 95% CI 1.09 to 1.47, and PFDR=5.26×10(-4), OR=1.30, 95% CI 1.14 to 1.48, respectively), but conditional logistic regression analysis showed that the association of IRAK1 rs1059702 with this subtype was explained by that of MECP2 rs17435. On the other hand, IRAK1 rs1059702 was consistently associated with presence of pulmonary fibrosis (PF), because statistical significance was observed when comparing SSc patients PF+ versus controls (PFDR=0.039, OR=1.30, 95% CI 1.07 to 1.58) and SSc patients PF+ versus SSc patients PF- (p=0.025, OR=1.26, 95% CI 1.03 to 1.55).

CONCLUSIONS

Our data clearly suggest the existence of two independent signals within the Xq28 region, one located in IRAK1 related to PF and another in MECP2 related to diffuse cutaneous SSc, indicating that both genes may have an impact on the clinical outcome of the disease.

A case-control study suggests that the CCR6 locus is not involved in the susceptibility to giant cell arteritis

OBJECTIVES

Polymorphisms of the CC chemokine receptor 6 (CCR6) gene have been recently reported to be associated with a number of autoimmune diseases. We aimed to investigate the possible influence of CCR6 rs3093024 gene variant in the susceptibility to and clinical expression of GCA.

METHODS

The CCR6 polymorphism rs3093024 was genotyped in a total of 463 Spanish patients diagnosed with biopsy-proven GCA and 920 healthy controls using a TaqMan® allelic discrimination assay. PLINK software was used for the statistical analyses.

RESULTS

No significant association between this CCR6 variant and GCA was observed (p=0.42, OR=0.94, CI95% 0.79-1.10). Similarly, when patients were stratified according to the specific clinical features of GCA such as polymyalgia rheumatica, visual ischaemic manifestations or irreversible occlusive disease, no statistical significant difference was detected either between the case subgroups and the control set or between GCA patients with and without the specific features of the disease.

CONCLUSIONS

Our results suggest that the CCR6 rs3093024 polymorphism may not play a relevant role in the GCA pathophysiology.

Evaluation of a shared autoimmune disease-associated polymorphism of TRAF6 in systemic sclerosis and giant cell arteritis

OBJECTIVE

We evaluated whether a single-nucleotide polymorphism (SNP) of the TRAF6 gene previously associated with systemic lupus erythematosus and rheumatoid arthritis may be a common risk factor for systemic sclerosis (SSc) and giant cell arteritis (GCA).

METHODS

A total of 1185 patients with SSc, 479 patients with biopsy-proven GCA, and 1442 unrelated healthy controls of white Spanish origin were genotyped for the rs540386 variant using a specifically designed TaqMan(©) allele discrimination assay.

RESULTS

No significant associations of this SNP with global SSc or GCA were found. This was also the case when the potential associations of the TRAF6 polymorphism with the main clinical phenotypes of the 2 diseases (e.g., limited cutaneous and diffuse cutaneous SSc, or presence of polymyalgia rheumatica and visual ischemic manifestations in GCA) were assessed.

CONCLUSION

Our data do not support a role of the rs540386 TRAF6 variant as a key component of the genetic network underlying SSc and GCA.

Identification of CSK as a systemic sclerosis genetic risk factor through Genome Wide Association Study follow-up

Systemic sclerosis (SSc) is complex autoimmune disease affecting the connective tissue; influenced by genetic and environmental components. Recently, we performed the first successful genome-wide association study (GWAS) of SSc. Here, we perform a large replication study to better dissect the genetic component of SSc. We selected 768 polymorphisms from the previous GWAS and genotyped them in seven replication cohorts from Europe. Overall significance was calculated for replicated significant SNPs by meta-analysis of the replication cohorts and replication-GWAS cohorts (3237 cases and 6097 controls). Six SNPs in regions not previously associated with SSc were selected for validation in another five independent cohorts, up to a total of 5270 SSc patients and 8326 controls. We found evidence for replication and overall genome-wide significance for one novel SSc genetic risk locus: CSK [P-value = 5.04 × 10(-12), odds ratio (OR) = 1.20]. Additionally, we found suggestive association in the loci PSD3 (P-value = 3.18 × 10(-7), OR = 1.36) and NFKB1 (P-value = 1.03 × 10(-6), OR = 1.14). Additionally, we strengthened the evidence for previously confirmed associations. This study significantly increases the number of known putative genetic risk factors for SSc, including the genes CSK, PSD3 and NFKB1, and further confirms six previously described ones.

Autoimmune disease-associated CD226 gene variants are not involved in giant cell arteritis susceptibility in the Spanish population

OBJECTIVES

CD226 genetic variants have been associated with a number of autoimmune diseases. The aim of this study was to investigate the potential implication of the CD226 loci in the susceptibility to and main clinical manifestations of giant cell arteritis (GCA).

METHODS

A Spanish Caucasian cohort of 455 patients diagnosed with biopsy-proven GCA and 1414 healthy controls were included in the study. Three CD226 polymorphisms, rs727088, rs34794968 and rs763361, were genotyped using the TaqMan® allelic discrimination technology. PLINK software was used for the statistical analyses.

RESULTS

No significant association between the CD226 polymorphisms and susceptibility to GCA was found (rs727088: p=0.92, OR=1.01, CI 95% 0.86-1.18; rs34794968: p=0.61, OR=1.04, CI 95% 0.89-1.22; rs763361: p=0.88, OR=0.99, CI 95% 0.84-1.16). Similarly, when patients were stratified according to the specific clinical features of GCA such as polymyalgia rheumatica, visual ischaemic manifestations or irreversible occlusive disease, no association was observed either between the case subgroups and the control set or between GCA patients with and without the specific features of the disease. Furthermore, the haplotype analysis revealed no significant association with the clinical manifestations of the disease.

CONCLUSIONS

Our results show that the three CD226 polymorphisms analysed do not play a relevant role in the susceptibility to GCA and clinical manifestations of this vasculitis.

Novel identification of the IRF7 region as an anticentromere autoantibody propensity locus in systemic sclerosis

OBJECTIVE

Systemic sclerosis (SSc) and systemic lupus erythematosus (SLE) are related chronic autoimmune diseases of complex aetiology in which the interferon (IFN) pathway plays a key role. Recent studies have reported an association between IRF7 and SLE which confers a risk to autoantibody production. A study was undertaken to investigate whether the IRF7 genomic region is also involved in susceptibility to SSc and the main clinical features.

METHODS

Two case-control sets of Caucasian origin from the USA and Spain, comprising a total of 2316 cases of SSc and 2347 healthy controls, were included in the study. Five single nucleotide polymorphisms (SNPs) in the PHRF1-IRF7-CDHR5 locus were genotyped using TaqMan allelic discrimination technology. A meta-analysis was performed to test the overall effect of these genetic variants on SSc.

RESULTS

Four out of five analysed SNPs were significantly associated with the presence of anticentromere autoantibodies (ACA) in the patients with SSc in the combined analysis (rs1131665: p(FDR)=6.14 × 10(-4), OR=0.78; rs4963128: p(FDR)=6.14 × 10(-4), OR=0.79; rs702966: p(FDR)=3.83 × 10(-3), OR=0.82; and rs2246614: p(FDR)=3.83 × 10(-3), OR=0.83). Significant p values were also obtained when the disease was tested globally; however, the statistical significance was lost when the ACA-positive patients were excluded from the study, suggesting that these associations rely on ACA positivity. Conditional logistic regression and allelic combination analyses suggested that the functional IRF7 SNP rs1131665 is the most likely causal variant.

CONCLUSIONS

The results show that variation in the IRF7 genomic region is associated with the presence of ACA in patients with SSc, supporting other evidence that this locus represents a common risk factor for autoantibody production in autoimmune diseases.

Development of the cornea of true moles (Talpidae): morphogenesis and expression of PAX6 and cytokeratins

Corneal development and structure were studied in the Iberian mole Talpa occidentalis, which has permanently closed eyelids, and the European mole Talpa europaea, in which the eyes are open. The vertebrate cornea typically maintains a three-layered structure - a stratified epithelium with protective and sensory function, an avascular, hypocellular, collagenous stroma, and an endothelium with both barrier and transport functions that regulates corneal hydration, hence maintaining transparency. Compared to mouse, both mole species had significant corneal specializations, but the Iberian mole had the most divergent phenotype, with no endothelium and a flattened monolayer epithelium. Nevertheless, normal epithelial cell junctions were observed and corneal transparency was maintained. Corneas of European moles have a dysmorphic phenotype that recapitulates the human disorder keratoconus for which no mouse model exists. Mole corneas are vascularized - a situation only previously observed in the manatee Trichechus- and have non-radial patterns of corneal innervation indicative of failure of corneal epithelial cell migration. The transcription factor Pax6 is required for corneal epithelial differentiation in mice, but was found to be dispensable in moles, which had mosaic patterns of PAX6 localization uniquely restricted, in European moles, to the apical epithelial cells. The apparently stalled or abnormal differentiation of corneas in adult moles is supported by their superficial similarity to the corneas of embryonic or neonatal mice, and their abnormal expression of cytokeratin-12 and cytokeratin-5. European moles seem to have maintained some barrier/protective function in their corneas. However, Iberian moles show a more significant corneal regression likely related to the permanent eyelid fusion. In this mole species, adaptation to the arid, harder, Southern European soils could have favoured the transfer of these functions to the permanently sealed eyelids.

Retinal development and function in a 'blind' mole

Animals adapted to dark ecotopes may experience selective pressure for retinal reduction. No previous studies have explicitly addressed the molecular basis of retinal development in any fossorial mammal. We studied retinal development and function in the Iberian mole Talpa occidentalis, which was presumed to be blind because of its permanently closed eyes. Prenatal retina development was relatively normal, with specification of all cell types and evidence of dorsoventral regionalization. Severe developmental defects occurred after birth, subsequent to lens abnormalities. 'Blind' Iberian moles had rods, cones and rod nuclear ultrastructure typical of diurnal mammals. DiI staining revealed only contralateral projections through the optic chiasm. Y-maze experiments demonstrated that moles retain a photoavoidance response. Over-representation of melanopsin-positive retinal ganglion cells that mediate photoperiodicity was observed. Hence, molecular pathways of eye development in Iberian moles retain the adaptive function of rod/cone primary vision and photoperiodicity, with no evidence that moles are likely to completely lose their eyes on an evolutionary time scale.

The molecular basis of defective lens development in the Iberian mole

Background

Fossorial mammals face natural selection pressures that differ from those acting on surface dwelling animals, and these may lead to reduced visual system development. We have studied eye development in a species of true mole, the Iberian mole Talpa occidentalis, and present the molecular basis of abnormal lens development. This is the first embryological developmental study of the eyes of any fossorial mammal at the molecular level.

Results

Lens fibre differentiation is not completed in the Iberian mole. Although eye development starts normally (similar to other model species), defects are seen after closure of the lens vesicle. PAX6 is not down-regulated in developing lens fibre nuclei, as it is in other species, and there is ectopic expression of FOXE3, a putative downstream effector of PAX6, in some, but not all lens fibres. FOXE3-positive lens fibres continue to proliferate within the posterior compartment of the embryonic lens, but unlike in the mouse, no proliferation was detected anywhere in the postnatal mole lens. The undifferentiated status of the anterior epithelial cells was compromised, and most of them undergo apoptosis. Furthermore, β-crystallin and PROX1 expression patterns are abnormal and our data suggest that genes encoding β-crystallins are not directly regulated by PAX6, c-MAF and PROX1 in the Iberian mole, as they are in other model vertebrates.

Conclusion

In other model vertebrates, genetic pathways controlling lens development robustly compartmentalise the lens into a simple, undifferentiated, proliferative anterior epithelium, and quiescent, anuclear, terminally differentiated posterior lens fibres. These pathways are not as robust in the mole, and lead to loss of the anterior epithelial phenotype and only partial differentiation of the lens fibres, which continue to express 'epithelial' genes. Paradigms of genetic regulatory networks developed in other vertebrates appear not to hold true for the Iberian mole.

Peritubular myoid cells are not the migrating population required for testis cord formation in the XY gonad

Cell migration is one of the earliest events required for development of the testis. Migration occurs only in XY gonads downstream of Sry expression and is required for the subsequent epithelialization of testis cords. Using organ culture experiments and tissue recombination, we and others speculated that peritubular myoid (PTM) cells were among the migratory cells and were likely the cell type required for cord formation. However, because no unique marker was found for PTM cells, their positive identification during or after migration remained unclear. alpha-Smooth Muscle Actin (alphaSma; approved gene symbol Acta2), a classic marker of adult PTM cells,is expressed broadly in testis interstitial cells at E12.5, and becomes highly enriched in PTM cells by E15.5-16.5. We used a novel transgenic line expressingEYFP under the control of an alphaSma promoter to determine whether alphaSma-EYFP positive cellsmigrate into the gonad. Surprisingly, mesonephroi expressing alphaSma-EYFP do not contribute any EYFP positive cells to XY gonads when used as donors in recombination cultures. These results indicate that alphaSma-EYFP cells do not migrate into the gonad during the critical window of sex determination and cannot be the migrating cell type required for testis cord formation. Our results suggest that PTM cells, and most other interstitial lineages, with the exception of endothelial cells, are induced within the gonad. These experiments suggest that endothelial cells are the migrating cell type required for epithelialization of testis cords.

Histone H3 lysine 9 acetylation pattern suggests that X and B chromosomes are silenced during entire male meiosis in a grasshopper

The facultative heterochromatic X chromosome in leptotene spermatocytes of the grasshopper Eyprepocnemis plorans showed marked hypoacetylation for lysine 9 in the H3 histone (H3-K9) with no sign of histone H2AX phosphorylation. Since H3-K9 hypoacetylation precedes the meiotic appearance of phosphorylated H2AX (gamma-H2AX), which marks the beginning of recombinational DNA double-strand breaks (DSBs), it seems that meiotic sex-chromosome inactivation (MSCI) in this grasshopper occurs prior to the beginning of recombination and hence synapsis (which in this species begins later than recombination). In addition, all constitutively heterochromatic chromosome regions harbouring a 180-bp tandem-repeat DNA and rDNA (B chromosomes and pericentromeric regions of A chromosomes) were H3-K9 hypoacetylated at early leptotene even though they will synapse at subsequent stages. This also suggests that meiotic silencing in this grasshopper might be independent of synapsis. The H3-K9 hypoacetylated state of facultative and constitutive heterochromatin persisted during subsequent meiotic stages and was even apparent in round spermatids. Finally, the fact that B chromosomes are differentially hypoacetylated in testis and embryo interphase cells suggests that they might be silenced early in development and remain this way for most (or all) life-cycle stages.

Histone H2AX phosphorylation is associated with most meiotic events in grasshopper

It is widely accepted that the H2AX histone in its phosphorylated form (gamma-H2AX) is related to the repair of DNA double-strand breaks (DSBs). In several organisms, gamma-H2AX presence has been demonstrated in meiotic processes such as recombination and sex chromosome inactivation during prophase I (from leptotene to pachytene). To test whether gamma-H2AX is present beyond pachytene, we have analysed the complete sequence of changes in H2AX phosphorylation during meiosis in grasshopper, a model organism for meiotic studies at the cytological level. We show the presence of phosphorylated H2AX during most of meiosis, with the exception only of diplotene and the end of each meiotic division. During the first meiotic division, gamma-H2AX is associated with i) recombination, as deduced from its presence in leptotene-zygotene over all chromosome length, ii) X chromosome inactivation, since at pachytene gamma-H2AX is present in the X chromosome only, and iii) chromosome segregation, as deduced from gamma-H2AX presence in centromere regions at first metaphase-anaphase. During second meiotic division, gamma-H2AX was very abundant at most chromosome lengths from metaphase to telophase, suggesting its possible association with the maintenance of chromosome condensation and segregation.