A global metagenomic map of urban microbiomes and antimicrobial resistance.

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

The Space Omics and Medical Atlas (SOMA) and international astronaut biobank.

Spaceflight induces molecular, cellular and physiological shifts in astronauts and poses myriad biomedical challenges to the human body, which are becoming increasingly relevant as more humans venture into space1-6. Yet current frameworks for aerospace medicine are nascent and lag far behind advancements in precision medicine on Earth, underscoring the need for rapid development of space medicine databases, tools and protocols. Here we present the Space Omics and Medical Atlas (SOMA), an integrated data and sample repository for clinical, cellular and multi-omic research profiles from a diverse range of missions, including the NASA Twins Study7, JAXA CFE study8,9, SpaceX Inspiration4 crew10-12, Axiom and Polaris. The SOMA resource represents a more than tenfold increase in publicly available human space omics data, with matched samples available from the Cornell Aerospace Medicine Biobank. The Atlas includes extensive molecular and physiological profiles encompassing genomics, epigenomics, transcriptomics, proteomics, metabolomics and microbiome datasets, which reveal some consistent features across missions, including cytokine shifts, telomere elongation and gene expression changes, as well as mission-specific molecular responses and links to orthologous, tissue-specific mouse datasets. Leveraging the datasets, tools and resources in SOMA can help to accelerate precision aerospace medicine, bringing needed health monitoring, risk mitigation and countermeasure data for upcoming lunar, Mars and exploration-class missions.

Implications for driving based on the risk of seizures after ischaemic stroke.

BACKGROUND: In addition to other stroke-related deficits, the risk of seizures may impact driving ability after stroke. METHODS: We analysed data from a multicentre international cohort, including 4452 adults with acute ischaemic stroke and no prior seizures. We calculated the Chance of Occurrence of Seizure in the next Year (COSY) according to the SeLECT2.0 prognostic model. We considered COSY<20% safe for private and <2% for professional driving, aligning with commonly used cut-offs. RESULTS: Seizure risks in the next year were mainly influenced by the baseline risk-stratified according to the SeLECT2.0 score and, to a lesser extent, by the poststroke seizure-free interval (SFI). Those without acute symptomatic seizures (SeLECT2.0 0-6 points) had low COSY (0.7%-11%) immediately after stroke, not requiring an SFI. In stroke survivors with acute symptomatic seizures (SeLECT2.0 3-13 points), COSY after a 3-month SFI ranged from 2% to 92%, showing substantial interindividual variability. Stroke survivors with acute symptomatic status epilepticus (SeLECT2.0 7-13 points) had the highest risk (14%-92%). CONCLUSIONS: Personalised prognostic models, such as SeLECT2.0, may offer better guidance for poststroke driving decisions than generic SFIs. Our findings provide practical tools, including a smartphone-based or web-based application, to assess seizure risks and determine appropriate SFIs for safe driving.

Telomeres and aging: on and off the planet!

Improving human healthspan in our rapidly aging population has never been more imperative. Telomeres, protective "caps" at the ends of linear chromosomes, are essential for maintaining genome stability of eukaryotic genomes. Due to their physical location and the "end-replication problem" first envisioned by Dr. Alexey Olovnikov, telomeres shorten with cell division, the implications of which are remarkably profound. Telomeres are hallmarks and molecular drivers of aging, as well as fundamental integrating components of the cumulative effects of genetic, lifestyle, and environmental factors that erode telomere length over time. Ongoing telomere attrition and the resulting limit to replicative potential imposed by cellular senescence serves a powerful tumor suppressor function, and also underlies aging and a spectrum of age-related degenerative pathologies, including reduced fertility, dementias, cardiovascular disease and cancer. However, very little data exists regarding the extraordinary stressors and exposures associated with long-duration space exploration and eventual habitation of other planets, nor how such missions will influence telomeres, reproduction, health, disease risk, and aging. Here, we briefly review our current understanding, which has advanced significantly in recent years as a result of the NASA Twins Study, the most comprehensive evaluation of human health effects associated with spaceflight ever conducted. Thus, the Twins Study is at the forefront of personalized space medicine approaches for astronauts and sets the stage for subsequent missions. We also extrapolate from current understanding to future missions, highlighting potential biological and biochemical strategies that may enable human survival, and consider the prospect of longevity in the extreme environment of space.

Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats.

BACKGROUND: Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated. RESULTS: Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh. CONCLUSIONS: The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.

Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations.

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

The landscape of long noncoding RNA expression in the goat brain.

The brain regulates multiple metabolic processes, such as food intake, energy expenditure, insulin secretion, hepatic glucose production, and glucose and fatty acid metabolism in adipose tissue, which are fundamental for the maintenance of energy and glucose homeostasis during lactation and pregnancy. In addition, brain expression has a fundamental impact on the development of maternal behavior. Although brain functions are partly regulated by long noncoding RNAs (lncRNAs), their expression profiles have not been characterized in depth in any ruminant species. We have sequenced the transcriptome of 12 brain tissues from 3 goats that were 1 mo pregnant and 4 nonpregnant goats to investigate their lncRNA expression patterns. Between 4,363 (adenohypophysis) and 4,604 (olfactory bulb) lncRNAs were expressed in brain tissues, leading us to establish a set of 794 already annotated lncRNAs and 5,098 novel lncRNA candidates. The detected lncRNAs shared features with those of other mammals, and tissue-specific lncRNAs were enriched in brain development-related terms. Differential expression analyses between goats that were 1 mo pregnant and nonpregnant goats showed that the lncRNA expression profiles of certain brain regions experience substantial changes associated with early pregnancy (238 lncRNAs are differentially expressed in the olfactory bulb), but others do not. Enrichment analysis showed that differentially expressed lncRNAs from the olfactory bulb are co-expressed with genes previously linked to behavioral changes related to pregnancy. These findings provide a first characterization of the landscape of lncRNA expression in the goat brain and provides valuable clues to understand the molecular events triggered by early pregnancy in the central nervous system.

Identification of nonsense variants in the genomes of 15 Murciano-Granadina bucks and analysis of their segregation in parent-offspring trios.

Nonsense variants can inactivate gene function by causing the synthesis of truncated proteins or by inducing nonsense mediated decay of messenger RNAs. The occurrence of such variants in the genomes of livestock species is modulated by multiple demographic and selective factors. Even though nonsense variants can have causal effects on embryo lethality, abortions, and disease, their genomic distribution and segregation in domestic goats have not been characterized in depth yet. In this work, we have sequenced the genomes of 15 Murciano-Granadina bucks with an average coverage of 32.92 × ± 1.45 × . Bioinformatic analysis revealed 947 nonsense variants consistently detected with SnpEff and Ensembl-VEP. These variants were especially abundant in the 3'end of the protein-coding regions. Genes related to olfactory perception, ATPase activity coupled to transmembrane movement of substances, defense to virus, hormonal response, and sensory perception of taste were particularly enriched in nonsense variants. Seventeen nonsense variants expected to have harmful effects on fitness were genotyped in parent-offspring trios. We observed that several nonsense variants predicted to be lethal based on mouse knockout data did not have such effect, a finding that could be explained by the existence of multiple mechanisms counteracting lethality. These findings demonstrate that predicting the effects of putative nonsense variants on fitness is extremely challenging. As a matter of fact, such a goal could only be achieved by generating a high quality telomere-to-telomere goat reference genome combined with carefully curated annotation and functional testing of promising candidate variants.

Effectiveness of a Second Dose of an mRNA Vaccine Against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Omicron Infection in Individuals Previously Infected by Other Variants.

BACKGROUND: Single-dose vaccination was widely recommended in the pre-Omicron era for persons with previous SARS-CoV-2 infection. The effectiveness of a second vaccine dose in this group in the Omicron era is unknown. METHODS: We linked nationwide population registries in Spain to identify community-dwelling individuals aged 18-64, with a positive SARS-CoV-2 test before single-dose mRNA vaccination (mRNA-1273 or BNT162b2). Every day between 3 January and 6 February 2022 we matched 1:1 individuals receiving a second mRNA vaccine dose and controls on sex, age, province, first dose type and time, month of primary infection, and number of previous tests. We then estimated Kaplan-Meier risks of confirmed SARS-CoV-2 reinfection. We performed a similar analysis in a Delta-dominant period, between 19 July and 30 November 2021. RESULTS: In the Omicron period, estimated effectiveness (95% CI) of a second dose was 62.2% (58.2-66.4%) 7-34 days after administration, similar across groups defined by age, sex, type of first vaccine, and time since the first dose. Estimated effectiveness was 65.4% (61.1-69.9%) for mRNA-1273 and 52.0% (41.8-63.1%) for BNT162b2. Estimated effectiveness was 78.5% (67.4-89.9%), 66.1% (54.9-77.5%), and 60.2% (55.5-64.8%) when primary infection had occurred in the Delta, Alpha, and pre-Alpha periods, respectively. In the Delta period, the estimated effectiveness of a second dose was 8.8% (-55.3% to 81.1%). CONCLUSIONS: Our results suggest that, over 1 month after administration, a second dose of mRNA vaccine increases protection against SARS-CoV-2 reinfection with the Omicron variant among individuals with single-dose vaccination and previously infected with another variant.

Near real-time surveillance of the SARS-CoV-2 epidemic with incomplete data.

When responding to infectious disease outbreaks, rapid and accurate estimation of the epidemic trajectory is critical. However, two common data collection problems affect the reliability of the epidemiological data in real time: missing information on the time of first symptoms, and retrospective revision of historical information, including right censoring. Here, we propose an approach to construct epidemic curves in near real time that addresses these two challenges by 1) imputation of dates of symptom onset for reported cases using a dynamically-estimated "backward" reporting delay conditional distribution, and 2) adjustment for right censoring using the NobBS software package to nowcast cases by date of symptom onset. This process allows us to obtain an approximation of the time-varying reproduction number (Rt) in real time. We apply this approach to characterize the early SARS-CoV-2 outbreak in two Spanish regions between March and April 2020. We evaluate how these real-time estimates compare with more complete epidemiological data that became available later. We explore the impact of the different assumptions on the estimates, and compare our estimates with those obtained from commonly used surveillance approaches. Our framework can help improve accuracy, quantify uncertainty, and evaluate frequently unstated assumptions when recovering the epidemic curves from limited data obtained from public health systems in other locations.

Estimates of mpox effective reproduction number in Spain, April-August 2022.

We analysed the transmission of the human mpox virus in Spain by estimating the effective reproduction number of the disease from official surveillance data. Our computations show that this decreased steadily after an initial burst phase, dropping below 1 on July 12, and thus the outbreak was expected to reduce in the following weeks. Differences in trends were found across geographical regions of the country and across MSM and heterosexual populations.

Univariate and multivariate genome-wide association studies for hematological traits in Murciano-Granadina goats.

Hematological traits are important indicators of health status, and they are frequently used as criteria for clinical diagnosis. In humans, the genomic architecture of blood traits has been investigated in depth and thousands of associations with genetic variants have been found. In contrast, the association between marker genotypes and the variation of hematological traits has not been investigated in goats yet. Herewith, we have recorded 12 hematological parameters in 882 Murciano-Granadina goats that were also genotyped with the Goat SNP50 BeadChip (Illumina). Performance of a univariate genome-wide association study (GWAS) made it possible to detect one genomic region on goat chromosome (CHI) 21 (19.2-19.5 Mb) associated, at the genome-wide level of significance, with 4 red blood cell traits. The three markers displaying the highest significances were rs268272996 (CHI21: 19225290 bp), rs268273004 (CHI21: 19565629 bp) and rs268239059 (CHI13: 9615190 bp). Consistently, a multivariate GWAS indicated that the rs268273004 marker on chromosome 21 is associated with seven blood cell traits. Interestingly, this marker maps close to the FA Complementation Group I (FANCI) gene (CHI21: 20021947-20077025 bp), which is functionally related to Fanconi anemia, a syndrome characterized by bone marrow failure, aplastic anemia, and congenital disorders. We have also uncovered additional chromosome-wide significant associations between genetic markers and erythrocyte and leukocyte traits in the univariate GWAS. These findings evidence that the phenotypic variation of hematological traits in goats is regulated, at least to some extent, by polygenic determinants distributed in multiple chromosomes.

A slightly different polyp: melanoma metastasis diagnosed by colonoscopy.

Melanoma metastases are rare in the colon. Its diagnosis is difficult because they do not usually produce symptoms. They can present through the endoscopic image of a non-pigmented polyp. This is the case of a 56-year-old woman diagnosed with melanoma metastasis through polypectomy of an unusual-looking polyp.

Heterogeneity in how women value risk-stratified breast screening.

PURPOSE: Risk-stratified screening has potential to improve the cost effectiveness of national breast cancer screening programs. This study aimed to inform a socially acceptable and equitable implementation framework by determining what influences a woman's decision to accept a personalized breast cancer risk assessment and what the relative impact of these key determinants is. METHODS: Multicriteria decision analysis was used to elicit the relative weights for 8 criteria that women reported influenced their decision. Preference heterogeneity was explored through cluster analysis. RESULTS: The 2 criteria valued most by the 347 participants related to program access, "Mode of invitation" and "Testing process". Both criteria significantly influenced participation (P < .001). A total of 73% preferred communication by letter/online. Almost all women preferred a multidisease risk assessment with potential for a familial high-risk result. Four preference-based subgroups were identified. Membership to the largest subgroup was predicted by lower educational attainment, and women in this subgroup were concerned with program access. Higher relative perceived breast cancer risk predicted membership to the smallest subgroup that was focused on test parameters, namely "Scope of test" and "Test specificity". CONCLUSION: Overall, Australian women would accept a personalized multidisease risk assessment, but when aligning with their preferences, it will necessitate a focus on program access and the development of online communication frameworks.

Endotracheal tube microbiome in hospitalized patients defined largely by hospital environment.

BACKGROUND: Studies of the respiratory tract microbiome primarily focus on airway and lung microbial diversity, but it is still unclear how these microbial communities may be affected by intubation and long periods in intensive care units (ICU), an aspect that today could aid in the understanding of COVID19 progression and disease severity. This study aimed to explore and characterize the endotracheal tube (ETT) microbiome by analyzing ETT-associated microbial communities. METHODS: This descriptive study was carried out on adult patients subjected to invasive mechanical ventilation from 2 to 21 days. ETT samples were obtained from 115 patients from ICU units in two hospitals. Bacteria isolated from endotracheal tubes belonging to the ESKAPE group were analyzed for biofilm formation using crystal violet quantification. Microbial profiles were obtained using Illumina sequencing of 16S rRNA gene. RESULTS: The ETT microbiome was mainly composed by the phyla Proteobacteria, Firmicutes and Bacteroidetes. Microbiome composition correlated with the ICU in which patients were hospitalized, while intubation time and diagnosis of ventilator-associated pneumonia (VAP) did not show any significant association. CONCLUSION: These results suggest that the ICU environment, or medical practices, could be a key to microbial colonization and have a direct influence on the ETT microbiomes of patients that require mechanical ventilation.

Direct and Indirect Effectiveness of mRNA Vaccination against Severe Acute Respiratory Syndrome Coronavirus 2 in Long-Term Care Facilities, Spain.

We conducted a registries-based cohort study of long-term care facility residents >65 years of age offered vaccination against severe acute respiratory syndrome coronavirus 2 before March 10, 2021, in Spain. Risk for infection in vaccinated and nonvaccinated persons was compared with risk in the same persons in a period before the vaccination campaign, adjusted by daily-varying incidence and reproduction number. We selected 299,209 persons; 99.0% had >1 dose, 92.6% had 2 doses, and 99.8% of vaccines were Pfizer/BioNTech (BNT162b2). For vaccinated persons with no previous infection, vaccine effectiveness was 81.8% (95% CI 81.0%-82.7%), and 11.6 (95% CI 11.3-11.9) cases were prevented per 10,000 vaccinated/day. In those with previous infection, effectiveness was 56.8% (95% CI 47.1%-67.7%). In nonvaccinated residents with no previous infection, risk decreased by up to 81.4% (95% CI 73.3%-90.3%). Our results confirm vaccine effectiveness in this population and suggest indirect protection in nonvaccinated persons.

Widespread Detection of Multiple Strains of Crimean-Congo Hemorrhagic Fever Virus in Ticks, Spain.

Human cases of Crimean-Congo hemorrhagic fever (CCHF) were first detected in Spain in 2016. National human and animal health authorities organized a large, multidisciplinary study focusing on ticks as sentinels to determine the nationwide distribution of ticks with CCHF virus. Ticks were collected from animals and vegetation, samples pooled (12,584 ticks; 4,556 pools), and molecular methods used to look for the virus. We detected the virus in 135 pools from most of the regions studied, indicating that it is widespread in Spain. We found sequences of CCHF virus genotypes I, III, and IV in the tick species collected, most commonly in Hyalomma lusitanicum, suggesting this tick has a prominent role in the virus's natural cycle. The red deer (Cervus elaphus) was the host that most frequently yielded positive ticks. Our study highlights the need for larger studies in Spain to ascertain the complete risk to public health.

Serial DaT-SPECT imaging in asymptomatic carriers of LRRK2 G2019S mutation: 8 years' follow-up.

BACKGROUND: Carriers of the G2019S mutation of LRRK2 provide a great opportunity to investigate the premotor stages of Parkinson's disease (PD). We have studied by serial clinical and dopamine transporter single photon emission computed tomography (DaT-SPECT) evaluations a cohort of asymptomatic carriers of the LRRK2-G2019S mutation in order to evaluate the usefulness of these tools as biomarkers. Here we report the results of the extended follow-up of this cohort at 8 years. METHODS: Seventeen participants, of the 25 available from the 4-year evaluation, completed the 8-year assessment. UPDRS-III, UPSIT test and DaT-SPECT imaging (123 I-ioflupane) were performed. We used repeated-measures linear mixed effects models to examine the changes in DaT binding over time. RESULTS: Three carriers had converted to PD at 4 years. One additional carrier converted at 8 years. PD-converters had lower striatal DaT binding at baseline than non-converters. There was a significant decline of DaT binding over time, with a mean annual rate of 3.5%, with somewhat inter-individual and intra-individual variability and comparable between PD-converters and non-converters. No carrier with DAT binding ratio above an undefined threshold between 0.5 and 0.8 developed PD symptoms. The age-adjusted UPSIT score did not change significantly over time. CONCLUSIONS: The rate of conversion to PD at 8 years in this cohort aged ~58 years at baseline was 16%. The observed decline of DaT binding over time and its association with the phenotype render DaT-SPECT a potentially useful tool for monitoring the premotor stage of the disease, although at the individual level its ability to predict phenoconversion is limited.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes.

BACKGROUND: Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs. RESULTS: By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits. CONCLUSIONS: Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

Occupation, occupational exposures and mammographic density in Spanish women.

INTRODUCTION: Mammographic density (MD), the proportion of radiologically dense breast tissue, is a strong risk factor for breast cancer. Our objective is to investigate the influence of occupations and occupational exposure to physical, chemical, and microbiological agents on MD in Spanish premenopausal women. METHODS: This is a cross-sectional study based on 1362 premenopausal workers, aged 39-50, who attended a gynecological screening in a breast radiodiagnosis unit of Madrid City Council. The work history was compiled through a personal interview. Exposure to occupational agents was evaluated using the Spanish job-exposure matrix MatEmESp. MD percentage was assessed using the validated semi-automated computer tool DM-Scan. The association between occupation, occupational exposures, and MD was quantified using multiple linear regression models, adjusted for age, educational level, body mass index, parity, previous breast biopsies, family history of breast cancer, energy intake, use of oral contraceptives, smoking, and alcohol consumption. RESULTS: Although no occupation was statistically significantly associated with MD, a borderline significant inverse association was mainly observed in orchard, greenhouse, nursery, and garden workers (ß = -6.60; 95% confidence interval (95%CI) = -14.27; 1.07) and information and communication technology technicians (ß = -7.27; 95%CI = -15.37; 0.84). On the contrary, a positive association was found among technicians in art galleries, museums, and libraries (ß = 8.47; 95%CI = -0.65; 17.60). Women occupationally exposed to fungicides, herbicides, and insecticides tended to have lower MD. The percentage of density decreased by almost 2% for every 5 years spent in occupations exposed to the mentioned agents. CONCLUSIONS: Although our findings point to a lack of association with the occupations and exposures analyzed, this study supports a deeper exploration of the role of certain occupational agents in MD, such as pesticides.