Two haplotype-resolved genomes reveal important flower traits in bigleaf hydrangea (Hydrangea macrophylla) and insights into Asterid evolution.

The Hydrangea genus belongs to the Hydrangeaceae family, in the Cornales order of flowering plants, which early diverged among the Asterids, and includes several species that are commonly used ornamental plants. Of them, Hydrangea macrophylla is one of the most valuable species in the nursery trade, yet few genomic resources are available for this crop or closely related Asterid species. Two high-quality haplotype-resolved reference genomes of hydrangea cultivars 'Veitchii' and 'Endless Summer' [highest quality at 2.22 gigabase pairs (Gb), 396 contigs, N50 22.8 megabase pairs (Mb)] were assembled and scaffolded into the expected 18 pseudochromosomes. Utilizing the newly developed high-quality reference genomes along with high-quality genomes of other related flowering plants, nuclear data were found to support a single divergence point in the Asterids clade where both the Cornales and Ericales diverged from the euasterids. Genetic mapping with an F1 hybrid population demonstrated the power of linkage mapping combined with the new genomic resources to identify the gene for inflorescence shape, CYP78A5 located on chromosome 4, and a novel gene, BAM3 located on chromosome 17, for causing double flower. Resources developed in this study will not only help to accelerate hydrangea genetic improvement but also contribute to understanding the largest group of flowering plants, the Asterids.

Complete chloroplast genome comparisons for Pityopsis (Asteraceae).

Pityopsis includes several regionally and one federally endangered species of herbaceous perennials. Four species are highly localized, including the federally endangered P. ruthii. The genus includes several ploidy levels and interesting ecological traits such as drought tolerance and fire-dependent flowering. Results from previous cladistic analyses of morphology and from initial DNA sequence studies did not agree with one another or with the infrageneric taxonomic classification, with the result that infrageneric relationships remain unresolved. We sequenced, assembled, and compared the chloroplast (cp) genomes of 12 species or varieties of Pityopsis to better understand generic evolution. A reference cp genome 152,569 bp in length was assembled de novo from P. falcata. Reads from other sampled species were then aligned to the P. falcata reference and individual chloroplast genomes were assembled for each, with manual gapfilling and polishing. After removing the duplicated second inverted region, a multiple sequence alignment of the cp genomes was used to construct a maximum likelihood (ML) phylogeny for the twelve cp genomes. Additionally, we constructed a ML phylogeny from the nuclear ribosomal repeat region after mapping reads to the Helianthus annuus region. The chloroplast phylogeny supported two clades. Previously proposed clades and taxonomic sections within the genus were largely unsupported by both nuclear and chloroplast phylogenies. Our results provide tools for exploring hybridity and examining the physiological and genetic basis for drought tolerance and fire-dependent flowering. This study will inform breeding and conservation practices, and general knowledge of evolutionary history, hybridization, and speciation within Pityopsis.

Low Genetic Diversity Suggests the Recent Introduction of Dogwood Powdery Mildew to North America.

Cornus florida (flowering dogwood) is a popular understory tree endemic to the eastern hardwood forests of the United States. In 1996, dogwood powdery mildew caused by Erysiphe pulchra, an obligate biotrophic fungus of large bracted dogwoods, reached epidemic levels throughout the C. florida growing region. In the late 1990s, both sexual and asexual stages of E. pulchra were regularly observed; thereafter, the sexual stage was found less frequently. We examined the genetic diversity and population structure of 167 E. pulchra samples on C. florida leaves using 15 microsatellite loci. Samples were organized into two separate collection zone data sets, separated as eight zones and two zones, for the subsequent analysis of microsatellite allele length data. Clone correction analysis reduced the sample size to 90 multilocus haplotypes. Our study indicated low genetic diversity, a lack of definitive population structure, low genetic distance among multilocus haplotypes, and significant linkage disequilibrium among zones. Evidence of a population bottleneck was also detected. The results of our study indicated a high probability that E. pulchra reproduces predominately via asexual conidia and lend support to the hypothesis that E. pulchra is an exotic pathogen to North America.[Formula: see text] Copyright © 2019 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Development of Genomic Resources for the Powdery Mildew, Erysiphe pulchra.

Powdery mildews (PMs) are important plant pathogens causing widespread damage. Here, we report the first draft genome of Erysiphe pulchra, the causative agent of PM of flowering dogwood, Cornus florida. The assembled genome was 63.5 Mbp and resulted in formation of 19,442 contigs (N50 = 11,686 bp) that contained an estimated 6,860 genes with a genome coverage of 62×. We found 102 candidate secreted effector proteins (CSEPs) in E. pulchra similar to E. necator genes that are potentially involved in disease development. This draft genome is an initial step for understanding the evolutionary history of the PMs and will also provide insight into evolutionary strategies that led to the wide host expansion and environmental adaptations so effectively employed by the PM lineages.

Comprehensive evaluation of RNA-seq analysis pipelines in diploid and polyploid species.

Background: The usual analysis of RNA sequencing (RNA-seq) reads is based on an existing reference genome and annotated gene models. However, when a reference for the sequenced species is not available, alternatives include using a reference genome from a related species or reconstructing transcript sequences with de novo assembly. In addition, researchers are faced with many options for RNA-seq data processing and limited information on how their decisions will impact the final outcome. Using both a diploid and polyploid species with a distant reference genome, we have tested the influence of different tools at various steps of a typical RNA-seq analysis workflow on the recovery of useful processed data available for downstream analysis. Findings: At the preprocessing step, we found error correction has a strong influence on de novo assembly but not on mapping results. After trimming, a greater percentage of reads could be used in downstream analysis by selecting gentle quality trimming performed with Skewer instead of strict quality trimming with Trimmomatic. This availability of reads correlated with size, quality, and completeness of de novo assemblies and with number of mapped reads. When selecting a reference genome from a related species to map reads, outcome was significantly improved when using mapping software tolerant of greater sequence divergence, such as Stampy or GSNAP. Conclusions: The selection of bioinformatic software tools for RNA-seq data analysis can maximize quality parameters on de novo assemblies and availability of reads in downstream analysis.

Population Structure and Genetic Diversity Within the Endangered Species Pityopsis ruthii (Asteraceae).

Pityopsis ruthii (Ruth's golden aster) is a federally endangered herbaceous perennial endemic to the Hiwassee and Ocoee Rivers in southeastern Tennessee, United States. Comprehensive genetic studies providing novel information to conservationists for preservation of the species are lacking. Genetic variation and gene flow were evaluated for 814 individuals from 33 discrete locations using polymorphic microsatellites: seven chloroplast and twelve nuclear. A total of 198 alleles were detected with the nuclear loci and 79 alleles with the chloroplast loci. Gene flow was estimated, with the Hiwassee River (Nm = 2.16; FST = 0.15) showing higher levels of gene flow and lower levels of population differentiation than the Ocoee River (Nm = 1.28; FST = 0.19). Population structure was examined using Bayesian cluster analyses. Nuclear and chloroplast analyses were incongruent. From the chloroplast microsatellites, three clusters were identified; all were present in sampling sites at both rivers, indicating a lack of allele fixation along rivers. Nuclear markers revealed two clusters and separated by river. When the Hiwassee River locations were analyzed, four clusters were identified for both the chloroplast and nuclear microsatellites, though the individuals clustered differently. Analysis of the Ocoee River revealed two clusters for the chloroplast microsatellites and three for the nuclear microsatellites. We recommend P. ruthii be managed as four populations for the Hiwassee River and three populations for the Ocoee River. Our results provide critical genetic information for P. ruthii that can be used for species management decisions to drive future population augmentation/reintroduction and ex situ conservation efforts.

Regulation of gene expression in roots of the pH-sensitive Vaccinium corymbosum and the pH-tolerant Vaccinium arboreum in response to near neutral pH stress using RNA-Seq.

BACKGROUND: Blueberries are one of the few horticultural crops adapted to grow in acidic soils. Neutral to basic soil pH is detrimental to all commonly cultivated blueberry species, including Vaccinium corymbosum (VC). In contrast, the wild species V. arboreum (VA) is able to tolerate a wider range of soil pH. To assess the molecular mechanisms involved in near neutral pH stress response, plants from pH-sensitive VC (tetraploid) and pH-tolerant VA (diploid) were grown at near neutral pH 6.5 and at the preferred pH of 4.5. RESULTS: Transcriptome sequencing of root RNA was performed for 4 biological replications per species x pH level interaction, for a total of 16 samples. Reads were mapped to the reference genome from diploid V. corymbosum, transforming ~55% of the reads to gene counts. A quasi-likelihood F test identified differential expression due to pH stress in 337 and 4867 genes in VA and VC, respectively. Both species shared regulation of genes involved in nutrient homeostasis and cell wall metabolism. VA and VC exhibited differential regulation of signaling pathways related to abiotic/biotic stress, cellulose and lignin biosynthesis, and nutrient uptake. CONCLUSIONS: The specific responses in VA likely facilitate tolerance to higher soil pH. In contrast, response in VC, despite affecting a greater number of genes, is not effective overcoming the stress induced by pH. Further inspection of those genes with differential expression that are specific in VA may provide insight on the mechanisms towards tolerance.

Variability of United States isolates of Macrophomina phaseolina based on simple sequence repeats and cross genus transferability to related genera within botryosphaeriaceae.

Twelve simple sequence repeat (SSRs) loci were used to evaluate genetic diversity of 109 isolates of Macrophomina phaseolina collected from different geographical regions and host species throughout the United States (US). Genetic diversity was assessed using Nei's minimum genetic distance, and the usefulness of each locus was determined by calculating the polymorphism information content (PIC). A total of 98 alleles were detected and of these 31 were unique to individual genotypes. Eight of twelve loci were highly informative with PIC values greater than 0.50. The majority of pairwise comparisons of genetic distance were greater than 0.60 indicating moderate to high genetic diversity. Dendrograms based on the genetic dissimilarities were created for the 109 isolates of which 79 were from soybean. Some clustering by host and geography was noted, but, the dendrograms generally grouped isolates independent of host or geography. Additionally, sequencing of the internal transcribed spacer region (ITS) for 10 isolates revealed that all of these isolates were 99% similar. Three SSR loci from M. phaseolina were cross amplified in other genera in the Botryosphaeriaceae. This was the first study of genotyping and assessing genetic diversity of M. phaseolina isolates collected from a widespread host and geographic range across the US with SSRs. With an additional 34 loci publically available for M. phaseolina, the results indicate that previously developed SSRs from one species can be used in future population, ecological, and genetic studies of M. phaseolina and other genera within the Botryosphaeriaceae.

Bright-Field and Fluorescence Microscopic Study of Development of Erysiphe polygoni in Susceptible and Resistant Bigleaf Hydrangea.

Temporal development of Erysiphe polygoni and responses of bigleaf hydrangeas (Hydrangea macrophylla) to the fungal attack were investigated using bright-field and fluorescence microscopy. Conidia germinated 2 h after inoculation (HAI) and formed primary appressoria at the tip of the primary germ tubes within 4 HAI. Secondary germ tubes were initiated from primary appressoria or other parts of conidia 12 HAI. Hyphae developed through elongation of secondary germ tubes, and paired lateral appressoria were formed along hyphae within 2 days after inoculation (DAI). Conidiophores and conidia were formed 5 DAI. In the susceptible cultivar Nikko Blue and the resistant cultivar Veitchii, the fungus established a parasitic relationship, which was indicated by the formation of haustoria under primary appressoria and development of secondary germ tubes at 1 DAI. A hypersensitive response (HR) and accumulation of callose were detected in both resistant and susceptible cultivars at 3 DAI. Resistance to powdery mildew in Veitchii was evident by manifestation of early accumulation of callose, relatively high percentage of necrotic infected cells, and restricted colony development compared to the susceptible cultivar Nikko Blue. Restricting hyphal growth and sporulation by early response of callose accumulation and HR are important resistance mechanisms that could be used in screening hydrangeas for resistance to powdery mildew.

Microsatellites from the charcoal rot fungus (Macrophomina phaseolina).

Microsatellite loci were identified from the charcoal rot fungus (Macrophomina phaseolina). Primer pairs for 46 loci were developed, and of these, 13 were optimized and screened using genomic DNA from 55 fungal isolates collected predominantly from two soybean fields in Mississippi. Twelve of the optimized loci were polymorphic and the number of alleles per locus ranged from 6 to 22. These microsatellites will be useful in population and pathogenicity studies to correspond with development of potential disease-resistant soybean and other susceptible crops.

Microsatellites from kousa dogwood (Cornus kousa).

Microsatellite loci were identified from Cornus kousa'National'. Primer pairs for 86 loci were developed and of these, eight were optimized and screened using genomic DNA from 22 kousa cultivars. All optimized loci were polymorphic and the number of alleles per locus ranged from three to 17. Observed heterozygosity ranged from 0 to 0.3 and expected heterozygosity ranged from 0.38 to 0.91. These microsatellites will be useful in population studies, and a breeding programme for cultivar development of Cornus species.

Rapid, high-throughput detection of azalea lace bug (Hemiptera: Tingidae) predation by Chrysoperla rufilabris (Neuroptera: Chrysopidae), using fluorescent-polymerase chain reaction primers.

Azalea lace bugs, Stephanitis pyrioides (Scott) (Hemiptera: Tingidae), are the most common pest of azaleas (Rhododendron spp.) in nursery production and the landscape. Although pesticides are commonly used to control lace bugs, natural enemies can be a significant source of lace bug mortality. Lacewings (Neuroptera: Chrysopidae) are natural enemies of lace bugs and easily consume them in laboratory studies. Field studies on lacewing biocontrol of azalea lace bugs are underway; however, monitoring lacewing predation in a nursery environment by direct observation is impractical. Here, we describe a fluorescent-polymerase chain reaction method to estimate S. pyrioides consumption based on the gut contents of lacewing predators. Lace bug DNA was detected in fed lacewings up to 32 h after ingestion. More than 80% of the ingested lace bugs were detected using our method with only one false positive result. The assay is both high-throughput and relatively inexpensive, making it a practical approach to documenting lace bug predation in the field.