Evaluation of Tomato Germplasm against Tomato Brown Rugose Fruit Virus and Identification of Resistance in Solanum pimpinellifolium

The tomato is one of the most important vegetable crops grown worldwide. Tomato brown rugose fruit virus (ToBRFV), a seed-borne tobamovirus, poses a serious threat to tomato production due to its ability to break the resistant genes (Tm-1, Tm-2, Tm-22) in tomatoes. The objective of this work was to identify new resistant source(s) of tomato germplasm against ToBRFV. To achieve this aim, a total of 476 accessions from 12 Solanum species were tested with the ToBRFV US isolate for their resistance and susceptibility. As a result, a total of 44 asymptomatic accessions were identified as resistant/tolerant, including thirty-one accessions of S. pimpinellifolium, one accession of S. corneliomulleri, four accessions of S. habrochaites, three accessions of S. peruvianum, and five accessions of S. subsection lycopersicon hybrid. Further analyses using serological tests identified four highly resistant S. pimpinellifolium lines, PI 390713, PI 390714, PI 390716, and PI 390717. The inheritance of resistance in the selected lines was verified in the next generation and confirmed using RT-qPCR. To our knowledge, this is a first report of high resistance to ToBRFV in S. pimpinellifolium. These new genetic resources will expand the genetic pool available for breeders to develop new resistant cultivars of tomato against ToBRFV.

Genome-Wide Association Study Reveals a QTL Region for Ashy Stem Blight Resistance in PRA154 Andean Common Bean

Ashy stem blight (ASB) caused by Macrophomina phaseolina (Tassi) Goidanich affects the common bean (Phaseolus vulgaris L.) at all growing stages. Higher levels of resistance were observed in Andean common beans, but specific resistant quantitative trait loci (QTLs) conferring resistance to this pathogen have not been reported in this gene pool. The objectives of this research were to: (i) conduct a genome-wide association study (GWAS) and QTL mapping for resistance in the Andean breeding line PRA154; and (ii) identify single nucleotide polymorphism (SNP) markers and candidate genes for ASB resistance. Phenotyping was conducted under greenhouse conditions by inoculating the 107 F6:7 recombinant inbred lines (RILs) derived from the cross between the susceptible cultivar 'Verano' and the partial-resistant breeding line PRA154 twice with the M. phaseolina isolate PRI21. Genotyping was performed with 109,040 SNPs distributed across all 11 P. vulgaris chromosomes. A novel major QTL was located between 28,761,668 and 31,263,845 bp, extending 2.5 Mbp on chromosome Pv07, and the highest significant SNP markers were Chr07_28761668_A_G, Chr07_29131720_G_A, and Chr07_31263845_C_T with the highest LOD (more than 10 in most of the cases) and R-squared values, explaining 40% of the phenotypic variance of the PRI21 isolate. The gene Phvul.007G173900 (methylcrotonyl-CoA carboxylase alpha chain and mitochondrial 3-methylcrotonyl-CoA carboxylase 1 [MCCA]) with a size of 10,891 bp, located between 29,131,591 and 29,142,481 bp on Pv07, was identified as one candidate for ASB resistance in PRA154, and it contained Chr07_29131720_G_A. The QTL and genetic marker information could be used to assist common bean breeders to develop germplasm and cultivars with ASB resistance through molecular breeding.

Genome-wide characterization and evolutionary analysis of the AP2/ERF gene family in lettuce (Lactuca sativa)

The APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) gene family plays vital roles in plants, serving as a key regulator in responses to abiotic stresses. Despite its significance, a comprehensive understanding of this family in lettuce remains incomplete. In this study, we performed a genome-wide search for the AP2/ERF family in lettuce and identified a total of 224 members. The duplication patterns provided evidence that both tandem and segmental duplications contributed to the expansion of this family. Ka/Ks ratio analysis demonstrated that, following duplication events, the genes have been subjected to purifying selection pressure, leading to selective constraints on their protein sequence. This selective pressure provides a dosage benefit against stresses in plants. Additionally, a transcriptome analysis indicated that some duplicated genes gained novel functions, emphasizing the contribution of both dosage effect and functional divergence to the family functionalities. Furthermore, an orthologous relationship study showed that 60% of genes descended from a common ancestor of Rosid and Asterid lineages, 28% from the Asterid ancestor, and 12% evolved in the lettuce lineage, suggesting lineage-specific roles in adaptive evolution. These results provide valuable insights into the evolutionary mechanisms of the AP2/ERF gene family in lettuce, with implications for enhancing abiotic stress tolerance, ultimately contributing to the genetic improvement of lettuce crop production.

Genome-Wide Association Study and Genomic Prediction of Fusarium Wilt Resistance in Common Bean Core Collection

The common bean (Phaseolus vulgaris L.) is a globally cultivated leguminous crop. Fusarium wilt (FW), caused by Fusarium oxysporum f. sp. phaseoli (Fop), is a significant disease leading to substantial yield loss in common beans. Disease-resistant cultivars are recommended to counteract this. The objective of this investigation was to identify single nucleotide polymorphism (SNP) markers associated with FW resistance and to pinpoint potential resistant common bean accessions within a core collection, utilizing a panel of 157 accessions through the Genome-wide association study (GWAS) approach with TASSEL 5 and GAPIT 3. Phenotypes for Fop race 1 and race 4 were matched with genotypic data from 4740 SNPs of BARCBean6K_3 Infinium Bea Chips. After ranking the 157-accession panel and revealing 21 Fusarium wilt-resistant accessions, the GWAS pinpointed 16 SNPs on chromosomes Pv04, Pv05, Pv07, Pv8, and Pv09 linked to Fop race 1 resistance, 23 SNPs on chromosomes Pv03, Pv04, Pv05, Pv07, Pv09, Pv10, and Pv11 associated with Fop race 4 resistance, and 7 SNPs on chromosomes Pv04 and Pv09 correlated with both Fop race 1 and race 4 resistances. Furthermore, within a 30 kb flanking region of these associated SNPs, a total of 17 candidate genes were identified. Some of these genes were annotated as classical disease resistance protein/enzymes, including NB-ARC domain proteins, Leucine-rich repeat protein kinase family proteins, zinc finger family proteins, P-loopcontaining nucleoside triphosphate hydrolase superfamily, etc. Genomic prediction (GP) accuracy for Fop race resistances ranged from 0.26 to 0.55. This study advanced common bean genetic enhancement through marker-assisted selection (MAS) and genomic selection (GS) strategies, paving the way for improved Fop resistance.

Genetic Architecture of Salt Tolerance in Cowpea (Vigna unguiculata (L.) Walp.) at Seedling Stage Using a Whole Genome Resequencing Approach

Cowpea (Vigna unguiculata (L.) Walp.) is a diploid legume crop used for human consumption, feed for livestock, and cover crops. Earlier reports have shown that salinity has been a growing threat to cowpea cultivation. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify SNP markers and to investigate candidate genes for salt tolerance in cowpea. A total of 331 cowpea genotypes were evaluated for salt tolerance by supplying a solution of 200 mM NaCl in our previous work. The cowpea panel was genotyped using a whole genome resequencing approach, generating 14,465,516 SNPs. Moreover, 5,884,299 SNPs were used after SNP filtering. GWAS was conducted on a total of 296 cowpea genotypes that have high-quality SNPs. BLINK was used for conducting GWAS. Results showed (1) a strong GWAS peak on an 890-bk region of chromosome 2 for leaf SPAD chlorophyll under salt stress in cowpea and harboring a significant cluster of nicotinamide adenine dinucleotide (NAD) dependent epimerase/dehydratase genes such as Vigun02g128900.1, Vigun02g129000.1, Vigun02g129100.1, Vigun02g129200.1, and Vigun02g129500.1; (2) two GWAS peaks associated with relative tolerance index for chlorophyll were identified on chromosomes 1 and 2. The peak on chromosome 1 was defined by a cluster of 10 significant SNPs mapped on a 5 kb region and was located in the vicinity of Vigun01g086000.1, encoding for a GATA transcription factor. The GWAS peak on chromosome 2 was defined by a cluster of 53 significant SNPs and mapped on a 68 bk region of chromosome 2, and (3) the highest GWAS peak was identified on chromosome 3, and this locus was associated with leaf score injury. This peak was within the structure of a potassium channel gene (Vigun03g144700.1). To the best of our knowledge, this is one the earliest reports on the salt tolerance study of cowpea using whole genome resequencing data.

High Dose Hyperfractionated Thoracic Radiotherapy vs. Standard Dose for Limited Stage Small-Cell Lung Cancer: A Multicenter, Open-Label Randomized, Phase 3 Trial

PURPOSE/OBJECTIVE(S)

Limited stage small-cell lung cancer (LS-SCLC) is associated with poor prognosis. We aimed to assess the efficacy and safety of high-dose, hyperfractionated thoracic radiotherapy of 54 Gy in 30 fractions compared with standard dose (45 Gy in 30 fractions) as a first-line treatment for LS-SCLC.

MATERIALS/METHODS

The study was an open-label, randomized, phase 3 trial, done at 16 public hospitals in China. Key inclusion criteria were patients aged 18-70 years, with previously histologically or cytologically confirmed LS-SCLC, previously untreated or received 1-2 courses of intravenous cisplatin (75 mg/m²of body-surface area, on day 1 or divided into two days of each cycle) or carboplatin (area under the curve of 5 mg/mL per min, day 1 of each cycle)and intravenous etoposide (100 mg/m²of body-surface area, on days 1-3 of each cycle), and an Eastern Cooperative Oncology Group (ECOG) performance status of 0-1.Eligible patients were randomly assigned (1:1) to receive volumetric-modulated arc radiotherapy (VMAT) of 45 Gy in 30 fractions or the simultaneous integrated boost VMAT (SIB-VMAT) of 54 Gy in 30 fractions to the primary lung tumor and lymph node metastases starting 0-42 days after the first chemotherapy course. Both groups of patients received thoracic radiotherapy twice per day and 10 fractions per week. Prophylactic cranial radiation (PCI, 25 Gy in 10 fractions) was implemented to patients with responsive disease. The primary endpoint was overall survival. Safety was analyzed in the as-treated population.

RESULTS

Between June 30, 2017, and April 6, 2021, 224 eligible patients were enrolled and randomly assigned to 54 Gy (n = 108) or 45 Gy (n = 116). Median follow-up for the primary analysis was 45 months (IQR 41-48). Median overall survival was significantly improved in the 54 Gy group (62.4 months) compared with the 45 Gy group (43.1 months; p = 0.001). Median progression-free survival was significantly improved in the 54 Gy group (30.5 months) compared with the 45 Gy group (16.7 months; p = 0.044). The most common grade 3-4 adverse events were neutropenia (30 [28%] of 108 patients in the 54 Gy group vs 27 [23%] of 116 patients in the 45 Gy group), neutropenic infections (6 [6%] vs 2 [2%]), thrombocytopenia (13 [12%] vs 12 [10%]), anemia (6 [6%] vs 4 [3%]), and esophagitis (1 [1%] vs 3 [3%]). Treatment-related serious adverse events occurred in 9 [8%] patients in the 54 Gy group and 16 [14%] patients in the 45 Gy group. There were one treatment-related deaths in 54 Gy group (myocardial infarction).

CONCLUSION

Compared with standard thoracic radiotherapy dose of 45 Gy, the high dose of 54 Gy improved overall survival and progression-free survival without increasing toxicities in patients with LS-SCLC, supporting twice-daily hyperfractionated thoracic radiotherapy of 54 Gy with concurrent chemotherapy is an alternative treatment option for LS-SCLC. This study is complete and registered with ClinicalTrials.gov, NCT03214003.

A Genome-Wide Association Study Reveals Region Associated with Seed Protein Content in Cowpea

Cowpea (Vigna unguiculata L. Walp., 2n = 2x = 22) is a protein-rich crop that complements staple cereals for humans and serves as fodder for livestock. It is widely grown in Africa and other developing countries as the primary source of protein in the diet; therefore, it is necessary to identify the protein-related loci to improve cowpea breeding. In the current study, we conducted a genome-wide association study (GWAS) on 161 cowpea accessions (151 USDA germplasm plus 10 Arkansas breeding lines) with a wide range of seed protein contents (21.8~28.9%) with 110,155 high-quality whole-genome single-nucleotide polymorphisms (SNPs) to identify markers associated with protein content, then performed genomic prediction (GP) for future breeding. A total of seven significant SNP markers were identified using five GWAS models (single-marker regression (SMR), the general linear model (GLM), Mixed Linear Model (MLM), Fixed and Random Model Circulating Probability Unification (FarmCPU), and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK), which are located at the same locus on chromosome 8 for seed protein content. This locus was associated with the gene Vigun08g039200, which was annotated as the protein of the thioredoxin superfamily, playing a critical function for protein content increase and nutritional quality improvement. In this study, a genomic prediction (GP) approach was employed to assess the accuracy of predicting seed protein content in cowpea. The GP was conducted using cross-prediction with five models, namely ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB), and Bayesian least absolute shrinkage and selection operator (BL), applied to seven random whole genome marker sets with different densities (10 k, 5 k, 2 k, 1 k, 500, 200, and 7), as well as significant markers identified through GWAS. The accuracies of the GP varied between 42.9% and 52.1% across the seven SNPs considered, depending on the model used. These findings not only have the potential to expedite the breeding cycle through early prediction of individual performance prior to phenotyping, but also offer practical implications for cowpea breeding programs striving to enhance seed protein content and nutritional quality.

Skim resequencing finely maps the downy mildew resistance loci RPF2 and RPF3 in spinach cultivars whale and Lazio

Commercial production of spinach (Spinacia oleracea L.) is centered in California and Arizona in the US, where downy mildew caused by Peronospora effusa is the most destructive disease. Nineteen typical races of P. effusa have been reported to infect spinach, with 16 identified after 1990. The regular appearance of new pathogen races breaks the resistance gene introgressed in spinach. We attempted to map and delineate the RPF2 locus at a finer resolution, identify linked single nucleotide polymorphism (SNP) markers, and report candidate downy mildew resistance (R) genes. Progeny populations segregating for RPF2 locus derived from resistant differential cultivar Lazio were infected using race 5 of P. effusa and were used to study for genetic transmission and mapping analysis in this study. Association analysis performed with low coverage whole genome resequencing-generated SNP markers mapped the RPF2 locus between 0.47 to 1.46 Mb of chromosome 3 with peak SNP (Chr3_1, 221, 009) showing a LOD value of 61.6 in the GLM model in TASSEL, which was within 1.08 Kb from Spo12821, a gene that encodes CC-NBS-LRR plant disease resistance protein. In addition, a combined analysis of progeny panels of Lazio and Whale segregating for RPF2 and RPF3 loci delineated the resistance section in chromosome 3 between 1.18-1.23 and 1.75-1.76 Mb. This study provides valuable information on the RPF2 resistance region in the spinach cultivar Lazio compared to RPF3 loci in the cultivar Whale. The RPF2 and RPF3 specific SNP markers, plus the resistant genes reported here, could add value to breeding efforts to develop downy mildew resistant cultivars in the future.

A genome-wide association study and genomic prediction for Phakopsora pachyrhizi resistance in soybean

Soybean brown rust (SBR), caused by Phakopsora pachyrhizi, is a devastating fungal disease that threatens global soybean production. This study conducted a genome-wide association study (GWAS) with seven models on a panel of 3,082 soybean accessions to identify the markers associated with SBR resistance by 30,314 high quality single nucleotide polymorphism (SNPs). Then five genomic selection (GS) models, including Ridge regression best linear unbiased predictor (rrBLUP), Genomic best linear unbiased predictor (gBLUP), Bayesian least absolute shrinkage and selection operator (Bayesian LASSO), Random Forest (RF), and Support vector machines (SVM), were used to predict breeding values of SBR resistance using whole genome SNP sets and GWAS-based marker sets. Four SNPs, namely Gm18_57,223,391 (LOD = 2.69), Gm16_29,491,946 (LOD = 3.86), Gm06_45,035,185 (LOD = 4.74), and Gm18_51,994,200 (LOD = 3.60), were located near the reported P. pachyrhizi R genes, Rpp1, Rpp2, Rpp3, and Rpp4, respectively. Other significant SNPs, including Gm02_7,235,181 (LOD = 7.91), Gm02_7234594 (LOD = 7.61), Gm03_38,913,029 (LOD = 6.85), Gm04_46,003,059 (LOD = 6.03), Gm09_1,951,644 (LOD = 10.07), Gm10_39,142,024 (LOD = 7.12), Gm12_28,136,735 (LOD = 7.03), Gm13_16,350,701(LOD = 5.63), Gm14_6,185,611 (LOD = 5.51), and Gm19_44,734,953 (LOD = 6.02), were associated with abundant disease resistance genes, such as Glyma.02G084100, Glyma.03G175300, Glyma.04g189500, Glyma.09G023800, Glyma.12G160400, Glyma.13G064500, Glyma.14g073300, and Glyma.19G190200. The annotations of these genes included but not limited to: LRR class gene, cytochrome 450, cell wall structure, RCC1, NAC, ABC transporter, F-box domain, etc. The GWAS based markers showed more accuracies in genomic prediction than the whole genome SNPs, and Bayesian LASSO model was the ideal model in SBR resistance prediction with 44.5% ~ 60.4% accuracies. This study aids breeders in predicting selection accuracy of complex traits such as disease resistance and can shorten the soybean breeding cycle by the identified markers.

A Genome-Wide Association Study and Genomic Prediction for Fiber and Sucrose Contents in a Mapping Population of LCP 85-384 Sugarcane

Sugarcane (Saccharum spp. hybrids) is an economically important crop for both sugar and biofuel industries. Fiber and sucrose contents are the two most critical quantitative traits in sugarcane breeding that require multiple-year and multiple-location evaluations. Marker-assisted selection (MAS) could significantly reduce the time and cost of developing new sugarcane varieties. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify DNA markers associated with fiber and sucrose contents and to perform genomic prediction (GP) for the two traits. Fiber and sucrose data were collected from 237 self-pollinated progenies of LCP 85-384, the most popular Louisiana sugarcane cultivar from 1999 to 2007. The GWAS was performed using 1310 polymorphic DNA marker alleles with three models of TASSEL 5, single marker regression (SMR), general linear model (GLM) and mixed linear model (MLM), and the fixed and random model circulating probability unification (FarmCPU) of R package. The results showed that 13 and 9 markers were associated with fiber and sucrose contents, respectively. The GP was performed by cross-prediction with five models, ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB) and Bayesian least absolute shrinkage and selection operator (BL). The accuracy of GP varied from 55.8% to 58.9% for fiber content and 54.6% to 57.2% for sucrose content. Upon validation, these markers can be applied in MAS and genomic selection (GS) to select superior sugarcane with good fiber and high sucrose contents.

Identification of a QTL region for ashy stem blight resistance using genome-wide association and linage analysis in common bean recombinant inbred lines derived from BAT 477 and NY6020-4

Ashy stem blight (ASB), caused by the fungus Macrophomina phaseolina (Tassi) Goidanich is an important disease of the common bean (Phaseolus vulgaris L.). It is important to identify quantitative trait loci (QTL) for ASB resistance and introgress into susceptible cultivars of the common bean. The objective of this research was to identify QTL and single nucleotide polymorphism (SNP) markers associated with ASB resistance in recombinant inbred lines (RIL) derived from a cross between BAT 477 and NY6020-4 common bean. One hundred and twenty-six F_6:7 RIL were phenotyped for ASB in the greenhouse. Disease severity was scored on a scale of 1-9. Genotyping was performed using whole genome resequencing with 2x common bean genome size coverage, and over six million SNPs were obtained. After being filtered, 72,017 SNPs distributed on 11 chromosomes were used to conduct the genome-wide association study (GWAS) and QTL mapping. A novel QTL region of ~4.28 Mbp from 35,546,329 bp to 39,826,434 bp on chromosome Pv03 was identified for ASB resistance. The two SNPs, Chr03_39824257 and Chr03_39824268 located at 39,824,257 bp and 39,824,268 bp on Pv03, respectively, were identified as the strongest markers associated with ASB resistance. The gene Phvul.003G175900 (drought sensitive, WD repeat-containing protein 76) located at 39,822,021 - 39,824,655 bp on Pv03 was recognized as one candidate for ASB resistance in the RIL, and the gene contained the two SNP markers. QTL and SNP markers may be used to select plants and lines for ASB resistance through marker-assisted selection (MAS) in common bean breeding.

Waterlogging during the reproductive growth stage causes physiological and biochemical modifications in the leaves of cowpea (Vigna unguiculata L.) genotypes with contrasting tolerance

Waterlogging causes various metabolic, physiological, and morphological changes in crops, resulting in yield loss of most legumes in rainfed and irrigated agriculture. However, research on cowpea genotypes using physiological and biochemical traits as a measure of tolerance to waterlogging stress is limited. We evaluated the impacts of 7 days of waterlogging (DOW) and 7 days of recovery (DOR) on the physiology and biochemistry of two cowpea (Vigna unguiculata (L.) Walp) genotypes (UCR 369 and EpicSelect.4) with contrasting waterlogging tolerance. Cowpea genotypes were grown in a controlled environment until the R2 stage and then subjected to 7 DOW. Later, the waterlogged plants were reoxygenated for an additional 7 DOR. Overall, cowpea genotypes had a contrasting response to waterlogging using different mechanisms. Compared to the control, the photosynthetic parameters of both cowpea genotypes were impaired under 7 DOW and could not recover at 7 DOR, with a larger decline in EpicSelect.4.7 DOW caused significant loss in the chlorophyll and carotenoid content of both genotypes. However, only waterlogged UCR 369 was not photo-inhibited and able to restore the levels of chlorophyll and carotenoids at 7 DOR. In addition, 7 DOW induced intense stress in UCR 369 with increased zeaxanthin, sucrose, and flavonoid content, while these metabolites were decreased in EpicSelect.4. On the other hand, glucose, fructose, and phenolic content were increased in EpicSelect.4 but decreased in UCR 369 at 7 DOR. In summary, compared to EpicSelect.4, UCR 369 restored their photosynthetic pigments and metabolites to the control levels at 7 DOR, indicating a likely tolerance to waterlogging stress.

Mapping and selection of downy mildew resistance in spinach cv. whale by low coverage whole genome sequencing

Spinach (Spinacia oleracea) is a popular leafy vegetable crop and commercial production is centered in California and Arizona in the US. The oomycete Peronospora effusa causes the most important disease in spinach, downy mildew. A total of nineteen races of P. effusa are known, with more than 15 documented in the last three decades, and the regular emergence of new races is continually overcoming the genetic resistance to the pathogen. This study aimed to finely map the downy mildew resistance locus RPF3 in spinach, identify single nucleotide polymorphism (SNP) markers associated with the resistance, refine the candidate genes responsible for the resistance, and evaluate the prediction performance using multiple machine learning genomic prediction (GP) methods. Segregating progeny population developed from a cross of resistant cultivar Whale and susceptible cultivar Viroflay to race 5 of P. effusa was inoculated under greenhouse conditions to determine downy mildew disease response across the panel. The progeny panel and the parents were resequenced at low coverage (1x) to identify genome wide SNP markers. Association analysis was performed using disease response phenotype data and SNP markers in TASSEL, GAPIT, and GENESIS programs and mapped the race 5 resistance loci (RPF3) to 1.25 and 2.73 Mb of Monoe-Viroflay chromosome 3 with the associated SNP in the 1.25 Mb region was 0.9 Kb from the NBS-LRR gene SOV3g001250. The RPF3 locus in the 1.22-1.23 Mb region of Sp75 chromosome 3 is 2.41-3.65 Kb from the gene Spo12821 annotated as NBS-LRR disease resistance protein. This study extended our understanding of the genetic basis of downy mildew resistance in spinach cultivar Whale and mapped the RPF3 resistance loci close to the NBS-LRR gene providing a target to pursue functional validation. Three SNP markers efficiently selected resistance based on multiple genomic selection (GS) models. The results from this study have added new genomic resources, generated an informed basis of the RPF3 locus resistant to spinach downy mildew pathogen, and developed markers and prediction methods to select resistant lines.

Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods

Downy mildew, commercially the most important disease of spinach, is caused by the obligate oomycete Peronospora effusa. In the past two decades, new pathogen races have repeatedly overcome the resistance used in newly released cultivars, urging the need for more durable resistance. Commercial spinach cultivars are bred with major R genes to impart resistance to downy mildew pathogens and are effective against some pathogen races/isolates. This work aimed to evaluate the worldwide USDA spinach germplasm collections and commercial cultivars for resistance to downy mildew pathogen in the field condition under natural inoculum pressure and conduct genome wide association analysis (GWAS) to identify resistance-associated genomic regions (alleles). Another objective was to evaluate the prediction accuracy (PA) using several genomic prediction (GP) methods to assess the potential implementation of genomic selection (GS) to improve spinach breeding for resistance to downy mildew pathogen. More than four hundred diverse spinach genotypes comprising USDA germplasm accessions and commercial cultivars were evaluated for resistance to downy mildew pathogen between 2017-2019 in Salinas Valley, California and Yuma, Arizona. GWAS was performed using single nucleotide polymorphism (SNP) markers identified via whole genome resequencing (WGR) in GAPIT and TASSEL programs; detected 14, 12, 5, and 10 significantly associated SNP markers with the resistance from four tested environments, respectively; and the QTL alleles were detected at the previously reported region of chromosome 3 in three of the four experiments. In parallel, PA was assessed using six GP models and seven unique marker datasets for field resistance to downy mildew pathogen across four tested environments. The results suggest the suitability of GS to improve field resistance to downy mildew pathogen. The QTL, SNP markers, and PA estimates provide new information in spinach breeding to select resistant plants and breeding lines through marker-assisted selection (MAS) and GS, eventually helping to accumulate beneficial alleles for durable disease resistance.

Growth and Photosynthetic Responses of Cowpea Genotypes under Waterlogging at the Reproductive Stage

Waterlogging is an important environmental stress limiting the productivity of crops worldwide. Cowpea (Vigna unguiculata L.) is particularly sensitive to waterlogging stress during the reproductive stage, with a consequent decline in pod formation and yield. However, little is known about the critical processes underlying cowpea’s responses to waterlogging during the reproductive stage. Thus, we investigated the key parameters influencing carbon fixation, including stomatal conductance (g_s), intercellular CO₂ concentration, chlorophyll content, and chlorophyll fluorescence, of two cowpea genotypes with contrasting waterlogging tolerance. These closely related genotypes have starkly contrasting responses to waterlogging during and after 7 days of waterlogging stress (DOW). In the intolerant genotype (‘EpicSelect.4’), waterlogging resulted in a gradual loss of pigment and decreased photosynthetic capacity as a consequent decline in shoot biomass. On the other hand, the waterlogging-tolerant genotype (‘UCR 369’) maintained CO₂ assimilation rate (A), stomatal conductance (g_s), biomass, and chlorophyll content until 5 DOW. Moreover, there was a highly specific downregulation of the mesophyll conductance (g_m), maximum rate of Rubisco (V_cmax), and photosynthetic electron transport rate (J_max) as non-stomatal limiting factors decreasing A in EpicSelect.4. Exposure of EpicSelect.4 to 2 DOW resulted in the loss of PSII photochemistry by downregulating the PSII quantum yield (F_v/F_m), photochemical efficiency (Φ_PSII), and photochemical quenching (qP). In contrast, we found no substantial change in the photosynthesis and chlorophyll fluorescence of UCR 369 in the first 5 DOW. Instead, UCR 369 maintained biomass accumulation, chlorophyll content, and Rubisco activity, enabling the genotype to maintain nutrient absorption and photosynthesis during the early period of waterlogging. However, compared to the control, both cowpea genotypes could not fully recover their photosynthetic capacity after 7 DOW, with a more significant decline in EpicSelect.4. Overall, our findings suggest that the tolerant UCR 369 genotype maintains higher photosynthesis under waterlogging stress attributable to higher photochemical efficiency, Rubisco activity, and less stomatal restriction. After recovery, the incomplete recovery of A can be attributed to the reduced g_s caused by severe waterlogging damage in both genotypes. Thus, promoting the rapid recovery of stomata from waterlogging stress may be crucial for the complete restoration of carbon fixation in cowpeas during the reproductive stage.

Genetic dissection of nitrogen induced changes in the shoot and root biomass of spinach

Efficient partitioning of above and below-ground biomass in response to nitrogen (N) is critical to the productivity of plants under sub-optimal conditions. It is particularly essential in vegetable crops like spinach with shallow root systems, a short growth cycle, and poor nitrogen use efficiency. In this study, we conducted a genome-wide association study (GWAS) to explore N-induced changes using spinach accessions with diverse genetic backgrounds. We evaluated phenotypic variations as percent changes in the shoot and root biomass in response to N using 201 spinach accessions grown in randomized complete blocks design in a soilless media under a controlled environment. A GWAS was performed for the percent changes in the shoot and root biomass in response to N in the 201 spinach accessions using 60,940 whole-genome resequencing generated SNPs. Three SNP markers, chr4_28292655, chr6_1531056, and chr6_37966006 on chromosomes 4 and 6, were significantly associated with %change in root weight, and two SNP markers, chr2_18480277 and chr4_47598760 on chromosomes 2 and 4, were significantly associated with % change shoot weight. The outcome of this study established a foundation for genetic studies needed to improve the partitioning of total biomass and provided a resource to identify molecular markers to enhance N uptake via marker-assisted selection or genomic selection in spinach breeding programs.

Genome-Wide Association Study and Genomic Prediction for Bacterial Wilt Resistance in Common Bean (Phaseolus vulgaris) Core Collection

Common bean (Phaseolus vulgaris) is one of the major legume crops cultivated worldwide. Bacterial wilt (BW) of common bean (Curtobacterium flaccumfaciens pv. flaccumfaciens), being a seed-borne disease, has been a challenge in common bean producing regions. A genome-wide association study (GWAS) was conducted to identify SNP markers associated with BW resistance in the USDA common bean core collection. A total of 168 accessions were evaluated for resistance against three different isolates of BW. Our study identified a total of 14 single nucleotide polymorphism (SNP) markers associated with the resistance to BW isolates 528, 557, and 597 using mixed linear models (MLMs) in BLINK, FarmCPU, GAPIT, and TASSEL 5. These SNPs were located on chromosomes Phaseolus vulgaris [Pv]02, Pv04, Pv08, and Pv09 for isolate 528; Pv07, Pv10, and Pv11 for isolate 557; and Pv04, Pv08, and Pv10 for isolate 597. The genomic prediction accuracy was assessed by utilizing seven GP models with 1) all the 4,568 SNPs and 2) the 14 SNP markers. The overall prediction accuracy (PA) ranged from 0.30 to 0.56 for resistance against the three BW isolates. A total of 14 candidate genes were discovered for BW resistance located on chromosomes Pv02, Pv04, Pv07, Pv08, and Pv09. This study revealed vital information for developing genetic resistance against the BW pathogen in common bean. Accordingly, the identified SNP markers and candidate genes can be utilized in common bean molecular breeding programs to develop novel resistant cultivars.

Insights into nitrogen metabolism in the wild and cultivated lettuce as revealed by transcriptome and weighted gene co-expression network analysis

Large amounts of nitrogen fertilizers applied during lettuce (Lactuca sativa L.) production are lost due to leaching or volatilization, causing severe environmental pollution and increased costs of production. Developing lettuce varieties with high nitrogen use efficiency (NUE) is the eco-friendly solution to reduce nitrogen pollution. Hence, in-depth knowledge of nitrogen metabolism and assimilation genes and their regulation is critical for developing high NUE varieties. In this study, we performed comparative transcriptomic analysis of the cultivated lettuce (L. sativa L.) and its wild progenitor (L. serriola) under high and low nitrogen conditions. A total of 2,704 differentially expressed genes were identified. Key enriched biological processes included photosynthesis, oxidation-reduction process, chlorophyll biosynthetic process, and cell redox homeostasis. The transcription factors (TFs) belonging to the ethylene responsive factor family and basic helix-loop-helix family were among the top differentially expressed TFs. Using weighted gene co-expression network analysis we constructed nine co-expression modules. Among these, two modules were further investigated because of their significant association with total nitrogen content and photosynthetic efficiency of photosystem II. Three highly correlated clusters were identified which included hub genes for nitrogen metabolism, secondary metabolites, and carbon assimilation, and were regulated by cluster specific TFs. We found that the expression of nitrogen transportation and assimilation genes varied significantly between the two lettuce species thereby providing the opportunity of introgressing wild alleles into the cultivated germplasm for developing lettuce cultivars with more efficient use of nitrogen.

Decades of Genetic Research on Soybean mosaic virus Resistance in Soybean

This review summarizes the history and current state of the known genetic basis for soybean resistance to Soybean mosaic virus (SMV), and examines how the integration of molecular markers has been utilized in breeding for crop improvement. SVM causes yield loss and seed quality reduction in soybean based on the SMV strain and the host genotype. Understanding the molecular underpinnings of SMV-soybean interactions and the genes conferring resistance to SMV has been a focus of intense research interest for decades. Soybean reactions are classified into three main responses: resistant, necrotic, or susceptible. Significant progress has been achieved that has greatly increased the understanding of soybean germplasm diversity, differential reactions to SMV strains, genotype-strain interactions, genes/alleles conferring specific reactions, and interactions among resistance genes and alleles. Many studies that aimed to uncover the physical position of resistance genes have been published in recent decades, collectively proposing different candidate genes. The studies on SMV resistance loci revealed that the resistance genes are mainly distributed on three chromosomes. Resistance has been pyramided in various combinations for durable resistance to SMV strains. The causative genes are still elusive despite early successes in identifying resistance alleles in soybean; however, a gene at the Rsv4 locus has been well validated.

Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm

White rust, caused by Albugo occidentalis, is one of the major yield-limiting diseases of spinach (Spinacia oleracea) in some major commercial production areas, particularly in southern Texas in the United States. The use of host resistance is the most economical and environment-friendly approach to managing white rust in spinach production. The objectives of this study were to conduct a genome-wide associating study (GWAS), to identify single nucleotide polymorphism (SNP) markers associated with white rust resistance in spinach, and to perform genomic prediction (GP) to estimate the prediction accuracy (PA). A GWAS panel of 346 USDA (US Dept. of Agriculture) germplasm accessions was phenotyped for white rust resistance under field conditions and GWAS was performed using 13 235 whole-genome resequencing (WGR) generated SNPs. Nine SNPs, chr2_53 049 132, chr3_58 479 501, chr3_95 114 909, chr4_9 176 069, chr4_17 807 168, chr4_83 938 338, chr4_87 601 768, chr6_1 877 096, and chr6_31 287 118, located on chromosomes 2, 3, 4, and 6 were associated with white rust resistance in this GWAS panel. Four scenarios were tested for PA using Pearson's correlation coefficient (r) between the genomic estimation breeding value (GEBV) and the observed values: (1) different ratios between the training set and testing set (fold), (2) different GP models, (3) different SNP numbers in three different SNP sets, and (4) the use of GWAS-derived significant SNP markers. The results indicated that a 2- to 10-fold difference in the various GP models had similar, although not identical, averaged r values in each SNP set; using GWAS-derived significant SNP markers would increase PA with a high r-value up to 0.84. The SNP markers and the high PA can provide valuable information for breeders to improve spinach by marker-assisted selection (MAS) and genomic selection (GS).

Genetic architecture of salt tolerance in a Multi-Parent Advanced Generation Inter-Cross (MAGIC) cowpea population

Background

Previous reports have shown that soil salinity is a growing threat to cowpea production, and thus the need for breeding salt-tolerant cowpea cultivars. A total of 234 Multi-Parent Advanced Generation Inter-Cross (MAGIC) lines along with their 8 founders were evaluated for salt tolerance under greenhouse conditions. The objectives of this study were to evaluate salt tolerance in a multi-parent advanced generation inter-cross (MAGIC) cowpea population, to identify single nucleotide polymorphism (SNP) markers associated with salt tolerance, and to assess the accuracy of genomic selection (GS) in predicting salt tolerance, and to explore possible epistatic interactions affecting salt tolerance in cowpea. Phenotyping was validated through the use of salt-tolerant and salt-susceptible controls that were previously reported. Genome-wide association study (GWAS) was conducted using a total of 32,047 filtered SNPs. The epistatic interaction analysis was conducted using the PLINK platform.

Results

Results indicated that: (1) large variation in traits evaluated for salt tolerance was identified among the MAGIC lines, (2) a total of 7, 2, 18, 18, 3, 2, 5, 1, and 23 were associated with number of dead plants, salt injury score, leaf SPAD chlorophyll under salt treatment, relative tolerance index for leaf SPAD chlorophyll, fresh leaf biomass under salt treatment, relative tolerance index for fresh leaf biomass, relative tolerance index for fresh stem biomass, relative tolerance index for the total above-ground fresh biomass, and relative tolerance index for plant height, respectively, with overlapping SNP markers between traits, (3) candidate genes encoding for proteins involved in ion transport such as Na+/Ca2+ K+ independent exchanger and H+/oligopeptide symporter were identified, and (4) epistatic interactions were identified.

Conclusions

These results will have direct applications in breeding programs aiming at improving salt tolerance in cowpea through marker-assisted selection. To the best of our knowledge, this study was one of the earliest reports using a MAGIC population to investigate the genetic architecture of salt tolerance in cowpea.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08332-y.

Screening of cowpea (Vigna unguiculata (L.) Walp.) genotypes for waterlogging tolerance using morpho-physiological traits at early growth stage

The majority of cowpea (Vigna unguiculata (L.) Walp.) produced in the U.S. is planted shortly after the summer rains and subsequently depends on rain or artificial irrigation. Therefore, excessive precipitation and poor soil drainage will cause cowpea plants to suffer temporary waterlogging, reducing the submerged tissue's oxygen level. Although cowpea is sensitive to waterlogging, excessive moisture can induce several morpho-physiological changes with adverse impacts on yield in its early stages of development. The current study subjected 30 cowpea genotypes to 10-days of waterlogging at the seedling stage under a controlled environment. The dynamic changes of 24 morpho-physiological parameters under waterlogging and optimal water conditions were analyzed to understand cowpea's response to waterlogging. Significant waterlogging treatment, cowpea genotypes, and their interactions (P < 0.001) were observed for most of the measured parameters. The results indicated that plant height (PH), leaf area (LA), fresh (FW) and dry weight (DW) of cowpea genotypes were significantly decreased under waterlogging compared to the control treatments. Similar results were obtained for net photosynthesis (P_n), stomatal conductance (g_s), intercellular CO₂ concentration (C_i), and transpiration rate (E). However, the water use efficiency (WUE) and adventitious roots (ARs) increased linearly under waterlogging conditions. Waterlogging also declined chlorophyll fluorescence parameters except non-photochemical quenching (qN), which increased with excess soil moisture. In addition, waterlogging tolerance coefficient (WTC) and multivariate analysis (MCA) methods were used to characterize cowpea genotypes for waterlogging tolerance. Accordingly, the cowpea genotype Dagupan Pangasinan, UCR 369, and Negro were classified as waterlogging tolerant, while EpicSelect.4 and ICARDA 140071, as the most waterlogging sensitive. The cowpea genotypes and morpho-physiological traits determined from this study may be useful for genetic engineering and breeding programs that integrate cowpea waterlogging tolerance.

Identification of Candidate Genes and Genomic Selection for Seed Protein in Soybean Breeding Pipeline

Soybean is a primary meal protein for human consumption, poultry, and livestock feed. In this study, quantitative trait locus (QTL) controlling protein content was explored via genome-wide association studies (GWAS) and linkage mapping approaches based on 284 soybean accessions and 180 recombinant inbred lines (RILs), respectively, which were evaluated for protein content for 4 years. A total of 22 single nucleotide polymorphisms (SNPs) associated with protein content were detected using mixed linear model (MLM) and general linear model (GLM) methods in Tassel and 5 QTLs using Bayesian interval mapping (IM), single-trait multiple interval mapping (SMIM), single-trait composite interval mapping maximum likelihood estimation (SMLE), and single marker regression (SMR) models in Q-Gene and IciMapping. Major QTLs were detected on chromosomes 6 and 20 in both populations. The new QTL genomic region on chromosome 6 (Chr6_18844283-19315351) included 7 candidate genes and the Hap.X ^AA at the Chr6_19172961 position was associated with high protein content. Genomic selection (GS) of protein content was performed using Bayesian Lasso (BL) and ridge regression best linear unbiased prediction (rrBULP) based on all the SNPs and the SNPs significantly associated with protein content resulted from GWAS. The results showed that BL and rrBLUP performed similarly; GS accuracy was dependent on the SNP set and training population size. GS efficiency was higher for the SNPs derived from GWAS than random SNPs and reached a plateau when the number of markers was >2,000. The SNP markers identified in this study and other information were essential in establishing an efficient marker-assisted selection (MAS) and GS pipelines for improving soybean protein content.

Genetic Diversity, Structure, and Selective Sweeps in Spinacia turkestanica Associated With the Domestication of Cultivated Spinach

Genotype-by-sequencing (GBS) was used to explore the genetic diversity and structure of Spinacia turkestanica, and the selective sweeps involved in domestication of cultivated spinach, S. oleracea, from S. turkestanica. A total 7,065 single nucleotide polymorphisms (SNPs) generated for 16 Spinacia oleracea and 76 S. turkestanica accessions placed the S. oleracea accessions in one group, Q1, and the 76 S. turkestanica accessions, which originated from Central Asia, in two distinct groups, Q2 and Q3. The Q2 group shared greater genetic identity with the S. oleracea accessions, Q1, than the Q3 S. turkestanica group. Likewise, the S. oleracea Q1 group had a smaller Fst (0.008) with the Q2 group than with the Q3 group (Fst = 0.012), and a greater gene flow (Nm = 30.13) with the Q2 group than with the Q3 group (Nm = 21.83). The Q2 accessions originated primarily from Uzbekistan while the Q3 accessions originated mostly from Tajikistan. The Zarafshan Mountain Range appears to have served as a physical barrier that largely separated members of the Q2 and Q3 groups of S. turkestanica. Accessions with admixtures of Q2 and Q3 were collected primarily from lower elevations at the southern end of the Zarafshan Mountain Range in Uzbekistan. Selective sweep regions identified at 32, 49, and 52 Mb on chromosomes 1, 2, and 3, respectively, appear to have played a vital role in the domestication of S. oleracea as they are correlated with important domestication traits, including day length sensitivity for bolting (flowering). High XP-CLR scores at the 52 Mb genomic region of chromosome three suggest that a selective sweep at this region was responsible for early differentiation of S. turkestanica into two groups in Central Asia.

Genome-wide association study and genomic selection for yield and related traits in soybean

Soybean [Glycine max (L.) Merr.] is a crop of great interest worldwide. Exploring molecular approaches to increase yield genetic gain has been one of the main challenges for soybean breeders and geneticists. Agronomic traits such as maturity, plant height, and seed weight have been found to contribute to yield. In this study, a total of 250 soybean accessions were genotyped with 10,259 high-quality SNPs postulated from genotyping by sequencing (GBS) and evaluated for grain yield, maturity, plant height, and seed weight over three years. A genome-wide association study (GWAS) was performed using a Bayesian Information and Linkage Disequilibrium Iteratively Nested Keyway (BLINK) model. Genomic selection (GS) was evaluated using a ridge regression best linear unbiased predictor (rrBLUP) model. The results revealed that 20, 31, 37, and 23 SNPs were significantly associated with maturity, plant height, seed weight, and yield, respectively; Many SNPs were mapped to previously described maturity and plant height loci (E2, E4, and Dt1) and a new plant height locus was mapped to chromosome 20. Candidate genes were found in the vicinity of the two SNPs with the highest significant levels associated with yield, maturity, plant height, seed weight, respectively. A 11.5-Mb region of chromosome 10 was associated with both yield and seed weight. Overall, the accuracy of GS was dependent on the trait, year, and population structure, and high accuracy indicates that these agronomic traits can be selected in molecular breeding through GS. The SNP markers identified in this study can be used to improve yield and agronomic traits through the marker-assisted selection and GS in breeding programs.

Population genetics and genome-wide association studies provide insights into the influence of selective breeding on genetic variation in lettuce

Genetic diversity is an important resource in crop breeding to improve cultivars with desirable traits. Selective breeding can lead to a reduction of genetic diversity. However, our understanding on this subject remains limited in lettuce (Lactuca sativa L.). Genotyping-by-sequencing (GBS) can provide a reduced version of the genome as a cost-effective method to identify genetic variants across the genome. We genotyped a diverse set of 441 lettuce accessions using the GBS method. Phylogenetic and population genetic analyses indicated substantial genetic divergence among four horticultural types of lettuce: butterhead, crisphead, leaf, and romaine. Genetic-diversity estimates between and within the four types indicated that the crisphead type was the most differentiated from other types, whereas its population was the most homogenous with the slowest linkage disequilibrium (LD) decay among the four types. These results suggested that crisphead lettuces had relatively less genetic variation across the genome as well as low gene flow from other types. We identified putative selective sweep regions that showed low genetic variation in the crisphead type. Genome-wide association study (GWAS) and quantitative trait loci (QTL) analyses provided evidence that these genomic regions were, in part, associated with delayed bolting, implicating the positive selection of delayed bolting in reducing variation. Our findings enhance the current understanding of genetic diversity and the impacts of selective breeding on patterning genetic variation in lettuce.

High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach

Background

Downy mildew, the most devastating disease of spinach (Spinacia oleracea L.), is caused by the oomycete Peronospora effusa [=P. farinosa f. sp. spinaciae]. The P. effusa shows race specificities to the resistant host and comprises 19 reported races and many novel isolates. Sixteen new P. effusa races were identified during the past three decades, and the new pathogen races are continually overcoming the genetic resistances used in commercial cultivars. A spinach breeding population derived from the cross between cultivars Whale and Lazio was inoculated with P. effusa race 16 in an environment-controlled facility; disease response was recorded and genotyped using genotyping by sequencing (GBS). The main objective of this study was to identify resistance-associated single nucleotide polymorphism (SNP) markers from the cultivar Whale against the P. effusa race 16.

Results

Association analysis conducted using GBS markers identified six significant SNPs (S3_658,306, S3_692697, S3_1050601, S3_1227787, S3_1227802, S3_1231197). The downy mildew resistance locus from cultivar Whale was mapped to a 0.57 Mb region on chromosome 3, including four disease resistance candidate genes (Spo12736, Spo12784, Spo12908, and Spo12821) within 2.69–11.28 Kb of the peak SNP.

Conclusions

Genomewide association analysis approach was used to map the P. effusa race 16 resistance loci and identify associated SNP markers and the candidate genes. The results from this study could be valuable in understanding the genetic basis of downy mildew resistance, and the SNP marker will be useful in spinach breeding to select resistant lines.

Genome-Wide Association Study and Genomic Prediction for Soybean Cyst Nematode Resistance in USDA Common Bean (Phaseolus vulgaris) Core Collection

Soybean cyst nematode (SCN, Heterodera glycines) has become the major yield-limiting biological factor in soybean production. Common bean is also a good host of SCN, and its production is challenged by this emerging pest in many regions such as the upper Midwest USA. The use of host genetic resistance has been the most effective and environmentally friendly method to manage SCN. The objectives of this study were to evaluate the SCN resistance in the USDA common bean core collection and conduct a genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers with SCN resistance. A total of 315 accessions of the USDA common bean core collection were evaluated for resistance to SCN HG Type 0 (race 6). The common bean core set was genotyped with the BARCBean6K_3 Infinium BeadChips, consisting of 4,654 SNPs. Results showed that 15 accessions were resistant to SCN with a Female Index (FI) at 4.8 to 9.4, and 62 accessions were moderately resistant (10 < FI < 30) to HG Type 0. The association study showed that 11 SNP markers, located on chromosomes Pv04, 07, 09, and 11, were strongly associated with resistance to HG Type 0. GWAS was also conducted for resistance to HG Type 2.5.7 and HG Type 1.2.3.5.6.7 based on the public dataset (N = 276), consisting of a diverse set of common bean accessions genotyped with the BARCBean6K_3 chip. Six SNPs associated with HG Type 2.5.7 resistance on Pv 01, 02, 03, and 07, and 12 SNPs with HG Type 1.2.3.5.6.7 resistance on Pv 01, 03, 06, 07, 09, 10, and 11 were detected. The accuracy of genomic prediction (GP) was 0.36 to 0.49 for resistance to the three SCN HG types, indicating that genomic selection (GS) of SCN resistance is feasible. This study provides basic information for developing SCN-resistant common bean cultivars, using the USDA core germ plasm accessions. The SNP markers can be used in molecular breeding in common beans through marker-assisted selection (MAS) and GS.

Learning from reproducing computational results: introducing three principles and the Reproduction Package

We carry out efforts to reproduce computational results for seven published articles and identify barriers to computational reproducibility. We then derive three principles to guide the practice and dissemination of reproducible computational research: (i) Provide transparency regarding how computational results are produced; (ii) When writing and releasing research software, aim for ease of (re-)executability; (iii) Make any code upon which the results rely as deterministic as possible. We then exemplify these three principles with 12 specific guidelines for their implementation in practice. We illustrate the three principles of reproducible research with a series of vignettes from our experimental reproducibility work. We define a novel Reproduction Package, a formalism that specifies a structured way to share computational research artifacts that implements the guidelines generated from our reproduction efforts to allow others to build, reproduce and extend computational science. We make our reproduction efforts in this paper publicly available as exemplar Reproduction Packages. This article is part of the theme issue 'Reliability and reproducibility in computational science: implementing verification, validation and uncertainty quantification in silico'.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

The novel Mechanical Ventilator Milano for the COVID-19 pandemic

This paper presents the Mechanical Ventilator Milano (MVM), a novel intensive therapy mechanical ventilator designed for rapid, large-scale, low-cost production for the COVID-19 pandemic. Free of moving mechanical parts and requiring only a source of compressed oxygen and medical air to operate, the MVM is designed to support the long-term invasive ventilation often required for COVID-19 patients and operates in pressure-regulated ventilation modes, which minimize the risk of furthering lung trauma. The MVM was extensively tested against ISO standards in the laboratory using a breathing simulator, with good agreement between input and measured breathing parameters and performing correctly in response to fault conditions and stability tests. The MVM has obtained Emergency Use Authorization by U.S. Food and Drug Administration (FDA) for use in healthcare settings during the COVID-19 pandemic and Health Canada Medical Device Authorization for Importation or Sale, under Interim Order for Use in Relation to COVID-19. Following these certifications, mass production is ongoing and distribution is under way in several countries. The MVM was designed, tested, prepared for certification, and mass produced in the space of a few months by a unique collaboration of respiratory healthcare professionals and experimental physicists, working with industrial partners, and is an excellent ventilator candidate for this pandemic anywhere in the world.

Loci discovery, network-guided approach, and genomic prediction for drought tolerance index in a multi-parent advanced generation intercross (MAGIC) cowpea population

Cowpea is a nutrient-dense legume that significantly contributes to the population's diet in sub-Saharan Africa and other regions of the world. Improving cowpea cultivars to be more resilient to abiotic stress such as drought would be of great importance. The use of a multi-parent advanced generation intercross (MAGIC) population has been shown to be efficient in increasing the frequency of rare alleles that could be associated with important agricultural traits. In addition, drought tolerance index has been reported to be a reliable parameter for assessing crop tolerance to water-deficit conditions. Therefore, the objectives of this study were to evaluate the drought tolerance index for plant growth habit, plant maturity, flowering time, 100-seed weight, and grain yield in a MAGIC cowpea population, to conduct genome-wide association study (GWAS) and identify single nucleotide polymorphism (SNP) markers associated with the drought tolerance indices, to investigate the potential relationship existing between the significant loci associated with the drought tolerance indices, and to conduct genomic selection (GS). These analyses were performed using the existing phenotypic and genotypic data published for the MAGIC population which consisted of 305 F8 recombinant inbred lines (RILs) developed at University of California, Riverside. The results indicated that: (1) large variation in drought tolerance indices existed among the cowpea genotypes, (2) a total of 14, 18, 5, 5, and 35 SNPs were associated with plant growth habit change due to drought stress, and drought tolerance indices for maturity, flowering time, 100-seed weight, and grain yield, respectively, (3) the network-guided approach revealed clear interactions between the loci associated with the drought tolerance traits, and (4) the GS accuracy varied from low to moderate. These results could be applied to improve drought tolerance in cowpea through marker-assisted selection (MAS) and genomic selection (GS). To the best of our knowledge, this is the first report on marker loci associated with drought tolerance indices in cowpea.

Genome Wide Association Studies in Multiple Spinach Breeding Populations Refine Downy Mildew Race 13 Resistance Genes

Downy mildew, caused by the oomycete Peronospora effusa, is the most economically important disease on spinach. Fourteen new races of P. effusa have been identified in the last three decades. The frequent emergence of new races of P. effusa continually overcome the genetic resistance to the pathogen. The objectives of this research were to more clearly map the downy mildew resistance locus RPF1 in spinach, to identify single nucleotide polymorphism (SNP) markers associated with the resistance, and to refine the candidate genes responsible for the resistance. Progeny from populations generated from crosses of cultivars resistant (due to RPF1) to race 13 of P. effusa (Swan, T-Bird, Squirrel, and Tonga) with race 13 susceptible cultivars (Whale and Polka) were inoculated and the downy mildew disease response determined. Association analysis was performed in TASSEL, GAPIT, PLINK, and GENESIS programs using SNP markers identified from genotyping by sequencing (GBS). Association analysis mapped the race 13 resistance loci (RPF1) to positions 0.39, 0.69, 0.94-0.98, and 1.2 Mb of chromosome 3. The associated SNPs were within 1-7 kb of the disease resistance genes Spo12784, Spo12719, Spo12905, and Spo12821, and 11-18 Kb from Spo12903. This study extended our understanding of the genetic basis of downy mildew resistance in spinach and provided the most promising candidate genes Spo12784 and Spo12903 near the RPF1 locus, to pursue functional validation. The SNP markers may be used to select for the resistant lines to improve genetic resistance against the downy mildew pathogen and in developing durably resistant cultivars.

Proteomics Analysis Indicates Greater Abundance of Proteins Involved in Major Metabolic Pathways in Lactuca sativa cv. Salinas than Lactuca serriola Accession US96UC23

Lettuce (Lactuca sativa), cultivated mainly for its edible leaves and stems, is an important vegetable crop worldwide. Genomes of cultivated lettuce (L. sativa cv. Salinas) and its wild relative L. serriola accession US96UC23 are sequenced, but a clear understanding of the genetic basis for divergence in phenotypes of the two species is lacking. Tandem mass tag (TMT) based mass spectrometry is used to quantitatively compare protein levels between these two species. Four-day old seedlings is transplanted into 500 mL pots filled with soil. Plants are grown for 8 weeks under 250 µmol m-2 sec-1 continuous light, 20 °C and relative humidity between 50-70%. Leaf discs (1 cm diameter) from three individuals per biological replicate are analyzed. A total of 3000 proteins are identified, of which the levels of 650 are significantly different between 'Salinas' and US96UC23. Pathway analysis indicated a higher flux of carbon in 'Salinas' than US96UC23. Many essential metabolic pathways such as tetrapyrrole metabolism and fatty acid biosynthesis are upregulated in 'Salinas' compared with US96UC23. This study provides a reference proteome for researchers interested in understanding lettuce biology and improving traits for cultivation.

Genome-wide association study and genomic selection for tolerance of soybean biomass to soybean cyst nematode infestation

Soybean cyst nematode (SCN), Heterodera glycines Ichinohe, is one of the most devastating pathogens affecting soybean production in the U.S. and worldwide. The use of SCN-resistant soybean cultivars is one of the most affordable strategies to cope with SCN infestation. Because of the limited sources of SCN resistance and changes in SCN virulence phenotypes, host resistance in current cultivars has increasingly been overcome by the pathogen. Host tolerance has been recognized as an additional tool to manage the SCN. The objectives of this study were to conduct a genome-wide association study (GWAS), to identify single nucleotide polymorphism (SNP) markers, and to perform a genomic selection (GS) study for SCN tolerance in soybean based on reduction in biomass. A total of 234 soybean genotypes (lines) were evaluated for their tolerance to SCN in greenhouse using four replicates. The tolerance index (TI = 100 × Biomass of a line in SCN infested / Biomass of the line without SCN) was used as phenotypic data of SCN tolerance. GWAS was conducted using a total of 3,782 high quality SNPs. GS was performed based upon the whole set of SNPs and the GWAS-derived SNPs, respectively. Results showed that (1) a large variation in soybean TI to SCN infection among the soybean genotypes was identified; (2) a total of 35, 21, and 6 SNPs were found to be associated with SCN tolerance using the models SMR, GLM (PCA), and MLM (PCA+K) with 6 SNPs overlapping between models; (3) GS accuracy was SNP set-, model-, and training population size-dependent; and (4) genes around Glyma.06G134900, Glyma.15G097500.1, Glyma.15G100900.3, Glyma.15G105400, Glyma.15G107200, and Glyma.19G121200.1 (Table 4). Glyma.06G134900, Glyma.15G097500.1, Glyma.15G100900.3, Glyma.15G105400, and Glyma.19G121200.1 are best candidates. To the best of our knowledge, this is the first report highlighting SNP markers associated with tolerance index based on biomass reduction under SCN infestation in soybean. This research opens a new approach to use SCN tolerance in soybean breeding and the SNP markers will provide a tool for breeders to select for SCN tolerance.

Transcriptional analyses of differential cultivars during resistant and susceptible interactions with Peronospora effusa, the causal agent of spinach downy mildew

Downy mildew of spinach is caused by the obligate oomycete pathogen, Peronospora effusa. The disease causes significant economic losses, especially in the organic sector of the industry where the use of synthetic fungicides is not permitted for disease control. New pathotypes of this pathogen are increasingly reported which are capable of breaking resistance. In this study, we took advantage of new spinach genome resources to conduct RNA-seq analyses of transcriptomic changes in leaf tissue of resistant and susceptible spinach cultivars Solomon and Viroflay, respectively, at an early stage of pathogen establishment (48 hours post inoculation, hpi) to a late stage of symptom expression and pathogen sporulation (168 hpi). Fold change differences in gene expression were recorded between the two cultivars to identify candidate genes for resistance. In Solomon, the hypersensitive inducible genes such as pathogenesis-related gene PR-1, glutathione-S-transferase, phospholipid hydroperoxide glutathione peroxidase and peroxidase were significantly up-regulated uniquely at 48 hpi and genes involved in zinc finger CCCH protein, glycosyltransferase, 1-aminocyclopropane-1-carboxylate oxidase homologs, receptor-like protein kinases were expressed at 48 hpi through 168 hpi. The types of genes significantly up-regulated in Solomon in response to the pathogen suggests that salicylic acid and ethylene signaling pathways mediate resistance. Furthermore, many genes involved in the flavonoid and phenylpropanoid pathways were highly expressed in Viroflay compared to Solomon at 168 hpi. As anticipated, an abundance of significantly down-regulated genes was apparent at 168 hpi, reflecting symptom development and sporulation in cultivar Viroflay, but not at 48 hpi. In the pathogen, genes encoding RxLR-type effectors were expressed during early colonization of cultivar Viroflay while crinkler-type effector genes were expressed at the late stage of the colonization. Our results provide insights on gene expression in resistant and susceptible spinach-P. effusa interactions, which can guide future studies to assess candidate genes necessary for downy mildew resistance in spinach.

Genome-wide identification and expression analysis of the CBF/DREB1 gene family in lettuce

The C-repeat binding factor (CBF)/dehydration-responsive element binding (DREB1) proteins play a prominent role in freezing tolerance and are highly conserved in higher plants. Here we performed a genome-wide search of the CBF/DREB1 gene family in lettuce (Lactuca sativa L.) and identified 14 members of the family with one member gene containing a non-sense mutation within the AP2 DNA-binding domain. A comprehensive phylogenetic analysis of the CBF/DREB1 family members in 20 plant species from the Asterid or Rosid clade provided evidence that tandem duplication played an important role in the expansion of the CBF/DREB1 family. Expression analysis showed that twelve of the lettuce CBF genes were responsive to low temperature (4 °C), and that three and six of them could also be responsive to salt and heat stresses, respectively. Unlike Arabidopsis thaliana whose members of the CBF/DREB1 family respond only to a particular stress, lettuce CBFs provide wider protection from combinations of abiotic stresses. A global transcriptome analysis revealed distinctive temporal expression patterns among the cold-regulated genes in lettuce plants exposed to low temperature. Genes induced throughout the cold treatment are enriched in functions associated with protection from UV and high-light intensity and the genes suppressed after 7 days of cold exposure are enriched in photosynthesis-associated functions. These results provide insight into the molecular evolutionary properties of the CBF/DREB1 gene family in lettuce and a reference for genetic improvement of the lettuce response to cold acclimation.

Prediction of venous thromboembolism using clinical and serum biomarker data from a military cohort of trauma patients

INTRODUCTION

Venous thromboembolism (VTE) is a frequent complication of trauma associated with high mortality and morbidity. Clinicians lack appropriate tools for stratifying trauma patients for VTE, thus have yet to be able to predict when to intervene. We aimed to compare random forest (RF) and logistic regression (LR) predictive modelling for VTE using (1) clinical measures alone, (2) serum biomarkers alone and (3) clinical measures plus serum biomarkers.

METHODS

Data were collected from 73 military casualties with at least one extremity wound and prospectively enrolled in an observational study between 2007 and 2012. Clinical and serum cytokine data were collected. Modelling was performed with RF and LR based on the presence or absence of deep vein thrombosis (DVT) and/or pulmonary embolism (PE). For comparison, LR was also performed on the final variables from the RF model. Sensitivity/specificity and area under the curve (AUC) were reported.

RESULTS

Of the 73 patients (median Injury Severity Score=16), nine (12.3%) developed VTE, four (5.5%) with DVT, four (5.5%) with PE, and one (1.4%) with both DVT and PE. In all sets of predictive models, RF outperformed LR. The best RF model generated with clinical and serum biomarkers included five variables (interleukin-15, monokine induced by gamma, vascular endothelial growth factor, total blood products at resuscitation and presence of soft tissue injury) and had an AUC of 0.946, sensitivity of 0.992 and specificity of 0.838.

CONCLUSIONS

VTE may be predicted by clinical and molecular biomarkers in trauma patients. This will allow the development of clinical decision support tools which can help inform the management of high-risk patients for VTE.

Transcript profiling for regulation of sweet potato skin color in Sushu8 and its mutant Zhengshu20

Sweet potato [Ipomoea batatas (L.) Lam.] (2n = 6x = 90) is an economic important autopolyploid species and its varieties differ regarding storage root skin and flesh colors. Two sweet potato genetic lines, Sushu8 (with red skin) and its mutant Zhengshu20, which produced different colored storage roots, were used in this study. The total flavonoid, carotenoid, and anthocyanin contents of the two lines were analyzed and revealed that anthocyanin was primarily responsible for the skin color difference. In addition, the early storage root expanding stage was the key period for anthocyanin accumulation in Sushu8. A total of 24 samples, including the skins of the fibrous root and the storage root at the early and middle expanding stages as well as the flesh of the storage root at the middle expanding stage, were analyzed based on differentially expressed genes identified by transcriptome sequencing and a weighted gene co-expression network analysis. Two gene modules highly related with the regulation of sweet potato skin color through stress responses as well as starch synthesis and glucose metabolism were identified. Furthermore, the WRKY75 transcription factor gene, fructose-bisphosphate aldolase 2 gene, and other DEGs highly related to the regulation of anthocyanin metabolism were enriched in the brown and green modules.

Genome-wide association study and genomic selection for soybean chlorophyll content associated with soybean cyst nematode tolerance

BACKGROUND

Soybean cyst nematode (SCN), Heterodera glycines Ichinohe, has been one of the most devastating pathogens affecting soybean production. In the United States alone, SCN damage accounted for more than $1 billion loss annually. With a narrow genetic background of the currently available SCN-resistant commercial cultivars, high risk of resistance breakdown can occur. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify QTL, SNP markers, and candidate genes associated with soybean leaf chlorophyll content tolerance to SCN infection, and to carry out a genomic selection (GS) study for the chlorophyll content tolerance.

RESULTS

A total of 172 soybean genotypes were evaluated for the effect of SCN HG Type 1.2.3.5.6.7 (race 4) on soybean leaf chlorophyll. The soybean lines were genotyped using a total of 4089 filtered and high-quality SNPs. Results showed that (1) a large variation in SCN tolerance based on leaf chlorophyll content indices (CCI); (2) a total of 22, 14, and 16 SNPs associated with CCI of non-SCN-infected plants, SCN-infected plants, and reduction of CCI SCN, respectively; (3) a new locus of chlorophyll content tolerance to SCN mapped on chromosome 3; (4) candidate genes encoding for Leucine-rich repeat protein, plant hormone signaling molecules, and biomolecule transporters; and (5) an average GS accuracy ranging from 0.31 to 0.46 with all SNPs and varying from 0.55 to 0.76 when GWAS-derived SNP markers were used across five models. This study demonstrated the potential of using genome-wide selection to breed chlorophyll-content-tolerant soybean for managing SCN.

CONCLUSIONS

In this study, soybean accessions with higher CCI under SCN infestation, and molecular markers associated with chlorophyll content related to SCN were identified. In addition, a total of 15 candidate genes associated with chlorophyll content tolerance to SCN in soybean were also identified. These candidate genes will lead to a better understanding of the molecular mechanisms that control chlorophyll content tolerance to SCN in soybean. Genomic selection analysis of chlorophyll content tolerance to SCN showed that using significant SNPs obtained from GWAS could provide better GS accuracy.

Genome Wide Association Study and Genomic Selection of Amino Acid Concentrations in Soybean Seeds

Soybean is a major source of protein for human consumption and animal feed. Releasing new cultivars with high nutritional value is one of the major goals in soybean breeding. To achieve this goal, genome-wide association studies of seed amino acid contents were conducted based on 249 soybean accessions from China, US, Japan, and South Korea. The accessions were evaluated for 15 amino acids and genotyped by sequencing. Significant genetic variation was observed for amino acids among the accessions. Among the 231 single nucleotide polymorphisms (SNPs) significantly associated with variations in amino acid contents, fifteen SNPs localized near 14 candidate genes involving in amino acid metabolism. The amino acids were classified into two groups with five in one group and seven amino acids in the other. Correlation coefficients among the amino acids within each group were high and positive, but the correlation coefficients of amino acids between the two groups were negative. Twenty-five SNP markers associated with multiple amino acids can be used to simultaneously improve multi-amino acid concentration in soybean. Genomic selection analysis of amino acid concentration showed that selection efficiency of amino acids based on the markers significantly associated with all 15 amino acids was higher than that based on random markers or markers only associated with individual amino acid. The identified markers could facilitate selection of soybean varieties with improved seed quality.

Genetic Diversity Analysis of Olive Germplasm (Olea europaea L.) With Genotyping-by-Sequencing Technology

Olive (Olea europaea L.) is a very important edible oil crop and has been cultivated for about 4,000 years in the Mediterranean area. Due to its nutritional and economic importance, researches on germplasm characterization received extensive attention. In this study, using the genotyping-by-sequencing (GBS) technology, we carried out genetic diversity analysis on 57 olive cultivars with different geographical origins. In total, 73,482 high-quality single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) > 5%, call rate > 50%, and heterozygosity rate < 10% were obtained at the whole genome level. Genetic structure and phylogenetic analysis showed that the 57 olive cultivars could be classified into two groups (Group I and Group II). No clear geographical distributions of cultivars were observed generally between the two groups. The average nucleotide diversities (π) specific for Group I and Group II were 0.317 and 0.305. The fixation index (F_ST) between Group I and Group II was 0.033. In Group II, cultivars could be further divided into two subgroups (Group IIa and Group IIb), which seem to be associated with their fruit sizes. The five Chinese-bred cultivars were all clustered in Group II, showing a closer genetic relationship with those from the central Mediterranean region and limited genetic background. It is therefore necessary for Chinese olive breeding programs to incorporate other genetic basis by utilizing germplasm from the other regions particularly from the east Mediterranean region as breeding parents. The results showed that GBS is an effective marker choice for cultivar characterization and genetic diversity analysis in olive and will help us better understand the genetic backgrounds of the crop.

Spinach Downy Mildew: Advances in Our Understanding of the Disease Cycle and Prospects for Disease Management

Downy mildew on spinach is caused by Peronospora effusa, an oomycete pathogen that poses a challenge to spinach production worldwide, especially in organic production. Following infection, P. effusa produces abundant amounts of asexual sporangia. Sporangia become windborne and initiate new infections locally or distantly, leading to widespread epidemics. Oospores produced from the union of opposite mating types have been observed within infected leaves and seeds and may remain viable for many years. Sexual reproduction increases the genetic diversity of P. effusa through sexual recombination, and thus, the movement of oospores on seed has likely fueled the rapid explosion of new pathotypes in different regions of the world over the past 20 years. This review summarizes recent advances in spinach downy mildew research, especially in light of the findings of oospores in contemporary commercial spinach seed lots as well as their germination. Knowledge of the role of the oospores and other aspects of the disease cycle can directly translate into new and effective disease management strategies.

CD56bright natural killer cells induce HBsAg reduction via cytolysis and cccDNA decay in long-term entecavir-treated patients switching to peginterferon alfa-2a

HBV surface antigen (HBsAg) reduction is well observed in chronic hepatitis B (CHB) patients treated with pegylated interferon alpha-2a (PegIFNα). However, the mechanism of HBsAg suppression has not been fully elucidated. Twenty-seven of 55 entecavir-treated CHB e antigen positive patients were switched to PegIFNα treatment (Group A) whereas 28 patients continued entecavir treatment (Group B). The percentage or absolute number of CD56^bright /CD56^dim NK cells, expression of receptors and cytokines were evaluated by flow cytometry for 48 weeks and correlated with treatment efficacy. In vitro, purified NK cells were co-cultured with HepAD38 cells for measurement of HBsAg, apoptosis and covalently closed circular DNA (cccDNA). In association with a reduction of HBsAg, the percentage and absolute number of CD56^bright NK cells was significantly elevated in patients in group A, especially in Virologic Responders (VRs, HBsAg decreased). Furthermore, the percentage of NKp30⁺ , NKp46⁺ , TRAIL⁺ , TNF-α⁺ and IFNγ⁺ CD56^bright NK cells were significantly expanded in Group A, which were positively correlated with the decline of HBsAg at week 48. In vitro, peripheral NK cells from Group A induced a decline of HBsAg in comparison with NK cells from Group B which was significantly inhibited by anti-TRAIL, anti-TNF-α and anti-IFNγ antibodies. Furthermore, apoptosis of HepAD38 cells and levels of cccDNA, were significantly reduced by TRAIL⁺ and TNF-α⁺ /IFNγ⁺ NK cells from Group A, respectively. A functional restoration of CD56^bright NK cells in entecavir-treated patients who were switched to PegIFNα contributes to HBsAg and cccDNA clearance through TRAIL-induced cytolysis and TNF-α/IFNγ-mediated noncytolytic pathways.

Effect of local scrotal heating on the expression of tight junction-associated molecule Occludin in boar testes

The aim of this study was to determine whether local scrotal heating (42°C, for 1 hr) had an effect on the expression of tight junction (TJ)-associated molecule Occludin in boar testes. Adult boars (Landrace, n = 6) were used and randomly divided into two groups (n = 3 each). Three boars were given local scrotal exposure to 42°C for approximately 1 h with a home-made electric blanket of controlled temperature as local scrotal heating group, the other three boars received no heat treatment and were left at standard room temperature as control group. After 6 hr, all boars were castrated and the testes were harvested. qRT-PCR, Western blotting and immunohistochemistry were used to explore the expression and localization of Occludin. qRT-PCR and Western blotting showed that the protein and mRNA levels of Occludin significantly decreased in local scrotal heating group as compared to the control. Furthermore, immunoreactivity staining of Occludin was localized at the sites of the blood-testis barrier (BTB) and formed an almost consecutive and strong immunoreactivity strand in the control, while Occludin was limited to Sertoli cells (SCs) and no obvious immunoreactivity strand was present in local scrotal heating group. These data indicated that local scrotal heating decreased the expression of TJ-associated molecule Occludin, which may be involved in heat-induced spermatogenesis damage.