Environmental Conditions Determine the Course and Outcome of Phytoplankton Chytridiomycosis

Chytrid fungi are highly potent parasites of phytoplankton. They are thought to force phytoplankton organisms into an evolutionary arms race with high population diversity as the outcome. The underlying selection regime is known as Red Queen dynamics. However, our study suggests a more complex picture for chytrid parasitism in the cyanobacterium Planktothrix. Laboratory experiments identified a “cold thermal refuge”, inside which Planktothrix can grow without chytrid infection. A field study in two Norwegian lakes underlined the ecological significance of this finding. The study utilized sediment DNA as a biological archive in combination with existing monitoring data. In one lake, temperature and light conditions forced Planktothrix outside the thermal refuge for most of the growing season. This probably resulted in Red Queen dynamics as suggested by a high parasitic pressure exerted by chytrids, an increase in Planktothrix genotype diversity over time, and a correlation between Planktothrix genotype diversity and duration of bloom events. In the second lake, a colder climate allowed Planktothrix to largely stay inside the thermal refuge. The parasitic pressure exerted by chytrids and Planktothrix genotype diversity remained low, indicating that Planktothrix successfully evaded the Red Queen dynamics. Episodic Planktothrix blooms were observed during spring and autumn circulation, in the metalimnion or under the ice. Interestingly, both lakes were dominated by the same or related Planktothrix genotypes. Taken together, our data suggest that, depending on environmental conditions, chytrid parasitism can impose distinct selection regimes on conspecific phytoplankton populations with similar genotype composition, causing these populations to behave and perhaps to evolve differently.     

SUBPOPULATION DIFFERENTIATION ASSOCIATED WITH NONRIBOSOMAL PEPTIDE SYNTHETASE GENE CLUSTER DYNAMICS IN THE CYANOBACTERIUM PLANKTOTHRIX SPP.1

The Plankthotrix Anagn. et Komárek population in the mesotrophic Lake Steinsfjorden has been intensively studied over several decades. This Planktothrix population produces a number of different classes of oligopetides. However, over the study period, only four main oligopeptide profiles (chemotypes) have been associated with the strains isolated from the lake. The chemotypes show distinct interactions with the environment, demonstrated by shifts in abundance along time series and vertical profiles. Here, we present genetic analysis of nonribosomal peptide synthetase (NRPS) gene regions in strains representing the four Planktothrix chemotypes in Lake Steinsfjorden. On the basis of phylogenetic analyses, we show that the NRPS genes for microcystin (mcy) and cyanopeptolin (oci) display the same clustering as do the chemotypes. Nucleotide diversity in mcy and oci was significantly higher between strains of different chemotypes than between strains of the same chemotype. K_a/K_s (nonsynonymous vs. synonymous mutations) values indicated positive selection in several polymorphic regions of the mcy and oci genes. Notably, incongruence between the phylogenetic trees for different gene segments and split decomposition analyses for segments of oci suggested horizontal gene transfer (HGT) events between strains showing different oligopeptide profiles. The oci HGT region encodes a module responsible for incorporating a variable amino acid in cyanopeptolin and is one of the regions suggested to be under positive selection. Taken together, our data suggest that there are four genetically distinct sympatric subpopulations—displayed as distinct chemotypes—in Lake Steinsfjorden. The diversification process of the chemotypes, and consequently the subpopulations, is driven by HGT and reinforced by positive selection of the corresponding NRPS gene regions.     

Frequent recombination shapes the epidemic population structure of Planktothrix (Cyanoprokaryota) in Italian subalpine lakes

The planktonic genus Planktothrix, as other cyanobacteria, shows signals of both homologous and nonhomologous recombination. However, the frequency of recombination and its effect on Planktothrix population structuring is unknown. We isolated 290 Planktothrix strains from seven neighboring lakes in the subalpine Italian region and analyzed these using multilocus sequence typing. Four of six loci analyzed were polymorphic, resulting in 20 distinct multilocus genotypes. Association indices among alleles at different loci were suggestive of an “epidemic population structure,” resulting from an explosive (and temporary) dominance of one genotype against a panmictic background. ClonalFrame analyses supported this view by detecting: (i) three major clades affected by three distinct recombination events, (ii) a recombination rate about equal to the mutation rate, and (iii) the fact that recombination had an impact on introducing molecular diversity more than double the mutation rate. Furthermore, analysis of molecular variance over an annual cycle in three of seven lakes revealed that both local clonal expansion and recombination processes affected among‐lake diversity. Our observations suggest that recombination affects microevolution of Planktothrix and that an epidemic structure can emerge in populations of this genus.     

Genetic variation of adenylation domains of the anabaenopeptin synthesis operon and evolution of substrate promiscuity

Anabaenopeptins (AP) are bioactive cyclic hexapeptides synthesized nonribosomally in cyanobacteria. APs are characterized by several conserved motifs, including the ureido bond, N-methylation in position 5, and d-Lys in position 2. All other positions of the AP molecule are variable, resulting in numerous structural variants. We have identified a nonribosomal peptide synthetase (NRPS) operon from Planktothrix agardhii strain CYA126/8 consisting of five genes (apnA to apnE) encoding six NRPS modules and have confirmed its role in AP synthesis by the generation of a mutant via insertional inactivation of apnC. In order to correlate the genetic diversity among adenylation domains (A domains) with AP structure variation, we sequenced the A domains of all six NRPS modules from seven Planktothrix strains differing in the production of AP congeners. It is remarkable that single strains coproduce APs bearing either of the chemically divergent amino acids Arg and Tyr in exocyclic position 1. Since the A domain of the initiation module (the ApnA A₁ domain) has been proposed to activate the amino acid incorporated into exocyclic position 1, we decided to analyze this domain both biochemically and phylogenetically. Only ApnA A₁ enzymes from strains producing AP molecules containing Arg or Tyr in position 1 were found to activate these two chemically divergent amino acids in vitro. Phylogenetic analysis of apn A domain sequences revealed that strains with a promiscuous ApnA A₁ domain are derived from an ancestor that activates only Arg. Surprisingly, positive selection appears to affect only three codons within the apnA A₁ gene, suggesting that this remarkable promiscuity has evolved from point mutations only.     

DISTRIBUTION OF OLIGOPEPTIDE CHEMOTYPES OF THE CYANOBACTERIUM PLANKTOTHRIX AND THEIR PERSISTENCE IN SELECTED LAKES IN FENNOSCANDIA1

Eighty‐seven Planktothrix (Anagnostidis and Komàrek 1988) strains isolated from 13 lakes in Scandinavia and Finland between 1964 and 2007 were screened for oligopeptides. Forty‐six individual compounds were detected in total, belonging to the structural classes anabaenopeptins (six variants), aeruginosins (six variants), cyanopeptolins (21 variants), microcystins (five variants), microginins (two variants), and microviridins (two variants). Oscillatorin was also found. Three additional compounds could not be assigned to known oligopeptide classes. Thirty oligopeptides have not been described in previous studies. Of these new compounds, five were aeruginosins and 20 cyanopeptolins. The number of oligopeptides per strain ranged from one to 13. No oligopeptide‐free strains were found, suggesting that oligopeptide production is vital for Planktothrix. On the basis of their oligopeptide patterns, the Planktothrix strains of the present study were assigned to 17 chemotypes. Three major chemotypes occurred in up to six lakes. One chemotype occurred in lakes around the city of Oslo (Norway), on the Finnish island Fasta Åland, which is situated in the Baltic Sea, and on the Finnish mainland. This wide distribution suggests that chemotypes can be subjects of recurrent dispersal and/or strong directional selection. Lake size, maximum depth, and nutrient availability appeared to be of minor importance for the ability of some chemotypes to colonize a water body successfully as long as the general requirements of Planktothrix were met. Four chemotypes were reisolated from the Oslo lake district over a period of 33–40 years, suggesting that they have been members of local Planktothrix populations for decades.     

Transposons Inactivate Biosynthesis of the Nonribosomal Peptide Microcystin in Naturally Occurring Planktothrix spp.

The filamentous cyanobacteria Planktothrix spp. occur in the temperate region of the Northern hemisphere. The red-pigmented Planktothrix rubescens bacteria occur in deep, physically stratified, and less eutrophic lakes. Planktothrix is a known producer of the toxic heptapeptide microcystin (MC), which is produced nonribosomally by a large enzyme complex consisting of peptide synthetases and polyketide synthases encoded by a total of nine genes (mcy genes). Planktothrix spp. differ in their cellular MC contents as well as the production of MC variants; however, the mechanisms favoring this diversity are not understood. Recently, the occurrence of Planktothrix strains containing all mcy genes but lacking MC has been reported. In this study, 29 such strains were analyzed to find out if mutations of the mcy genes lead to the inability to synthesize MC. Two deletions, spanning 400 bp (in mcyB; one strain) and 1,869 bp (in mcyHA; three strains), and three insertions (IS), spanning 1,429 bp (in mcyD; eight strains), 1,433 bp (in mcyEG; one strain), and 1,433 bp (in mcyA; one strain), were identified. Though found in different genes and different isolates and transcribed in opposite directions, IS were found to be identical and contained conserved domains assigned to transposable elements. Using mutation-specific primers, two insertions (in mcyD and mcyA) and one deletion (in mcyHA) were found regularly in populations of P. rubescens in different lakes. The results demonstrate for the first time that different mutations resulting in inactivation of MC synthesis do occur frequently and make up a stable proportion of the mcy gene pool in Planktothrix populations over several years.     

Multiple barriers to recombination between divergent HIV-1 variants revealed by a dual-marker recombination assay

Recombination is a major force for generating human immunodeficiency virus type 1 (HIV-1) diversity and produces numerous recombinants circulating in the human population. We previously established a cell-based system using green fluorescent protein gene (gfp) as a reporter to study the mechanisms of HIV-1 recombination. We now report an improved system capable of detecting recombination using authentic viral sequences. Frameshift mutations were introduced into the gag gene so that parental viruses do not express full-length Gag; however, recombination can generate a progeny virus that expresses a functional Gag. We demonstrate that this Gag reconstitution assay can be used to detect recombination between two group M HIV-1 variants of the same or of different subtypes. Using both gfp and gag assays, we found that, similar to group M viruses, group O viruses also recombine frequently. When recombination between a group M virus and a group O virus was examined, we found three distinct barriers for intergroup recombination. First, similar to recombination within group M viruses, intergroup recombination is affected by the identity of the dimerization initiation signal (DIS); variants with the same DIS recombined at a higher rate than those with different DIS. Second, using the gfp recombination assay, we showed that intergroup recombination occurs much less frequently than intragroup recombination, even though the gfp target sequence is identical in all viruses. Finally, Gag reconstitution between variants from different groups is further reduced compared with green fluorescent protein, indicating that sequence divergence interferes with recombination efficiency in the gag gene. Compared with identical sequences, we estimate that recombination rates are reduced by 3-fold and by 10- to 13-fold when the target regions in gag contain 91% and 72-73% sequence identities, respectively. These results show that there are at least three distinct mechanisms preventing exchange of genetic information between divergent HIV-1 variants from different groups.     

Comparison of Cyanopeptolin Genes in Planktothrix, Microcystis, and Anabaena Strains: Evidence for Independent Evolution within Each Genus

The major cyclic peptide cyanopeptolin 1138, produced by Planktothrix strain NIVA CYA 116, was characterized and shown to be structurally very close to the earlier-characterized oscillapeptin E. A cyanopeptolin gene cluster likely to encode the corresponding peptide synthetase was sequenced from the same strain. The 30-kb oci gene cluster contains two novel domains previously not detected in nonribosomal peptide synthetase gene clusters (a putative glyceric acid-activating domain and a sulfotransferase domain), in addition to seven nonribosomal peptide synthetase modules. Unlike in two previously described cyanopeptolin gene clusters from Anabaena and Microcystis, a halogenase gene is not present. The three cyanopeptolin gene clusters show similar gene and domain arrangements, while the binding pocket signatures deduced from the adenylation domain sequences and the additional tailoring domains vary. This suggests loss and gain of tailoring domains within each genus, after the diversification of the three clades, as major events leading to the present diversity. The ABC transporter genes associated with the cyanopeptolin gene clusters form a monophyletic clade and accordingly are likely to have evolved as part of the functional unit. Phylogenetic analyses of adenylation and condensation domains, including domains from cyanopeptolins and microcystins, show a closer similarity between the Planktothrix and Microcystis cyanopeptolin domains than between these and the Anabaena domain. No clear evidence of recombination between cyanopeptolins and microcystins could be detected. There were no strong indications of horizontal gene transfer of cyanopeptolin gene sequences across the three genera, supporting independent evolution within each genus.     

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes,virulence genes,and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.     

Genetic diversity and population structure of Leptosphaeria maculans isolates in Western Canada

Leptosphaeria maculans is a serious concern for canola production worldwide. For effective disease management, knowledge of the pathogen's genetic variability and population structure is a prerequisite. In this study, whole-genome sequencing was performed for 162 of 1590 L. maculans isolates collected in the years 2007–2008 and 2012–2014 in Western Canada. DNA variants in genome-wide and specific regions including avirulence (Avr) genes were characterized. A total of 31,870 high-quality polymorphic DNA variants were used to study L. maculans genetic diversity and population structure. Cluster analysis showed that 150 isolates were clustered into 2 main groups and 4 subgroups by DNA variants located in either Avr or small secreted protein-encoding genes and into 2 main groups and 6 subgroups by genome-wide variants. The analysis of nucleotide diversity and differentiation also confirmed genetic variation within a population and among populations. Principal component analysis with genome-wide variants showed that the isolates collected in 2012–2014 were more genetically diverse than those collected in 2007–2008. Population structure analysis discovered three distinct sub-populations. Although isolates from Saskatchewan and Alberta were of similar genetic composition, Manitoba isolates were highly diverse. Genome-wide association study detected DNA variants in genes AvrLm4-7, Lema_T86300, and Lema_T86310 associated with the years of collection.     

Species Persistence with Hybridization in Toad-Headed Lizards Driven by Divergent Selection and Low Recombination

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.     

Phylogenetic and Structural Diversity in the Feline Leukemia Virus Env Gene

Feline leukemia virus (FeLV) belongs to the genus Gammaretrovirus, and causes a variety of neoplastic and non-neoplastic diseases in cats. Alteration of viral env sequences is thought to be associated with disease specificity, but the way in which genetic diversity of FeLV contributes to the generation of such variants in nature is poorly understood. We isolated FeLV env genes from naturally infected cats in Japan and analyzed the evolutionary dynamics of these genes. Phylogenetic reconstructions separated our FeLV samples into three distinct genetic clusters, termed Genotypes I, II, and III. Genotype I is a major genetic cluster and can be further classified into Clades 1–7 in Japan. Genotypes were correlated with geographical distribution; Genotypes I and II were distributed within Japan, whilst FeLV samples from outside Japan belonged to Genotype III. These results may be due to geographical isolation of FeLVs in Japan. The observed structural diversity of the FeLV env gene appears to be caused primarily by mutation, deletion, insertion and recombination, and these variants may be generated de novo in individual cats. FeLV interference assay revealed that FeLV genotypes did not correlate with known FeLV receptor subgroups. We have identified the genotypes which we consider to be reliable for evaluating phylogenetic relationships of FeLV, which embrace the high structural diversity observed in our sample. Overall, these findings extend our understanding of Gammaretrovirus evolutionary patterns in the field, and may provide a useful basis for assessing the emergence of novel strains and understanding the molecular mechanisms of FeLV transmission in cats.     

Population Genetic Structure of a Microalgal Species under Expansion

Biological invasions often cause major perturbations in the environment and are well studied among macroorganisms. Less is known about invasion by free-living microbes. Gonyostomum semen (Raphidophyceae) is a freshwater phytoplankton species that has increased in abundance in Northern Europe since the 1980''s and has expanded its habitat range. In this study, we aimed to determine the genetic population structure of G. semen in Northern Europe and to what extent it reflects the species'' recent expansion. We sampled lakes from 12 locations (11 lakes) in Norway, Sweden and Finland. Multiple strains from each location were genotyped using Amplified Fragment Length Polymorphism (AFLP). We found low differentiation between locations, and low gene diversity within each location. Moreover, there was an absence of genetic isolation with distance (Mantel test, p = 0.50). According to a Bayesian clustering method all the isolates belonged to the same genetic population. Together our data suggest the presence of one metapopulation and an overall low diversity, which is coherent with a recent expansion of G. semen.     

Genetic diversity of NBS–LRR class disease-resistance gene analogs in cultivated and wild eggplants

The genes encoding the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) motifs constitute a large gene family in plants and have attracted much interest, because most of the plant disease-resistance genes that have been cloned are from this gene family. In this study, degenerate oligonucleotide primers, designed on the basis of conserved regions of the NBS domains from known plant resistance genes, were used to isolate resistance gene analogs (RGAs) from cultivated and wild eggplants, i.e., S. melongena, S. aethiopicum gr. Gilo, S. linnaeanum, S. integrifolium, S. sisymbriifolium, and S. khasianum. Sequence analysis indicated that the cloned eggplant RGAs belong to the non-TIR–NBS–LRR type, which are very similar to the R genes or the RGAs identified in other plant species, especially Solanaceae plants, suggesting the existence of common ancestors. Wide genetic diversity of eggplant RGAs was observed both in interspecific and intraspecific sequences, and eight distinct families of eggplant RGAs were identified. Further studies revealed a high average ratio of synonymous to non-synonymous substitution and a low level of recombination. These results suggest that NBS-encoding sequences of RGAs in cultivated and wild eggplants are subject to gradual accumulation of mutations leading to purifying selection. This is the first report of NBS–LRR class RGAs in eggplants.     

Population Structure and Antimicrobial Resistance Profiles of Streptococcus suis Serotype 2 Sequence Type 25 Strains

Strains of serotype 2 Streptococcus suis are responsible for swine and human infections. Different serotype 2 genetic backgrounds have been defined using multilocus sequence typing (MLST). However, little is known about the genetic diversity within each MLST sequence type (ST). Here, we used whole-genome sequencing to test the hypothesis that S. suis serotype 2 strains of the ST25 lineage are genetically heterogeneous. We evaluated 51 serotype 2 ST25 S. suis strains isolated from diseased pigs and humans in Canada, the United States of America, and Thailand. Whole-genome sequencing revealed numerous large-scale rearrangements in the ST25 genome, compared to the genomes of ST1 and ST28 S. suis strains, which result, among other changes, in disruption of a pilus island locus. We report that recombination and lateral gene transfer contribute to ST25 genetic diversity. Phylogenetic analysis identified two main and distinct Thai and North American clades grouping most strains investigated. These clades also possessed distinct patterns of antimicrobial resistance genes, which correlated with acquisition of different integrative and conjugative elements (ICEs). Some of these ICEs were found to be integrated at a recombination hot spot, previously identified as the site of integration of the 89K pathogenicity island in serotype 2 ST7 S. suis strains. Our results highlight the limitations of MLST for phylogenetic analysis of S. suis, and the importance of lateral gene transfer and recombination as drivers of diversity in this swine pathogen and zoonotic agent.     

Prevalence of Viral Photosynthetic and Capsid Protein Genes from Cyanophages in Two Large and Deep Perialpine Lakes

Cyanophages are important components of aquatic ecosystems, but their genetic diversity has been little investigated in freshwaters. A yearlong survey was conducted in surface waters of the two largest natural perialpine lakes in France (Lake Annecy and Lake Bourget) to investigate part of this cyanophage diversity through the analysis of both structural (e.g., g20) and functional (e.g., psbA) genes. We found that these cyanophage signature genes were prevalent throughout the year but that the community compositions of g20 cyanomyoviruses were significantly different between the two lakes. In contrast, psbA-containing cyanophages seemed to be more similar between the two ecosystems. We also found that a large proportion of g20 sequences grouped with cyanomyophage isolates. psbA sequences, belonging to phages of Synechococcus spp., were characterized by distinct triplet motifs (with a novel viral triplet motif, EFE). Thus, our results show that cyanophages (i) are a diverse viral community in alpine lakes and (ii) are clearly distinct from some other freshwater and marine environments, suggesting the influence of unique biogeographic factors.     

Insights into an Unusual Nonribosomal Peptide Synthetase Biosynthesis: IDENTIFICATION AND CHARACTERIZATION OF THE GE81112 BIOSYNTHETIC GENE CLUSTER*

The GE81112 tetrapeptides (1–3) represent a structurally unique class of antibiotics, acting as specific inhibitors of prokaryotic protein synthesis. Here we report the cloning and sequencing of the GE81112 biosynthetic gene cluster from Streptomyces sp. L-49973 and the development of a genetic manipulation system for Streptomyces sp. L-49973. The biosynthetic gene cluster for the tetrapeptide antibiotic GE81112 (getA-N) was identified within a 61.7-kb region comprising 29 open reading frames (open reading frames), 14 of which were assigned to the biosynthetic gene cluster. Sequence analysis revealed the GE81112 cluster to consist of six nonribosomal peptide synthetase (NRPS) genes encoding incomplete di-domain NRPS modules and a single free standing NRPS domain as well as genes encoding other biosynthetic and modifying proteins. The involvement of the cloned gene cluster in GE81112 biosynthesis was confirmed by inactivating the NRPS gene getE resulting in a GE81112 production abolished mutant. In addition, we characterized the NRPS A-domains from the pathway by expression in Escherichia coli and in vitro enzymatic assays. The previously unknown stereochemistry of most chiral centers in GE81112 was established from a combined chemical and biosynthetic approach. Taken together, these findings have allowed us to propose a rational model for GE81112 biosynthesis. The results further open the door to developing new derivatives of these promising antibiotic compounds by genetic engineering.     

Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

Analyses of sequence polymorphism and haplotype diversity of LEAFY genes revealed post-domestication selection in the Chinese elite maize inbred lines

Post-domestication selection refers to the artificial selection on the loci controlling important agronomic traits during the process of genetic improvement in a population. The maize genes Zfl1 and Zfl2, duplicate orthologs of Arabidopsis LEAFY, are key regulators in plant branching, inflorescence and flower development, and reproduction. In this study, the full gene sequences of Zfl1 and Zfl2 from 62 Chinese elite inbred lines were amplified to evaluate their nucleotide polymorphisms and haplotype diversities. A total of 254 and 192 variants that included SNPs and indels were identified from the full sequences of Zfl1 and Zfl2, respectively. Although most of the variants were found to be located in the non-coding regions, the polymorphisms of CDS sequences classified Zfl1 into 16 haplotypes encoding 16 different proteins and Zfl2 into 18 haplotypes encoding eight different proteins. The population of Huangzaosi and its derived lines showed statistically significant signals of post-domestication selection on the Zfl1 CDS sequences, as well as lower nucleotide polymorphism and haplotype diversity than the whole set. However, the Zfl2 locus was only selected for in the heterotic group Reid. Further evidence revealed that at least 17 recombination events contributed to the genetic and haplotype diversities at the Zfl1 locus and 16 recombination events at the Zfl2 locus.     

Pervasive prophage recombination occurs during evolution of spore-forming Bacilli

Phages are the main source of within-species bacterial diversity and drivers of horizontal gene transfer, but we know little about the mechanisms that drive genetic diversity of these mobile genetic elements (MGEs). Recently, we showed that a sporulation selection regime promotes evolutionary changes within SPβ prophage of Bacillus subtilis, leading to direct antagonistic interactions within the population. Herein, we reveal that under a sporulation selection regime, SPβ recombines with low copy number phi3Ts phage DNA present within the B. subtilis population. Recombination results in a new prophage occupying a different integration site, as well as the spontaneous release of virulent phage hybrids. Analysis of Bacillus sp. strains suggests that SPβ and phi3T belong to a distinct cluster of unusually large phages inserted into sporulation-related genes that are equipped with a spore-related genetic arsenal. Comparison of Bacillus sp. genomes indicates that similar diversification of SPβ-like phages takes place in nature. Our work is a stepping stone toward empirical studies on phage evolution, and understanding the eco-evolutionary relationships between bacteria and their phages. By capturing the first steps of new phage evolution, we reveal striking relationship between survival strategy of bacteria and evolution of their phages.Subject terms: Bacterial genetics, Evolution