Chromosome-scale genome assembly of Poa trivialis and population genomics reveal widespread gene flow in a cool-season grass seed production system.

Poa trivialis (L.) is a cool-season grass species found in various environments worldwide. In addition to being a desired turfgrass species, it is a common weed of agricultural systems and natural areas. As a weed, it is an important contaminant of commercial cool-season grass seed lots, resulting in widespread gene flow facilitated by human activities and causing significant economic losses to farmers. To better understand and manage infestations, we assembled and annotated a haploid genome of P. trivialis and studied troublesome field populations from Oregon, the largest cool-season grass seed producing region in the United States. The genome assembly resulted in 1.35 Gb of DNA sequence distributed among seven chromosome-scale scaffolds, revealing a high content of transposable elements, conserved synteny with Poa annua, and a close relationship with other C3 grasses. A reduced-representation sequencing analysis of field populations revealed limited genetic diversity and suggested potential gene flow and human-assisted dispersal in the region. The genetic resources and insights into P. trivialis provided by this study will improve weed management strategies and enable the development of molecular detection tests for contaminated seed lots to limit seed-mediated gene flow. These resources should also be beneficial for turfgrass breeders seeking to improve desirable traits of commercial P. trivialis varieties and help to guide breeding efforts in other crops to enhance the resiliency of agricultural ecosystems under climate change. Significance Statement: The chromosome-scale assembly of Poa trivialis and population genomic analyses provide crucial insights into the gene flow of weedy populations in agricultural systems and contribute a valuable genomic resource for the plant science community.

Chinook salmon depth distributions on the continental shelf are shaped by interactions between location, season, and individual condition.

BACKGROUND: Ecological and physical conditions vary with depth in aquatic ecosystems, resulting in gradients of habitat suitability. Although variation in vertical distributions among individuals provides evidence of habitat selection, it has been challenging to disentangle how processes at multiple spatio-temporal scales shape behaviour. METHODS: We collected thousands of observations of depth from > 300 acoustically tagged adult Chinook salmon Oncorhynchus tshawytscha, spanning multiple seasons and years. We used these data to parameterize a machine-learning model to disentangle the influence of spatial, temporal, and dynamic oceanographic variables while accounting for differences in individual condition and maturation stage. RESULTS: The top performing machine learning model used bathymetric depth ratio (i.e., individual depth relative to seafloor depth) as a response. We found that bathymetry, season, maturation stage, and spatial location most strongly influenced Chinook salmon depth. Chinook salmon bathymetric depth ratios were deepest in shallow water, during winter, and for immature individuals. We also identified non-linear interactions among covariates, resulting in spatially-varying effects of zooplankton concentration, lunar cycle, temperature and oxygen concentration. CONCLUSIONS: Our results suggest Chinook salmon vertical habitat use is a function of ecological interactions, not physiological constraints. Temporal and spatial variation in depth distributions could be used to guide management decisions intended to reduce fishery impacts on Chinook salmon. More generally, our findings demonstrate how complex interactions among bathymetry, seasonality, location, and life history stage regulate vertical habitat selection.

Marine heatwaves disrupt ecosystem structure and function via altered food webs and energy flux.

The prevalence and intensity of marine heatwaves is increasing globally, disrupting local environmental conditions. The individual and population-level impacts of prolonged heatwaves on marine species have recently been demonstrated, yet whole-ecosystem consequences remain unexplored. We leveraged time series abundance data of 361 taxa, grouped into 86 functional groups, from six long-term surveys, diet information from a new diet database, and previous modeling efforts, to build two food web networks using an extension of the popular Ecopath ecosystem modeling framework, Ecotran. We compare ecosystem models parameterized before and after the onset of recent marine heatwaves to evaluate the cascading effects on ecosystem structure and function in the Northeast Pacific Ocean. While the ecosystem-level contribution (prey) and demand (predators) of most functional groups changed following the heatwaves, gelatinous taxa experienced the largest transformations, underscored by the arrival of northward-expanding pyrosomes. We show altered trophic relationships and energy flux have potentially profound consequences for ecosystem structure and function, and raise concerns for populations of threatened and harvested species.

An updated end-to-end ecosystem model of the Northern California Current reflecting ecosystem changes due to recent marine heatwaves.

The Northern California Current is a highly productive marine upwelling ecosystem that is economically and ecologically important. It is home to both commercially harvested species and those that are federally listed under the U.S. Endangered Species Act. Recently, there has been a global shift from single-species fisheries management to ecosystem-based fisheries management, which acknowledges that more complex dynamics can reverberate through a food web. Here, we have integrated new research into an end-to-end ecosystem model (i.e., physics to fisheries) using data from long-term ocean surveys, phytoplankton satellite imagery paired with a vertically generalized production model, a recently assembled diet database, fishery catch information, species distribution models, and existing literature. This spatially-explicit model includes 90 living and detrital functional groups ranging from phytoplankton, krill, and forage fish to salmon, seabirds, and marine mammals, and nine fisheries that occur off the coast of Washington, Oregon, and Northern California. This model was updated from previous regional models to account for more recent changes in the Northern California Current (e.g., increases in market squid and some gelatinous zooplankton such as pyrosomes and salps), to expand the previous domain to increase the spatial resolution, to include data from previously unincorporated surveys, and to add improved characterization of endangered species, such as Chinook salmon (Oncorhynchus tshawytscha) and southern resident killer whales (Orcinus orca). Our model is mass-balanced, ecologically plausible, without extinctions, and stable over 150-year simulations. Ammonium and nitrate availability, total primary production rates, and model-derived phytoplankton time series are within realistic ranges. As we move towards holistic ecosystem-based fisheries management, we must continue to openly and collaboratively integrate our disparate datasets and collective knowledge to solve the intricate problems we face. As a tool for future research, we provide the data and code to use our ecosystem model.

The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss.

To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.

Homoeologous evolution of the allotetraploid genome of Poa annua L.

BACKGROUND: Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty. RESULTS: We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis. CONCLUSIONS: The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Chromosome-Scale Genome Assembly and Annotation of Allotetraploid Annual Bluegrass (Poa annua L.).

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Accumulation of six PFAS compounds by woody and herbaceous plants: potential for phytoextraction.

Per and polyfluoroalkyl substances (PFAS) consist of a large group of compounds used to make products more resistant to stains, grease, and water and for fire suppression. They have been widely detected in the environment and exposure has been linked to adverse human health effects. Phytoremediation could be used to remediate PFAS-impacted sites, but there is little information on herbaceous and woody plant species uptake of PFAS compounds from soil. A greenhouse study evaluated the potential for eight herbaceous and seven woody plant species to absorb PFAS compounds. Six PFAS compounds: PFPeA, PFHxA, PFOA, PFBS, PFHxS, and PFOS were added weekly to irrigation water, and the plants grown for up to 14 weeks after an initial establishment period. Significant accumulation of all PFAS compounds occurred in at least one plant species. Mass recovery in above-ground tissue by the best performing plant ranged from a low of 3.8% for PFOS by Festuca rubra to a high of 42% for PFPeA by Schedonorus arundinaceus. Hyperaccumulation, defined as tissue/soil concentrations >10/1, was observed for all six PFAS compounds in at least one plant species. These results demonstrate the potential use of phytoremediation as a tool for remediating PFAS-contaminated sites.

Climate vulnerability assessment for Pacific salmon and steelhead in the California Current Large Marine Ecosystem.

Major ecological realignments are already occurring in response to climate change. To be successful, conservation strategies now need to account for geographical patterns in traits sensitive to climate change, as well as climate threats to species-level diversity. As part of an effort to provide such information, we conducted a climate vulnerability assessment that included all anadromous Pacific salmon and steelhead (Oncorhynchus spp.) population units listed under the U.S. Endangered Species Act. Using an expert-based scoring system, we ranked 20 attributes for the 28 listed units and 5 additional units. Attributes captured biological sensitivity, or the strength of linkages between each listing unit and the present climate; climate exposure, or the magnitude of projected change in local environmental conditions; and adaptive capacity, or the ability to modify phenotypes to cope with new climatic conditions. Each listing unit was then assigned one of four vulnerability categories. Units ranked most vulnerable overall were Chinook (O. tshawytscha) in the California Central Valley, coho (O. kisutch) in California and southern Oregon, sockeye (O. nerka) in the Snake River Basin, and spring-run Chinook in the interior Columbia and Willamette River Basins. We identified units with similar vulnerability profiles using a hierarchical cluster analysis. Life history characteristics, especially freshwater and estuary residence times, interplayed with gradations in exposure from south to north and from coastal to interior regions to generate landscape-level patterns within each species. Nearly all listing units faced high exposures to projected increases in stream temperature, sea surface temperature, and ocean acidification, but other aspects of exposure peaked in particular regions. Anthropogenic factors, especially migration barriers, habitat degradation, and hatchery influence, have reduced the adaptive capacity of most steelhead and salmon populations. Enhancing adaptive capacity is essential to mitigate for the increasing threat of climate change. Collectively, these results provide a framework to support recovery planning that considers climate impacts on the majority of West Coast anadromous salmonids.

Draft Genome Sequence of Salmacisia buchloëana (Basidiomycota), Which Induces Hermaphroditism in Dioecious Buffalograss.

Here, we report the draft genome of Salmacisia buchloëana type strain OK1, a monotypic species of fungi that induces ovary development in genetic male plants and increases reproductive biomass allocation in its host buffalograss. This research will enhance our understanding of fungal manipulation of host development at the genomic level.

Past and estimated future impact of invasive alien mammals on insular threatened vertebrate populations.

Invasive mammals on islands pose severe, ongoing threats to global biodiversity. However, the severity of threats from different mammals, and the role of interacting biotic and abiotic factors in driving extinctions, remain poorly understood at a global scale. Here we model global extirpation patterns for island populations of threatened and extinct vertebrates. Extirpations are driven by interacting factors including invasive rats, cats, pigs, mustelids and mongooses, native species taxonomic class and volancy, island size, precipitation and human presence. We show that controlling or eradicating the relevant invasive mammals could prevent 41-75% of predicted future extirpations. The magnitude of benefits varies across species and environments; for example, managing invasive mammals on small, dry islands could halve the extirpation risk for highly threatened birds and mammals, while doing so on large, wet islands may have little benefit. Our results provide quantitative estimates of conservation benefits and, when combined with costs in a return-on-investment framework, can guide efficient conservation strategies.

A spatially distinct history of the development of california groundfish fisheries.

During the past century, commercial fisheries have expanded from small vessels fishing in shallow, coastal habitats to a broad suite of vessels and gears that fish virtually every marine habitat on the globe. Understanding how fisheries have developed in space and time is critical for interpreting and managing the response of ecosystems to the effects of fishing, however time series of spatially explicit data are typically rare. Recently, the 1933-1968 portion of the commercial catch dataset from the California Department of Fish and Wildlife was recovered and digitized, completing the full historical series for both commercial and recreational datasets from 1933-2010. These unique datasets include landing estimates at a coarse 10 by 10 minute "grid-block" spatial resolution and extends the entire length of coastal California up to 180 kilometers from shore. In this study, we focus on the catch history of groundfish which were mapped for each grid-block using the year at 50% cumulative catch and total historical catch per habitat area. We then constructed generalized linear models to quantify the relationship between spatiotemporal trends in groundfish catches, distance from ports, depth, percentage of days with wind speed over 15 knots, SST and ocean productivity. Our results indicate that over the history of these fisheries, catches have taken place in increasingly deeper habitat, at a greater distance from ports, and in increasingly inclement weather conditions. Understanding spatial development of groundfish fisheries and catches in California are critical for improving population models and for evaluating whether implicit stock assessment model assumptions of relative homogeneity of fisheries removals over time and space are reasonable. This newly reconstructed catch dataset and analysis provides a comprehensive appreciation for the development of groundfish fisheries with respect to commonly assumed trends of global fisheries patterns that are typically constrained by a lack of long-term spatial datasets.

Pistil Smut Infection Increases Ovary Production, Seed Yield Components, and Pseudosexual Reproductive Allocation in Buffalograss.

Sex expression of dioecious buffalograss [Bouteloua dactyloides Columbus (syn. Buchloë dactyloides (Nutt.) Engelm.)] is known to be environmentally stable with approximate 1:1, male to female, sex ratios. Here we show that infection by the pistil smut fungus [Salmacisia buchloëana Huff & Chandra (syn. Tilletia buchloëana Kellerman and Swingle)] shifts sex ratios of buffalograss to be nearly 100% phenotypically hermaphroditic. In addition, pistil smut infection decreased vegetative reproductive allocation, increased most seed yield components, and increased pseudosexual reproductive allocation in both sex forms compared to uninfected clones. In female sex forms, pistil smut infection resulted in a 26 fold increase in ovary production and a 35 fold increase in potential harvest index. In male sex forms, pistil smut infection resulted in 2.37 fold increase in floret number and over 95% of these florets contained a well-developed pistil. Although all ovaries of infected plants are filled with fungal teliospores and hence reproductively sterile, an average male-female pair of infected plants exhibited an 87 fold increase in potential harvest index compared to their uninfected clones. Acquiring an ability to mimic the effects of pistil smut infection would enhance our understanding of the flowering process in grasses and our efforts to increase seed yield of buffalograss and perhaps other grasses.

Green sturgeon distribution in the Pacific Ocean estimated from modeled oceanographic features and migration behavior.

The green sturgeon (Acipenser medirostris), which is found in the eastern Pacific Ocean from Baja California to the Bering Sea, tends to be highly migratory, moving long distances among estuaries, spawning rivers, and distant coastal regions. Factors that determine the oceanic distribution of green sturgeon are unclear, but broad-scale physical conditions interacting with migration behavior may play an important role. We estimated the distribution of green sturgeon by modeling species-environment relationships using oceanographic and migration behavior covariates with maximum entropy modeling (MaxEnt) of species geographic distributions. The primary concentration of green sturgeon was estimated from approximately 41-51.5° N latitude in the coastal waters of Washington, Oregon, and Vancouver Island and in the vicinity of San Francisco and Monterey Bays from 36-37° N latitude. Unsuitably cold water temperatures in the far north and energetic efficiencies associated with prevailing water currents may provide the best explanation for the range-wide marine distribution of green sturgeon. Independent trawl records, fisheries observer records, and tagging studies corroborated our findings. However, our model also delineated patchily distributed habitat south of Monterey Bay, though there are few records of green sturgeon from this region. Green sturgeon are likely influenced by countervailing pressures governing their dispersal. They are behaviorally directed to revisit natal freshwater spawning rivers and persistent overwintering grounds in coastal marine habitats, yet they are likely physiologically bounded by abiotic and biotic environmental features. Impacts of human activities on green sturgeon or their habitat in coastal waters, such as bottom-disturbing trawl fisheries, may be minimized through marine spatial planning that makes use of high-quality species distribution information.

Habitat-related benthic macrofaunal assemblages of bays and estuaries of the western United States.

Data from 7 coastwide and regional benthic surveys were combined and used to assess the number and distribution of estuarine benthic macrofaunal assemblages of the western United States. Q-mode cluster analysis was applied to 714 samples and site groupings were tested for differences in 4 habitat factors (latitude, salinity, sediment grain size, and depth). Eight macrofaunal assemblages, structured primarily by latitude, salinity, and sediment grain size, were identified: (A) Puget Sound fine sediment, (B) Puget Sound coarse sediment, (C) southern California marine bays, (D) polyhaline central San Francisco Bay, (E) shallow estuaries and wetlands, (F) saline very coarse sediment, (G) mesohaline San Francisco Bay, and (H) limnetic and oligohaline. The Puget Sound, southern California, and San Francisco Bay assemblages were geographically distinct, while Assemblages E, F and H were distributed widely along the entire coast. A second Q-mode cluster analysis was conducted after adding replicate samples that were available from some of the sites and temporal replicates that were available for sites that were sampled in successive years. Variabilities due to small spatial scale habitat heterogeneity and temporal change were both low in Puget Sound, but temporal variability was high in the San Francisco estuary where large fluctuations in freshwater inputs and salinity among years leads to spatial relocation of the assemblages.

Green sturgeon physical habitat use in the coastal Pacific Ocean.

The green sturgeon (Acipenser medirostris) is a highly migratory, oceanic, anadromous species with a complex life history that makes it vulnerable to species-wide threats in both freshwater and at sea. Green sturgeon population declines have preceded legal protection and curtailment of activities in marine environments deemed to increase its extinction risk. Yet, its marine habitat is poorly understood. We built a statistical model to characterize green sturgeon marine habitat using data from a coastal tracking array located along the Siletz Reef near Newport, Oregon, USA that recorded the passage of 37 acoustically tagged green sturgeon. We classified seafloor physical habitat features with high-resolution bathymetric and backscatter data. We then described the distribution of habitat components and their relationship to green sturgeon presence using ordination and subsequently used generalized linear model selection to identify important habitat components. Finally, we summarized depth and temperature recordings from seven green sturgeon present off the Oregon coast that were fitted with pop-off archival geolocation tags. Our analyses indicated that green sturgeon, on average, spent a longer duration in areas with high seafloor complexity, especially where a greater proportion of the substrate consists of boulders. Green sturgeon in marine habitats are primarily found at depths of 20-60 meters and from 9.5-16.0°C. Many sturgeon in this study were likely migrating in a northward direction, moving deeper, and may have been using complex seafloor habitat because it coincides with the distribution of benthic prey taxa or provides refuge from predators. Identifying important green sturgeon marine habitat is an essential step towards accurately defining the conditions that are necessary for its survival and will eventually yield range-wide, spatially explicit predictions of green sturgeon distribution.

Mixed-source reintroductions lead to outbreeding depression in second-generation descendents of a native North American fish.

Reintroductions are commonly employed to preserve intraspecific biodiversity in fragmented landscapes. However, reintroduced populations are frequently smaller and more geographically isolated than native populations. Mixing genetically, divergent sources are often proposed to attenuate potentially low genetic diversity in reintroduced populations that may result from small effective population sizes. However, a possible negative tradeoff for mixing sources is outbreeding depression in hybrid offspring. We examined the consequences of mixed-source reintroductions on several fitness surrogates at nine slimy sculpin (Cottus cognatus) reintroduction sites in south-east Minnesota. We inferred the relative fitness of each crosstype in the reintroduced populations by comparing their growth rate, length, weight, body condition and persistence in reintroduced populations. Pure strain descendents from a single source population persisted in a greater proportion than expected in the reintroduced populations, whereas all other crosstypes occurred in a lesser proportion. Length, weight and growth rate were lower for second-generation intra-population hybrid descendents than for pure strain and first-generation hybrids. In the predominant pure strain, young-of the-year size was significantly greater than any other crosstype. Our results suggested that differences in fitness surrogates among crosstypes were consistent with disrupted co-adapted gene complexes associated with beneficial adaptations in these reintroduced populations. Future reintroductions may be improved by evaluating the potential for local adaptation in source populations or by avoiding the use of mixed sources by default when information on local adaptations or other genetic characteristics is lacking.

Genetic variability in populations of the southern chinch bug, Blissus insularis, assessed using AFLP analysis.

Southern chinch bug, Blissus insularis Barber (Heteroptera: Blissidae), is the most destructive insect pest of St. Augustine grass, Stenotaphrum secundatum Waltz (Kuntze), in the southern United States. The present study is focused on assessing genetic variability in five populations of B. insularis collected from Texas and Florida where St. Augustine grass is widely grown. The amplified fragment length polymorphism technique was used to DNA fingerprint individuals from each population (a total of 46 individuals) using five primer combinations (EcoRI/MSeI). Analysis of molecular variance results show no evidence to support significant genetic variability among Texas and Florida populations of B. insularis. Nearly all genetic variation was found to reside within populations (95%), with only approximately 3% residing among populations between regions. Low G(ST) values obtained from POPGENE and low F(ST) values obtained from the analysis of molecular variance both support the conclusion for high levels of gene flow resulting from interbreeding and/or migratory events among the populations. A Mantel test of Euclidean squared distances showed no correlation between the genetic distance and geographic distance matrices of tested populations of B. insularis. The results of the present study suggests that gene flow is occurring among populations of B. insularis and, therefore, breeders need to be aware that this resistance will most likely not remain localized, and it has the potential to spread as a result of migratory events.

A fungal parasite regulates a putative female-suppressor gene homologous to maize tasselseed2 and causes induced hermaphroditism in male buffalograss.

Parasitically induced hermaphroditism is a fascinating illustration of floral sex organ modification; however, knowledge of how parasites induce hermaphroditism in plants is limited. Here, we show the fungal parasite pistil smut induces development of female sex organs (pistils) in flowers of male buffalograss, potentially by downregulating a putative female-suppressor gene, BdTs2, homologous to maize Tasselseed2 (ZmTs2). Full-length BdTs2, isolated using rapid amplification of cDNA ends, exhibits 89% nucleotide sequence similarity with ZmTs2 and 85% amino acid sequence homology with ZmTs2 protein. Scanning electron micrographs demonstrate that unisexual buffalograss flowers develop through a process of selective abortion of opposite sex organs within hermaphroditic floral primordia. Quantitative real-time polymerase chain reaction showed that high expression levels of BdTs2 within male inflorescences correlate with the selective abortion of gynoecium, leading to the development of unisexual male flowers. RNA in situ hybridization confirmed the expression of BdTs2 precisely within vestigial gynoeciums of male flowers and not in other floral organs of the inflorescence. Furthermore, we show that BdTs2 expression is downregulated by pistil smut infection, which corresponds to the presence of pistils in flowers otherwise destined to become unisexual male. This study provides a potential molecular basis for pistil smut-induced hermaphroditism in male buffalograss.