Search

In 2013–2014, an outbreak involving 14 patients infected by an extensively drug-resistant strain of Pseudomonas aeruginosa was detected in a hospital in Madrid, Spain. Our objective was to evaluate an alternative strategy for investigating the outbreak in depth by means of molecular and genomic approaches. Pulsed-field gel electrophoresis (PFGE) was applied as a first-line approach, followed by a more refined whole genome sequencing analysis. Single nucleotide polymorphisms identified by whole genome sequencing were used to design a specific polymerase chain reaction (PCR) for screening unsuspected cases infected by the outbreak strain. Whole genome sequencing alerted us to the existence of greater genetic diversity than was initially assumed, splitting the PFGE-associated outbreak isolates into 4 groups, 2 of which represented coincidental transmission unrelated to the outbreak. A multiplex allele-specific PCR targeting outbreak-specific single nucleotide polymorphisms was applied to 290 isolates, which allowed us to identify 25 additional cases related to the outbreak during 2011–2017. Whole genome sequencing coupled with an outbreak-strain-specific PCR enabled us to markedly redefine the initial picture of the outbreak by 1) ruling out initially suspected cases, 2) defining likely independent coincidental transmission events, 3) predating the starting point of the outbreak, 4) capturing new unsuspected cases, and 5) revealing that the outbreak was still active.

The SARS-CoV-2 variant of concern (VOC) Delta (B.617.2 lineage) displaced the predominant VOC Alpha (B.1.1.7 lineage) in the United Kingdom. In Madrid, recent start of the decline of predominant VOC Alpha suggested an equivalent phenomenon. However, 11 different variants, none overrepresented in frequency, occupied progressively over a period of 7 weeks the niche previously dominated by VOC Alpha. Only after these 7 weeks, VOC Delta started to emerge. Viral competition due to the entry of VOC Delta is not the major force driving the start of VOC Alpha decline in Madrid. IMPORTANCE Our data indicate that the dynamics of SARS-CoV-2 VOCs turnover in our setting differ from those proposed for other countries. A systematic genomic analysis, updated on a weekly basis, of representative randomly selected samples of SARS-CoV-2 circulating variants allowed us to define a lapse of 7 weeks between the start of VOC Alpha decline and the final emergence of VOC Delta. During this period, VOC Alpha showed a sustained decline, while 11 VOCs, variants of interest (VOIs), and other identified variants, none overrepresented, occupied the niche left by VOC Alpha. Only after these 7 weeks, emergence of VOC Delta occurred, indicating that viral competition involving VOC Delta was not the exclusive direct driving force behind the starting of VOC Alpha decline.

12 páginas, 2 figuras, 1 tabla., A successful Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variant, B.1.1.7, has recently been reported in the UK, causing global alarm. Most likely, the new variant emerged in a persistently infected patient, justifying a special focus on these cases. Our aim in this study was to explore certain clinical profiles involving severe immunosuppression that may help explain the prolonged persistence of viable viruses. We present three severely immunosuppressed cases (A, B, and C) with a history of lymphoma and prolonged SARS-CoV-2 shedding (2, 4, and 6 months), two of whom finally died. Whole-genome sequencing of 9 and 10 specimens from Cases A and B revealed extensive within-patient acquisition of diversity, 12 and 28 new single nucleotide polymorphisms, respectively, which suggests ongoing SARS-CoV-2 replication. This diversity was not observed for Case C after analysing 5 sequential nasopharyngeal specimens and one plasma specimen, and was only observed in one bronchoaspirate specimen, although viral viability was still considered based on constant low Ct values throughout the disease and recovery of the virus in cell cultures. The acquired viral diversity in Cases A and B followed different dynamics. For Case A, new single nucleotide polymorphisms were quickly fixed (13-15 days) after emerging as minority variants, while for Case B, higher diversity was observed at a slower emergence: fixation pace (1-2 months). Slower SARS-CoV-2 evolutionary pace was observed for Case A following the administration of hyperimmune plasma. This work adds knowledge on SARS-CoV-2 prolonged shedding in severely immunocompromised patients and demonstrates viral viability, noteworthy acquired intra-patient diversity, and different SARS-CoV-2 evolutionary dynamics in persistent cases., This research was funded by Instituto de Salud Carlos III (Ref COV20/00140: SeqCOVID—Consorcio para la epidemiología genómica de SARS-CoV-2 en España) and by Consejo Superior de Investigaciones Científicas (CSIC) (PTI Salud Global). Miguel Servet Contracts CP15/00075 and CPII20/00001) to L.P.-L.

páginas, 4 figuras, 3 tablas, Background: Growing international migration has increased the complexity of tuberculosis transmission patterns. Italy's decision to close its borders in 2018 made of Spain the new European porte entrée for migration from the Horn of Africa (HA). In one of the first rescues of migrants from this region at the end of 2018, tuberculosis was diagnosed in eight subjects, mainly unaccompanied minors. Methods: Mycobacterium tuberculosis isolates from these recently arrived migrants were analysed by Mycobacterial Interspersed Repetitive-Unit/Variable-Number of Tandem Repeat (MIRU-VNTR) and subsequent whole genome sequencing (WGS) analysis. Data were compared with those from collections from other European countries receiving migrants from the HA and a strain-specific PCR was applied for a fast searching of common strains. Infections in a cellular model were performed to assess strain virulence. Results: MIRU-VNTR analysis allowed identifying an epidemiological cluster involving three of the eight cases from Somalia (0 single-nucleotide polymorphisms between isolates, HA cluster). Following detailed interviews revealed that two of these cases had shared the same migratory route in most of the trip and had spent a long time at a detention camp in Libya. To confirm potential en route transmission for the three cases, we searched the same strain in collections from other European countries receiving migrants from the HA. MIRU-VNTR, WGS and a strain-specific PCR for the HA strain were applied. The same strain was identified in 12 cases from Eritrea diagnosed soon after their arrival in 2018 to the Netherlands, Belgium and Italy. Intracellular replication rate of the strain did not reveal abnormal virulence. Conclusions: Our study suggests a potential en route transmission of a pan-susceptible strain, which caused at least 15 tuberculosis cases in Somalian and Eritrean migrants diagnosed in four different European countries., This work was supported by ERANet-LAC [ELAC2015/T08-0664] and Instituto de Salud Carlos III [AC16/00057, FIS15/01554, FIS13/01207, CP15/00075, PI19/00331] and cofounded by the European Regional Development Funds from the European Commission: ‘A way of making Europe’ and by Junta de Andalucia PI0488-2017.

INTRODUCTION: A newly identified SARS-CoV-2 variant, VOC202012/01 originating lineage B.1.1.7, recently emerged in the United Kingdom. The rapid spread in the UK of this new variant has caused other countries to be vigilant. MATERIAL AND METHODS: We based our initial screening of B.1.1.7 on the dropout of the S gene signal in the TaqPath assay, caused by the 69/70 deletion. Subsequently, we confirmed the B.1.1.7 candidates by whole genome sequencing. RESULTS: We describe the first three imported cases of this variant from London to Madrid, subsequent post-arrival household transmission to three relatives, and the two first cases without epidemiological links to UK. One case required hospitalization. In all cases, drop-out of gene S was correctly associated to the B.1.1.7 variant, as all the corresponding sequences carried the 17 lineage-marker mutations. CONCLUSION: The first identifications of the SARS-CoV-2 B.1.1.7 variant in Spain indicate the role of independent introductions from the UK coexisting with post-arrival transmission in the community, since the early steps of this new variant in our country.

Objectives: The aim of this this study was to describe the presence of different variants of the fosA gene in fosfomycin-resistant Escherichia coli strains in Madrid, Spain. Methods: fos genes were searched for in 55 E. coli strains collected from seven representative hospitals located in Madrid. A phenotypic screening test was performed following the disk diffusion method with sodium phosphonoformate added as described by Nakamura et al. Additionally, a molecular study based on PCR was used to confirm the screening results. Positive strains for fos genes were further subjected to whole-genome sequencing (WGS). Results: Phenotypic screening was positive in 9/55 strains (16.4%), although genotypic detection was positive in only 3 (fosA3, fosA4 and fosA6). Thus, the prevalence of fos genes in Madrid was 5.5% (3/55). WGS data were not available for the fosA6-positive strain. One isolate with fosA3 (ST69) carried a blaCTX-M-55 gene and seven virulence genes (air, eilA, iha, iss, lpfA, sat and senB). The fosA4-positive isolate (ST4038) carried the virulence genes iss, lpfA, iroN and mchF. Both fos genes were located between two IS26 mobile elements of a plasmid. Conclusion: We detected the presence of different variants of plasmid-mediated fosA genes in fosfomycin-resistant E. coli strains in Madrid, Spain. Despite the few reports in Europe, it would be of interest to monitor the spread of these acquired resistance genes.

Control of infections caused by carbapenem-resistant Klebsiella pneumoniae continues to be challenging. The success of this pathogen is favored by its ability to acquire antimicrobial resistance and to spread and persist in both the environment and in humans.The emergence of clinically important clones, such as sequence types 11, 15, 101, and 258, has been reported worldwide. However,the mechanisms promoting the dissemination of such highrisk clones are unknown. Unraveling the factors that play a role in the pathobiology and epidemicity of K. pneumoniae is therefore important for managing infections. To address this issue, we studied a carbapenem-resistant ST-15 K. pneumoniae isolate (Kp3380) that displayed a remarkable adherent phenotype with abundant pilus-like structures. Genome sequencing enabled us to identify a chaperone-usher pili system (Kpi) in Kp3380. Analysis of a large K. pneumoniae population from 32 European countries showed that the Kpi system is associated with the ST-15 clone. Phylogenetic analysis of the operon revealed that Kpi belongs to the littlecharacterized γ2-fimbrial clade. We demonstrate that Kpi contributes positively to the ability of K.pneumoniae to form biofilms and adhere to different host tissues. Moreover, the in vivo intestinal colonizing capacity of the Kpi-defective mutant was significantly reduced, as was its ability to infect Galleria mellonella. The findings provide information about the pathobiology and epidemicity of Kpi+ K. pneumoniae and indicate that the presence of Kpi may explain the success of the ST-15 clone. Disrupting bacterial adherence to the intestinal surface could potentially target gastrointestinal colonization., This research was supported by Projects p-01216A and IJCI-2016-29524 (to A.P.), funded by the Spanish Society of Infectious Diseases and Clinical Microbiology (SEIMC) and the Minestry of Economy and Competetiveness(MINECO), respectively. It was also supported by Projects PI11/01034 (to M.P.), PI14/00059 and PI17/1482 (to M.P. and A.B.), and PI18/00501 (to G.B.),included in the National Plan for Scientific Research, Development and Technological Innovation 2013-2016 and funded by the Instituto de Salud Carlos III (ISCIII) and Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía, Industria y Competitividad, Spanish Network for Research in Infectious Diseases (REIPI RD16/0016/006) cofinanced by European Development Regional Fund “A way to achieve Europe” and operative program Intelligent Growth 2014-2020. Grant BFU2016-77835-R of the MINECO (to A.R.) also supported this research. E.G. was financially supported by the SEIMC project. J.C.V.-U. was financially supported by the PFIS (Contratos Predoctorales de Formación en Investigación en Salud) program (F18/00315);J.A.V. was financially supported by IN607A 2016/22; M.M.-G. was financially supported by a Clara Roy grant (SEIMC); A.B. was financially supported by the Miguel Servet program (ISCIII, Spain); B.K.R.-J. was financially supported by Marie S. Curie Action SaPhaDe project (MSCA-IF-GF-836754); and A.P. was financially supported by the Juan de la Cierva program (MINECO, IJCI-2016-29524).

Carbapenemase-producing (CP) Pseudomonas aeruginosa is rare compared with mutation-driven carbapenem-resistance, but this situation may be changing. A collection of CP P. aeruginosa isolates was characterized in this study. In 2016, 232 unduplicated carbapenem-resistant P. aeruginosa isolates, of which 71 (30.6%) carried carbapenemase genes, were submitted to the Spanish antibiotic reference laboratory and were further analysed by whole-genome sequencing (WGS). Of the 71 CP P. aeruginosa, 39 (54.9%) carried blaVIM-2, 14 (19.7%) blaVIM-1, 8 (11.3%) blaIMP-8, 6 (8.5%) blaVIM-20, 2 (2.8%) blaVIM-2 plus blaKPC-2, one (1.4%) blaIMP-13 and one (1.4%) blaVIM-1 plus blaIMP-18. Four sequence types (ST175, ST244, ST815 and ST155) encompassed 83.1% of the 71 CP P. aeruginosa; ST175 was detected in hospitals from seven provinces. Using core genome multilocus sequence typing (cgMLST), four clusters were detected: Cluster 1 included nine ST815/VIM-2 isolates; Cluster 2 included five ST175/VIM-2 isolates; Cluster 3 included seven ST244 isolates (five VIM-2 and two VIM-2 plus KPC-2); and Cluster 4 included 11 ST175 isolates (seven VIM-2 and four IMP-8). The average number of acquired resistance genes was significantly higher in the blaVIM-1-carying isolates (7.1 ± 0.94) than in the blaVIM-2-carrying isolates (4.5 ± 0.20). CP P. aeruginosa isolates are spreading in Spain, mainly due to the dissemination of high-risk clones such as ST175 and ST244 producing VIM and IMP carbapenemases. Emergence of CP P. aeruginosa is a cause of clinical and epidemiological concern.

The emergence of linezolid-resistant Enterococcus spp. (LRE) due to transferable resistance determinants is a matter of concern. To understand the contribution of the plasmid-encoded optrA and poxtA genes to the emergence of LRE, clinical isolates from different Spanish hospitals submitted to the Spanish Reference Laboratory from 2015-2018 were analysed. Linezolid resistance mechanisms were screened in all isolates by PCR and sequencing. Genetic relatedness of Enterococcus spp. carrying optrA and poxtA was studied by PFGE and MLST. Antimicrobial susceptibility was tested by broth microdilution using EUCAST standards. A total of 97 LRE isolates were studied, in 94 (96.9%) of which at least one resistance determinant was detected; 84/97 isolates (86.6%) presented a single resistance mechanism as follows: 45/84 (53.6%) carried the optrA gene, 38/84 (45.2%) carried the G2576T mutation and 1/84 (1.2%) carried the poxtA gene. In addition, 5/97 isolates (5.2%) carried both optrA and the G2576T mutation and 5/97 (5.2%) carried both optrA and poxtA. The optrA gene was more frequent in Enterococcus faecalis (83.6%) than Enterococcus faecium (11.1%) and was mainly associated with community-acquired urinary tract infections. Carriage of the poxtA gene was more frequent in E. faecium (13.9%) than E. faecalis (1.6%). Among the optrA-positive E. faecalis isolates, two main clusters were detected by PFGE. These two clusters belonged to ST585 and ST480 and were distributed throughout 11 and 6 Spanish provinces, respectively. This is the first description of LRE carrying the poxtA gene in Spain, including the co-existence of optrA and poxtA in five isolates.

Aspergillus fumigatus is considered a common causative agent of human fungal infections. A restricted number of virulence factors have been described, and none of them lead to a differentiation in the virulence level among different strains. Variations in the virulence phenotype depending on the isolate origin, measured as survival percentage in animal infection models, have been previously reported. In this study, we analyzed the whole-genome sequence of A. fumigatus isolates from clinical and environmental origins to determine their virulence genetic content. The sample included four isolates sequenced at the University Medical Center Groningen (UMCG), three clinical (two of them isolated from the same patient) and the experimental strain B5233, and the draft genomes of one reference strain, two environmental and two clinical isolates obtained from a public database. The fungal genomes were screened for the presence of virulence-related genes (VRGs) using an in-house database of 244 genes related to thermotolerance, resistance to immune responses, cell wall formation, nutrient uptake, signaling and regulation, and production of toxins and secondary metabolites and allergens. In addition, we performed a variant calling analysis to compare the isolates sequenced at the UMCG and investigated their genetic relatedness using the TRESP (Tandem Repeats located within Exons of Surface Protein coding genes) genotyping method. We neither observed a difference in the virulence genetic content between the clinical isolates causing an invasive infection and a colonizing clinical isolate nor between isolates from the clinical and environmental origin. The four novel A. fumigatus sequences had a different TRESP genotype and a total number of genetic variants ranging from 48,590 to 68,352. In addition, a comparative genomics analysis showed the presence of single nucleotide polymorphisms in VRGs and repetitive genetic elements located next to VRG groups, which could influence the regulation of these genes. In conclusion, our genomic analysis revealed a high genetic diversity between environmental and clinical A. fumigatus isolates, as well as between clinical isolates from the same patient, indicating an infection with a mixed-population in the latter case. However, all isolates had a similar virulence genetic content, demonstrating their pathogenic potential at least at the genomic level. This work was partly supported by the INTERREG VA (202085) funded project EurHealth-1Health, part of a Dutch-German cross-border network supported by the European Commission, the Dutch Ministry of Health, Welfare and Sport (VWS), the Ministry of Economy, Innovation, Digitalisation and Energy of the German Federal State of North Rhine-Westphalia, and the German Federal State of Lower Saxony. FP-B was supported by the Erasmus Mundus Joint Master Degree (EMJMD) scholarship of the Erasmus + EU-Programme awarded under the International Master in Innovative Medicine (IMIM) programme. Sí

Aspergillus fumigatus is an opportunistic airborne pathogen and one of the most common causative agents of human fungal infections. A restricted number of virulence factors have been described but none of them lead to a differentiation of the virulence level among different strains. In this study, we analyzed the whole-genome sequence of a set of A. fumigatus isolates from clinical and environmental origin to compare their genomes and to determine their virulence profiles. For this purpose, a database containing 244 genes known to be associated with virulence was built. The genes were classified according to their biological function into factors involved in thermotolerance, resistance to immune responses, cell wall structure, toxins and secondary metabolites, allergens, nutrient uptake and signaling and regulation. No difference in virulence profiles was found between clinical isolates causing an infection and a colonizing clinical isolate, nor between isolates from clinical and environmental origin. We observed the presence of genetic repetitive elements located next to virulence related gene groups, which could potentially influence their regulation. In conclusion, our genomic analysis reveals that A. fumigatus, independently of their source of isolation, are potentially pathogenic at the genomic level, which may lead to fatal infections in vulnerable patients. However, other determinants such as genetic variations in virulence related genes and host-pathogen interactions most likely influence A. fumigatus pathogenicity and further studies should be performed.ImportanceAspergillus spp. infections are among the most clinically relevant fungal infections also presenting treatment difficulties due to increasing antifungal resistance. The lack of key virulence factors and a broad genomic diversity complicates the development of targeted diagnosis and novel treatment strategies. A widely spread variability in virulence has been reported for experimental, clinical and environmental isolates. Here we provide supporting evidence that members of this species are fully capable of establishing an infection in immunosuppressed hosts according to their virulence content at the genomic level. Due to the possible clinical complications, studies are urgently required linking strain’s virulent phenotype with the genotype to better understand the virulence activation of this important fungal pathogen.

Objectives NDM carbapenemases have spread worldwide. However, little information exists about the impact of NDM-producing Enterobacteriaceae in Spain. By WGS, we sought to elucidate the population structure of NDM-like-producing Klebsiella pneumoniae and Escherichia coli in Spain and to determine the plasmids harbouring blaNDM-like genes. Methods High-resolution SNP typing, core-genome MLST and plasmid reconstruction (PlasmidID) were performed on 59 NDM-like-producing K. pneumoniae and 8 NDM-like-producing E. coli isolated over an 8 year period in Spain. Results Five major epidemic clones of NDM-producing K. pneumoniae caused five important nationwide outbreaks: ST437/NDM-7, ST437/NDM-1, ST147/NDM-1, ST11/NDM-1 and ST101/NDM-1; in contrast, the spread of NDM-producing E. coli was polyclonal. Three blaNDM types were identified: blaNDM-1, 61.2%; blaNDM-7, 32.8%; and blaNDM-5, 6%. Five K. pneumoniae isolates co-produced other carbapenemases (three blaOXA-48 and two blaVIM-1). The average number of acquired resistance genes was higher in K. pneumoniae than in E. coli. The plasmids encoding blaNDM-like genes belonged to IncFII, IncFIB, IncX3, IncR, IncN and IncC types, of which IncF, IncR and IncC were associated with MDR. The genetic surroundings of blaNDM-like genes showed a highly variable region upstream of ISAba125. Conclusions In recent years NDM-producing K. pneumoniae and E. coli have emerged in Spain; the spread of a few high-risk K. pneumoniae clones such as ST437/NDM-7, ST437/NDM-1, ST147/NDM-1, ST11/NDM-1 and ST101/NDM-1 have caused several interregional outbreaks. In contrast, the spread of NDM-producing E. coli has been polyclonal. Plasmid types IncFII, IncFIB, IncX3, IncR, IncN and IncC carried blaNDM, and the same IncX3 plasmid was detected in K. pneumoniae and E. coli.

a b s t r a c t A new Pistacia variety, VIGROS, has been characterized here using nuclear and chloroplastidial DNA sequences, which are suitable tools for genetic improvement and cultivar identification. The nuclear DNA sequences (ribosomal ITS1, 5.8S, ITS2) determined that VIGROS is an interspecific hybrid between Pistacia vera L. and Pistacia atlantica Desf. The chloroplastidial DNA sequences (trnC-D, trnL-F regions) demonstrated that VIGROS bears a haplotype 100% coincident with that of P. vera, which would then be regarded as the female parent, while P. atlantica would be the male one. Due to its vigorous phenotype VIGROS has the potential to be used as a rootstock. This characterization may benefit producers and contribute to improved grafting of pistachios. Also, this paper demonstrates the validity of this set of molecular markers to characterize interspecific Pistacia hybrids.

The bud–rootstock compatibility is decisive in the success of grafting pistachios. Climate conditions and appropriate rootstock choice are the two key factors influencing grafting success. Despite producers’ efforts, fluctuations in these factors result in frequent differences in annual yield. To minimize production losses by increasing the percentage of successful grafting would greatly benefit the introduction of this crop in new areas. Here, we analyze the viability of Pistacia terebinthus L. differentially treated with mycorrhizae or phytohormones and used as rootstock for Pistacia vera buds. Our results, on an experimental plot of 12,905 plants, demonstrate that mycorrhiza-treated plants reached ∼80% positive grafts, while the phytohormone-treated plants and controls had 32.3% and 38.4% success, respectively. The increase in grafting success could be explained by more efficient nutrient uptake in mycorrhiza-treated plants. An analysis of chemical element accumulation and assimilation in leaves reveals that mycorrhiza-treated plants selectively accumulated Ca, Fe, Mg, N, Al, S, Sr, Ti, V, Mn, and Tl, but lacked K with respect to the rest of the plants. Mycorrhiza-treated plants were shorter, but no significant differences were found in trunk diameter and circumference. We propose the use of mycorrhizae to increase bud–rootstock compatibility for this tandem of species.

Two new centromeric satellite DNAs in flatfish (Order Pleuronectiformes) have been characterized. The SacI-family from Hippoglossus hippoglossus, restricted to this species, had a monomeric size of 334 base pair (bp) and was located in most of the centromeres of its karyotype. The PvuII-family, with a monomeric size of 177 bp, was initially isolated from the genome of Solea senegalensis, and fluorescent in situ hybridization (FISH) localized the repeat to centromeres of most of the chromosomes. This family could only be amplified in 2 other species of the genus Solea (Solea solea and Solea lascaris). Molecular features and chromosomal location indicated a possible structural and/or functional role of these sequence repeats. The presence of species-specific satellite-DNA families in the centromeres and their possible role in the speciation processes in this group of fishes is discussed.

This paper represents the first molecular cytogenetic characterization of the strictly dioecious pistachio tree (Pistacia vera L.). The karyotype was characterized by fluorescent in situ hybridization (FISH) with probes for 5S and 45S rDNAs, and the pistachio specific satellite DNAs PIVE-40, and PIVE-180, together with DAPI-staining. PIVE-180 has a monomeric unit of 176-178 bp and high sequence homology between family members; PIVE-40 has a 43 bp consensus monomeric unit, and is most likely arranged in higher order repeats (HORs) of two units. The P. vera genome is highly heterochromatic, and prominent DAPI positive blocks are detected in most chromosomes. Despite the difficulty in classifying chromosomes according to morphology, 10 out of 15 pairs (2n = 30) could be distinguished by their unique banding patterns using a combination of FISH probes. Significantly, the largest pair, designated HC1, is strongly heteropycnotic, shows differential condensation, and has massive enrichment in PIVE-40 repeats. There are two types of HC1 chromosomes (type-I and type-II) with differing PIVE-40 hybridization signal. Only type-I/II heterozygotes and type-I homozygotes individuals were found. We speculate that the differentiation between the two HC1 chromosomes is due to suppression of homologous recombination at meiosis, reinforced by the presence of PIVE-40 HORs and differences in PIVE-40 abundance. This would be compatible with a ZW sex-determination system in the pistachio tree.