A comprehensive review on the documented characteristics of four Reticulitermes termites (Rhinotermitidae, Blattodea) of China / Uma revisão abrangente sobre as características documentadas de quatro cupins Reticulitermes da China

Termites are known as social insects worldwide. Presently in China 473 species, 44 genera and 4 families of termites have been reported. Of them, 111 Reticulitermes species are widely spread in different zones of China. The dispersion flight season of these Chinese Reticulitermes species are usually started from February to June, but in some regions different species are distributed, sharing their boundaries and having overlapping flight seasons. These reasons become important sources of hybridization between two different heterospecific populations of termites. It was confirmed that the fertilized eggs and unfertilized eggs of some Reticulitermes termites have the capacity of cleavage. While the unfertilized eggs of R. aculabialis, R. chinensis and R. labralis cleaved normally and the only R. aculabialis unfertilized eggs develop in embryos. While, the R. flaviceps and R. chinensis were observed with their abnormal embryonic development, and not hatching of eggs parthenogenetically. They were reported more threatening to Chinese resources as they propagate with parthenogenesis, hybridization and sexual reproduction. Eggshell and macrophiles of eggs play important roles in species identification and control. Although, they are severe pests and cause a wide range of damages to wooden structures and products in homes, buildings, building materials, trees, crops, and forests in China's Mainland.

Os cupins são conhecidos como insetos sociais em todo o mundo. Atualmente na China foram relatadas 473 espécies, 44 gêneros e 4 famílias de cupins. Destas, 111 espécies de Reticulitermes estão amplamente distribuídas em diferentes zonas da China. A temporada de voo de dispersão dessas espécies chinesas de Reticulitermes geralmente começa de fevereiro a junho, mas em algumas regiões diferentes espécies são distribuídas, compartilhando seus limites e tendo temporadas de voo sobrepostas. Essas razões tornam-se importantes fontes de hibridização entre duas populações heteroespecíficas de cupins. Foi confirmado que os ovos fertilizados e não fertilizados de alguns cupins Reticulitermes possuem capacidade de clivagem. Já os ovos não fertilizados de R. aculabialis, R. chinensis e R. labralis clivaram normalmente, e os únicos ovos não fertilizados de R. aculabialis se desenvolvem em embriões. R. flaviceps e R. chinensis foram observados com desenvolvimento embrionário anormal, e não eclosão de ovos por partenogênese. Eles foram relatados como mais ameaçadores para os recursos chineses à medida que se propagam com partenogênese, hibridização e reprodução sexual. Casca de ovo e macrófilos de ovos desempenham papéis importantes na identificação e controle de espécies, embora sejam pragas graves e causem uma ampla gama de danos a estruturas e produtos de madeira em residências, edifícios, materiais de construção, árvores, plantações e florestas na China continental.

Genomic variability and population structure of six Colombian cattle breeds.

Analyses of the genetic diversity of indigenous cattle are essential for implementing conservation programs, promoting their sustainable use and maintaining productive advantages offered by these breeds in local conditions. The aim of this study was to investigate the genetic diversity and population structure of six Colombian cattle breeds: Blanco Orejinegro (BON), Costeño con Cuernos (CCC), Romosinuano (ROM), Sanmartinero (SAM), Casanareño (CAS), and Hartón del Valle (HDV). Two additional breed groups were included for comparison: Zebu (CEB) and a crossbreed of Colombian cattle breeds × Zebu. Genetic diversity within breeds was analyzed using expected heterozygosity (He), inbreeding coefficient (f), and runs of homozygosity (ROH). Population structure was assessed using model-based clustering (ADMIXTURE) and principal components analysis (PCA). Zebu cattle showed the lowest genetic diversity (He = 0.240). Breeds with the highest genetic diversity level were HDV and BON (He = 0.350 and 0.340, respectively). Inbreeding was lower for Colombian cattle breeds ranging between 0.005 and 0.045. Overall, the largest average genetic distance was found among the group of Colombian cattle breeds and Zebu, while the smallest was found between ROM and CCC. Model-based clustering revealed some level of admixture among HDV and CAS cattle which is consistent with their recent history. The results of the present study provide a useful insight on the genetic structure of Colombian cattle breeds.

Identification of fertility restorers for A2 cytoplasm-based CGMS lines and development of short duration hybrids in pigeonpea (Cajanus cajan(L.) Millsp.).

Pigeonpea is the second most important legume crop grown in India after chickpea. India is the largest producer of pigeonpea in the world. However, the productivity of pigeonpea in India remains stagnant over the years. The productivity of pigeonpea can be improved through exploitation of heterosis. The cytoplasmic genetic male sterility is the predominant method employed in hybrid development in pigeonpea during the recent days owing to the advantages involved. The present study involved the identification of fertility restorers for three Cajanus scarabaeoides(A2) based short duration (120-130 days) male sterile lines, namely CORG 990047A, CORG 990052A and CORG 7A. A total of 77 inbreds were involved in the hybridization programme. The pollen fertility of the 186 hybrids ranged from 0.00 to 94.89%. The independent confirmation of fertility restoration based on pollen fertility and pod set by selfing showed that, the hybrids, namely CORG 990047A 9 AK 261322, CORG 990052A 9 AK 261322 and CORG 7A 9 AK 261322 were identified as fertile. The inbred AK 261322 was the potential restorer of fertility in A2 male sterile lines. The hybrids, namely CORG 990047A 9 AK 261322 (35.19%), CORG 990052A 9 AK 261322 (12.75%) and CORG 7A 9 AK 261322 (19.77%) showed high heterosis for single plant yield over CO(Rg)7, a commercial check variety. The hybrids identified in the present study can be exploited for commercial cultivation after evaluation under various yield trials to estimate its performance. The polymorphic SSR markers identified in the present study can be utilized in future to assess the genetic purity of the hybrids.

A community genomics approach to natural hybridization.

Hybridization is a complicated, oft-misunderstood process. Once deemed unnatural and uncommon, hybridization is now recognized as ubiquitous among species. But hybridization rates within and among communities are poorly understood despite the relevance to ecology, evolution and conservation. To clarify, we examined hybridization across 75 freshwater fish communities within the Ozarks of the North American Interior Highlands (USA) by single nucleotide polymorphism (SNP) genotyping 33 species (N = 2865 individuals; double-digest restriction site-associated DNA sequencing (ddRAD)). We found evidence of hybridization (70 putative hybrids; 2.4% of individuals) among 18 species-pairs involving 73% (24/33) of study species, with the majority being concentrated within one family (Leuciscidae/minnows; 15 species; 66 hybrids). Interspecific genetic exchange-or introgression-was evident from 24 backcrossed individuals (10/18 species-pairs). Hybrids occurred within 42 of 75 communities (56%). Four selected environmental variables (species richness, protected area extent, precipitation (May and annually)) exhibited 73-78% accuracy in predicting hybrid occurrence via random forest classification. Our community-level assessment identified hybridization as spatially widespread and environmentally dependent (albeit predominantly within one diverse, omnipresent family). Our approach provides a more holistic survey of natural hybridization by testing a wide range of species-pairs, thus contrasting with more conventional evaluations.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator.

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis.

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Disomic 1M (1B) Triticum aestivum-Aegilops comosa Substitution Line with Favorable Protein Properties.

Aegilops comosa (2n = 2x = 14, MM) contains many excellent genes/traits for wheat breeding. Wheat-Ae. comosa introgression lines have potential value in the genetic improvement of wheat quality. A disomic 1M (1B) Triticum aestivum-Ae. comosa substitution line NAL-35 was identified by fluorescence in situ hybridization and genomic in situ hybridization analysis from a hybridization cross between a disomic 1M (1D) substitution line NB 4-8-5-9 with CS N1BT1D. The observation of pollen mother cells showed that NAL-35 had normal chromosome pairing, suggesting that NAL-35 could be used for the quality test. NAL-35 with alien Mx and My subunits showed positive effects on some protein-related parameters including high protein content and high ratios of high-molecular-weight glutenin subunits (HMW-GSs)/glutenin and HMW-GS/low-molecular-weight glutenin subunits. The changes in gluten composition improved the rheological properties of the dough of NAL-35, resulting in a tighter and more uniform microstructure. NAL-35 is a potential material for wheat quality improvement that transferred quality-related genes from Ae. comosa.

A comparison between the FlowCam 8100, microscopy, and sandwich hybridization assay for quantifying abundances of the saxitoxin-producing dinoflagellate, Alexandrium catenella.

Light microscopy, FlowCam, and sandwich hybridization assay (SHA) are three approaches that facilitate the monitoring of harmful algal bloom (HAB) forming phytoplankton. Yet, cross-comparisons among these techniques have not been conducted. This study addressed that gap using the saxitoxin-producing 'red tide' dinoflagellate Alexandrium catenella, a species responsible for blooms and paralytic shellfish poisoning worldwide. To achieve this goal, the dynamic ranges of each technique were compared using A. catenella cultures spanning low (pre-bloom), moderate (bloom), and high (dense bloom) levels. To assess field detection, water samples containing very low (<3 cells mL-1) A. catenella levels were collected from Long Island Sound, USA (Jun-Aug 2021) and evaluated using each method. Field samples were also spiked with A. catenella to high (160 cells mL-1) or low (40 cells mL-1) concentrations. In general, microscopy, FlowCam, and SHA returned comparable A. catenella cell concentrations for all tests. Mean cell concentrations from laboratory intercalibration experiments were not significantly different for any method or concentration (ANOVA, p > 0.05). However, relative to microscopy at times SHA produced non-detect signals <2 cells mL-1 in field samples and the FlowCam slightly underestimated cell concentrations when A. catenella abundances were high in laboratory and field samples. Mean cell concentrations of spike experiments were not significantly different for any test date, sampling location, or method, despite variability among methods within the high concentration treatment (ANOVA, p > 0.05 for all treatments). Findings are relevant to HAB researchers, managers, and public health officials because they help reconcile disparate cell abundance datasets that inform numerical models and enhance HAB monitoring and prediction. Results are also likely broadly applicable to several HAB species.

Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization.

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

The Genomic Variation and Differentially Expressed Genes on the 6P Chromosomes in Wheat-Agropyron cristatum Addition Lines 5113 and II-30-5 Confer Different Desirable Traits.

Wild relatives of wheat are essential gene pools for broadening the genetic basis of wheat. Chromosome rearrangements and genomic variation in alien chromosomes are widespread. Knowledge of the genetic variation between alien homologous chromosomes is valuable for discovering and utilizing alien genes. In this study, we found that 5113 and II-30-5, two wheat-A. cristatum 6P addition lines, exhibited considerable differences in heading date, grain number per spike, and grain weight. Genome resequencing and transcriptome analysis revealed significant differences in the 6P chromosomes of the two addition lines, including 143,511 single-nucleotide polymorphisms, 62,103 insertion/deletion polymorphisms, and 757 differentially expressed genes. Intriguingly, genomic variations were mainly distributed in the middle of the chromosome arms and the proximal centromere region. GO and KEGG analyses of the variant genes and differentially expressed genes showed the enrichment of genes involved in the circadian rhythm, carbon metabolism, carbon fixation, and lipid metabolism, suggesting that the differential genes on the 6P chromosome are closely related to the phenotypic differences. For example, the photosynthesis-related genes PsbA, PsbT, and YCF48 were upregulated in II-30-5 compared with 5113. ACS and FabG are related to carbon fixation and fatty acid biosynthesis, respectively, and both carried modification variations and were upregulated in 5113 relative to II-30-5. Therefore, this study provides important guidance for cloning desirable genes from alien homologous chromosomes and for their effective utilization in wheat improvement.

Molecular Cytological Analysis and Specific Marker Development in Wheat-Psathyrostachys huashanica Keng 3Ns Additional Line with Elongated Glume.

Psathyrostachys huashanica Keng (2n = 2x = 14, NsNs) is an excellent gene resource for wheat breeding, which is characterized by early maturity, low plant height, and disease resistance. The wheat-P. huashanica derivatives were created by the elite genes of P. huashanica and permeate into common wheat through hybridization. Among them, a long-glume material 20JH1155 was identified, with larger grains and longer spike than its parents. In the present study, the methods of cytological observation, GISH, and sequential FISH analysis showed that 20JH1155 contained 21 pairs of wheat chromosomes and a pair of P. huashanica. There were some differences in 5A and 7B chromosomes between 20JH1155 and parental wheat 7182. Molecular marker, FISH, and sequence cloning indicated 20JH1155 alien chromosomes were 3Ns of P. huashanica. In addition, differentially expressed genes during immature spikelet development of 20JH1155 and 7182 and predicted transcription factors were obtained by transcriptome sequencing. Moreover, a total of 7 makers derived from Ph#3Ns were developed from transcriptome data. Taken together, the wheat-P. huashanica derived line 20JH1155 provides a new horizon on distant hybridization of wheat and accelerates the utilization of genes of P. huashanica.

Regulatory Divergence as a Mechanism for X-Autosome Incompatibilities in Caenorhabditis Nematodes.

The world's astounding biodiversity results from speciation, the process of formation of distinct species. Hybrids between species often have reduced fitness due to negative epistatic interactions between divergent genetic factors, as each lineage accumulated substitutions independently in their evolutionary history. Such negative genetic interactions can manifest as gene misexpression due to divergence in gene regulatory controls from mutations in cis-regulatory elements and trans-acting factors. Gene misexpression due to differences in regulatory controls can ultimately contribute to incompatibility within hybrids through developmental defects such as sterility and inviability. We sought to quantify the contributions of regulatory divergence to postzygotic reproductive isolation using sterile interspecies hybrids of two Caenorhabditis nematodes: Caenorhabditis briggsae and Caenorhabditis nigoni. We analyzed previous transcriptome profiles for two introgression lines with distinct homozygous X-linked fragments from C. briggsae in a C. nigoni genomic background that confers male sterility, owing to defects in spermatogenesis (Li R, et al. 2016. Specific down-regulation of spermatogenesis genes targeted by 22G RNAs in hybrid sterile males associated with an X-chromosome introgression. Genome Res. 26:1219-1232). Our analysis identified hundreds of genes that show distinct classes of nonadditive expression inheritance and regulatory divergence. We find that these nonoverlapping introgressions affect many of the same genes in the same way and demonstrate that the preponderance of transgressive gene expression is due to regulatory divergence involving compensatory and joint effects of cis- and trans-acting factors. The similar transcriptomic responses to nonoverlapping genetic perturbations of the X-chromosome implicate multiway incompatibilities as an important feature contributing to hybrid male sterility in this system.

Genetic factors predict hybrid formation in the British flora.

Natural hybridization can have a profound evolutionary impact, with consequences ranging from the extinction of rare taxa to the origin of new species. Natural hybridization is particularly common in plants; however, our understanding of the general factors that promote or prevent hybridization is hampered by the highly variable outcomes in different lineages. Here, we quantify the influence of different predictors on hybrid formation across species from an entire flora. We combine estimates of hybridization with ecological attributes and a new species-level phylogeny for over 1,100 UK flowering plant species. Our results show that genetic factors, particularly parental genetic distance, as well as phylogenetic position and ploidy, are key determinants of hybrid formation, whereas many other factors such as range overlap and genus size explain much less variation in hybrid formation. Overall, intrinsic genetic factors shape the evolutionary and ecological consequences of natural hybridization across species in a flora.

Selection and adaptive introgression guided the complex evolutionary history of the European common bean.

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro's expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

Genomic prediction of rice mesocotyl length indicative of directing seeding suitability using a half-sib hybrid population.

Direct seeding has been widely adopted as an economical and labor-saving technique in rice production, though problems such as low seedling emergence rate, emergence irregularity and poor lodging resistance are existing. These problems are currently partially overcome by increasing seeding rate, however it is not acceptable for hybrid rice due to the high seed cost. Improving direct seeding by breeding is seen as the ultimate solution to these problems. For hybrid breeding, identifying superior hybrids among a massive number of hybrids from crossings between male and female parental populations by phenotypic evaluation is tedious and costly. Contrastingly, genomic selection/prediction (GS/GP) could efficiently detect the superior hybrids capitalizing on genomic data, which holds a great potential in plant hybrids breeding. In this study, we utilized 402 rice inbred varieties and 401 hybrids to investigate the effectiveness of GS on rice mesocotyl length, a representative indicative trait of direct seeding suitability. Several GP methods and training set designs were studied to seek the optimal scenario of hybrid prediction. It was shown that using half-sib hybrids as training set with the phenotypes of all parental lines being fitted as a covariate could optimally predict mesocotyl length. Partitioning the molecular markers into trait-associated and -unassociated groups based on genome-wide association study using all parental lines and hybrids could further improve the prediction accuracy. This study indicates that GS could be an effective and efficient method for hybrid breeding for rice direct seeding.

A new hypothesis for the origin of the lager yeast Saccharomyces pastorianus.

Saccharomyces pastorianus, which is responsible for the production of bottom-fermented lager beer, is a hybrid species that arose from the mating of the top-fermenting ale yeast Saccharomyces cerevisiae and the cold-tolerant Saccharomyces eubayanus around the start of the 17th century. Based on detailed analysis of Central European brewing records, we propose that the critical event for the hybridization was the introduction of top-fermenting S. cerevisiae into an environment where S. eubayanus was present, rather than the other way around. Bottom fermentation in parts of Bavaria preceded the proposed hybridization date by a couple of hundred years and we suggest that this was carried out by mixtures of yeasts, which may have included S. eubayanus. A plausible case can be made that the S. cerevisiae parent came either from the Schwarzach wheat brewery or the city of Einbeck, and the formation of S. pastorianus happened in the Munich Hofbräuhaus between 1602 and 1615 when both wheat beer and lager were brewed contemporaneously. We also describe how the distribution of strains from the Munich Spaten brewery, and the development by Hansen and Linder of methods for producing pure starter cultures, facilitated the global spread of the Bavarian S. pastorianus lineages.

Meiotic Recognition of Evolutionarily Diverged Homologs: Chromosomal Hybrid Sterility Revisited.

Hybrid sterility (HS) is an early postzygotic reproductive isolation mechanism observed in all sexually reproducing species. Infertility of hybrids prevents gene flow between incipient species and leads to speciation. While Drosophila studies have focused almost exclusively on the genic control of HS, two other model species, Mus musculus and budding yeast, provided the first experimental evidence of hybrid sterility governed by the nongenic effects of DNA sequence divergence. Here, we propose that the nongenic effect of increasing DNA divergence between closely related species may impair mutual recognition of homologous chromosomes and disrupt their synapsis. Unsynapsed or mispaired homologs can induce early meiotic arrest, or their random segregation can cause aneuploidy of spermatids and sperm cells. Impaired recognition of homologs may thus act as a universal chromosomal checkpoint contributing to the complexity of genetic control of HS. Chromosomal HS controlled by the Prdm9 gene in mice and HS driven by the mismatch repair machinery in yeast are currently the most advanced examples of chromosomal homology search-based HS. More focus on the cellular and molecular phenotypes of meiosis will be needed to further validate the role of homolog recognition in hybrid sterility and speciation.

Multiple colonizations and hybridization of a freshwater fish group on a satellite island of Sulawesi.

Repeated colonizations and resultant hybridization may increase lineage diversity on an island if introgression occurs only in a portion of the indigenous island lineage. Therefore, to precisely understand how island biodiversity was shaped, it is essential to reconstruct the history of secondary colonization and resultant hybridization both in time and space. In this study, we reconstructed the history of multiple colonizations of the Oryzias woworae species group, a freshwater fish group of the family Adrianichthyidae, from Sulawesi Island to its southeast satellite island, Muna Island. Phylogenetic and species tree analyses using genome-wide single-nucleotide polymorphisms revealed that all local populations on Muna Island were monophyletic, but that there were several genetically distinct lineages within the island. Population structure and phylogenetic network analyses demonstrated that colonization of this island occurred more than once, and that secondary colonization and resultant introgressive hybridization occurred only in one local population on the island. The spatially heterogeneous introgression induced by the multiple colonizations were also supported by differential admixture analyses. In addition, the differential admixture analyses detected reverse colonization from Muna Island to the Sulawesi mainland. Coalescence-based demographic inference estimated that these mutual colonizations occurred during the middle to late Quaternary period, during which sea level repeatedly declined; this indicates that the colonizations occurred via land bridges. We conclude that these mutual colonizations between Muna Island and the Sulawesi mainland, and the resultant spatially heterogeneous introgression shaped the current biodiversity of this species group in this area.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.).

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.