Meiotic recombination profiling of interspecific hybrid F1 tomato pollen by linked read sequencing

Genome wide screening of pooled pollen samples from a single interspecific F1 hybrid obtained from a cross between tomato, Solanum lycopersicum and its wild relative, Solanum pimpinellifolium using linked read sequencing of the haploid nuclei, allowed profiling of the crossover (CO) and gene conversion (GC) landscape. We observed a striking overlap between cold regions of CO in the male gametes and our previously established F6 recombinant inbred lines (RILs) population. COs were overrepresented in non‐coding regions in the gene promoter and 5′UTR regions of genes. Poly‐A/T and AT rich motifs were found enriched in 1 kb promoter regions flanking the CO sites. Non‐crossover associated allelic and ectopic GCs were detected in most chromosomes, confirming that besides CO, GC represents also a source for genetic diversity and genome plasticity in tomato. Furthermore, we identified processed break junctions pointing at the involvement of both homology directed and non‐homology directed repair pathways, suggesting a recombination machinery in tomato that is more complex than currently anticipated.     

Effects of XRCC2 and RAD51B mutations on somatic and meiotic recombination in Arabidopsis thaliana

Homologous recombination is key to the maintenance of genome integrity and the creation of genetic diversity. At the mechanistic level, recombination involves the invasion of a homologous DNA template by broken DNA ends, repair of the break and exchange of genetic information between the two DNA molecules. Invasion of the template in eukaryotic cells is catalysed by the RAD51 and DMC1 recombinases, assisted by a number of accessory proteins, including the RAD51 paralogues. Eukaryotic genomes encode a variable number of RAD51 paralogues, ranging from two in yeast to five in animals and plants. The RAD51 paralogues form at least two distinct protein complexes, believed to play roles in the assembly and stabilization of the RAD51‐DNA nucleofilament. Somatic recombination assays and immunocytology confirm that the three ‘non‐meiotic’ paralogues of Arabidopsis, RAD51B, RAD51D and XRCC2, are involved in somatic homologous recombination, and that they are not required for the formation of radioinduced RAD51 foci. Given the presence of all five proteins in meiotic cells, the apparent absence of a meiotic role for RAD51B, RAD51D and XRCC2 is surprising, and perhaps simply the result of a more subtle meiotic phenotype in the mutants. Analysis of meiotic recombination confirms this, showing that the absence of XRCC2, and to a lesser extent RAD51B, but not RAD51D, increases rates of meiotic crossing over. The roles of RAD51B and XRCC2 in recombination are thus not limited to mitotic cells.     

Cryptic ecology among host generalist Campylobacter jejuni in domestic animals

Homologous recombination between bacterial strains is theoretically capable of preventing the separation of daughter clusters, and producing cohesive clouds of genotypes in sequence space. However, numerous barriers to recombination are known. Barriers may be essential such as adaptive incompatibility, or ecological, which is associated with the opportunities for recombination in the natural habitat. Campylobacter jejuni is a gut colonizer of numerous animal species and a major human enteric pathogen. We demonstrate that the two major generalist lineages of C. jejuni do not show evidence of recombination with each other in nature, despite having a high degree of host niche overlap and recombining extensively with specialist lineages. However, transformation experiments show that the generalist lineages readily recombine with one another in vitro. This suggests ecological rather than essential barriers to recombination, caused by a cryptic niche structure within the hosts.     

Comparative analyses of clinical and environmental populations of Cryptococcus neoformans in Botswana

Cryptococcus neoformans var. grubii (Cng) is the most common cause of fungal meningitis, and its prevalence is highest in sub‐Saharan Africa. Patients become infected by inhaling airborne spores or desiccated yeast cells from the environment, where the fungus thrives in avian droppings, trees and soil. To investigate the prevalence and population structure of Cng in southern Africa, we analysed isolates from 77 environmental samples and 64 patients. We detected significant genetic diversity among isolates and strong evidence of geographic structure at the local level. High proportions of isolates with the rare MAT a allele were observed in both clinical and environmental isolates; however, the mating‐type alleles were unevenly distributed among different subpopulations. Nearly equal proportions of the MAT a and MATα mating types were observed among all clinical isolates and in one environmental subpopulation from the eastern part of Botswana. As previously reported, there was evidence of both clonality and recombination in different geographic areas. These results provide a foundation for subsequent genomewide association studies to identify genes and genotypes linked to pathogenicity in humans.     

Reanalysis suggests that genomic islands of speciation are due to reduced diversity,not reduced gene flow

The metaphor of ‘genomic islands of speciation’ was first used to describe heterogeneous differentiation among loci between the genomes of closely related species. The biological model proposed to explain these differences was that the regions showing high levels of differentiation were resistant to gene flow between species, while the remainder of the genome was being homogenized by gene flow and consequently showed lower levels of differentiation. However, the conditions under which such differentiation can occur at multiple unlinked loci are restrictive; additionally, essentially, all previous analyses have been carried out using relative measures of divergence, which can be misleading when regions with different levels of recombination are compared. Here, we test the model of differential gene flow by asking whether absolute divergence is also higher in the previously identified ‘islands’. Using five species pairs for which full sequence data are available, we find that absolute measures of divergence are not higher in genomic islands. Instead, in all cases examined, we find reduced diversity in these regions, a consequence of which is that relative measures of divergence are abnormally high. These data therefore do not support a model of differential gene flow among loci, although islands of relative divergence may represent loci involved in local adaptation. Simulations using the program IMa2 further suggest that inferences of any gene flow may be incorrect in many comparisons. We instead present an alternative explanation for heterogeneous patterns of differentiation, one in which postspeciation selection generates patterns consistent with multiple aspects of the data.     

BLM and SLX4 play opposing roles in recombination‐dependent replication at human telomeres

Alternative lengthening of telomeres (ALT) is a telomere lengthening pathway that predominates in aggressive tumors of mesenchymal origin; however, the underlying mechanism of telomere synthesis is not fully understood. Here, we show that the BLM–TOP3A–RMI (BTR) dissolvase complex is required for ALT‐mediated telomere synthesis. We propose that recombination intermediates formed during strand invasion are processed by the BTR complex, initiating rapid and extensive POLD3‐dependent telomere synthesis followed by dissolution, with no overall exchange of telomeric DNA. This process is counteracted by the SLX4–SLX1–ERCC4 complex, which promotes resolution of the recombination intermediate, resulting in telomere exchange in the absence of telomere extension. Our data are consistent with ALT being a conservative DNA replication process, analogous to break‐induced replication, which is dependent on BTR and counteracted by SLX4 complex‐mediated resolution events.     

Molecular Characterization of Whole Genomic RNA2 From Iranian Isolates of Grapevine Fanleaf Virus

Grapevine fanleaf virus (GFLV) is the major causal agent of the grapevine degeneration disease. To characterize the genomic RNA2 segment from Iranian isolates of GFLV, leaf samples were collected from infected vineyards in different locations with a long history of vine cultivation. Four isolates were selected for cloning and sequencing on the basis of the restriction profiles of RT‐PCR products. The sequencing data revealed that the RNA2 of the Iranian GFLV isolates were the shortest compared with that of all previously described GFLV isolates. The sizes were 3730 nucleotides (nt) for Shir‐Amin and Urmia isolates and 3749 nt for Takestan and Bonab isolates (excluding the poly (A) tail), due to deletion events in both 5′ and 3′ non‐coding regions. In the phylogenetic tree based on the full‐length nucleotide sequences of GFLV RNA2, all the GFLV isolates clustered into two groups with the exception of the Hungarian isolate (GHu). The Iranian isolates grouped as a distinct cluster. Recombination analyses showed that GFLV‐NW (Germany), GFLV‐F13 (reference isolate), GFLV isolate Shir‐Amin (Iran) and Arabis mosaic virus isolate Lv were recombinant isolates and one of their parents belonged to the same lineage as the Iranian isolates. These findings suggest that these isolates originated from a common ancestor.     

Interspecific plastidial recombination in the diatom genus Pseudo‐nitzschia

Plastids are usually uni‐parentally inherited and genetic recombination between these organelles is seldom observed. The genus Pseudo‐nitzschia, a globally relevant marine diatom, features bi‐parental plastid inheritance in the course of sexual reproduction. This observation inspired the recombination detection we pursued in this paper over a ~1,400‐nucleotide‐long region of the plastidial rbcL, a marker used in both molecular taxonomy and phylogenetic studies in diatoms. Among all the rbcL‐sequences available in web‐databases for Pseudo‐nitzschia, 42 haplotypes were identified and grouped in five clusters by Bayesian phylogeny. Signs of hybridization were evident in four of five clusters, at both intra‐ and interspecific levels, suggesting that, in diatoms, (i) plastidial recombination is not absent and (ii) hybridization can play a role in speciation of Pseudo‐nitzschia spp.     

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.     

First evidence of the occurrence of Turnip mosaic virus in Ukraine and molecular characterization of its isolate

A total of 54 samples of Brassicaceae crops showing symptoms of mosaic, mottling, vein banding and/or leaf deformation were collected in Kyiv region (northern central part of Ukraine) in 2014–2015. A half of collected samples was found to be infected with Turnip mosaic virus (TuMV), and TuMV was detected in samples from Brassica oleracea var. capitata (cabbage), Raphanus sativus, Brassica juncea, Raphanus sp., Sinapis alba, Camelina sativa and Bunias orientalis (weed). The full‐length sequence of the genomic RNA of a Ukrainian isolate (UKR9), which was isolated from cabbage, was determined. Recombination analysis of UKR9 isolate showed that this isolate was an interlineage recombinant of world‐Brassica and Asian‐Brassica/Raphanus phylogenetic groups. This study shows for the first time the occurrence of TuMV in Ukraine.     

Recombination Signals In The rpoC1 Gene Indicate Gene‐Flow Between Planktothrix (Cyanoprokaryota) Species

The delineation of species boundaries in the potentially harmful cyanobacterium Planktothrix Anagnostidis et Komárek 1988 is particularly tangled. Genetic recombination has been invoked to explain the occurrence of overlapping biological traits among recognized species. Although horizontal gene transfer is shown as a driver of diversification in this genus, clear evidence for homologous recombination at the single gene level is still lacking. Several Planktothrix strains (n = 244) sampled in eight fresh water lakes in north Italy were characterized by sequencing the rpoC1 gene, a molecular marker previously proposed to discriminate between species. Six haplotypes were detected, four of which are newly described. A relevant number of rpoC1 sequences (n = 54) showed evidence of homologous recombination. By comparing the sequences produced in the work presented here to those available on GenBank for the genus, multiple recombination events were tracked between haplotypes associated to P. rubescens, P. suspensa and P. agardhii, the latter being a species not found in our survey. Recombination signals were in form of (i) a vast mosaic structure present in the alignment of rpoC1 haplotypes, (ii) multiple and statistically significant paths in the split decomposition network connecting these haplotypes and (iii) many individual crossing‐over events detected by means of recombination detection tests. Data suggest that the molecular evolution of the rpoC1 gene in the genus Planktothrix appears as strongly influenced by homologous recombination. In addition, rpoC1 diversity effectively tracks recombinational processes among species in the complex made up by P. rubescens, P. agardhii and P. suspensa, which are not isolated in terms of gene‐flow.     

Extensive paternal mtDNA leakage in natural populations of Drosophila melanogaster

Strict maternal inheritance is considered a hallmark of animal mtDNA. Although recent reports suggest that paternal leakage occurs in a broad range of species, it is still considered an exceptionally rare event. To evaluate the impact of paternal leakage on the evolution of mtDNA, it is essential to reliably estimate the frequency of paternal leakage in natural populations. Using allele‐specific real‐time quantitative PCR (RT‐qPCR), we show that heteroplasmy is common in natural populations with at least 14% of the individuals carrying multiple mitochondrial haplotypes. However, the average frequency of the minor mtDNA haplotype is low (0.8%), which suggests that this pervasive heteroplasmy has not been noticed before due to a lack of power in sequencing surveys. Based on the distribution of mtDNA haplotypes in the offspring of heteroplasmic mothers, we found no evidence for strong selection against one of the haplotypes. We estimated that the rate of paternal leakage is 6% and that at least 100 generations are required for complete sorting of mtDNA haplotypes. Despite the high proportion of heteroplasmic individuals in natural populations, we found no evidence for recombination between mtDNA molecules, suggesting that either recombination is rare or recombinant haplotypes are counter‐selected. Our results indicate that evolutionary studies using mtDNA as a marker might be biased by paternal leakage in this species.     

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole‐genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non‐UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.     

The condensin complexes play distinct roles to ensure normal chromosome morphogenesis during meiotic division in Arabidopsis

Meiosis is a specialized cell division essential for sexual reproduction. During meiosis the chromosomes are highly organized, and correct chromosome architecture is required for faithful segregation of chromosomes at anaphase I and II. Condensin is involved in chromosome organization during meiotic and mitotic cell divisions. Three condensin subunits, AtSMC4 and the condensin I and II specific subunits AtCAP‐D2 and AtCAP‐D3, respectively, have been studied for their role in meiosis. This has revealed that both the condensin I and condensin II complexes are required to maintain normal structural integrity of the meiotic chromosomes during the two nuclear divisions. Their roles appear functionally distinct in that condensin I is required to maintain normal compaction of the centromeric repeats and 45S rDNA, whereas loss of condensin II was associated with extensive interchromosome connections at metaphase I. Depletion of condensin is also associated with a slight reduction in crossover formation, suggesting a role during meiotic prophase I.     

The adaptive significance of chromosomal inversion polymorphisms in Drosophila melanogaster

Chromosomal inversions, structural mutations that reverse a segment of a chromosome, cause suppression of recombination in the heterozygous state. Several studies have shown that inversion polymorphisms can form clines or fluctuate predictably in frequency over seasonal time spans. These observations prompted the hypothesis that chromosomal rearrangements might be subject to spatially and/or temporally varying selection. Here, we review what has been learned about the adaptive significance of inversion polymorphisms in the vinegar fly Drosophila melanogaster, the species in which they were first discovered by Sturtevant in 1917. A large body of work provides compelling evidence that several inversions in this system are adaptive; however, the precise selective mechanisms that maintain them polymorphic in natural populations remain poorly understood. Recent advances in population genomics, modelling and functional genetics promise to greatly improve our understanding of this long‐standing and fundamental problem in the near future.     

Genome of an iconic Australian bird: High‐quality assembly and linkage map of the superb fairy‐wren (Malurus cyaneus)

The superb fairy‐wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high‐quality reference. To rectify that, we present the first high‐quality genome assembly of a Meliphagides representative: the superb fairy‐wren. We combined Illumina shotgun and mate‐pair sequences, PacBio long‐reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07‐Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high‐quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree‐based bird genetic maps, we find that the fairy‐wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high‐quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy‐wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.     

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora,Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1–3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52–57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.     

Chromosome arm‐specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector,Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub‐Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine‐scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of F_ST, absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.