Meta-Transcriptomic Analysis Uncovers the Presence of Four Novel Viruses and Multiple Known Virus Genera in a Single Hibiscus rosa-sinensis Plant in Colombia.

Hibiscus is not native to Colombia but well suited to its arid soil and dry climates. A single hibiscus plant from Risaralda, showing black spots on upper and lower sides of its leaves, was collected for virome analysis using meta-transcriptomic high-throughput sequencing technology. Bioinformatic analysis identified 12.5% of the total reads in the Ribo-Zero cDNA library which mapped to viral genomes. BLAST searches revealed the presence of carlavirus, potexvirus, and of known members of the genera Betacarmovirus, Cilevirus, Nepovirus, and Tobamovirus in the sample; confirmed by RT-PCR with virus-specific primers followed by amplicon sequencing. Furthermore, in silico analysis suggested the possibility of a novel soymovirus, and a new hibiscus strain of citrus leprosis virus C2 in the mixed infection. Both RNA dependent RNA polymerase and coat protein gene sequences of the potex and carla viruses shared less than 72% nucleotide and 80% amino acid identities with any alphaflexi- and betaflexi-virus sequences available in GenBank, identifying three novel carlavirus and one potexvirus species in the Hibiscus rosa-sinensis plant. The detection of physalis vein necrosis nepovirus and passion fruit green spot cilevirus in hibiscus are also new reports from Colombia. Overall, the meta-transcriptome analysis identified the complex virome associated with the black spot symptoms on hibiscus leaves and demonstrated the diversity of virus genera tolerated in the mixed infection of a single H. rosa-sinensis plant.

Go west: Population genomics reveals unexpected population fluctuations and little gene flow in Western hemisphere populations of the predatory lady beetle, Hippodamia convergens.

Hippodamia convergens-the convergent lady beetle, has been used extensively in augmentative biological control of aphids, thrips, and whiteflies across its native range in North America, and was introduced into South America in the 1950s. Overwintering H. convergens populations from its native western range in the United States are commercially collected and released across its current range in the eastern USA, with little knowledge of the effectiveness of its augmentative biological control. Here we use a novel ddRADseq-based SNP/haplotype discovery approach to estimate its range-wide population diversity, differentiation, and recent evolutionary history. Our results indicate (1) significant population differentiation among eastern USA, western USA, and South American populations of H. convergens, with (2) little to no detectable recent admixture between them, despite repeated population augmentation, and (3) continued recent population size expansion across its range. These results contradict previous findings using microsatellite markers. In light of these new findings, the implications for the effectiveness of augmentative biological control using H. convergens are discussed. Additionally, because quantifying the non-target effects of augmentative biological control is a difficult problem in migratory beetles, our results could serve as a cornerstone in improving and predicting the efficacy of future releases of H. convergens across its range.

Genetic Analysis of the Emerging Citrus Yellow Vein Clearing Virus Reveals a Divergent Virus Population in American Isolates.

Citrus yellow vein clearing virus is a previously reported citrus virus from Asia with widespread distribution in China. In 2022, the California Department of Food and Agriculture conducted a multipest citrus survey targeting multiple citrus pathogens including citrus yellow vein clearing virus (CYVCV). In March 2022, a lemon tree with symptoms of vein clearing, chlorosis, and mottling in a private garden in the city of Tulare, California, tested positive for CYVCV, which triggered an intensive survey in the surrounding areas. A total of 3,019 plant samples, including citrus and noncitrus species, were collected and tested for CYVCV using conventional reverse transcription polymerase chain reaction, reverse transcription quantitative polymerase chain reaction, and Sanger sequencing. Five hundred eighty-six citrus trees tested positive for CYVCV, including eight citrus species not previously recorded infected under field conditions. Comparative genomic studies were conducted using 17 complete viral genomes. Sequence analysis revealed two major phylogenetic groups. Known Asian isolates and five California isolates from this study made up the first group, whereas all other CYVCV isolates from California formed a second group, distinct from all worldwide isolates. Overall, the CYVCV population shows rapid expansion and high differentiation indicating a population bottleneck typical of a recent introduction into a new geographic area.

Comparative Analysis of Tomato Brown Rugose Fruit Virus Isolates Shows Limited Genetic Diversity.

Tomato is an important vegetable in the United States and around the world. Recently, tomato brown rugose fruit virus (ToBRFV), an emerging tobamovirus, has impacted tomato crops worldwide and can result in fruit loss. ToBRFV causes severe symptoms, such as mosaic, puckering, and necrotic lesions on leaves; other symptoms include brown rugose and marbling on fruits. More importantly, ToBRFV can overcome resistance in tomato cultivars carrying the Tm-22 locus. In this study, we recovered ToBRFV sequences from tomato seeds, leaves, and fruits from the U.S., Mexico, and Peru. Samples were pre-screened using a real-time RT-PCR assay prior to high-throughput sequencing. Virus draft genomes from 22 samples were assembled and analyzed against more than 120 publicly available genomes. Overall, most sequenced isolates were similar to each other and did not form a distinct population. Phylogenetic analysis revealed three clades within the ToBRFV population. Most of the isolates (95%) clustered in clade 3. Genetic analysis revealed differentiation between the three clades indicating minor divergence occurring. Overall, pairwise identity showed limited genetic diversity among the isolates in this study with worldwide isolates, with a pairwise identity ranging from 99.36% and 99.97%. The overall population is undergoing high gene flow and population expansion with strong negative selection pressure at all ToBRFV genes. Based on the results of this study, it is likely that the limited ToBRFV diversity is associated with the rapid movement and eradication of ToBRFV-infected material between countries.

Complete Genome Sequence of a Tomato Brown Rugose Fruit Virus Isolated in the United States.

The complete genome sequence of a U.S. isolate of a Tomato brown rugose fruit virus (ToBRFV) (CA18-01) was obtained through Illumina and MinION sequencing. The U.S. ToBRFV isolate shared a high nucleic acid sequence identity (>99%) with known ToBRFV isolates. Phylogenetic analysis revealed a tight cluster for ToBRFV isolates throughout the world, suggesting a short evolutionary history.

Mutations in a Novel Cadherin Gene Associated with Bt Resistance in Helicoverpa zea.

Transgenic corn and cotton produce crystalline (Cry) proteins derived from the soil bacterium Bacillus thuringiensis (Bt) that are toxic to lepidopteran larvae. Helicoverpa zea, a key pest of corn and cotton in the U.S., has evolved widespread resistance to these proteins produced in Bt corn and cotton. While the genomic targets of Cry selection and the mutations that produce resistant phenotypes are known in other lepidopteran species, little is known about how selection by Cry proteins shape the genome of H. zea We scanned the genomes of Cry1Ac-selected and unselected H. zea lines, and identified twelve genes on five scaffolds that differed between lines, including cadherin-86C (cad-86C), a gene from a family that is involved in Cry1A resistance in other lepidopterans. Although this gene was expressed in the H. zea larval midgut, the protein it encodes has only 17 to 22% identity with cadherin proteins from other species previously reported to be involved in Bt resistance. An analysis of midgut-expressed cDNAs showed significant between-line differences in the frequencies of putative nonsynonymous substitutions (both SNPs and indels). Our results indicate that cad-86C is a likely target of Cry1Ac selection in H. zea It remains unclear, however, whether genomic changes at this locus directly disrupt midgut binding of Cry1Ac and cause Bt resistance, or indirectly enhance fitness of H. zea in the presence of Cry1Ac by some other mechanism. Future work should investigate phenotypic effects of these nonsynonymous substitutions and their impact on fitness of H. zea larvae that ingest Cry1Ac.

Estimation of contemporary effective population size and population declines using RAD sequence data.

Large genomic data sets generated with restriction site-associated DNA sequencing (RADseq), in combination with demographic inference methods, are improving our ability to gain insights into the population history of species. We used a simulation approach to examine the potential for RADseq data sets to accurately estimate effective population size (N e) over the course of stable and declining population trends, and we compare the ability of two methods of analysis to accurately distinguish stable from steadily declining populations over a contemporary time scale (20 generations). Using a linkage disequilibrium-based analysis, individual sampling (i.e., n ≥ 30) had the greatest effect on N e estimation and the detection of population size declines, with declines reliably detected across scenarios ~10 generations after they began. Coalescent-based inference required fewer sampled individuals (i.e., n = 15), and instead was most influenced by the size of the SNP data set, with 25,000-50,000 SNPs required for accurate detection of population trends and at least 20 generations after decline began. The number of samples available and targeted number of RADseq loci are important criteria when choosing between these methods. Neither method suffered any apparent bias due to the effects of allele dropout typical of RAD data. With an understanding of the limitations and biases of these approaches, researchers can make more informed decisions when designing their sampling and analyses. Overall, our results reveal that demographic inference using RADseq data can be successfully applied to infer recent population size change and may be an important tool for population monitoring and conservation biology.

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species.

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census-based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark-recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent-based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark-recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome-based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.

Patterns of amphibian infection prevalence across wetlands on the Savannah River Site, South Carolina, USA.

Amphibian diseases, such as chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd) and ranaviral disease caused by ranaviruses, are often linked to global amphibian population declines, yet the ecological dynamics of both pathogens are poorly understood. The goal of our study was to determine the baseline prevalence, pathogen loads, and co-infection rate of Bd and ranavirus across the Savannah River Site (SRS) in South Carolina, USA, a region with rich amphibian diversity and a history of amphibian-based research. We tested over 1000 individuals, encompassing 21 amphibian species from 11 wetlands for both Bd and ranavirus. The prevalence of Bd across individuals was 9.7%. Using wetland means, the mean (±SE) Bd prevalence was 7.9 ± 2.9%. Among toad species, Anaxyrus terrestris had 95 and 380% greater odds of being infected with Bd than Scaphiopus holbrookii and Gastrophryne carolinensis, respectively. Odds of Bd infection in adult A. terrestris and Lithobates sphenocephalus were 75 to 77% greater in metal-contaminated sites. The prevalence of ranavirus infections across all individuals was 37.4%. Mean wetland ranavirus prevalence was 29.8 ± 8.8% and was higher in post-metamorphic individuals than in aquatic larvae. Ambystoma tigrinum had 83 to 85% higher odds of ranavirus infection than A. opacum and A. talpoideum. We detected a 4.8% co-infection rate, with individuals positive for ranavirus having a 5% higher occurrence of Bd. In adult Anaxyrus terrestris, odds of Bd infection were 13% higher in ranavirus-positive animals and odds of co-infection were 23% higher in contaminated wetlands. Overall, we found the pathogen prevalence varied by wetland, species, and life stage.

Development and characterization of microsatellite loci for the endangered scrub lupine, Lupinus aridorum (Fabaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed in scrub lupine (Lupinus aridorum, Fabaceae), an endemic species to Florida that is listed as endangered in the United States, to assess connectivity among populations, identify hybrids, and examine genetic diversity. METHODS AND RESULTS: We isolated and characterized 12 microsatellite loci polymorphic in scrub lupine or in closely related species (i.e., sky-blue lupine [L. diffusus] and Gulf Coast lupine [L. westianus]). Loci showed low to moderate polymorphism, ranging from two to 14 alleles per locus and 0.01 to 0.86 observed heterozygosity. CONCLUSIONS: These loci are the first developed for Florida species of lupine and will be used to determine differentiation among species and to aid in conservation of the endangered scrub lupine.

32 species validation of a new Illumina paired-end approach for the development of microsatellites.

Development and optimization of novel species-specific microsatellites, or simple sequence repeats (SSRs) remains an important step for studies in ecology, evolution, and behavior. Numerous approaches exist for identifying new SSRs that vary widely in terms of both time and cost investments. A recent approach of using paired-end Illumina sequence data in conjunction with the bioinformatics pipeline, PAL_FINDER, has the potential to substantially reduce the cost and labor investment while also improving efficiency. However, it does not appear that the approach has been widely adopted, perhaps due to concerns over its broad applicability across taxa. Therefore, to validate the utility of the approach we developed SSRs for 32 species representing 30 families, 25 orders, 11 classes, and six phyla and optimized SSRs for 13 of the species. Overall the IPE method worked extremely well and we identified 1000s of SSRs for all species (mean = 128,485), with 17% of loci being potentially amplifiable loci, and 25% of these met our most stringent criteria designed to that avoid SSRs associated with repetitive elements. Approximately 61% of screened primers yielded strong amplification of a single locus.

Characterization of 42 polymorphic microsatellite loci in Mimulus ringens (Phrymaceae) using Illumina sequencing.

PREMISE OF THE STUDY: Microsatellite markers were isolated and characterized in Mimulus ringens (Phrymaceae), a herbaceous wetland perennial, to facilitate studies of mating patterns and population genetic structure. • METHODS AND RESULTS: A total of 42 polymorphic loci were identified from a sample of 24 individuals from a single population in Ohio, USA. The number of alleles per locus ranged from two to nine, and median observed heterozygosity was 0.435. • CONCLUSIONS: This large number of polymorphic loci will enable researchers to quantify male fitness, patterns of multiple paternity, selfing, and biparental inbreeding in large natural populations of this species. These markers will also permit detailed study of fine-scale patterns of genetic structure.

Development and characterization of microsatellite markers for Berberis thunbergii (Berberidaceae).

PREMISE OF THE STUDY: Microsatellite markers were isolated and characterized in Berberis thunbergii, an invasive and ornamental shrub in the eastern United States, to assess genetic diversity among populations and potentially identify horticultural cultivars. METHODS AND RESULTS: A total of 12 loci were identified for the species. Eight of the loci were polymorphic and were screened in 24 individuals from two native (Tochigi and Ibaraki prefectures, Japan) and one invasive (Connecticut, USA) population and 21 horticultural cultivars. The number of alleles per locus ranged from three to seven, and observed heterozygosity ranged from 0.048 to 0.636. CONCLUSIONS: These new markers will provide tools for examining genetic relatedness of B. thunbergii plants in the native and invasive range, including phylogeographic studies and assessment of rapid evolution in the invasive range. These markers may also provide tools for examining hybridization with other related species in the invasive range.