Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages

BACKGROUND

Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation.

METHODS

We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination.

RESULTS

In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties.

CONCLUSIONS

We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

Controlling Recombination to Evolve Bacteriophages

Recombination among different phages sometimes facilitates their ability to grow on new hosts. Protocols to direct the evolution of phage host range, as might be used in the application of phage therapy, would then benefit from including steps to enable recombination. Applying mathematical and computational models, in addition to experiments using phages T3 and T7, we consider ways that a protocol may influence recombination levels. We first address coinfection, which is the first step to enabling recombination. The multiplicity of infection (MOI, the ratio of phage to cell concentration) is insufficient for predicting (co)infection levels. The force of infection (the rate at which cells are infected) is also critical but is more challenging to measure. Using both a high force of infection and high MOI (>1) for the different phages ensures high levels of coinfection. We also apply a four-genetic-locus model to study protocol effects on recombinant levels. Recombinants accumulate over multiple generations of phage growth, less so if one phage outgrows the other. Supplementing the phage pool with the low-fitness phage recovers some of this 'lost' recombination. Overall, fine tuning of phage recombination rates will not be practical with wild phages, but qualitative enhancement can be attained with some basic procedures.

Methylomes as key features for predicting recombination in some plant species

Knowing how chromosome recombination works is essential for plant breeding. It enables the design of crosses between different varieties to combine desirable traits and create new ones. This is because the meiotic crossovers between homologous chromatids are not purely random, and various strategies have been developed to describe and predict such exchange events. Recent studies have used methylation data to predict chromosomal recombination in rice using machine learning models. This approach proved successful due to the presence of a positive correlation between the CHH context cytosine methylation and recombination rates in rice chromosomes. This paper assesses the question if methylation can be used to predict recombination in four plant species: Arabidopsis, maize, sorghum, and tomato. The results indicate a positive association between CHH context methylation and recombination rates in certain plant species, with varying degrees of strength in their relationships. The CG and CHG methylation contexts show negative correlation with recombination. Methylation data was key effectively in predicting recombination in sorghum and tomato, with a mean determination coefficient of 0.65 ± 0.11 and 0.76 ± 0.05, respectively. In addition, the mean correlation values between predicted and experimental recombination rates were 0.83 ± 0.06 for sorghum and 0.90 ± 0.05 for tomato, confirming the significance of methylomes in both monocotyledonous and dicotyledonous species. The predictions for Arabidopsis and maize were not as accurate, likely due to the comparatively weaker relationships between methylation contexts and recombination, in contrast to sorghum and tomato, where stronger associations were observed. To enhance the accuracy of predictions, further evaluations using data sets closely related to each other might prove beneficial. In general, this methylome-based method holds great potential as a reliable strategy for predicting recombination rates in various plant species, offering valuable insights to breeders in their quest to develop novel and improved varieties.

Chromosome recombination and modification by LoxP-mediated evolution in Vibrio natriegens using CRISPR-associated transposases

Chromosome rearrangement by LoxP-mediated evolution has emerged as a powerful approach to studying how chromosome architecture impacts phenotypes. However, it relies on the in vitro synthesis of artificial chromosomes. The recently reported CRISPR-associated transposases (CASTs) held great promise for the efficient insertion of abundant LoxP sites directly onto the genome of wild-type strains. In this study, with the fastest-growing bacterium Vibrio natrigens (V. natriegens) as an object, a multiplex genome integration tool derived from CASTs was employed to achieve the insertion of cargo genes at eight specific genomic loci within 2 days. Next, we introduced 30 LoxP sites onto chromosome 2 (Chr2) of V. natriegens. Rigorously induced Cre recombinase was used to demonstrate Chromosome Rearrangement and Modification by LoxP-mediated Evolution (CRaMbLE). Growth characterization and genome sequencing showed that the ~358 kb fragment on Chr2 was accountable for the rapid growth of V. natriegens. The enabling tools we developed can help identify genomic regions that influence the rapid growth of V. natriegens without a prior understanding of genome mechanisms. This groundbreaking demonstration may also be extended to other organisms such as Escherichia coli, Pseudomonas putida, Bacillus subtilis, and so on.

Molecular underpinnings and environmental drivers of loss of heterozygosity in Drosophila intestinal stem cells

During development and aging, genome mutation leading to loss of heterozygosity (LOH) can uncover recessive phenotypes within tissue compartments. This phenomenon occurs in normal human tissues and is prevalent in pathological genetic conditions and cancers. While studies in yeast have defined DNA repair mechanisms that can promote LOH, the predominant pathways and environmental triggers in somatic tissues of multicellular organisms are not well understood. Here, we investigate mechanisms underlying LOH in intestinal stem cells in Drosophila. Infection with the pathogenic bacteria, Erwinia carotovora carotovora 15, but not Pseudomonas entomophila, increases LOH frequency. Using whole genome sequencing of somatic LOH events, we demonstrate that they arise primarily via mitotic recombination. Molecular features and genetic evidence argue against a break-induced replication mechanism and instead support cross-over via double Holliday junction-based repair. This study provides a mechanistic understanding of mitotic recombination, an important mediator of LOH, and its effects on stem cells in vivo.

The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication

Panax notoginseng is one of the most valuable medicinal species. However, its mitochondrial genome has not been reported yet. We aimed to determine the mitogenome sequence of P. notoginseng. We de novo assembled the mitogenome with Illumina short reads and Nanopore long reads. The mitochondrial genome of P. notoginseng has a multipartite structure consisting of interconversion between a "master circle" and numerous "subgenomic circles" through recombinations mediated by 64 pairs of repetitive sequences. Among the multipartite structure, seven subgenomic circles were best supported. Six of the seven subgenomic circles shared an 852 bp conserved fragment. The complete mitogenome of P. notoginseng was 662,479 bp long including 34 mitochondrial protein-coding genes (PCGs), three rRNA, and 19 tRNA genes. We identified 166 microsatellite repeats and 26 long-tandem repeats. Phylogenetic analysis resolved a tree that was mostly congruent with the phylogeny of Apiales species described in the APG IV system and the tree built with the chloroplast genome sequences. A total of 12 mitochondrial plastid DNA fragments were identified. Lastly, we predicted 591C-to-U RNA editing sites in the coding regions of mitochondrial PCGs. The mitochondrial genome will lay the foundation for understanding the evolution of Panax species.

The proper interplay between the expression of Spo11 splice isoforms and the structure of the pseudoautosomal region promotes XY chromosomes recombination

XY chromosome missegregation is relatively common in humans and can lead to sterility or the generation of aneuploid spermatozoa. A leading cause of XY missegregation in mammals is the lack of formation of double-strand breaks (DSBs) in the pseudoautosomal region (PAR), a defect that may occur in mice due to faulty expression of Spo11 splice isoforms. Using a knock-in (ki) mouse that expresses only the single Spo11β splice isoform, here we demonstrate that by varying the genetic background of mice, the length of chromatin loops extending from the PAR axis and the XY recombination proficiency varies. In spermatocytes of C57Spo11βki/- mice, in which loops are relatively short, recombination/synapsis between XY is fairly normal. In contrast, in cells of C57/129Spo11βki/- males where PAR loops are relatively long, formation of DSBs in the PAR (more frequently the Y-PAR) and XY synapsis fails at a high rate, and mice produce sperm with sex-chromosomal aneuploidy. However, if the entire set of Spo11 splicing isoforms is expressed by a wild type allele in the C57/129 background, XY recombination and synapsis is recovered. By generating a Spo11αki mouse model, we prove that concomitant expression of SPO11β and SPO11α isoforms, boosts DSB formation in the PAR. Based on these findings, we propose that SPO11 splice isoforms cooperate functionally in promoting recombination in the PAR, constraining XY asynapsis defects that may arise due to differences in the conformation of the PAR between mouse strains.

Recombination between heterologous human acrocentric chromosomes

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

Rhometa: Population recombination rate estimation from metagenomic read datasets

Prokaryotic evolution is influenced by the exchange of genetic information between species through a process referred to as recombination. The rate of recombination is a useful measure for the adaptive capacity of a prokaryotic population. We introduce Rhometa (https://github.com/sid-krish/Rhometa), a new software package to determine recombination rates from shotgun sequencing reads of metagenomes. It extends the composite likelihood approach for population recombination rate estimation and enables the analysis of modern short-read datasets. We evaluated Rhometa over a broad range of sequencing depths and complexities, using simulated and real experimental short-read data aligned to external reference genomes. Rhometa offers a comprehensive solution for determining population recombination rates from contemporary metagenomic read datasets. Rhometa extends the capabilities of conventional sequence-based composite likelihood population recombination rate estimators to include modern aligned metagenomic read datasets with diverse sequencing depths, thereby enabling the effective application of these techniques and their high accuracy rates to the field of metagenomics. Using simulated datasets, we show that our method performs well, with its accuracy improving with increasing numbers of genomes. Rhometa was validated on a real S. pneumoniae transformation experiment, where we show that it obtains plausible estimates of the rate of recombination. Finally, the program was also run on ocean surface water metagenomic datasets, through which we demonstrate that the program works on uncultured metagenomic datasets.

A Cre-LoxP-based approach for combinatorial chromosome rearrangements in human HAP1 cells

Alterations of human karyotype caused by chromosomal rearrangements are often associated with considerable phenotypic effects. Studying molecular mechanisms underlying these effects requires an efficient and scalable experimental model. Here, we propose a Cre-LoxP-based approach for the generation of combinatorial diversity of chromosomal rearrangements. We demonstrate that using the developed system, both intra- and inter-chromosomal rearrangements can be induced in the human haploid HAP1 cells, although the latter is significantly less effective. The obtained genetically modified HAP1 cell line can be used to dissect genomic effects associated with intra-chromosomal structural variations.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Rad51-mediated interhomolog recombination during budding yeast meiosis is promoted by the meiotic recombination checkpoint and the conserved Pif1 helicase

During meiosis, recombination between homologous chromosomes (homologs) generates crossovers that promote proper segregation at the first meiotic division. Recombination is initiated by Spo11-catalyzed DNA double strand breaks (DSBs). 5' end resection of the DSBs creates 3' single strand tails that two recombinases, Rad51 and Dmc1, bind to form presynaptic filaments that search for homology, mediate strand invasion and generate displacement loops (D-loops). D-loop processing then forms crossover and non-crossover recombinants. Meiotic recombination occurs in two temporally distinct phases. During Phase 1, Rad51 is inhibited and Dmc1 mediates the interhomolog recombination that promotes homolog synapsis. In Phase 2, Rad51 becomes active and functions with Rad54 to repair residual DSBs, making increasing use of sister chromatids. The transition from Phase 1 to Phase 2 is controlled by the meiotic recombination checkpoint through the meiosis-specific effector kinase Mek1. This work shows that constitutive activation of Rad51 in Phase 1 results in a subset of DSBs being repaired by a Rad51-mediated interhomolog recombination pathway that is distinct from that of Dmc1. Strand invasion intermediates generated by Rad51 require more time to be processed into recombinants, resulting in a meiotic recombination checkpoint delay in prophase I. Without the checkpoint, Rad51-generated intermediates are more likely to involve a sister chromatid, thereby increasing Meiosis I chromosome nondisjunction. This Rad51 interhomolog recombination pathway is specifically promoted by the conserved 5'-3' helicase PIF1 and its paralog, RRM3 and requires Pif1 helicase activity and its interaction with PCNA. This work demonstrates that (1) inhibition of Rad51 during Phase 1 is important to prevent competition with Dmc1 for DSB repair, (2) Rad51-mediated meiotic recombination intermediates are initially processed differently than those made by Dmc1, and (3) the meiotic recombination checkpoint provides time during prophase 1 for processing of Rad51-generated recombination intermediates.

Efficient assembly of long DNA fragments and multiple genes with improved nickase-based cloning and Cre/loxP recombination

Functional genomics, synthetic biology and metabolic engineering require efficient tools to deliver long DNA fragments or multiple gene constructs. Although numerous DNA assembly methods exist, most are complicated, time-consuming and expensive. Here, we developed a simple and flexible strategy, unique nucleotide sequence-guided nicking endonuclease (UNiE)-mediated DNA assembly (UNiEDA), for efficient cloning of long DNAs and multigene stacking. In this system, a set of unique 15-nt 3' single-strand overhangs were designed and produced by nicking endonucleases (nickases) in vectors and insert sequences. We introduced UNiEDA into our modified Cre/loxP recombination-mediated TransGene Stacking II (TGSII) system to generate an improved multigene stacking system we call TGSII-UNiE. Using TGSII-UNiE, we achieved efficient cloning of long DNA fragments of different sizes and assembly of multiple gene cassettes. Finally, we engineered and validated the biosynthesis of betanin in wild tobacco (Nicotiana benthamiana) leaves and transgenic rice (Oryza sativa) using multigene stacking constructs based on TGSII-UNiE. In conclusion, UNiEDA is an efficient, convenient and low-cost method for DNA cloning and multigene stacking, and the TGSII-UNiE system has important application prospects for plant functional genomics, genetic engineering and synthetic biology research.

Construction of a High-Density Recombination Bin-Based Genetic Map Facilitates High-Resolution Mapping of a Major QTL Underlying Anthocyanin Pigmentation in Eggplant

High-density genetic maps can significantly improve the resolution of QTL mapping. We constructed a high-density recombination bin-based genetic map of eggplant based on 200 F₂ plants from an interspecific cross (Solanum melongena × S. incanum) using the whole genome resequencing strategy. The map was 2022.8 cM long, covering near 99% of the eggplant genome. The map contained 3776 bins, with 3644 (96.5%) being effective (position non-redundant) ones, giving a nominal average distance of 0.54 cM and an effective average distance of 0.56 cM between adjacent bins, respectively. Using this map and 172 F_2:3 lines, a major QTL with pleiotropic effects on two anthocyanin pigmentation-related traits, leaf vein color (LVC) and fruit pericarp color (FPC), was steadily detected in a bin interval of 2.28 cM (or 1.68 Mb) on chromosome E10 in two cropping seasons, explaining ~65% and 55% of the phenotypic variation in LVC and FPC, respectively. Genome-wide association analysis in this population validated the QTL and demonstrated the correctness of mapping two bins of chromosome E02 onto E10. Bioinformatics analysis suggested that a WDR protein gene inside the bin interval with reliable effective variation between the two parents could be a possible candidate gene of the QTL.

Pandemic-scale phylogenomics reveals the SARS-CoV-2 recombination landscape

Accurate and timely detection of recombinant lineages is crucial for interpreting genetic variation, reconstructing epidemic spread, identifying selection and variants of interest, and accurately performing phylogenetic analyses1-4. During the SARS-CoV-2 pandemic, genomic data generation has exceeded the capacities of existing analysis platforms, thereby crippling real-time analysis of viral evolution5. Here, we use a new phylogenomic method to search a nearly comprehensive SARS-CoV-2 phylogeny for recombinant lineages. In a 1.6 million sample tree from May 2021, we identify 589 recombination events, which indicate that around 2.7% of sequenced SARS-CoV-2 genomes have detectable recombinant ancestry. Recombination breakpoints are inferred to occur disproportionately in the 3' portion of the genome that contains the spike protein. Our results highlight the need for timely analyses of recombination for pinpointing the emergence of recombinant lineages with the potential to increase transmissibility or virulence of the virus. We anticipate that this approach will empower comprehensive real-time tracking of viral recombination during the SARS-CoV-2 pandemic and beyond.

Neural networks for self-adjusting mutation rate estimation when the recombination rate is unknown

Estimating the mutation rate, or equivalently effective population size, is a common task in population genetics. If recombination is low or high, optimal linear estimation methods are known and well understood. For intermediate recombination rates, the calculation of optimal estimators is more challenging. As an alternative to model-based estimation, neural networks and other machine learning tools could help to develop good estimators in these involved scenarios. However, if no benchmark is available it is difficult to assess how well suited these tools are for different applications in population genetics.

Here we investigate feedforward neural networks for the estimation of the mutation rate based on the site frequency spectrum and compare their performance with model-based estimators. For this we use the model-based estimators introduced by Fu, Futschik et al., and Watterson that minimize the variance or mean squared error for no and free recombination. We find that neural networks reproduce these estimators if provided with the appropriate features and training sets. Remarkably, using the model-based estimators to adjust the weights of the training data, only one hidden layer is necessary to obtain a single estimator that performs almost as well as model-based estimators for low and high recombination rates, and at the same time provides a superior estimation method for intermediate recombination rates. We apply the method to simulated data based on the human chromosome 2 recombination map, highlighting its robustness in a realistic setting where local recombination rates vary and/or are unknown.

single-layer feedforward neural networks learn the established model-based linear estimators for high and low recombination rates

neural networks learn good estimators for intermediate recombination rates where computation of model-based optimal estimators is hardly possible

a single neural network estimator can automatically adapt to a variable recombination rate and performs close to optimal

this is advantageous when recombination rates vary along the chromosome according to a recombination map

using the known estimators as a benchmark to adapt the training error function improves the estimates of the neural networks

Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes

Many organisms have sex chromosomes with large nonrecombining regions that have expanded stepwise, generating "evolutionary strata" of differentiation. The reasons for this remain poorly understood, but the principal hypotheses proposed to date are based on antagonistic selection due to differences between sexes. However, it has proved difficult to obtain empirical evidence of a role for sexually antagonistic selection in extending recombination suppression, and antagonistic selection has been shown to be unlikely to account for the evolutionary strata observed on fungal mating-type chromosomes. We show here, by mathematical modeling and stochastic simulation, that recombination suppression on sex chromosomes and around supergenes can expand under a wide range of parameter values simply because it shelters recessive deleterious mutations, which are ubiquitous in genomes. Permanently heterozygous alleles, such as the male-determining allele in XY systems, protect linked chromosomal inversions against the expression of their recessive mutation load, leading to the successive accumulation of inversions around these alleles without antagonistic selection. Similar results were obtained with models assuming recombination-suppressing mechanisms other than chromosomal inversions and for supergenes other than sex chromosomes, including those without XY-like asymmetry, such as fungal mating-type chromosomes. However, inversions capturing a permanently heterozygous allele were found to be less likely to spread when the mutation load segregating in populations was lower (e.g., under large effective population sizes or low mutation rates). This may explain why sex chromosomes remain homomorphic in some organisms but are highly divergent in others. Here, we model a simple and testable hypothesis explaining the stepwise extensions of recombination suppression on sex chromosomes, mating-type chromosomes, and supergenes in general.

Bayesian inference of ancestral recombination graphs

We present a novel algorithm, implemented in the software ARGinfer, for probabilistic inference of the Ancestral Recombination Graph under the Coalescent with Recombination. Our Markov Chain Monte Carlo algorithm takes advantage of the Succinct Tree Sequence data structure that has allowed great advances in simulation and point estimation, but not yet probabilistic inference. Unlike previous methods, which employ the Sequentially Markov Coalescent approximation, ARGinfer uses the Coalescent with Recombination, allowing more accurate inference of key evolutionary parameters. We show using simulations that ARGinfer can accurately estimate many properties of the evolutionary history of the sample, including the topology and branch lengths of the genealogical tree at each sequence site, and the times and locations of mutation and recombination events. ARGinfer approximates posterior probability distributions for these and other quantities, providing interpretable assessments of uncertainty that we show to be well calibrated. ARGinfer is currently limited to tens of DNA sequences of several hundreds of kilobases, but has scope for further computational improvements to increase its applicability.

RecPD: A Recombination-aware measure of phylogenetic diversity

A critical step in studying biological features (e.g., genetic variants, gene families, metabolic capabilities, or taxa) is assessing their diversity and distribution among a sample of individuals. Accurate assessments of these patterns are essential for linking features to traits or outcomes of interest and understanding their functional impact. Consequently, it is of crucial importance that the measures employed for quantifying feature diversity can perform robustly under any evolutionary scenario. However, the standard measures used for quantifying and comparing the distribution of features, such as prevalence, phylogenetic diversity, and related approaches, either do not take into consideration evolutionary history, or assume strictly vertical patterns of inheritance. Consequently, these approaches cannot accurately assess diversity for features that have undergone recombination or horizontal transfer. To address this issue, we have devised RecPD, a novel recombination-aware phylogenetic-diversity statistic for measuring the distribution and diversity of features under all evolutionary scenarios. RecPD utilizes ancestral-state reconstruction to map the presence / absence of features onto ancestral nodes in a species tree, and then identifies potential recombination events in the evolutionary history of the feature. We also derive several related measures from RecPD that can be used to assess and quantify evolutionary dynamics and correlation of feature evolutionary histories. We used simulation studies to show that RecPD reliably reconstructs feature evolutionary histories under diverse recombination and loss scenarios. We then applied RecPD in two diverse real-world scenarios including a preliminary study type III effector protein families secreted by the plant pathogenic bacterium Pseudomonas syringae and growth phenotypes of the Pseudomonas genus and demonstrate that prevalence is an inadequate measure that obscures the potential impact of recombination. We believe RecPD will have broad utility for revealing and quantifying complex evolutionary processes for features at any biological level.

Phylogenetic diversity is an important concept utilized in evolutionary ecology which has extensive applications in population genetics to help us understand how evolutionary processes have distributed genetic variation among individuals of a species, and how this impacts phenotypic diversification over time. However, existing approaches for studying phylogenetic diversity largely assume that the genetic features follow vertical inheritance, which is frequently violated in the case of microbial genomes due to horizontal transfer. To address this shortcoming, we present RecPD, a recombination-aware phylogenetic diversity measure, which incorporates ancestral state reconstruction to quantify the phylogenetic diversity of genetic features mapped onto a species phylogeny. Through simulation experiments we show that RecPD robustly reconstructs the evolutionary histories of features evolving under various scenarios of recombination and loss. When applied to a real-world example of type III secreted effector protein families from the plant pathogenic bacterium Pseudomonas syringae, RecPD reveals that horizontal transfer has played an important role in shaping the phylogenetic distributions of a substantial proportion of families across the P. syringae species complex. Furthermore, we demonstrate that the traditional measures of feature prevalence are unsuitable as a measure for comparing feature diversity. We also provide a R package implementation of RecPD for public use: https://github.com/cedatorma/recpd.

Mechanism of Bloom syndrome complex assembly required for double Holliday junction dissolution and genome stability

The RecQ-like helicase BLM cooperates with topoisomerase IIIα, RMI1, and RMI2 in a heterotetrameric complex (the "Bloom syndrome complex") for dissolution of double Holliday junctions, key intermediates in homologous recombination. Mutations in any component of the Bloom syndrome complex can cause genome instability and a highly cancer-prone disorder called Bloom syndrome. Some heterozygous carriers are also predisposed to breast cancer. To understand how the activities of BLM helicase and topoisomerase IIIα are coupled, we purified the active four-subunit complex. Chemical cross-linking and mass spectrometry revealed a unique architecture that links the helicase and topoisomerase domains. Using biochemical experiments, we demonstrated dimerization mediated by the N terminus of BLM with a 2:2:2:2 stoichiometry within the Bloom syndrome complex. We identified mutations that independently abrogate dimerization or association of BLM with RMI1, and we show that both are dysfunctional for dissolution using in vitro assays and cause genome instability and synthetic lethal interactions with GEN1/MUS81 in cells. Truncated BLM can also inhibit the activity of full-length BLM in mixed dimers, suggesting a putative mechanism of dominant-negative action in carriers of BLM truncation alleles. Our results identify critical molecular determinants of Bloom syndrome complex assembly required for double Holliday junction dissolution and maintenance of genome stability.

Resolving the structure of phage-bacteria interactions in the context of natural diversity

Microbial communities are shaped by viral predators. Yet, resolving which viruses (phages) and bacteria are interacting is a major challenge in the context of natural levels of microbial diversity. Thus, fundamental features of how phage-bacteria interactions are structured and evolve in the wild remain poorly resolved. Here we use large-scale isolation of environmental marine Vibrio bacteria and their phages to obtain estimates of strain-level phage predator loads, and use all-by-all host range assays to discover how phage and host genomic diversity shape interactions. We show that lytic interactions in environmental interaction networks (as observed in agar overlay) are sparse-with phage predator loads being low for most bacterial strains, and phages being host-strain-specific. Paradoxically, we also find that although overlap in killing is generally rare between tailed phages, recombination is common. Together, these results suggest that recombination during cryptic co-infections is an important mode of phage evolution in microbial communities. In the development of phages for bioengineering and therapeutics it is important to consider that nucleic acids of introduced phages may spread into local phage populations through recombination, and that the likelihood of transfer is not predictable based on lytic host range.

Assembly of the Complete Mitochondrial Genome of Chinese Plum (Prunus salicina): Characterization of Genome Recombination and RNA Editing Sites

Despite the significant progress that has been made in the genome sequencing of Prunus, this area of research has been lacking a systematic description of the mitochondrial genome of this genus for a long time. In this study, we assembled the mitochondrial genome of the Chinese plum (Prunus salicina) using Illumina and Oxford Nanopore sequencing data. The mitochondrial genome size of P. salicina was found to be 508,035 base pair (bp), which is the largest reported in the Rosaceae family to date, and P. salicina was shown to be 63,453 bp longer than sweet cherry (P. avium). The P. salicina mitochondrial genome contained 37 protein-coding genes (PCGs), 3 ribosomal RNA (rRNA) genes, and 16 transfer RNA (tRNA) genes. Two plastid-derived tRNA were identified. We also found two short repeats that captured the nad3 and nad6 genes and resulted in two copies. In addition, nine pairs of repeat sequences were identified as being involved in the mediation of genome recombination. This is crucial for the formation of subgenomic configurations. To characterize RNA editing sites, transcriptome data were used, and we identified 480 RNA editing sites in protein-coding sequences. Among them, the initiation codon of the nad1 gene confirmed that an RNA editing event occurred, and the genomic encoded ACG was edited as AUG in the transcript. Combined with previous reports on the chloroplast genome, our data complemented our understanding of the last part of the organelle genome of plum, which will facilitate our understanding of the evolution of organelle genomes.

De novo deletions and duplications at recombination hotspots in mouse germlines

Numerous DNA double-strand breaks (DSBs) arise during meiosis to initiate homologous recombination. These DSBs are usually repaired faithfully, but here, we uncover a distinct type of mutational event in which deletions form via joining of ends from two closely spaced DSBs (double cuts) within a single hotspot or at adjacent hotspots on the same or different chromatids. Deletions occur in normal meiosis but are much more frequent when DSB formation is dysregulated in the absence of the ATM kinase. Events between chromosome homologs point to multi-chromatid damage and aborted gap repair. Some deletions contain DNA from other hotspots, indicating that double cutting at distant sites creates substrates for insertional mutagenesis. End joining at double cuts can also yield tandem duplications or extrachromosomal circles. Our findings highlight the importance of DSB regulation and reveal a previously hidden potential for meiotic mutagenesis that is likely to affect human health and genome evolution.

MicroRNA-directed pathway discovery elucidates an miR-221/222-mediated regulatory circuit in class switch recombination

MicroRNAs (miRNAs, miRs) regulate cell fate decisions by post-transcriptionally tuning networks of mRNA targets. We used miRNA-directed pathway discovery to reveal a regulatory circuit that influences Ig class switch recombination (CSR). We developed a system to deplete mature, activated B cells of miRNAs, and performed a rescue screen that identified the miR-221/222 family as a positive regulator of CSR. Endogenous miR-221/222 regulated B cell CSR to IgE and IgG1 in vitro, and miR-221/222-deficient mice exhibited defective IgE production in allergic airway challenge and polyclonal B cell activation models in vivo. We combined comparative Ago2-HITS-CLIP and gene expression analyses to identify mRNAs bound and regulated by miR-221/222 in primary B cells. Interrogation of these putative direct targets uncovered functionally relevant downstream genes. Genetic depletion or pharmacological inhibition of Foxp1 and Arid1a confirmed their roles as key modulators of CSR to IgE and IgG1.

Fundamental evolution of all Orthocoronavirinae including three deadly lineages descendent from Chiroptera-hosted coronaviruses: SARS-CoV, MERS-CoV and SARS-CoV-2

The severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in humans in 2002. Despite reports showing Chiroptera as the original animal reservoir of SARS-CoV, many argue that Carnivora-hosted viruses are the most likely origin. The emergence of the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 also involves Chiroptera-hosted lineages. However, factors such as the lack of comprehensive phylogenies hamper our understanding of host shifts once MERS-CoV emerged in humans and Artiodactyla. Since 2019, the origin of SARS-CoV-2, causative agent of coronavirus disease 2019 (COVID-19), added to this episodic history of zoonotic transmission events. Here we introduce a phylogenetic analysis of 2006 unique and complete genomes of different lineages of Orthocoronavirinae. We used gene annotations to align orthologous sequences for total evidence analysis under the parsimony optimality criterion. Deltacoronavirus and Gammacoronavirus were set as outgroups to understand spillovers of Alphacoronavirus and Betacoronavirus among ten orders of animals. We corroborated that Chiroptera-hosted viruses are the sister group of SARS-CoV, SARS-CoV-2 and MERS-related viruses. Other zoonotic events were qualified and quantified to provide a comprehensive picture of the risk of coronavirus emergence among humans. Finally, we used a 250 SARS-CoV-2 genomes dataset to elucidate the phylogenetic relationship between SARS-CoV-2 and Chiroptera-hosted coronaviruses.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

Directed evolution of a family of AAV capsid variants enabling potent muscle-directed gene delivery across species

Replacing or editing disease-causing mutations holds great promise for treating many human diseases. Yet, delivering therapeutic genetic modifiers to specific cells in vivo has been challenging, particularly in large, anatomically distributed tissues such as skeletal muscle. Here, we establish an in vivo strategy to evolve and stringently select capsid variants of adeno-associated viruses (AAVs) that enable potent delivery to desired tissues. Using this method, we identify a class of RGD motif-containing capsids that transduces muscle with superior efficiency and selectivity after intravenous injection in mice and non-human primates. We demonstrate substantially enhanced potency and therapeutic efficacy of these engineered vectors compared to naturally occurring AAV capsids in two mouse models of genetic muscle disease. The top capsid variants from our selection approach show conserved potency for delivery across a variety of inbred mouse strains, and in cynomolgus macaques and human primary myotubes, with transduction dependent on target cell expressed integrin heterodimers.

Pericentromeric chromatin reorganisation follows the initiation of recombination and coincides with early events of synapsis in cereals

The reciprocal exchange of genetic information between homologous chromosomes during meiotic recombination is essential to secure balanced chromosome segregation and to promote genetic diversity. The chromosomal position and frequency of reciprocal genetic exchange shapes the efficiency of breeding programmes and influences crop improvement under a changing climate. In large genome cereals, such as wheat and barley, crossovers are consistently restricted to subtelomeric chromosomal regions, thus preventing favourable allele combinations being formed within a considerable proportion of the genome, including interstitial and pericentromeric chromatin. Understanding the key elements driving crossover designation is therefore essential to broaden the regions available for crossovers. Here, we followed early meiotic chromatin dynamism in cereals through the visualisation of a homologous barley chromosome arm pair stably transferred into the wheat genetic background. By capturing the dynamics of a single chromosome arm at the same time as detecting the undergoing events of meiotic recombination and synapsis, we showed that subtelomeric chromatin of homologues synchronously transitions to an open chromatin structure during recombination initiation. By contrast, pericentromeric and interstitial regions preserved their closed chromatin organisation and become unpackaged only later, concomitant with initiation of recombinatorial repair and the initial assembly of the synaptonemal complex. Our results raise the possibility that the closed pericentromeric chromatin structure in cereals may influence the fate decision during recombination initiation, as well as the spatial development of synapsis, and may also explain the suppression of crossover events in the proximity of the centromeres.

Identification of stable heat tolerance QTLs using inter-specific recombinant inbred line population derived from GPF 2 and ILWC 292

Heat stress during reproductive stages has been leading to significant yield losses in chickpea (Cicer arietinum L.). With an aim of identifying the genomic regions or QTLs responsible for heat tolerance, 187 F8 recombinant inbred lines (RILs) derived from the cross GPF 2 (heat tolerant) × ILWC 292 (heat sensitive) were evaluated under late-sown irrigated (January-May) and timely-sown irrigated environments (November-April) at Ludhiana and Faridkot in Punjab, India for 13 heat tolerance related traits. The pooled ANOVA for both locations for the traits namely days to germination (DG), days to flowering initiation (DFI), days to 50% flowering (DFF), days to 100% flowering (DHF), plant height (PH), pods per plant (NPP), biomass (BIO), grain yield (YLD), 100-seed weight (HSW), harvest index (HI), membrane permeability index (MPI), relative leaf water content (RLWC) and pollen viability (PV)) showed a highly significant difference in RILs. The phenotyping data coupled with the genetic map comprising of 1365 ddRAD-Seq based SNP markers were used for identifying the QTLs for heat tolerance. Composite interval mapping provided a total of 28 and 23 QTLs, respectively at Ludhiana and Faridkot locations. Of these, 13 consensus QTLs for DG, DFI, DFF, DHF, PH, YLD, and MPI have been identified at both locations. Four QTL clusters containing QTLs for multiple traits were identified on the same genomic region at both locations. Stable QTLs for days to flowering can be one of the major factors for providing heat tolerance as early flowering has an advantage of more seed setting due to a comparatively longer reproductive period. Identified QTLs can be used in genomics-assisted breeding to develop heat stress-tolerant high yielding chickpea cultivars.

Germinal center-dependent and -independent memory B cells produced throughout the immune response

Memory B cells comprise a heterogenous group of cells that differ in origin and phenotype. During the early phases of the immune response, activated B cells can differentiate into IgM-expressing memory cells, short-lived plasma cells, or seed germinal centers (GCs). The memory compartment is subsequently enriched by B cells that have been through several rounds of division and selection in the GC. Here, we report on the use of an unbiased lineage-tracking approach to explore the origins and properties of memory B cell subsets in mice with an intact immune system. We find that activated B cells continue to differentiate into memory B cells throughout the immune response. When defined on the basis of their origins, the memory B cells originating from activated B cells or GCs differ in isotype and overall gene expression, somatic hypermutation, and their affinity for antigen.

Effective RNAi-Mediated Silencing of the Mismatch Repair MSH2 Gene Induces Sterility of Tomato Plants but Not an Increase in Meiotic Recombination

In plant breeding, the ability to manipulate meiotic recombination aids in the efficient construction of new allelic compositions of chromosomes and facilitates gene transfer from wild relatives of crop plants. The DNA mismatch repair system antagonizes meiotic recombination. In this research, a trial was conducted to evaluate transgenic tomato plants carrying an RNA interference (RNAi) construct designed to inhibit the expression of the mismatch repair MSH2 gene. To drive the RNAi construct, we used either a pro-SmAMP2 promoter from Stellaria media ANTIMICROBIAL PEPTIDE2 or a Cauliflower mosaic virus 35S promoter (CaMV35S). The results of real-time PCR showed that, with a 16 h light/8 h dark photoperiod, MSH2-RNAi tomato transgenic plants exhibited MSH2 gene transcript contents ranging from 0% to 3% in the leaves, relative to untransformed controls. However, with this lighting mode, the MSH2-RNAi transgenic plants grew slowly, flowered poorly, and did not form seed sets. During cultivation with a 12 h light/12 h dark photoperiod, MSH2-RNAi transgenic plants exhibited MSH2 gene transcript contents ranging from 3% to 42%, relative to untransformed controls. Under these conditions, F1 hybrid seed sets formed for most of the MSH2-RNAi transgenic plants with the RNAi construct driven by the CaMV35S promoter, and for one transformant with the RNAi construct driven by the pro-SmAMP2 promoter. Under conditions of a 12 h light/12 h dark photoperiod, most of the F1 transgenic hybrids showed MSH2 gene transcript contents ranging from 3% to 34% and formed F2 offspring sets, which made it possible to assess the meiotic recombination frequency. We showed that the effective inhibition of MSH2 in MSH2-RNAi tomato transgenic plants is not associated with an increase in meiotic recombination compared to the control, but it stimulates the sterility of plants. It was established that the expression of the MSH2 gene in tomato plants is about 50 times higher with a 12 h light/12 h dark than with a 16 h light/8 h dark photoperiod. It is discussed that, in Solanum lycopersicum tomato plants, which are not sensitive to the day length for flowering, changing the lighting time may be a means of controlling the meiotic recombination frequency within certain limits.

A novel mouse line with epididymal initial segment-specific expression of Cre recombinase driven by the endogenous Lcn9 promoter

Spermatozoa released from testes undergo a maturation process and acquire the capacity to fertilize ova through epididymal transit. The epididymis is divided into four regions, each with unique morphology, gene profile, luminal microenvironment and distinct function. To study the functions of relevant genes in the epididymal initial segment (IS), a novel IS-specific mouse model, Lcn9-Cre knock-in (KI) mouse line was generated via CRISPR/Cas9 technology. The TAG stop codon was replaced by a 2A-NLS-Cre cassette, resulting in the co-expression of Lcn9 and Cre recombinase. IS-specific Cre expression was first observed from postnatal day 17. Using the Rosa26^tdTomato reporter mice, the Cre-mediated DNA recombination was detected exclusively in principal cells. The epididymal IS-specific Cre activity in vivo was further confirmed using Lcn9-Cre mice crossed with a mouse strain carrying Tsc1 floxed alleles (Tsc1^flox/+). Cre expression did not affect either normal development or male fecundity. Different from any epididymis-specific Cre mice reported previously, the novel Lcn9-Cre mouse line can be used to introduce entire IS-specific conditional gene editing for gene functional study.

Assembly of the complete mitochondrial genome of an endemic plant, Scutellaria tsinyunensis, revealed the existence of two conformations generated by a repeat-mediated recombination

We assembled the complete mitochondrial genome of Scutellaria tsinyunensis in this study. Repeat-mediated recombination resulted in the formation of two conformations of the mitochondrial genome in S. tsinyunensis. Scutellaria tsinyunensis belongs to the family Lamiaceae, distributed only in the Jinyun Mountain, Chongqing, China. As a valuable endemic and small population species, it is regarded as a natural resource potentially with significant economic and ecological importance. In this study, we assembled a complete and gap-free mitochondrial genome of S. tsinyunensis. This genome had a length of 354,073 bp and the base composition of the genome was A (27.44%), T (27.30%), C (22.58%), and G (22.68%). This genome encodes 59 genes, including 32 protein-coding genes, 24 tRNA genes, and 3 rRNA genes. The Sanger sequencing and Oxford Nanopore sequencing confirmed a pair of direct repeats had mediated genome recombination, resulting in the formation of two conformations. The gene conversation between plastome and mitochondrial genome was also observed in S. tsinyunensis by detecting gene migration, including six tRNA genes (namely, trnW-CCA, trnI-CAU, trnH-UUU, trnD-GUC, trnN-GUU, and trnM-CAU), five protein-coding gene fragments, and the fragments from 2 rRNA genes. Moreover, the dN/dS analysis revealed the atp9 gene had undergone strong negative selection, and four genes (atp4, mttB, ccmFc, and ccmB) probably had undergone positive selection during evolution in Lamiales. This work reported the first mitochondrial genome of S. tsinyunensis, which could be used as a reference genome for the important medicinal plants of the genus Scutellaria, and also provide much-desired information for molecular breeding.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Mutational Analysis of Redβ Single Strand Annealing Protein: Roles of the 14 Lysine Residues in DNA Binding and Recombination In Vivo

Redβ is a 261 amino acid protein from bacteriophage λ that promotes a single-strand annealing (SSA) reaction for repair of double-stranded DNA (dsDNA) breaks. While there is currently no high-resolution structure available for Redβ, models of its DNA binding domain (residues 1-188) have been proposed based on homology with human Rad52, and a crystal structure of its C-terminal domain (CTD, residues 193-261), which binds to λ exonuclease and E. coli single-stranded DNA binding protein (SSB), has been determined. To evaluate these models, the 14 lysine residues of Redβ were mutated to alanine, and the variants tested for recombination in vivo and DNA binding and annealing in vitro. Most of the lysines within the DNA binding domain, including K36, K61, K111, K132, K148, K154, and K172, were found to be critical for DNA binding in vitro and recombination in vivo. By contrast, none of the lysines within the CTD, including K214, K245, K251, K253, and K258 were required for DNA binding in vitro, but two, K214 and K253, were critical for recombination in vivo, likely due to their involvement in binding to SSB. K61 was identified as a residue that is critical for DNA annealing, but not for initial ssDNA binding, suggesting a role in binding to the second strand of DNA incorporated into the complex. The K148A variant, which has previously been shown to be defective in oligomer formation, had the lowest affinity for ssDNA, and was the only variant that was completely non-cooperative, suggesting that ssDNA binding is coupled to oligomerization.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

ATM controls meiotic DNA double-strand break formation and recombination and affects synaptonemal complex organization in plants

Meiosis is a specialized cell division that gives rise to genetically distinct gametic cells. Meiosis relies on the tightly controlled formation of DNA double-strand breaks (DSBs) and their repair via homologous recombination for correct chromosome segregation. Like all forms of DNA damage, meiotic DSBs are potentially harmful and their formation activates an elaborate response to inhibit excessive DNA break formation and ensure successful repair. Previous studies established the protein kinase ATM as a DSB sensor and meiotic regulator in several organisms. Here we show that Arabidopsis ATM acts at multiple steps during DSB formation and processing, as well as crossover (CO) formation and synaptonemal complex (SC) organization, all vital for the successful completion of meiosis. We developed a single-molecule approach to quantify meiotic breaks and determined that ATM is essential to limit the number of meiotic DSBs. Local and genome-wide recombination screens showed that ATM restricts the number of interference-insensitive COs, while super-resolution STED nanoscopy of meiotic chromosomes revealed that the kinase affects chromatin loop size and SC length and width. Our study extends our understanding of how ATM functions during plant meiosis and establishes it as an integral factor of the meiotic program.

Genetic mapping and QTL analysis for peanut smut resistance

BACKGROUND

Peanut smut is a disease caused by the fungus Thecaphora frezii Carranza & Lindquist to which most commercial cultivars in South America are highly susceptible. It is responsible for severely decreased yield and no effective chemical treatment is available to date. However, smut resistance has been identified in wild Arachis species and further transferred to peanut elite cultivars. To identify the genome regions conferring smut resistance within a tetraploid genetic background, this study evaluated a RIL population {susceptible Arachis hypogaea subsp. hypogaea (JS17304-7-B) × resistant synthetic amphidiploid (JS1806) [A. correntina (K 11905) × A. cardenasii (KSSc 36015)] × A. batizocoi (K 9484)^4×} segregating for the trait.

RESULTS

A SNP based genetic map arranged into 21 linkage groups belonging to the 20 peanut chromosomes was constructed with 1819 markers, spanning a genetic distance of 2531.81 cM. Two consistent quantitative trait loci (QTLs) were identified qSmIA08 and qSmIA02/B02, located on chromosome A08 and A02/B02, respectively. The QTL qSmIA08 at 15.20 cM/5.03 Mbp explained 17.53% of the phenotypic variance, while qSmIA02/B02 at 4.0 cM/3.56 Mbp explained 9.06% of the phenotypic variance. The combined genotypic effects of both QTLs reduced smut incidence by 57% and were stable over the 3 years of evaluation. The genome regions containing the QTLs are rich in genes encoding proteins involved in plant defense, providing new insights into the genetic architecture of peanut smut resistance.

CONCLUSIONS

A major QTL and a minor QTL identified in this study provide new insights into the genetic architecture of peanut smut resistance that may aid in breeding new varieties resistant to peanut smut.

High-throughput measuring of meiotic recombination rates in barley pollen nuclei using Crystal Digital PCRTM

Breeding exploits novel allelic combinations assured by meiotic recombination. Barley (Hordeum vulgare) single pollen nucleus genotyping enables measurement of meiotic recombination rates in gametes before fertilization without the need for segregating populations. However, so far, established methods rely on whole-genome amplification of every single pollen nucleus due to their limited DNA content, thus restricting the number of analyzed samples. In this study, high-throughput measurements of meiotic recombination rates in barley pollen nuclei without whole-genome amplification were performed through a Crystal Digital PCR^TM -based genotyping assay. Meiotic recombination rates within two centromeric and two distal chromosomal intervals were measured in hybrid plants by genotyping a total of >42 000 individual pollen nuclei (up to 4900 nuclei analyzed per plant). Determined recombination frequencies in pollen nuclei were similar to frequencies in segregating populations. We improved the efficiency of the genotyping by pretreating the pollen nuclei with a thermostable restriction enzyme. Additional opportunities for a higher sample throughput and a further increase of the genotyping efficiency are presented and discussed. Taken together, single barley pollen nucleus genotyping based on Crystal Digital PCR^TM enables reliable, rapid and high-throughput meiotic recombination measurements within defined chromosomal intervals of intraspecific hybrid plants. The successful encapsulation of nuclei from a range of species with different nuclear and genome sizes suggests that the proposed method is broadly applicable to genotyping single nuclei.

DNA methylation occurring in Cre-expressing cells inhibits loxP recombination and silences loxP-sandwiched genes

The low DNA recombination efficiency of site-specific recombinase systems in plants limits their application; however, the underlying mechanism is unknown. We evaluate the gene deletion performance of four recombinase systems (Cre/loxP, Flp/FRT, KD/KDRT and B3/B3RT) in tobacco where the recombinases are under the control of germline-specific promoters. We find that the expression of these recombinases results mostly in gene silencing rather than gene deletion. Using the Cre/loxP system as a model, we reveal that the region flanked by loxP sites (floxed) is hypermethylated, which prevents floxed genes from deletion while silencing the expression of the genes. We further show CG methylation alone in the recombinase binding element of the loxP site is unable to impede gene deletion; instead, CHH methylation in the crossover region is required to inhibit loxP recombination. Our study illustrates the important role of recombinase-induced DNA methylation in the inhibition of site-specific DNA recombination and uncovers the mechanism underlying recombinase-associated gene silence in plants.

A genome-wide library of MADM mice for single-cell genetic mosaic analysis

Mosaic analysis with double markers (MADM) offers one approach to visualize and concomitantly manipulate genetically defined cells in mice with single-cell resolution. MADM applications include the analysis of lineage, single-cell morphology and physiology, genomic imprinting phenotypes, and dissection of cell-autonomous gene functions in vivo in health and disease. Yet, MADM can only be applied to <25% of all mouse genes on select chromosomes to date. To overcome this limitation, we generate transgenic mice with knocked-in MADM cassettes near the centromeres of all 19 autosomes and validate their use across organs. With this resource, >96% of the entire mouse genome can now be subjected to single-cell genetic mosaic analysis. Beyond a proof of principle, we apply our MADM library to systematically trace sister chromatid segregation in distinct mitotic cell lineages. We find striking chromosome-specific biases in segregation patterns, reflecting a putative mechanism for the asymmetric segregation of genetic determinants in somatic stem cell division.

Herpes Simplex Virus 1 Coinfection Modifies Adeno-associated Virus Genome End Recombination

Wild-type adeno-associated virus (AAV) can only replicate in the presence of helper factors, which can be provided by coinfecting helper viruses such as adenoviruses and herpesviruses. The AAV genome consists of a linear, single-stranded DNA (ssDNA), which is converted into different molecular structures within the host cell. Using high-throughput sequencing, we found that herpes simplex virus 1 (HSV-1) coinfection leads to a shift in the type of AAV genome end recombination. In particular, open-end inverted terminal repeat (ITR) recombination was enhanced, whereas open-closed ITR recombination was reduced in the presence of HSV-1. We demonstrate that the HSV-1 protein ICP8 plays an essential role in HSV-1-mediated interference with AAV genome end recombination, indicating that the previously described ICP8-driven mechanism of HSV-1 genome recombination may be underlying the observed changes. We also provide evidence that additional factors, such as products of true late genes, are involved. Although HSV-1 coinfection significantly changed the type of AAV genome end recombination, no significant change in the amount of circular AAV genomes was identified. IMPORTANCE Adeno-associated virus (AAV)-mediated gene therapy represents one of the most promising approaches for the treatment of genetic diseases. Currently, various GMP-compatible production methods can be applied to manufacture clinical-grade vector, including methods that employ helper factors derived from herpes simplex virus 1 (HSV-1). Yet, to date, we do not fully understand how HSV-1 interacts with AAV. We observed that HSV-1 modulates AAV genome ends similarly to the genome recombination events observed during HSV-1 replication and postulate that further improvements of the HSV-1 production platform may enhance packaging of the recombinant AAV particles.

Tousled-Like Kinases Suppress Innate Immune Signaling Triggered by Alternative Lengthening of Telomeres

The Tousled-like kinases 1 and 2 (TLK1/2) control histone deposition through the ASF1 histone chaperone and influence cell cycle progression and genome maintenance, yet the mechanisms underlying TLK-mediated genome stability remain uncertain. Here, we show that TLK loss results in severe chromatin decompaction and altered genome accessibility, particularly affecting heterochromatic regions. Failure to maintain heterochromatin increases spurious transcription of repetitive elements and induces features of alternative lengthening of telomeres (ALT). TLK depletion culminates in a cGAS-STING-TBK1-mediated innate immune response that is independent of replication-stress signaling and attenuated by the depletion of factors required to produce extra-telomeric DNA. Analysis of human cancers reveals that chromosomal instability correlates with high TLK2 and low STING levels in many cohorts. Based on these findings, we propose that high TLK levels contribute to immune evasion in chromosomally unstable and ALT+ cancers.

•

TLK-deficient cells have increased accessibility at heterochromatin regions

•

TLK1/2 suppress spurious transcription and telomere hyper-recombination

•

Extra-telomeric DNA generated upon TLK loss promotes innate immune signaling

•

cGAS-STING-TBK1 signaling in TLK-deficient cells is independent of replication stress

Inhibition of glioma by adenovirus KGHV500 encoding anti-p21Ras scFv and carried by cytokine-induced killer cells

Ras gene mutation or overexpression can lead to tumorigenesis in multiple kinds of cancer, including glioma. However, no drugs targeting Ras or its expression products have been approved for clinical application thus far. Adenoviral gene therapy is a promising method for the treatment of glioma. In this study, the human glioma cell line U251 was co-cultured with recombinant adenovirus KGHV500, and the anti-tumor effects of KGHV500 were determined by MTT, scratch test, Transwell invasion, and apoptosis assays. Then, KGHV500 was delivered via the intravenous injection of CIK cells into glioma xenografts. Tumor volume, ki67 proliferation index, apoptosis levels, and anti-p21Ras scFv expression were tested to evaluate targeting ability, anti-tumor efficacy, and safety. We found that the KGHV500 exhibited anti-tumor activity in U251 cells and increased the intracellular expression of anti-p21Ras scFv compared with that in the control groups. CIK cells delivered KGHV500 to U251 glioma cell xenografts and enhanced anti-tumor activity against glioma xenografts compared to that produced by the control treatment. In conclusion, targeting Ras is a useful therapeutic strategy for gliomas and other Ras-driven cancers, and the delivery of anti-p21Ras scFv by recombinant adenovirus and CIK cells may play an essential role in the therapy of Ras-driven cancers.

Complete genome sequences of three newly discovered cacao mild mosaic virus isolates from Theobroma cacao L. in Brazil and Puerto Rico and evidence for recombination

To analyze the DNA virome associated with cacao (Theobroma cacao L.) trees showing virus-like symptoms in Brazil (BR) and Puerto Rico (PR) during 2018-2019, total DNA was isolated from symptomatic leaves and subjected to high-throughput Illumina sequencing. The assembled complete badnaviral genome sequences were verified by PCR amplification, cloning, and DNA sequencing. Based on pairwise distances and phylogenetic analysis, three badnaviral genomes were identified, and these viruses were found to be isolates of the previously described cacao mild mosaic virus (CaMMV). The three genomes were 7,520, 7,524, and 7,514 bp in size for the isolates CaMMV-BR321, CaMMV-BR322, and CaMMV-PR3, respectively. Each genome contained four predicted open reading frames: ORFs 1-3 and ORFY. The CaMMV-PR3 isolate was identified as a probable recombinant, with a CaMMV-BR-like virus as the major parent.

Evolutionary dynamics of the human pseudoautosomal regions

Recombination between the X and Y human sex chromosomes is limited to the two pseudoautosomal regions (PARs) that present quite distinct evolutionary origins. Despite the crucial importance for male meiosis, genetic diversity patterns and evolutionary dynamics of these regions are poorly understood. In the present study, we analyzed and compared the genetic diversity of the PAR regions using publicly available genomic sequences encompassing both PAR1 and PAR2. Comparisons were performed through allele diversities, linkage disequilibrium status and recombination frequencies within and between X and Y chromosomes. In agreement with previous studies, we confirmed the role of PAR1 as a male-specific recombination hotspot, but also observed similar characteristic patterns of diversity in both regions although male recombination occurs at PAR2 to a much lower extent (at least one recombination event at PAR1 and in ≈1% in normal male meioses at PAR2). Furthermore, we demonstrate that both PARs harbor significantly different allele frequencies between X and Y chromosomes, which could support that recombination is not sufficient to homogenize the pseudoautosomal gene pool or is counterbalanced by other evolutionary forces. Nevertheless, the observed patterns of diversity are not entirely explainable by sexually antagonistic selection. A better understanding of such processes requires new data from intergenerational transmission studies of PARs, which would be decisive on the elucidation of PARs evolution and their role in male-driven heterosomal aneuploidies.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

Recombinant Nontypeable Genotype II Human Noroviruses in the Americas

We report multiple nontypeable genotype II noroviruses circulating in South America; nucleotides differed by >25% from those of other genotypes. These viruses have been circulating in the Americas for ≈20 years and show recombination with other genotypes. Clues to norovirus natural history can guide development of treatment and prevention plans.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.