A Simulation Study to Reveal the Epidemiology and Aerosol Transmission Characteristics of Botrytis cinerea in Grape Greenhouses.

Most previously studies had considered that plant fungal disease spread widely and quickly by airborne fungi spore. However, little is known about the release dynamics, aerodynamic diameter, and pathogenicity threshold of fungi spore in air of the greenhouse environment. Grape gray mold is caused by Botrytis cinerea; the disease spreads in greenhouses by spores in the air and the spore attaches to the leaf and infects plant through the orifice. In this study, 120 µmol/L propidium monoazide (PMA) were suitable for treatment and quantitation viable spore by quantitative real-time PCR, with a limit detection of 8 spores/mL in spore suspension. In total, 93 strains of B. cinerea with high pathogenicity were isolated and identified from the air samples of grapevines greenhouses by a portable sampler. The particle size of B. cinerea aerosol ranged predominately from 0.65-3.3 µm, accounting for 71.77% of the total amount. The B. cinerea spore aerosols were infective to healthy grape plants, with the lowest concentration that could cause disease being 42 spores/m3. Botrytis cinerea spores collected form six greenhouse in Shandong Province were quantified by PMA-qPCR, with a higher concentration (1182.89 spores/m3) in May and June and a lower concentration in July and August (6.30 spores/m3). This study suggested that spore dispersal in aerosol is an important route for the epidemiology of plant fungal disease, and these data will contribute to the development of new strategies for the effective alleviation and control of plant diseases.

First report of Coniella vitis causing white rot on Virginia creeper (Parthenocissus quinquefolia [L.] Planch.) in China.

Virginia creeper (Parthenocissus quinquefolia [L.] Planch.) belongs to the genus of Parthenocissus and Vitaceae family, which is very common in vineyards and where wild grape occurs (Bergh et al., 2011). In September of 2021, a severe white rot disease was observed on Virginia creeper around the vineyard of wine grapevine (Cabernet Sauvignon) located in Penglai city (37º 75'38" N, 120º 84'28" E), Shandong province of China. The disease incidence was about 75%, and infected leaf of Virginia creeper exhibited irregular necrotic lesion with brown center, and most lesion occurred on leaf margin, black pycnidia were also observed on the infected leaf at the late stage of infection. To determine the causal agent, symptomatic leaves with typical lesions were cut into small pieces (5 mm × 3 mm), surface sterilized with 75% ethanol for 1 min, followed by three times rinsed in sterile water. Leaf sections were plated onto potato dextrose agar (PDA) medium and incubated at 28°C for 3 days. Totally, five isolates (referred to as JD01, JD07, JD09, JD12 and JD16) were collected and transferred on to fresh PDA medium for incubation at 28°C. Seven days later, colonies on PDA plates had crenulated edges with concentric rings, the upper surface of colonies was mostly flat and white with many pycnidia. The conidia were hyaline at immature and became brown later, spherical or ellipsoid, aseptate, and 7.92 ± 1.20 µm × 5.18 ± 0.61 µm (n=50), length : width ratio is nearly 2. Morphologically, the isolates were identified as Coniella vitis (Chethana et al., 2017). Further to confirm the fungal species, the internal transcribed spacer region (ITS) of the ribosomal RNA gene, large subunit rRNA gene (LSU) and the translation elongation factor 1-alpaha gene (TEF1-α) were amplified using primers ITS1/ ITS4, LR7/ LROR, and TEF1- 728F/ TEF1- 986R (Chethana et al., 2017; Raudabaugh et al., 2018). The amplification products were sequenced and deposited in GenBank database. The sequences were compared to type sequences in GenBank. The results showed that ITS (GenBank accession numbers ON329769, ON329770, ON329771, ON329772 and ON329773), LSU (ON358423，ON358424, ON358425, ON358426 and ON358427) and TEF (ON297671, ON229071, ON229072, ON229073 and ON297672) sequences of the five isolates were 99.66%, 96.90% and 98.79% identical with the sequences data from C. vitis isolates in GeneBank (MFLUCC 18-0093, JZB3700020 and MFLUCC 18-0093, respectively). Furthermore, concatenated sequences of the three genes (ITS, LSU and TEF) were used to conduct a phylogenetic tree using maximum likehood MEGA-X (Raudabaugh et al., 2018). The phylogenetic analysis showed that the five isolates (JD01, JD07, JD09, JD12 and JD16) belong to C. vitis clade among the 41strains of Coniella spp. In the pathogenicity tests, detached leaves of Virginia creeper (1-year-old) were inoculated with mycelia plugs (5 mm diameter) (form 3-day-old of isolate JD07 culture), and control were inoculated with PDA plugs (5 mm diameter). Virginia creeper live plants (1-year-old) were inoculated with conidial suspension (2.5×106 spores/ml) of the isolate JD07 of one week old, and control plants were inoculated with sterile water. All treated Virginia creeper plants (detached leaves) were placed in a greenhouse maintained at 28°C and 95% relative humidity. Virginia creeper plants (detached leaves) inoculated with the conidial suspension (fungal mycelia) had brown lesion on leaves, the disease symptoms were similar to those observed in field. No such symptoms were observed on control plants (detached leaves). The pathogen was reisolated from inoculated Virginia creeper plants and re-identified, thus fulfilling Koch's postulates. C. vitis had been reported to cause grape white rot in China (Chethana et al., 2017). Virginia creeper, as an excellent host of C. vitis, will increase the transmission risk of the pathogens. To our knowledge, this is the first report of C. vitis causing white rot on Virginia creeper, and this finding will provide useful information for developing effective control strategies for white rot disease.

Insights into the Genomic Regions and Candidate Genes of Senescence-Related Traits in Upland Cotton via GWAS.

Senescence is the last stage of plant development and is controlled by both internal and external factors. Premature senescence significantly affects the yield and quality of cotton. However, the genetic architecture underlying cotton senescence remains unclear. In this study, genome-wide association studies (GWAS) were performed based on 3,015,002 high-quality SNP markers from the resequencing data of 355 upland cotton accessions to detect genomic regions for cotton senescence. A total of 977 candidate genes within 55 senescence-related genomic regions (SGRs), SGR1-SGR55, were predicted. Gene ontology (GO) analysis of candidate genes revealed that a set of biological processes was enriched, such as salt stress, ethylene processes, and leaf senescence. Furthermore, in the leaf senescence GO term, one candidate gene was focused on: Gohir.A12G270900 (GhMKK9), located in SGR36, which encodes a protein of the MAP kinase kinase family. Quantitative real-time PCR (qRT-PCR) analysis showed that GhMKK9 was up-regulated in old cotton leaves. Overexpression of GhMKK9 in Arabidopsis accelerated natural leaf senescence. Virus-induced gene silencing (VIGS) of GhMKK9 in cotton increased drought tolerance. These results suggest that GhMKK9 is a positive regulator and might be involved in drought-induced senescence in cotton. The results provide new insights into the genetic basis of cotton senescence and will be useful for improving cotton breeding in the future.

Characterization and Functional Analysis of GhNAC82, A NAM Domain Gene, Coordinates the Leaf Senescence in Upland Cotton (Gossypium hirsutum L.).

In the process of growth and development, cotton exhibits premature senescence under various abiotic stresses, impairing yield and fiber quality. NAC (NAM, ATAF1,2, and CUC2) protein widely distributed in land plants, play the critical role in responding to abiotic stress and regulating leaf senescence. We have identified and functional analyzed a NAM domain gene GhNAC82 in upland cotton, it was located on the A11 chromosome 4,921,702 to 4,922,748 bp, only containing one exon. The spatio-temporal expression pattern analysis revealed that it was highly expressed in root, torus, ovule and fiber development stage. The results of qRT-PCR validated that GhNAC82 negatively regulated by salt stress, drought stress, H2O2 stress, IAA treatment, and ethylene treatment, positively regulated by the ABA and MeJA treatment. Moreover, heterologous overexpression of GhNAC82 results in leaf premature senescence and delays root system development in Arabidopsis thaliana. The phenotype of delayed-senescence was performed after silencing GhNAC82 by VIGS in premature cotton. Taken together, GhNAC82 was involved in different abiotic stress pathways and play important roles in negatively regulating leaf premature senescence.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton.

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

Uncovering Novel Genomic Regions and Candidate Genes for Senescence-Related Traits by Genome-Wide Association Studies in Upland Cotton (Gossypium hirsutum L.).

Senescence in plants is a complex trait, which is controlled by both genetic and environmental factors and can affect the yield and quality of cotton. However, the genetic basis of cotton senescence remains relatively unknown. In this study, we reported genome-wide association studies (GWAS) based on 185 accessions of upland cotton and 26,999 high-quality single-nucleotide polymorphisms (SNPs) to reveal the genetic basis of cotton senescence. To determine cotton senescence, we evaluated eight traits/indices. Our results revealed a high positive correlation (r>0.5) among SPAD value 20 days after topping (SPAD20d), relative difference of SPAD (RSPAD), nodes above white flower on topping day (NAWF0d), nodes above white flower 7 days after topping (NAWF7d), and number of open bolls on the upper four branches (NB), and genetic analysis revealed that all traits had medium or high heritability ranging from 0.53 to 0.86. Based on a multi-locus method (FASTmrMLM), a total of 63 stable and significant quantitative trait nucleotides (QTNs) were detected, which represented 50 genomic regions (GWAS risk loci) associated with cotton senescence. We observed three reliable loci located on chromosomes A02 (A02_105891088_107196428), D03 (D03_37952328_38393621) and D13 (D13_59408561_60730103) because of their high repeatability. One candidate gene (Ghir_D03G011060) was found in the locus D03_37952328_38393621, and its Arabidopsis thaliana homologous gene (AT5G23040) encodes a cell growth defect factor-like protein (CDF1), which might be involved in chlorophyll synthesis and cell death. Moreover, qRT-PCR showed that the transcript level of Ghir_D03G011060 was down-regulated in old cotton leaves, and virus-induced gene silencing (VIGS) indicated that silencing of Ghir_D03G011060 resulted in leaf chlorosis and promoted leaf senescence. In addition, two candidate genes (Ghir_A02G017660 and Ghir_D13G021720) were identified in loci A02_105891088_107196428 and D13_59408561_60730103, respectively. These results provide new insights into the genetic basis of cotton senescence and will serve as an important reference for the development and implementation of strategies to prevent premature senescence in cotton breeding programs.

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.).

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses.

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F 2 population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton.

KEY MESSAGE: Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton. A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

High-density genetic linkage map construction by F2 populations and QTL analysis of early-maturity traits in upland cotton (Gossypium hirsutum L.).

Due to China's rapidly increasing population, the total arable land area has dramatically decreased; as a consequence, the competition for farming land allocated for grain and cotton production has become fierce. Therefore, to overcome the existing contradiction between cotton grain and fiber production and the limited farming land, development of early-maturing cultivars is necessary. In this research, a high-density linkage map of upland cotton was constructed using genotyping by sequencing (GBS) to discover single nucleotide polymorphism (SNP) markers associated with early maturity in 170 F2 individuals derived from a cross between LU28 and ZHONG213. The high-density genetic map, which was composed of 3978 SNP markers across the 26 cotton chromosomes, spanned 2480 cM with an average genetic distance of 0.62 cM. Collinearity analysis showed that the genetic map was of high quality and accurate and agreed well with the Gossypium hirsutum reference genome. Based on this high-density linkage map, QTL analysis was performed on cotton early-maturity traits, including FT, FBP, WGP, NFFB, HNFFB and PH. A total 47 QTLs for the six traits were detected; each of these QTLs explained between 2.61% and 32.57% of the observed phenotypic variation. A major region controlling early-maturity traits in Gossypium hirsutum was identified for FT, FBP, WGP, NFFB and HNFFB on chromosome D03. QTL analyses revealed that phenotypic variation explained (PVE) ranged from 10.42% to 32.57%. Two potential candidate genes, Gh_D03G0885 and Gh_D03G0922, were predicted in a stable QTL region and had higher expression levels in the early-maturity variety ZHONG213 than in the late-maturity variety LU28. However, further evidence is required for functional validation. This study could provide useful information for the dissection of early-maturity traits and guide valuable genetic loci for molecular-assisted selection (MAS) in cotton breeding.