Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep.

BACKGROUND: Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs). RESULTS: Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms. CONCLUSIONS: In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Pedigree diversity and implications for genetic selection of Katahdin sheep.

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (Ne) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the Ne is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current Ne suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes.

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Key genetic variants associated with variation of milk oligosaccharides from diverse human populations.

Human milk oligosaccharides (HMO), the third most abundant component of human milk, are thought to be important contributors to infant health. Studies have provided evidence that geography, stage of lactation, and Lewis and secretor blood groups are associated with HMO profile. However, little is known about how variation across the genome may influence HMO composition among women in various populations. In this study, we performed genome-wide association analyses of 395 women from 8 countries to identify genetic regions associated with 19 different HMO. Our data support FUT2 as the most significantly associated (P < 4.23-9 to P < 4.5-70) gene with seven HMO and provide evidence of balancing selection for FUT2. Although polymorphisms in FUT3 were also associated with variation in lacto-N-fucopentaose II and difucosyllacto-N-tetrose, we found little evidence of selection on FUT3. To our knowledge, this is the first report of the use of genome-wide association analyses on HMO.

A novel understanding of global DNA methylation in bobcat (Lynx rufus).

Epigenetic mechanisms may provide a novel prospective of bobcat (Lynx rufus) adaptation to habitat loss/fragmentation. Previous research has focused on bobcat behavior and genetics, but epigenetics has not been studied in bobcat. The aim of this study was to determine the quantity of global DNA methylation in the liver of 30 bobcats. DNA was extracted from liver samples obtained from the Vermont Fish and Wildlife Department. The percent of global DNA methylation was quantified and calculated using the MethylFlashTM Methylated DNA 5-mC Quantification Kit from Epigentek (Farmingdale, NY, USA). Age, sex, and carcass weight data were collected at sampling and analyzed with percent of global DNA methylation. Global DNA methylation was found to range from 0.46% to 2.76%. Age ranged from <1 to 12 years old and weight ranged from 3.18 to 13.61 kg. Further analysis of differential methylation may provide insight into novel means of bobcat conservation within different regions of Vermont. These results reinforce the need for genome-wide epigenetic studies in conservation biology.

Meiotic Recombination Differences in Rams from Three Breeds of Sheep in the US.

Meiotic recombination is an important contributor to genetic variation and ensures proper chromosome segregation during gametogenesis. Previous studies suggest that at least 1 crossover (CO) per chromosome arm is important to avoid mis-segregation. While the total number of COs per spermatocyte is known to differ in mice, this is only beginning to be evaluated in sheep. This study used a cytogenetic approach to quantify and compare the number of COs per spermatocyte in rams from 3 breeds of sheep: Suffolk, Icelandic, and Targhee. In total, 2,758 spermatocytes and over 170,000 COs were examined. Suffolk rams exhibited the lowest mean number of COs (61.1 ± 0.15) compared to Icelandic (63.5 ± 0.27) and Targhee (65.9 ± 0.26) rams. Significant differences in the number of COs per spermatocyte were observed between Suffolk, Icelandic, and Targhee breeds as well as within each breed. Additionally, the number and location of COs were characterized for homologous chromosomes in a subset of spermatocytes for each ram. A positive correlation was identified between the number of COs and the length of the homologous chromosome pair. Suffolk and Icelandic rams exhibited up to 7 COs per chromosome, while Targhee rams exhibited up to 9. Further, distinct CO location preferences on homologous chromosome pairs with 1, 2, 3, and 4 COs were observed in all 3 breeds. These data in sheep will aid in elucidating the mechanism of mammalian meiotic recombination, an important contributor to genetic diversity.

Altered cohesin gene dosage affects Mammalian meiotic chromosome structure and behavior.

Based on studies in mice and humans, cohesin loss from chromosomes during the period of protracted meiotic arrest appears to play a major role in chromosome segregation errors during female meiosis. In mice, mutations in meiosis-specific cohesin genes cause meiotic disturbances and infertility. However, the more clinically relevant situation, heterozygosity for mutations in these genes, has not been evaluated. We report here evidence from the mouse that partial loss of gene function for either Smc1b or Rec8 causes perturbations in the formation of the synaptonemal complex (SC) and affects both synapsis and recombination between homologs during meiotic prophase. Importantly, these defects increase the frequency of chromosomally abnormal eggs in the adult female. These findings have important implications for humans: they suggest that women who carry mutations or variants that affect cohesin function have an elevated risk of aneuploid pregnancies and may even be at increased risk of transmitting structural chromosome abnormalities.

Bisphenol A alters early oogenesis and follicle formation in the fetal ovary of the rhesus monkey.

Widespread use of the endocrine disrupting chemical bisphenol A (BPA) in consumer products has resulted in nearly continuous human exposure. In rodents, low-dose exposures have been reported to adversely affect two distinct stages of oogenesis in the developing ovary: the events of prophase at the onset of meiosis in the fetal ovary and the formation of follicles in the perinatal ovary. Because these effects could influence the reproductive longevity and success of the exposed individual, we conducted studies in the rhesus monkey to determine whether BPA induces similar disturbances in the developing primate ovary. The routes and levels of human exposure are unclear; hence, two different exposure protocols were used: single daily oral doses and continuous exposure via subdermal implant. Our analyses of second trimester fetuses exposed at the time of meiotic onset suggest that, as in mice, BPA induces subtle disturbances in the prophase events that set the stage for chromosome segregation at the first meiotic division. Our analyses of third-trimester fetuses exposed to single daily oral doses during the time of follicle formation revealed an increase in multioocyte follicles analogous to that reported in rodents. However, two unique phenotypes were evident in continuously exposed animals: persistent unenclosed oocytes in the medullary region and small, nongrowing oocytes in secondary and antral follicles. Because effects on both stages of oogenesis were elicited using doses that yield circulating levels of BPA analogous to those reported in humans, these findings raise concerns for human reproductive health.

Corrigendum: Genetic association of wool quality characteristics in United States Rambouillet sheep.

[This corrects the article DOI: 10.3389/fgene.2022.1081175.].

Genome-wide associations with longevity and reproductive traits in U.S. rangeland ewes.

Introduction: Improving ewe longevity is an important breeding and management goal, as death loss and early culling of mature ewes are economic burdens in the sheep industry. Ewe longevity can be improved by selecting for positive reproductive outcomes. However, the breeding approaches for accomplishing this come with the challenge of recording a lifetime trait. Characterizing genetic factors underpinning ewe longevity and related traits could result in the development of genomic selection strategies to improve the stayability of sheep through early, informed selection of replacement ewes. Methods: Towards this aim, a genome-wide association study (GWAS) was performed to identify genetic markers associated with ewe longevity, reproductive, and production traits. Traits evaluated included longevity (i.e., length of time in the flock), parity and the lifetime number of lambs born, lambs born alive, lambs weaned, and weight of lambs weaned. Ewe records from previous studies were used. Specifically, Rambouillet (n = 480), Polypay (n = 404), Suffolk (n = 182), and Columbia (n = 64) breed ewes (N = 1,130) were analyzed against 503,617 SNPs in across-breed and within-breed GWAS conducted with the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) model in R. Results: The across-breed GWAS identified 25 significant SNPs and the within-breed GWAS for Rambouillet, Polypay, and Suffolk ewes identified an additional 19 significant SNPs. The most significant markers were rs411309094 (13:22,467,143) associated with longevity in across-breed GWAS (p-value = 8.3E-13) and rs429525276 (2:148,398,336) associated with both longevity (p-value = 6.4E-15) and parity (p-value = 4.8E-15) in Rambouillet GWAS. Significant SNPs were identified within or in proximity (±50 kb) of genes with known or proposed roles in reproduction, dentition, and the immune system. These genes include ALPL, ANOS1, ARHGEF26, ASIC2, ASTN2, ATP8A2, CAMK2D, CEP89, DISC1, ITGB6, KCNH8, MBNL3, MINDY4, MTSS1, PLEKHA7, PRIM2, RNF43, ROBO2, SLCO1A2, TMEM266, TNFRSF21, and ZNF804B. Discussion: This study proposes multiple SNPs as candidates for use in selection indices and suggests genes for further research towards improving understanding of the genetic factors contributing to longevity, reproductive, and production traits of ewes.

Myokines Produced by Cultured Bovine Satellite Cells Harvested from 3- and 11-Month-Old Angus Steers.

The myokines interleukin 6 (IL-6), interleukin 15 (IL-15), myonectin (CTRP15), fibronectin type III domain containing protein 5/irisin (FNDC5), and brain-derived neurotrophic factor (BDNF) are associated with skeletal muscle cell proliferation, differentiation, and muscle hypertrophy in biomedical model species. This study evaluated whether these myokines are produced by cultured bovine satellite cells (BSCs) harvested from 3- and 11-month-old commercial black Angus steers and if the expression and secretion of these targets change across 0, 12, 24, and 48 h in vitro. IL-6, IL-15, FNDC5, and BDNF expression were greater (p ≤ 0.05) in the differentiated vs. undifferentiated BSCs at 0, 12, 24, and 48 h. CTRP15 expression was greater (p ≤ 0.03) in the undifferentiated vs. differentiated BSCs at 24 and 48 h. IL-6 and CTRP15 protein from culture media were greater (p ≤ 0.04) in undifferentiated vs. differentiated BSCs at 0, 12, 24, and 48 h. BDNF protein was greater in the media of differentiated vs. undifferentiated BSCs at 0, 12, 24, and 48 h. IL-6, 1L-15, FNDC5, and BDNF are expressed in association with BSC differentiation, and CTRP15 appears to be expressed in association with BSC proliferation. This study also confirms IL-6, IL-15, CTRP15, and BDNF proteins present in media collected from primary cultures of BSCs.

Reducing fecal egg count through selective breeding alters dorper lamb response to Haemonchus contortus in an artificial challenge trial.

Infection by gastrointestinal nematodes (GIN), particularly Haemonchus contortus, can be detrimental to sheep health and performance. Genetic susceptibility to GIN varies between breeds, with those lacking high levels of natural resistance often requiring frequent anthelmintic treatment when facing parasitic challenge. Genetic technology can serve as a tool to decrease GIN susceptibility via selection for sheep with reduced fecal egg count (FEC) estimated breeding values (EBVs). However, the physiological changes that result from implementation of this strategy are not well described. Additionally, there is a need for comparison of animals from recent selective breeding against breeds with inherent GIN resistance. In this study we administered a challenge of H. contortus to Dorper x White Dorper (DWD; n = 92) lambs that have been genetically selected for either low (DWD-) or high (DWD+) FEC EBVs and Barbados Blackbelly x Mouflon (BBM; n = 19) lambs from a genetically resistant breed backgrounds. Lamb FEC, packed-cell volume (PCV) and serum IgG were measured at intermittent levels over 5 weeks. At day 21 and day 35, the selectively bred DWD- had a lower mean FEC compared to DWD+, but were higher than BBM. Reductions in both PCV and serum IgG from initial day 0 levels were observed in DWD lambs, but not in BBM. Furthermore, from a subset of lambs (n = 24) harvested at day 21, DWD- only tended (p = 0.056) to have lower mean worm counts than DWD+, with BBM having the lowest mean worm count. Differentially expressed genes (DEGs) identified via RNA-sequencing of abomasal tissue at day 21 indicate a more pronounced Th2 immune response and more rapid worm expulsion occurred in iBBM than iDWD- and iDWD+ lambs. However, gene expression in DWD- suggests an association between reduced FEC EBV and gastric acid secretion and the ability to limit worm fecundity. Ultimately, selection of Dorper sheep for low FEC EBV can reduce susceptibility to GIN, but it will likely require multiple generations with this trait as a breeding priority before presenting a similar resistance level to Caribbean breeds.

Expression and secretion of SPARC, FGF-21 and DCN in bovine muscle cells: Effects of age and differentiation.

Skeletal muscle growth is an economically important trait in the cattle industry. Secreted muscle-derived proteins, referred to as myokines, have important roles in regulating the growth, metabolism, and health of skeletal muscle in human and biomedical research models. Accumulating evidence supports the importance of myokines in skeletal muscle and whole-body health, though little is known about the potential presence and functional significance of these proteins in cattle. This study evaluates and confirms that secreted proteins acidic and rich in cysteine (SPARC), fibroblast growth factor 21 (FGF-21), myostatin (MSTN), and decorin (DCN) are expressed and SPARC, FGF-21, and DCN are secreted by primary bovine satellite cells from 3- (BSC3; n = 3) and 11- (BSC11; n = 3) month -old commercial angus steers. Cells were cultured and collected at zero, 12, 24, and 48 hours to characterize temporal expression and secretion from undifferentiated and differentiated cells. The expression of SPARC was higher in the undifferentiated (p = 0.04) and differentiated (p = 0.07) BSC11 than BSC3. The same was observed with protein secretion from undifferentiated (p <0.0001) BSC11 compared to BSC3. Protein secretion of FGF-21 was higher in undifferentiated BSC11 (p < 0.0001) vs. BSC3. DCN expression was higher in differentiated BSC11 (p = 0.006) vs. BSC3. Comparing undifferentiated vs. differentiated BSC, MSTN expression was higher in differentiated BSC3 (p ≤ 0.001) for 0, 12, and 24 hours and in BSC11 (p ≤ 0.03) for 0, 12, 24, and 48 hours. There is also a change over time for SPARC expression (p ≤ 0.03) in undifferentiated and differentiated BSC and protein secretion (p < 0.0001) in undifferentiated BSC, as well as FGF-21 expression (p = 0.007) in differentiated BSC. This study confirms SPARC, FGF-21, and DCN are secreted, and SPARC, FGF-21, MSTN, and DCN are expressed in primary bovine muscle cells with age and temporal differences.

Identifying Regions of the Genome Associated with Conception Rate to the First Service in Holstein Heifers Bred by Artificial Insemination and as Embryo Transfer Recipients.

Heifer conception rate to the first service (HCR1) is defined as the number of heifers that become pregnant to the first breeding service compared to the heifers bred. This study aimed to identify loci associated and gene sets enriched for HCR1 for heifers that were bred by artificial insemination (AI, n = 2829) or were embryo transfer (ET, n = 2086) recipients, by completing a genome-wide association analysis and gene set enrichment analysis using SNP data (GSEA-SNP). Three unique loci, containing four positional candidate genes, were associated (p < 1 × 10-5) with HCR1 for ET recipients, while the GSEA-SNP identified four gene sets (NES ≥ 3) and sixty-two leading edge genes (LEGs) enriched for HCR1. While no loci were associated with HCR1 bred by AI, one gene set and twelve LEGs were enriched (NES ≥ 3) for HCR1 with the GSEA-SNP. This included one gene (PKD2) shared between HCR1 AI and ET services. Identifying loci associated or enriched for HCR1 provides an opportunity to use them as genomic selection tools to facilitate the selection of cattle with higher reproductive efficiency, and to better understand embryonic loss.

The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content.

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

The Ruminant Telomere-to-Telomere (RT2T) Consortium.

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Human milk oligosaccharides are associated with maternal genetics and respiratory health of human milk-fed children.

Breastfeeding provides many health benefits, but its impact on respiratory health remains unclear. This study addresses the complex and dynamic nature of the mother-milk-infant triad by investigating maternal genomic factors regulating human milk oligosaccharides (HMOs), and their associations with respiratory health among human milk-fed infants. Nineteen HMOs are quantified from 980 mothers of the CHILD Cohort Study. Genome-wide association studies identify HMO-associated loci on chromosome 19p13.3 and 19q13.33 (lowest P = 2.4e-118), spanning several fucosyltransferase (FUT) genes. We identify novel associations on chromosome 3q27.3 for 6'-sialyllactose (P = 2.2e-9) in the sialyltransferase (ST6GAL1) gene. These, plus additional associations on chromosomes 7q21.32, 7q31.32 and 13q33.3, are replicated in the independent INSPIRE Cohort. Moreover, gene-environment interaction analyses suggest that fucosylated HMOs may modulate overall risk of recurrent wheeze among preschoolers with variable genetic risk scores (P < 0.01). Thus, we report novel genetic factors associated with HMOs, some of which may protect the respiratory health of children.