Molecular characterization of emerging Echovirus 11 (E11) shed light on the recombinant origin of a variant associated with severe hepatitis in neonates.

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. Due to the limited data available, the World Health Organization (WHO) considers public health risk to the general population to be low. The present study investigated the genetic variation and molecular evolution of E11 genomes collected from May to December 2023. Whole genome sequencing (WGS) was performed for 16 E11 strains. Phylogenetic analysis on WG showed how all Italian strains belonged to genogroup D5, similarly to other E11 strains recently reported in France and Germany all together aggregated into separate clusters. A cluster-specific recombination pattern was also identified using phylogenetic analysis of different genome regions. Echovirus 6 was identified as the major recombinant virus in 3Cpro and 3Dpol regions. The molecular clock analysis revealed that the recombination event probably occurred in June 2018 (95% HPD interval: Jan 2016-Jan 2020). Shannon entropy analyses, within P1 region, showed how 11 amino acids exhibited relatively high entropy. Five of them were exposed on the canyon region which is responsible for receptor binding with the neonatal Fc receptor. The present study showed the recombinant origin of a new lineage of E11 associated with severe neonatal infections.

Genetic recombination-mediated evolutionary interactions between phages of potential industrial importance and prophages of their hosts within or across the domains of Escherichia, Listeria, Salmonella, Campylobacter, and Staphylococcus.

BACKGROUND: The in-depth understanding of the role of lateral genetic transfer (LGT) in phage-prophage interactions is essential to rationalizing phage applications for human and animal therapy, as well as for food and environmental safety. This in silico study aimed to detect LGT between phages of potential industrial importance and their hosts. METHODS: A large array of genetic recombination detection algorithms, implemented in SplitsTree and RDP4, was applied to detect LGT between various Escherichia, Listeria, Salmonella, Campylobacter, Staphylococcus, Pseudomonas, and Vibrio phages and their hosts. PHASTER and RAST were employed respectively to identify prophages across the host genome and to annotate LGT-affected genes with unknown functions. PhageAI was used to gain deeper insights into the life cycle history of recombined phages. RESULTS: The split decomposition inferences (bootstrap values: 91.3-100; fit: 91.433-100), coupled with the Phi (0.0-2.836E-12) and RDP4 (P being well below 0.05) statistics, provided strong evidence for LGT between certain Escherichia, Listeria, Salmonella, and Campylobacter virulent phages and prophages of their hosts. The LGT events entailed mainly the phage genes encoding for hypothetical proteins, while some of these genetic loci appeared to have been affected even by intergeneric recombination in specific E. coli and S. enterica virulent phages when interacting with their host prophages. Moreover, it is shown that certain L. monocytogenes virulent phages could serve at least as the donors of the gene loci, involved in encoding for the basal promoter specificity factor, for L. monocytogenes. In contrast, the large genetic clusters were determined to have been simultaneously exchanged by many S. aureus prophages and some Staphylococcus temperate phages proposed earlier as potential therapeutic candidates (in their native or modified state). The above genetic clusters were found to encompass multiple genes encoding for various proteins, such as e.g., phage tail proteins, the capsid and scaffold proteins, holins, and transcriptional terminator proteins. CONCLUSIONS: It is suggested that phage-prophage interactions, mediated by LGT (including intergeneric recombination), can have a far-reaching impact on the co-evolutionary trajectories of industrial phages and their hosts especially when excessively present across microbially rich environments.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure.

BACKGROUND: Shallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon's true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India. RESULTS: Shallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia. CONCLUSION: This study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Tandem gene duplication selected by activation of horizontally transferred gene in bacteria.

Horizontal gene transfer occurs frequently in bacteria, but the mechanism driving activation and optimization of the expression of horizontally transferred genes (HTGs) in new recipient strains is not clear. Our previous study found that spontaneous tandem DNA duplication resulted in rapid activation of HTGs. Here, we took advantage of this finding to develop a novel technique for tandem gene duplication, named tandem gene duplication selected by activation of horizontally transferred gene in bacteria (TDAH), in which tandem duplication was selected by the activation of horizontally transferred selectable marker gene. TDAH construction does not contain any reported functional elements based on homologous or site-specific recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection and 9-bp-microhomology border sequences for precise illegitimate recombination. One transformation and 3 days were enough to produce a high-copy strain, so its procedure is simple and fast. Without subsequent knockout of the endogenous recombination system, TDAH could also generate the relatively stable high-copy tandem duplication for plasmid-carried and genome-integrated DNA. TDAH also showed an excellent capacity for increase gene expression and worked well in different industrial bacteria. We also applied TDAH to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5 times and remained stable even after 12 subcultures. TDAH is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. KEY POINTS: • We develop a novel and efficient technique (TDAH) for tandem gene duplication in bacterium. TDAH is based on the mechanism of HTG rapid activation. TDAH does not contain any reported functional elements based on homologous recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection, so its construction and procedure are very simple and fast. • TDAH is the first reported selected and stable tandem-gene-duplication technique in which the selected high-copy plasmid-carried and genome-integrated DNA could remain stable without the subsequent knockout of recombination system. • TDAH showed an excellent capacity for regulating gene expression and worked well in different industrial bacteria, indicating it is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. • TDAH was applied to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5-fold and remained stable even after 12 subcultures.

Far-East Asian Toxoplasma isolates share ancestry with North and South/Central American recombinant lineages.

Toxoplasma gondii is a global protozoan pathogen. Clonal lineages predominate in Europe, North America, Africa, and China, whereas highly recombinant parasites are endemic in South/Central America. Far East Asian T. gondii isolates are not included in current global population genetic structure analyses at WGS resolution. Here we report a genome-wide population study that compared eight Japanese and two Chinese isolates against representative worldwide T. gondii genomes using POPSICLE, a novel population structure analyzing software. Also included were 7 genomes resurrected from non-viable isolates by target enrichment sequencing. Visualization of the genome structure by POPSICLE shows a mixture of Chinese haplogroup (HG) 13 haploblocks introgressed within the genomes of Japanese HG2 and North American HG12. Furthermore, two ancestral lineages were identified in the Japanese strains; one lineage shares a common ancestor with HG11 found in both Japanese strains and North American HG12. The other ancestral lineage, found in T. gondii isolates from a small island in Japan, is admixed with genetically diversified South/Central American strains. Taken together, this study suggests multiple ancestral links between Far East Asian and American T. gondii strains and provides insight into the transmission history of this cosmopolitan organism.

Predicting class switch recombination in B-cells from antibody repertoire data.

Statistical and machine learning methods have proved useful in many areas of immunology. In this paper, we address for the first time the problem of predicting the occurrence of class switch recombination (CSR) in B-cells, a problem of interest in understanding antibody response under immunological challenges. We propose a framework to analyze antibody repertoire data, based on clonal (CG) group representation in a way that allows us to predict CSR events using CG level features as input. We assess and compare the performance of several predicting models (logistic regression, LASSO logistic regression, random forest, and support vector machine) in carrying out this task. The proposed approach can obtain an unweighted average recall of 71 % $71\%$ with models based on variable region descriptors and measures of CG diversity during an immune challenge and, most notably, before an immune challenge.

Screening Key Sites of Class 2 Integron Integrase that Impact Recombination Efficiency.

By capturing and expressing exogenous resistance gene cassettes through site-specific recombination, integrons play important roles in the horizontal transfer of antimicrobial resistant genes among bacteria. The characteristics of integron integrase make it to be a potential gene editing tool enzyme. In this study, a random mutation library using error-prone PCR was constructed, and amino acid residues mutants that impact on attI2 × attC or attC × attC recombination efficiency were screened and analyzed. Thirteen amino acid mutations were identified to be critical impacted on site-specific recombination of IntI2, including the predicted catalyzed site Y301. Nine of 13 mutated amino acid residues that have critically impacted on IntI2 activity were relative concentrated and near the predicted catalyzed site Y301 in the predicted three-dimensional structure indicated the importance of this area in maintain the activity of IntI2. No mutant with obviously increased recombination activity (more than four-fold as high as that of wild IntI2) was found in library screening, except P95S, R100K slightly increased (within two-fold) the excision activity of IntI2, and S243T slightly increased (within two-fold) both excision and integration activity of IntI2. These findings will provide clues for further specific modification of integron integrase to be a tool enzyme as well as establishing a new gene editing system and applied practically.

High prevalence of PRDM9-independent recombination hotspots in placental mammals.

In many mammals, recombination events are concentrated in hotspots directed by a sequence-specific DNA-binding protein named PRDM9. Intriguingly, PRDM9 has been lost several times in vertebrates, and notably among mammals, it has been pseudogenized in the ancestor of canids. In the absence of PRDM9, recombination hotspots tend to occur in promoter-like features such as CpG islands. It has thus been proposed that one role of PRDM9 could be to direct recombination away from PRDM9-independent hotspots. However, the ability of PRDM9 to direct recombination hotspots has been assessed in only a handful of species, and a clear picture of how much recombination occurs outside of PRDM9-directed hotspots in mammals is still lacking. In this study, we derived an estimator of past recombination activity based on signatures of GC-biased gene conversion in substitution patterns. We quantified recombination activity in PRDM9-independent hotspots in 52 species of boreoeutherian mammals. We observe a wide range of recombination rates at these loci: several species (such as mice, humans, some felids, or cetaceans) show a deficit of recombination, while a majority of mammals display a clear peak of recombination. Our results demonstrate that PRDM9-directed and PRDM9-independent hotspots can coexist in mammals and that their coexistence appears to be the rule rather than the exception. Additionally, we show that the location of PRDM9-independent hotspots is relatively more stable than that of PRDM9-directed hotspots, but that PRDM9-independent hotspots nevertheless evolve slowly in concert with DNA hypomethylation.

Lytic and Latent Genetic Diversity of the Epstein-Barr Virus Reveals Raji-Related Variants from Southeastern Brazil Associated with Recombination Markers.

Epstein-Barr virus (EBV) is a ubiquitous gammaherpesvirus etiologically associated with benign and malignant diseases. Since the pathogenic mechanisms of EBV are not fully understood, understanding EBV genetic diversity is an ongoing goal. Therefore, the present work describes the genetic diversity of the lytic gene BZLF1 in a sampling of 70 EBV-positive cases from southeastern Brazil. Additionally, together with the genetic regions previously characterized, the aim of the present study was to determine the impact of viral genetic factors that may influence EBV genetic diversity. Accordingly, the phylogenetic analysis of the BZLF1 indicated two main clades with high support, BZ-A and BZ-B (PP > 0.85). Thus, the BZ-A clade was the most diverse clade associated with the main polymorphisms investigated, including the haplotype Type 1 + V3 (p < 0.001). Furthermore, the multigene phylogenetic analysis (MLA) between BZLF1 and the oncogene LMP1 showed specific clusters, revealing haplotypic segregation that previous single-gene phylogenies from both genes failed to demonstrate. Surprisingly, the LMP1 Raji-related variant clusters were shown to be more diverse, associated with BZ-A/B and the Type 2/1 + V3 haplotypes. Finally, due to the high haplotypic diversity of the Raji-related variants, the number of DNA recombination-inducing motifs (DRIMs) was evaluated within the different clusters defined by the MLA. Similarly, the haplotype BZ-A + Raji was shown to harbor a greater number of DRIMs (p < 0.001). These results call attention to the high haplotype diversity of EBV in southeast Brazil and strengthen the hypothesis of the recombinant potential of South American Raji-related variants via the LMP1 oncogene.

BRCA1 safeguards genome integrity by activating chromosome asynapsis checkpoint to eliminate recombination-defective oocytes.

In the meiotic prophase, programmed DNA double-strand breaks are repaired by meiotic recombination. Recombination-defective meiocytes are eliminated to preserve genome integrity in gametes. BRCA1 is a critical protein in somatic homologous recombination, but studies have suggested that BRCA1 is dispensable for meiotic recombination. Here we show that BRCA1 is essential for meiotic recombination. Interestingly, BRCA1 also has a function in eliminating recombination-defective oocytes. Brca1 knockout (KO) rescues the survival of Dmc1 KO oocytes far more efficiently than removing CHK2, a vital component of the DNA damage checkpoint in oocytes. Mechanistically, BRCA1 activates chromosome asynapsis checkpoint by promoting ATR activity at unsynapsed chromosome axes in Dmc1 KO oocytes. Moreover, Brca1 KO also rescues the survival of asynaptic Spo11 KO oocytes. Collectively, our study not only unveils an unappreciated role of chromosome asynapsis in eliminating recombination-defective oocytes but also reveals the dual functions of BRCA1 in safeguarding oocyte genome integrity.

The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird.

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

Bridging the gap between the evolutionary dynamics and the molecular mechanisms of meiosis: A model based exploration of the PRDM9 intra-genomic Red Queen.

Molecular dissection of meiotic recombination in mammals, combined with population-genetic and comparative studies, have revealed a complex evolutionary dynamic characterized by short-lived recombination hotspots. Hotspots are chromosome positions containing DNA sequences where the protein PRDM9 can bind and cause crossing-over. To explain these fast evolutionary dynamic, a so-called intra-genomic Red Queen model has been proposed, based on the interplay between two antagonistic forces: biased gene conversion, mediated by double-strand breaks, resulting in hotspot extinction (the hotspot conversion paradox), followed by positive selection favoring mutant PRDM9 alleles recognizing new sequence motifs. Although this model predicts many empirical observations, the exact causes of the positive selection acting on new PRDM9 alleles is still not well understood. In this direction, experiment on mouse hybrids have suggested that, in addition to targeting double strand breaks, PRDM9 has another role during meiosis. Specifically, PRDM9 symmetric binding (simultaneous binding at the same site on both homologues) would facilitate homology search and, as a result, the pairing of the homologues. Although discovered in hybrids, this second function of PRDM9 could also be involved in the evolutionary dynamic observed within populations. To address this point, here, we present a theoretical model of the evolutionary dynamic of meiotic recombination integrating current knowledge about the molecular function of PRDM9. Our modeling work gives important insights into the selective forces driving the turnover of recombination hotspots. Specifically, the reduced symmetrical binding of PRDM9 caused by the loss of high affinity binding sites induces a net positive selection eliciting new PRDM9 alleles recognizing new targets. The model also offers new insights about the influence of the gene dosage of PRDM9, which can paradoxically result in negative selection on new PRDM9 alleles entering the population, driving their eviction and thus reducing standing variation at this locus.

Expanding the Scope of Adenoviral Vectors by Utilizing Novel Tools for Recombination and Vector Rescue.

Recombinant adenoviruses are widely used in clinical and laboratory applications. Despite the wide variety of available sero- and genotypes, only a fraction is utilized in vivo. As adenoviruses are a large group of viruses, displaying many different tropisms, immune epitopes, and replication characteristics, the merits of translating these natural benefits into vector applications are apparent. This translation, however, proves difficult, since while research has investigated the application of these viruses, there are no universally applicable rules in vector design for non-classical adenovirus types. In this paper, we describe a generalized workflow that allows vectorization, rescue, and cloning of all adenoviral species to enable the rapid development of new vector variants. We show this using human and simian adenoviruses, further modifying a selection of them to investigate their gene transfer potential and build potential vector candidates for future applications.

Development of a Cell Culture Model for Inducible SARS-CoV-2 Replication.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) induces direct cytopathic effects, complicating the establishment of low-cytotoxicity cell culture models for studying its replication. We initially developed a DNA vector-based replicon system utilizing the CMV promoter to generate a recombinant viral genome bearing reporter genes. However, this system frequently resulted in drug resistance and cytotoxicity, impeding model establishment. Herein, we present a novel cell culture model with SARS-CoV-2 replication induced by Cre/LoxP-mediated DNA recombination. An engineered SARS-CoV-2 transcription unit was subcloned into a bacterial artificial chromosome (BAC) vector. To enhance biosafety, the viral spike protein gene was deleted, and the nucleocapsid gene was replaced with a reporter gene. An exogenous sequence was inserted within NSP1 as a modulatory cassette that is removable after Cre/LoxP-mediated DNA recombination and subsequent RNA splicing. Using the PiggyBac transposon strategy, the transcription unit was integrated into host cell chromatin, yielding a stable cell line capable of inducing recombinant SARS-CoV-2 RNA replication. The model exhibited sensitivity to the potential antivirals forsythoside A and verteporfin. An innovative inducible SARS-CoV-2 replicon cell model was introduced to further explore the replication and pathogenesis of the virus and facilitate screening and assessment of anti-SARS-CoV-2 therapeutics.

Identification of the novel HLA-A*02:829 allele by sequencing-based typing.

The novel HLA-A*02:829 allele is likely generated from the recombination of both HLA-A*02:05:01:01 and HLA-A*32:01:01:01 alleles.

Genome comparison of long-circulating field CnmeGV isolates from the same region.

Cnaphalocrocis medinalis granulovirus (CnmeGV), belonging to Betabaculovirus cnamedinalis, can infect the rice pest, the rice leaf roller. In 1979, a CnmeGV isolate, CnmeGV-EP, was collected from Enping County, China. In 2014, we collected another CnmeGV isolate, CnmeGV-EPDH3, at the same location and obtained the complete virus genome sequence using Illumina and ONT sequencing technologies. By combining these two virus isolates, we updated the genome annotation of CnmeGV and conducted an in-depth analysis of its genome features. CnmeGV genome contains abundant tandem repeat sequences, and the repeating units in the homologous regions (hrs) exhibit overlapping and nested patterns. The genetic variations within EPDH3 population show the high stability of CnmeGV genome, and tandem repeats are the only region of high genetic variation in CnmeGV genome replication. Some defective viral genomes formed by recombination were found within the population. Comparison analysis of the two virus isolates collected from Enping showed that the proteins encoded by the CnmeGV-specific genes were less conserved relative to the baculovirus core genes. At the genomic level, there are a large number of SNPs and InDels between the two virus isolates, especially in and around the bro genes and hrs. Additionally, we discovered that CnmeGV acquired a segment of non-ORF sequence from its host, which does not provide any new proteins but rather serves as redundant genetic material integrated into the viral genome. Furthermore, we observed that the host's transposon piggyBac has inserted into some virus genes. Together, dsDNA viruses could acquire non-coding genetic material from their hosts to expand the size of their genomes. These findings provide new insights into the evolution of dsDNA viruses.

Fidelity Characterization of Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus and NADC30-like Strain.

The porcine reproductive and respiratory syndrome virus (PRRSV) has significantly impacted the global pork industry for over three decades. Its high mutation rates and frequent recombination greatly intensifies its epidemic and threat. To explore the fidelity characterization of Chinese highly pathogenic PRRSV JXwn06 and the NADC30-like strain CHsx1401, self-recombination and mutation in PAMs, MARC-145 cells, and pigs were assessed. In vitro, CHsx1401 displayed a higher frequency of recombination junctions and a greater diversity of junction types than JXwn06. In vivo, CHsx1401 exhibited fewer junction types yet maintained a higher junction frequency. Notably, JXwn06 showed more accumulation of mutations. To pinpoint the genomic regions influencing their fidelity, chimeric viruses were constructed, with the exchanged nsp9-10 regions between JXwn06 and CHsx1401. The SJn9n10 strain, which incorporates JXwn06's nsp9-10 into the CHsx1401 genome, demonstrated reduced sensitivity to nucleotide analogs compared to CHsx1401. Conversely, compared with JXwn06, the JSn9n10 strain showed increased sensitivity to these inhibitors. The swapped nsp9-10 also influences the junction frequency and accumulated mutations as their donor strains. The results indicate a propensity for different types of genetic variations between these two strains and further highlight the nsp9-10 region as a critical determinant of their fidelity.

Evolutionary Profile of Mayaro Virus in the Americas: An Update into Genome Variability.

The Mayaro virus (MAYV) is an arbovirus with emerging potential, though with a limited understanding of its epidemiology and evolution due to the lack of studies and surveillance. Here, we investigated 71 MAYV genome sequences from the Americas available at GenBank and characterized the phylogenetic relationship among virus strains. A phylogenetic analysis showed that sequences were grouped according to the genotypes L, D, and N. Genotype D sequences were closely related to sequences collected in adjacent years and from their respective countries, suggesting that isolates may have originated from circulating lineages. The coalescent analysis demonstrated similar results, indicating the continuous circulation of the virus between countries as well. An unidentified sequence from the USA was grouped with genotype D, suggesting the insertion of this genotype in the country. Furthermore, the recombination analysis detected homologous and three heterologous hybrids which presented an insertion into the nsP3 protein. Amino acid substitutions among sequences indicated selective pressure sites, suggesting viral adaptability. This also impacted the binding affinity between the E1-E2 protein complex and the Mxra8 receptor, associated with MAYV entry into human cells. These results provide information for a better understanding of genotypes circulating in the Americas.

Tools for Cre-Mediated Conditional Deletion of Floxed Alleles from Developing Cerebellar Purkinje Cells.

The Cre-lox system is an indispensable tool in neuroscience research for targeting gene deletions to specific cellular populations. Here we assess the utility of several transgenic Cre lines, along with a viral approach, for targeting cerebellar Purkinje cells (PCs) in mice. Using a combination of a fluorescent reporter line (Ai14) to indicate Cre-mediated recombination and a floxed Dystroglycan line (Dag1flox ), we show that reporter expression does not always align precisely with loss of protein. The commonly used Pcp2Cre line exhibits a gradual mosaic pattern of Cre recombination in PCs from Postnatal Day 7 (P7) to P14, while loss of Dag1 protein is not complete until P30. Ptf1aCre drives recombination in precursor cells that give rise to GABAergic neurons in the embryonic cerebellum, including PCs and molecular layer interneurons. However, due to its transient expression in precursors, Ptf1aCre results in stochastic loss of Dag1 protein in these neurons. NestinCre , which is often described as a "pan-neuronal" Cre line for the central nervous system, does not drive Cre-mediated recombination in PCs. We identify a Calb1Cre line that drives efficient and complete recombination in embryonic PCs, resulting in loss of Dag1 protein before the period of synaptogenesis. AAV8-mediated delivery of Cre at P0 results in gradual transduction of PCs during the second postnatal week, with loss of Dag1 protein not reaching appreciable levels until P35. These results characterize several tools for targeting conditional deletions in cerebellar PCs at different developmental stages and illustrate the importance of validating the loss of protein following recombination.

Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra.

Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.