Comparative transcriptomics of soybean genotypes with partial resistance towards Phytophthora sojae, Conrad, and M92-220, to moderately susceptible fast neutron mutant soybeans and Sloan.

The breeding of disease-resistant soybeans cultivars to manage Phytophthora root and stem rot caused by the pathogen Phytophthora sojae involves combining quantitative disease resistance (QDR) and Rps gene-mediated resistance. To identify and confirm potential mechanisms of QDR towards P. sojae, we conducted a time course study comparing changes in gene expression among Conrad and M92-220 with high QDR to susceptible genotypes, Sloan and 3 mutants derived from fast neutron (FN) irradiation of M92-220. Differentially expressed genes from Conrad and M92-220 indicated several shared defense-related pathways at the transcriptomic level, but also defense pathways unique to each cultivar such as stilbenoid, diarylheptanoid and gingerol biosynthesis, and monobactam biosynthesis. Gene Ontology pathway analysis showed that the susceptible FN mutants lacked enrichment of three terpenoid related-pathways and two cell wall-related pathways at either one or both timepoints, in contrast to M92-220. The susceptible mutants also lacked enrichment of potentially important KEGG pathways at either one or both timepoints, including sesquiterpenoid and triterpenoid biosynthesis, thiamine metabolism, arachidonic acid, stilbenoid, diarylheptanoid and gingerol biosynthesis, and monobactam biosynthesis. Additionally, thirty-one genes which were differentially expressed in M92-220 following P. sojae infection were not expressed in the mutants. These 31 genes have annotations related to unknown proteins, valine, leucine, and isoleucine biosynthesis and protein and lipid metabolic processes. The results of this study confirm previously proposed mechanisms of QDR, provide evidence for potential novel QDR pathways in M92-220, and furthers our understanding of the complex network associated with QDR mechanisms in soybean towards P. sojae.

In vitro, in planta, and comparative genomic analyses of Pseudomonas syringae pv. syringae strains of pepper (Capsicum annuum var. annuum).

Pseudomonas syringae pv. syringae (Pss) is an emerging phytopathogen that causes Pseudomonas leaf spot (PLS) disease in pepper plants. Pss can cause serious economic damage to pepper production, yet very little is known about the virulence factors carried by Pss that cause disease in pepper seedlings. In this study, Pss strains isolated from pepper plants showing PLS symptoms in Ohio between 2013 and 2021 (n = 16) showed varying degrees of virulence (Pss populations and disease symptoms on leaves) on 6-week-old pepper seedlings. In vitro studies assessing growth in nutrient-limited conditions, biofilm production, and motility also showed varying degrees of virulence, but in vitro and in planta variation in virulence between Pss strains did not correlate. Comparative whole-genome sequencing studies identified notable virulence genes including 30 biofilm genes, 87 motility genes, and 106 secretion system genes. Additionally, a total of 27 antimicrobial resistance genes were found. A multivariate correlation analysis and Scoary analysis based on variation in gene content (n = 812 variable genes) and single nucleotide polymorphisms within virulence genes identified no significant correlations with disease severity, likely due to our limited sample size. In summary, our study explored the virulence and antimicrobial gene content of Pss in pepper seedlings as a first step toward understanding the virulence and pathogenicity of Pss in pepper seedlings. Further studies with additional pepper Pss strains will facilitate defining genes in Pss that correlate with its virulence in pepper seedlings, which can facilitate the development of effective measures to control Pss in pepper and other related P. syringae pathovars. IMPORTANCE: Pseudomonas leaf spot (PLS) caused by Pseudomonas syringae pv. syringae (Pss) causes significant losses to the pepper industry. Highly virulent Pss strains under optimal environmental conditions (cool-moderate temperatures, high moisture) can cause severe necrotic lesions on pepper leaves that consequently can decrease pepper yield if the disease persists. Hence, it is important to understand the virulence mechanisms of Pss to be able to effectively control PLS in peppers. In our study, in vitro, in planta, and whole-genome sequence analyses were conducted to better understand the virulence and pathogenicity characteristics of Pss strains in peppers. Our findings fill a knowledge gap regarding potential virulence and pathogenicity characteristics of Pss in peppers, including virulence and antimicrobial gene content. Our study helps pave a path to further identify the role of specific virulence genes in causing disease in peppers, which can have implications in developing strategies to effectively control PLS in peppers.

Impact of salinity and time on structure and functional potential of wastewater treatment biofilms in intermittent sand bioreactors.

High salt wastewater is produced in industries, including seafood and pickling processing. The salinity in such wastewaters has been shown to negatively impact biological treatment efficacy. Little is known about the changes in the microbial community structure in the mature biological 2 treatment systems, the impacts of salinity on community composition, and the shifts over time during operation. This study aimed to identify the changes in the microbial community due to both salt and days of operation through 16s rRNA sequencing and KEGG functional predictions. Intermittent sand bioreactors (ISBs) with a focus on ammonia treatment were utilized. Results showed that the overall community structure and diversity were distinct as wastewater salinity varied from 0%-1.3%. At 1.3% salinity Zoogloea, a common genus in wastewater treatment plants, was not present and Aequorovita, Thauera and Dokdonella became the dominant genera. Nitrosomonas, an important ammonia oxidizing bacteria, increased in abundance with days of operation but was not significantly impacted by an increase in salinity. This finding was further supported by an increase in predicted nitrification potential with time of operation within all intermittent sand bioreactors tested. These results provide a deeper understanding of the impacts of salinity on microbial community development in biological treatment systems and elucidate the shifts in community structure occurring during early operations and into system maturity.

Exploring the effects of transport duration on the fecal microbial communities of surplus dairy calves.

Transportation significantly affects the health and welfare of surplus dairy calves, largely due to the various stressors and pathogen exposures encountered during the process. Concurrently, an animal's microbiome is known to correlate with its health status, with stress-induced alterations in the microbiota potentially precipitating various diseases. This study aimed to compare the effects of transportation durations of 6, 12, or 16 h on the fecal microbiota in young surplus dairy calves. We used a randomized controlled design in which surplus dairy calves aged 1 to 19 d from 5 commercial dairy farms in Ontario were allocated into 1 of 3 transportation groups (6, 12, and 16 h of continuous transportation). Health assessments were conducted before, immediately after, and for 2 wk following transportation. Fecal samples were collected before, immediately after, and at 24 and 72 h after transportation and subjected to 16S rRNA sequencing. Alpha diversity metrics showed no significant differences between the 3 transportation groups at any of the sampling time points. Although ß diversity metrics revealed no clustering by transportation groups, they indicated significant differences across sampling time points within each group. The overall analysis revealed a total of 22 phyla and 353 genera, with Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, and Fusobacteria being the most abundant phyla. Bacteroides, Escherichia/Shigella, Lactobacillus, Collinsella, and Bifidobacterium were the most abundant genera. The reduction in Fusobacteria abundance before and after transport was significantly larger in the 16-h transportation group when compared with the 6-h transportation group. We also identified several genus-level and amplicon sequence variation-level taxa that displayed significant differences in their abundances across various transportation groups, observed at all sampling time points investigated. This research identifies microbiota changes due to varying transportation durations in surplus dairy calves, providing a broad understanding of the microbial shifts in surplus dairy calves after transportation across varying durations. Although these variations may not directly correlate with overall calf health or indicate dysbiosis, these results emphasize the importance of further investigating transportation practices to enhance calf health and well-being. Further studies are warranted to elucidate the relationship between microbiota and calf health.

Genomic Diversity, Antimicrobial Resistance, Plasmidome, and Virulence Profiles of Salmonella Isolated from Small Specialty Crop Farms Revealed by Whole-Genome Sequencing.

Salmonella is the leading cause of death associated with foodborne illnesses in the USA. Difficulty in treating human salmonellosis is attributed to the development of antimicrobial resistance and the pathogenicity of Salmonella strains. Therefore, it is important to study the genetic landscape of Salmonella, such as the diversity, plasmids, and presence antimicrobial resistance genes (AMRs) and virulence genes. To this end, we isolated Salmonella from environmental samples from small specialty crop farms (SSCFs) in Northeast Ohio from 2016 to 2021; 80 Salmonella isolates from 29 Salmonella-positive samples were subjected to whole-genome sequencing (WGS). In silico serotyping revealed the presence of 15 serotypes. AMR genes were detected in 15% of the samples, with 75% exhibiting phenotypic and genotypic multidrug resistance (MDR). Plasmid analysis demonstrated the presence of nine different types of plasmids, and 75% of AMR genes were located on plasmids. Interestingly, five Salmonella Newport isolates and one Salmonella Dublin isolate carried the ACSSuT gene cassette on a plasmid, which confers resistance to ampicillin, chloramphenicol, streptomycin, sulfonamide, and tetracycline. Overall, our results show that SSCFs are a potential reservoir of Salmonella with MDR genes. Thus, regular monitoring is needed to prevent the transmission of MDR Salmonella from SSCFs to humans.

Whole genome characterization of thermophilic Campylobacter species isolated from dairy manure in small specialty crop farms of Northeast Ohio.

Introduction: With more public interest in consuming locally grown produce, small specialty crop farms (SSCF) are a viable and growing segment of the food production chain in the United States. Methods: The goal of this study was to investigate the genomic diversity of Campylobacter isolated from dairy manure (n = 69) collected from 10 SSCF in Northeast Ohio between 2018 and 2020. Results: A total of 56 C. jejuni and 13 C. coli isolates were sequenced. Multi-locus sequence typing (MLST) identified 22 sequence types (STs), with ST-922 (18%) and ST-61 (13%) predominant in C. jejuni and ST-829 (62%) and ST-1068 (38%) predominant in C. coli. Interestingly, isolates with similar genomic and gene contents were detected within and between SSCF over time, suggesting that Campylobacter could be transmitted between farms and may persist in a given SSCF over time. Virulence-associated genes (n = 35) involved in the uptake and utilization of potassium and organic compounds (succinate, gluconate, oxoglutarate, and malate) were detected only in the C. jejuni isolates, while 45 genes associated with increased resistance to environmental stresses (capsule production, cell envelope integrity, and iron uptake) were detected only in the C. coli isolates. Campylobacter coli isolates were also sub-divided into two distinct clusters based on the presence of unique prophages (n = 21) or IncQ conjugative plasmid/type-IV secretion system genes (n = 15). Campylobacter coli isolates harbored genes associated with resistance to streptomycin (aadE-Cc; 54%) and quinolone (gyrA-T86I; 77%), while C. jejuni had resistance genes for kanamycin (aph3'-IIIa; 20%). Both species harbored resistance genes associated with ß-lactam (especially, blaOXA-193; up to 100%) and tetracycline (tetO; up to 59%). Discussion/Conclusion: Our study demonstrated that Campylobacter genome plasticity associated with conjugative transfer might provide resistance to certain antimicrobials and viral infections via the acquisition of protein-encoding genes involved in mechanisms such as ribosomal protection and capsule modification.

Estimation of species divergence times in presence of cross-species gene flow.

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent].

Dietary Impacts on Changes in Diversity and Abundance of the Murine Microbiome during Progression and Treatment of Cancer.

The intestinal microbial population is recognized for its impact on cancer treatment outcomes. Little research has reported microbiome changes during cancer progression or the interplay of disease progression, dietary sugar/fat intake, and the microbiome through surgery and chemotherapy. In this study, the murine gut microbiome was used as a model system, and changes in microbiome diversity, richness, and evenness over the progression of the cancer and treatment were analyzed. Mice were categorized into four diet cohorts, combinations of either high or low sucrose and high or low omega-3 fatty acids, and two treatment cohorts, saline vehicle or chemotherapy, for a total of eight groups. Fecal samples were collected at specific timepoints to assess changes due to diet implementation, onset of cancer, lumpectomy, and chemotherapy. Akkermansia muciniphila abundance was very high in some samples and negatively correlated with overall Amplicon Sequence Variant (ASV) richness (r(64) = -0.55, p = 3 × 10-8). Throughout the disease progression, ASV richness significantly decreased and was impacted by diet and treatment. Alpha-diversity and differential microbial abundance were significantly affected by disease progression, diet, treatment, and their interactions. These findings help establish a baseline for understanding how cancer progression, dietary macronutrients, and specific treatments impact the murine microbiome, which may influence outcomes.

Invasive Aedes japonicus Mosquitoes Dominate the Aedes Fauna Collected with Gravid Traps in Wooster, Northeastern Ohio, USA.

Aedes japonicus (Diptera: Culicidae), or the Asian rock pool mosquito, is an invasive mosquito in Europe and America. It was first detected outside of Asia in 1990 in Oceania. It has since expanded to North America and Europe in 1998 and 2000, respectively. Even though it is classified as a secondary vector of pathogens, it is competent to several arboviruses and filarial worms, and it is contributing to the transmission of La Crosse virus (LACV) and West Nile virus (WNV). In this study, CDC light, BG-sentinel, and gravid traps were used to collect mosquitoes between June and October 2021, in Wooster, Northeastern Ohio, USA. Morphological identification or/and Sanger sequencing were performed to identify the collected mosquitoes. Our results revealed that (adult) Ae. japonicus mosquitoes were the most abundant mosquito species collected with gravid traps in Wooster in 2021, confirming its establishment in Ohio. Molecular analyses of Ae. japonicus showed 100% nucleotide similarity with Ae. japonicus collected in Iowa (USA) and Canada, suggesting multiple introductions. Its presence may increase the risk of future arbovirus outbreaks in Wooster, Ohio. This study stresses the importance of actively monitoring the density and distribution of all members of the Ae. japonicus complex.

RADseq data reveal a lack of admixture in a mouse lemur contact zone contrary to previous microsatellite results.

Microsatellites have been a workhorse of evolutionary genetic studies for decades and are still commonly in use for estimating signatures of genetic diversity at the population and species level across a multitude of taxa. Yet, the very high mutation rate of these loci is a double-edged sword, conferring great sensitivity at shallow levels of analysis (e.g. paternity analysis) but yielding considerable uncertainty for deeper evolutionary comparisons. For the present study, we used reduced representation genome-wide data (restriction site-associated DNA sequencing (RADseq)) to test for patterns of interspecific hybridization previously characterized using microsatellite data in a contact zone between two closely related mouse lemur species in Madagascar (Microcebus murinus and Microcebus griseorufus). We revisit this system by examining populations in, near, and far from the contact zone, including many of the same individuals that had previously been identified as hybrids with microsatellite data. Surprisingly, we find no evidence for admixed nuclear ancestry. Instead, re-analyses of microsatellite data and simulations suggest that previously inferred hybrids were false positives and that the program NewHybrids can be particularly sensitive to erroneously inferring hybrid ancestry. Combined with results from coalescent-based analyses and evidence for local syntopic co-occurrence, we conclude that the two mouse lemur species are in fact completely reproductively isolated, thus providing a new understanding of the evolutionary rate whereby reproductive isolation can be achieved in a primate.

Pedigree-based and phylogenetic methods support surprising patterns of mutation rate and spectrum in the gray mouse lemur.

Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10-8 (95% credible interval: 1.28 × 10-8-1.78 × 10-8) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.

A vertebrate adaptive radiation is assembled from an ancient and disjunct spatiotemporal landscape.

To investigate the origins and stages of vertebrate adaptive radiation, we reconstructed the spatial and temporal histories of adaptive alleles underlying major phenotypic axes of diversification from the genomes of 202 Caribbean pupfishes. On a single Bahamian island, ancient standing variation from disjunct geographic sources was reassembled into new combinations under strong directional selection for adaptation to the novel trophic niches of scale-eating and molluscivory. We found evidence for two longstanding hypotheses of adaptive radiation: hybrid swarm origins and temporal stages of adaptation. Using a combination of population genomics, transcriptomics, and genome-wide association mapping, we demonstrate that this microendemic adaptive radiation of novel trophic specialists on San Salvador Island, Bahamas experienced twice as much adaptive introgression as generalist populations on neighboring islands and that adaptive divergence occurred in stages. First, standing regulatory variation in genes associated with feeding behavior (prlh, cfap20, and rmi1) were swept to fixation by selection, then standing regulatory variation in genes associated with craniofacial and muscular development (itga5, ext1, cyp26b1, and galr2) and finally the only de novo nonsynonymous substitution in an osteogenic transcription factor and oncogene (twist1) swept to fixation most recently. Our results demonstrate how ancient alleles maintained in distinct environmental refugia can be assembled into new adaptive combinations and provide a framework for reconstructing the spatiotemporal landscape of adaptation and speciation.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Ecology and morphology of mouse lemurs (Microcebus spp.) in a hotspot of microendemism in northeastern Madagascar, with the description of a new species.

Delimitation of cryptic species is increasingly based on genetic analyses but the integration of distributional, morphological, behavioral, and ecological data offers unique complementary insights into species diversification. We surveyed communities of nocturnal mouse lemurs (Microcebus spp.) in five different sites of northeastern Madagascar, measuring a variety of morphological parameters and assessing reproductive states for 123 individuals belonging to five different lineages. We documented two different non-sister lineages occurring in sympatry in two areas. In both cases, sympatric species pairs consisted of a locally restricted (M. macarthurii or M. sp. #3) and a more widespread lineage (M. mittermeieri or M. lehilahytsara). Estimated Extents of Occurrence (EOO) of these lineages differed remarkably with 560 and 1,500 km2 versus 9,250 and 50,700 km2 , respectively. Morphometric analyses distinguished unambiguously between sympatric species and detected more subtle but significant differences among sister lineages. Tail length and body size were most informative in this regard. Reproductive schedules were highly variable among lineages, most likely impacted by phylogenetic relatedness and environmental variables. While sympatric species pairs differed in their reproductive timing (M. sp. #3/M. lehilahytsara and M. macarthurii/M. mittermeieri), warmer lowland rainforests were associated with a less seasonal reproductive schedule for M. mittermeieri and M. lehilahytsara compared with populations occurring in montane forests. Distributional, morphological, and ecological data gathered in this study support the results of genomic species delimitation analyses conducted in a companion study, which identified one lineage, M. sp. #3, as meriting formal description as a new species. Consequently, a formal species description is included. Worryingly, our data also show that geographically restricted populations of M. sp. #3 and its sister species (M. macarthurii) are at high risk of local and perhaps permanent extinction from both deforestation and habitat fragmentation.

Molecular Clocks without Rocks: New Solutions for Old Problems.

Molecular data have been used to date species divergences ever since they were described as documents of evolutionary history in the 1960s. Yet, an inadequate fossil record and discordance between gene trees and species trees are persistently problematic. We examine how, by accommodating gene tree discordance and by scaling branch lengths to absolute time using mutation rate and generation time, multispecies coalescent (MSC) methods can potentially overcome these challenges. We find that time estimates can differ - in some cases, substantially - depending on whether MSC methods or traditional phylogenetic methods that apply concatenation are used, and whether the tree is calibrated with pedigree-based mutation rates or with fossils. We discuss the advantages and shortcomings of both approaches and provide practical guidance for data analysis when using these methods.

Conservation genomic analysis reveals ancient introgression and declining levels of genetic diversity in Madagascar's hibernating dwarf lemurs.

Madagascar's biodiversity is notoriously threatened by deforestation and climate change. Many of these organisms are rare, cryptic, and severely threatened, making population-level sampling unrealistic. Such is the case with Madagascar's dwarf lemurs (genus Cheirogaleus), the only obligate hibernating primate. We here apply comparative genomic approaches to generate the first genome-wide estimates of genetic diversity within dwarf lemurs. We generate a reference genome for the fat-tailed dwarf lemur, Cheirogaleus medius, and use this resource to facilitate analyses of high-coverage (~30×) genome sequences for wild-caught individuals representing species: C. sp. cf. medius, C. major, C. crossleyi, and C. sibreei. This study represents the largest contribution to date of novel genomic resources for Madagascar's lemurs. We find concordant phylogenetic relationships among the four lineages of Cheirogaleus across most of the genome, and yet detect a number of discordant genomic regions consistent with ancient admixture. We hypothesized that these regions could have resulted from adaptive introgression related to hibernation, indeed finding that genes associated with hibernation are present, though most significantly, that gene ontology categories relating to transcription are over-represented. We estimate levels of heterozygosity and find particularly low levels in an individual sampled from an isolated population of C. medius that we refer to as C. sp. cf. medius. Results are consistent with a recent decline in effective population size, which is evident across species. Our study highlights the power of comparative genomic analysis for identifying species and populations of conservation concern, as well as for illuminating possible mechanisms of adaptive phenotypic evolution.

Purifying Selection in Corvids Is Less Efficient on Islands.

Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.

Epistatic mutations under divergent selection govern phenotypic variation in the crow hybrid zone.

The evolution of genetic barriers opposing interspecific gene flow is key to the origin of new species. Drawing from information on over 400 admixed genomes sourced from replicate transects across the European hybrid zone between all-black carrion crows and grey-coated hooded crows, we decipher the interplay between phenotypic divergence and selection at the molecular level. Over 68% of plumage variation was explained by epistasis between the gene NDP and a ~2.8-megabase region on chromosome 18 with suppressed recombination. Both pigmentation loci showed evidence for divergent selection resisting introgression. This study reveals how few, large-effect loci can govern prezygotic isolation and shield phenotypic divergence from gene flow.

Don't throw out the sympatric speciation with the crater lake water: fine-scale investigation of introgression provides equivocal support for causal role of secondary gene flow in one of the clearest examples of sympatric speciation.

Genomic data has revealed complex histories of colonization and repeated gene flow previously unrecognized in some of the most celebrated examples of sympatric speciation and radiation. However, much of the evidence for secondary gene flow into these radiations comes from summary statistics calculated from sparse genomic sampling without knowledge of which specific genomic regions introgressed. This tells us little about how gene flow potentially influenced sympatric diversification. Here, we investigated whole genomes of Barombi Mbo crater lake cichlids for fine-scale patterns of introgression with neighboring riverine cichlid populations. We found evidence of secondary gene flow into the radiation scattered across <0.24% of the genome; however, from our analyses, it is not clear if the functional diversity in these regions contributed to the ecological, sexual, and morphological diversity found in the lake. Unlike similar studies, we found no obvious candidate genes for adaptive introgression and we cannot rule out that secondary gene flow was predominantly neutral with respect to the diversification process. We also found evidence for differential assortment of ancestral polymorphisms found in riverine populations between sympatric sister species, suggesting the presence of an ancestral hybrid swarm. Although the history of gene flow and colonization is more complicated than previously assumed, the lack of compelling evidence for secondary gene flow's role in species diversification suggests that we should not yet rule out one of the most celebrated examples of sympatric speciation in nature without a more thorough investigation of the timing and functional role of each introgressed region.

Speciation in sympatry with ongoing secondary gene flow and a potential olfactory trigger in a radiation of Cameroon cichlids.

The process of sympatric speciation in nature remains a fundamental unsolved problem. Cameroon crater lake cichlid radiations were long regarded as one of the most compelling examples; however, recent work showed that their origins were more complex than a single colonization event followed by isolation. Here, we performed a detailed investigation of the speciation history of a radiation of Coptodon cichlids from Lake Ejagham, Cameroon, using whole-genome sequencing data. The existence of the Lake Ejagham Coptodon radiation is remarkable as this 0.5 km2 lake offers limited scope for divergence across a shallow depth gradient, disruptive selection is currently weak, and the species are sexually monochromatic. We infer that Lake Ejagham was colonized by Coptodon cichlids soon after its formation 9,000 years ago, yet speciation occurred only in the last 1,000-2,000 years. We show that secondary gene flow from a nearby riverine species has been ongoing, into ancestral as well as extant lineages, and we identify and date river-to-lake admixture blocks. One block contains a cluster of olfactory receptor genes that introgressed near the time of the first speciation event and coincides with a higher overall rate of admixture. Olfactory signalling is a key component of mate choice and species recognition in cichlids. A functional role for this introgression event is consistent with previous findings that sexual isolation appears much stronger than ecological isolation in Ejagham Coptodon. We conclude that speciation in this radiation took place in sympatry, yet may have benefited from ongoing riverine gene flow.