Mitotic gene conversion, reciprocal recombination and gene replacement at the benA, beta-tubulin, locus of Aspergillus nidulans

Summary We have developed a procedure for determining the rates of mitotic recombination of an interrupted duplication created by integration of transforming plasmid sequences at the benA, beta-tubulin, locus of Aspergillus nidulans. Transformation of a strain carrying a benomyl-resistant benA allele with plasmid AIpGM4, which carries the wild-type benA allele and the pyr4 (orotidine-5-phosphate decarboxylase) gene of Neurospora crassa, creates an interrupted duplication with plasmid sequences flanked by two benA alleles, one wild type and one benomyl resdistant. Such transformants will not grow in the presence of high levels of benomyl. Mitotic recombination causes the loss of the wild-type benA allele or conversion of the wild-type to the mutant allele resulting in nuclei carrying only the benomylresistant allele. Conidia containing such nuclei can be selected on media with high benomyl allowing easy quantitation of mitotic recombination. We found that the rate of recombination giving rise to benomyl-resistant conidia was 4.6×10^-4. Reciprocal recombination leading to benomyl-resistant conidia lacking plasmid sequences occurred at a rate of 2.0×10^-4 and gene conversion leading to benomylresistant conidia occurred at a rate of 2.6×10^-4. We selected for reciprocal recombination leading to loss of pyr4 sequences on 5-fluoro-orotic acid and used this selection for two-step gene replacement of a mutant benA allele with the wild-type allele.     

rpsL+: a dominant selectable marker for gene replacement in mycobacteria

Molecular genetic manipulations in mycobacteria would benefit from procedures which efficiently select for double-crossover events by homologous recombination. Here we describe a vector-host system for gene replacement in mycobacteria, the utility of which was investigated using functional inactivation of the pyrF gene in Mycobacterium smegmatis as a model. This system is based on the expression of the wild-type rpsL gene coding for ribosomal protein S12 in a streptomycin-resistant host. Owing to the absence of a mycobacterial origin the plasmids are unable to replicate autonomously in mycobacteria. The first selection for maintenance of cloned sequences is conferred by the kanamycin-resistance gene. The second simultaneous selection by streptomycin is against maintenance of cloned sequences which contain the gene encoding the streptomycin-sensitive allele of the rpsL gene. By placing the gene for positive selection and that used for negative selection within and outside the target gene of interest, respectively, gene replacement is obtained. A one-step double selection procedure provides a means to distinguish strictly between gene replacement by double crossover versus homologous recombination by single crossover events. The system should have considerable potential for genera or species where single-crossover events or even illegitimate recombination are the predominant recombination mechanisms. It should also be of wide use for the construction of mutants without a selectable phenotype.     

Large-scale DNA polymorphism study of <Emphasis Type="Italic">Oryza sativa</Emphasis> and <Emphasis Type="Italic">O. rufipogon</Emphasis> reveals the origin and divergence of Asian rice

Polymorphism over ∼26 kb of DNA sequence spanning 22 loci and one region distributed on chromosomes 1, 2, 3 and 4 was studied in 30 accessions of cultivated rice, Oryza sativa, and its wild relatives. Phylogenetic analysis using all the DNA sequences suggested that O. sativa ssp. indica and ssp. japonica were independently domesticated from a wild species O. rufipogon. O. sativa ssp. indica contained substantial genetic diversity (π = 0.0024), whereas ssp. japonica exhibited extremely low nucleotide diversity (π = 0.0001) suggesting the origin of the latter from a small number of founders. O. sativa ssp. japonica contained a larger number of derived and fixed non-synonymous substitutions as compared to ssp. indica. Nucleotide diversity and genealogical history substantially varied across the 22 loci. A locus, RLD15 on chromosome 2, showed a distinct genealogy with ssp. japonica sequences distantly separated from those of O. rufipogon and O. sativa ssp. indica. Linkage disequilibrium (LD) was analyzed in two different regions. LD in O. rufipogon decays within 5 kb, whereas it extends to ∼50 kb in O. sativa ssp. indica. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Recombination, balancing selection and phylogenies in MHC and self-incompatibility genes.

Using a coalescent model of multiallelic balancing selection with recombination, the genealogical process as a function of recombinational distance from a site under selection is investigated. We find that the shape of the phylogenetic tree is independent of the distance to the site under selection. Only the timescale changes from the value predicted by Takahata's allelic genealogy at the site under selection, converging with increasing recombination to the timescale of the neutral coalescent. However, if nucleotide sequences are simulated over a recombining region containing a site under balancing selection, a phylogenetic tree constructed while ignoring such recombination is strongly affected. This is true even for small rates of recombination. Published studies of multiallelic balancing selection, i.e., the major histocompatibility complex (MHC) of vertebrates, gametophytic and sporophytic self-incompatibility of plants, and incompatibility of fungi, all observe allelic genealogies with unexpected shapes. We conclude that small absolute levels of recombination are compatible with these observed distortions of the shape of the allelic genealogy, suggesting a possible cause of these observations. Furthermore, we illustrate that the variance in the coalescent with recombination process makes it difficult to locate sites under selection and to estimate the selection coefficient from levels of variability.     

<Emphasis Type="Italic">Mhc</Emphasis> class II diversity and balancing selection in greater prairie-chickens

The major histocompatibility complex (Mhc) of domestic chickens has been characterized as small and relatively simple compared with that of mammals. However, there is growing evidence that the Mhc of many bird lineages may be more complex, even within the Order Galliformes. In this study, we measured genetic variation and balancing selection at Mhc loci in another galliform, the greater prairie-chicken. We cloned and sequenced a 239 bp fragment of Mhc Class II β-chain (BLB) exon 2 in 14 individuals. There was a total of 10 unique sequences and a minimum of four BLB loci. The d _N/d _S ratio at peptide-binding codons was significantly greater than one, suggesting balancing selection is acting on the BLB. We also recovered two YLB sequences, which clustered tightly with YLB sequences from three other species: domestic chicken, black grouse and common quail. The relatively large number of loci revealed in our study suggests that even closely related galliforms differ in the level of Mhc variation and structure.     

Structure,diversity and evolutionary patterns of expressed MHC class II<Emphasis Type="Italic">B</Emphasis> genes in chub (<Emphasis Type="Italic">Squalius cephalus</Emphasis>), a cyprinid fish species from Europe

The polymorphism of exon 2 of the DAB genes (major histocompatibility complex [MHC] class IIB) was investigated for the first time in the freshwater cyprinid fish species, Squalius cephalus, in the wide range of its distribution in Europe. We identified 111 different MHC class IIB variants in 15 chub populations distributed from Finland to Spain. The sequence analysis showed that many structurally important amino acid sites that were conserved among tetrapods were also conserved in chub. The analysis of recombination indicated that it does not play an important role in producing and maintaining the variation of DAB genes analyzed in the present study. The exon 2 was shown to be subjected to intense positive selection. Phylogenetic analysis and sequence identities suggest the presence of two class IIB loci (DAB1-like and DAB3-like) in chub. Nevertheless, the presence of three DAB3-like sequence variants in several individuals indicates the duplication of the DAB3 gene. A contrasting selection pattern was found in DAB1-like and DAB3-like genes, which suggests the potential functional differences between these genes. Some DAB sequence variants were shared among the populations of different mtDNA lineages. The phylogenetic analyses did not confirm any biogeographical pattern of the genetic structure of MHC IIB in chub, which is in line with balancing selection and trans-species polymorphism in MHC genes. Nevertheless, cluster analysis based on the presence/absence of DAB sequence variants in the populations showed the phylogeophraphical pattern corresponding to the mtDNA lineages, which indicates that neutral selection can partially explain the MHC IIB evolution in chub.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Evidence for Recombination in the Microcystin Synthetase (<Emphasis Type="Italic">mcy</Emphasis>) Genes ofToxic Cyanobacteria <Emphasis Type="Italic">Microcystis</Emphasis><Emphasis Type="Bold">spp</Emphasis>

Recombination has been suggested to be an important factor for the genetic variation of bacterial genes, but few studies have dealt with intragenic recombination between the same or closely related species of cyanobacteria. Here we provide strong evidence for recombination in the microcystin synthetase (mcy) gene cluster of the toxic cyanobacteria Microcystis spp. This gene cluster contains 10 genes (mcyA to J) that encode a mixed polyketide synthase (PKS)/nonribosomal peptide synthetase (NRPS) complex. mcy gene sequences were determined for four selected regions (within mcyA, D, G, and J) within the mcy gene cluster from 1 Canadian and 10 Asian toxic Microcystis and compared with previously published mcy sequences. Split decomposition analysis indicated a reticulate phylogeny of mcyA, and several potential recombination tracts of mcyA were identified by the RDP analysis and a runs test implemented in GENECONV. In contrast, no recombination was detected in the mcyD, G, and J sequences. However, discrepancies among the four mcy gene genealogies were evident from the results of independent split decomposition analyses, which were further supported by incongruence length difference (ILD) tests. Taken together, these findings suggest that both intragenic and intergenic recombination within the mcy gene cluster contributes to the genetic diversity of the mcy genes of Microcystis spp.This article contains online supplementary material.     

Duplicate Gene Evolution Toward Multiple Fates at the <Emphasis Type="Italic">Drosophila melanogaster HIP</Emphasis>/<Emphasis Type="Italic">HIP-Replacement</Emphasis> Locus

Hsc/Hsp70-interacting protein (HIP) is a rapidly evolving Hsp70 cofactor. Analyses of multiple Drosophila species indicate that the HIP gene is duplicated only in D. melanogaster. The HIP region, in fact, contains seven distinctly evolving duplicated genes. The regional duplication occurred in two steps, fixed rapidly, and illustrates multiple modes of duplicate gene evolution. HIP and its duplicate HIP-R are adaptively evolving in a manner unique to the region: they exhibit elevated divergence from other drosophilids and low polymorphism within D. melanogaster. HIP and HIP-R are virtually identical, share polymorphisms, and are subject to gene conversion. In contrast, two other duplicate genes in the region, CG33221 and GP-CG32779, are pseudogenes, and the chimeric gene Crg1 is subject to balancing selection. HIP and HIP-R are evolving rapidly and adaptively; however, positive selection is not sufficient to explain the molecular evolution of the region as a whole. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Contrasting Roles of Interallelic Recombination at the Hla-a and Hla-B Loci

A statistical study of DNA sequences of alleles at the highly polymorphic class I MHC loci of humans, HLA-A and HLA-B, showed evidence of both large-scale recombination events (involving recombination of exons 1-2 of one allele with exons 3-8 of another) and small-scale recombination events (involving apparent exchange of short DNA segments). The latter events occurred disproportionately in the region of the gene encoding the antigen recognition site (ARS) of the class I molecule. Furthermore, they involved the ARS codons which are under the strongest selection favoring allelic diversity at the amino acid level. Thus, the frequency of recombinant alleles appears to have been increased by some form of balancing selection (such as overdominant selection) favoring heterozygosity in the ARS. These analyses also revealed a striking difference between the A and B loci. Recombination events appear to have occurred about twice as frequently at the B locus, and recombinants at the B locus were significantly more likely to affect polymorphic sites in the ARS. At the A locus, there are well-defined allelic lineages that have persisted since prior to the human-chimpanzee divergence; but at the B locus, there is no evidence for such long-lasting allelic lineages. Thus, relatively frequent interallelic recombination has apparently been a feature of the long-term evolution of the B locus but not of the A locus.     

Selection-Marker-Free Modification of the Murine β-Casein Gene Using a lox2722 Site

Gene targeting and site-specific recombination strategies allow the precise modification of the eukaryotic genome. Many of the recombination strategies currently used, however, will introduce a selection marker gene at the modified site. DNA sequences of prokaryotic origin like vector sequences, selection marker, and reporter genes have been shown to markedly influence the regulation of the modified genomic loci. In order to avoid the insertion of excess sequences, a biphasic recombination strategy involving homologous recombination and Cre-recombinase-mediated cassette exchange (RMCE) was devised and used to insert a foreign gene into the β-casein gene in murine embryonic stem cells. The incompatibility of the heterospecific lox sites used for the recombinase-mediated cassette exchange was found to be critical for the success of the strategy. The frequently used mutant site lox511, which differs from the natural loxP site by a single point mutation, proved unsuitable for this approach. A mutant lox site carrying two point mutations, however, was highly effective and 90% of the selected cell clones carried the desired modification. This biphasic recombination strategy allows for the efficient and precise modification of gene loci without the concomitant introduction of a selectable marker gene.     

KRN4 Controls Quantitative Variation in Maize Kernel Row Number

Kernel row number (KRN) is an important component of yield during the domestication and improvement of maize and controlled by quantitative trait loci (QTL). Here, we fine-mapped a major KRN QTL, KRN4, which can enhance grain productivity by increasing KRN per ear. We found that a ~3-Kb intergenic region about 60 Kb downstream from the SBP-box gene Unbranched3 (UB3) was responsible for quantitative variation in KRN by regulating the level of UB3 expression. Within the 3-Kb region, the 1.2-Kb Presence-Absence variant was found to be strongly associated with quantitative variation in KRN in diverse maize inbred lines, and our results suggest that this 1.2-Kb transposon-containing insertion is likely responsible for increased KRN. A previously identified A/G SNP (S35, also known as Ser220Asn) in UB3 was also found to be significantly associated with KRN in our association-mapping panel. Although no visible genetic effect of S35 alone could be detected in our linkage mapping population, it was found to genetically interact with the 1.2-Kb PAV to modulate KRN. The KRN4 was under strong selection during maize domestication and the favorable allele for the 1.2-Kb PAV and S35 has been significantly enriched in modern maize improvement process. The favorable haplotype (Hap1) of 1.2-Kb-PAV-S35 was selected during temperate maize improvement, but is still rare in tropical and subtropical maize germplasm. The dissection of the KRN4 locus improves our understanding of the genetic basis of quantitative variation in complex traits in maize.     

The genealogy of sequences containing multiple sites subject to strong selection in a subdivided population

A stochastic model for the genealogy of a sample of recombining sequences containing one or more sites subject to selection in a subdivided population is described. Selection is incorporated by dividing the population into allelic classes and then conditioning on the past sizes of these classes. The past allele frequencies at the selected sites are thus treated as parameters rather than as random variables. The purpose of the model is not to investigate the dynamics of selection, but to investigate effects of linkage to the selected sites on the genealogy of the surrounding chromosomal region. This approach is useful for modeling strong selection, when it is natural to parameterize the past allele frequencies at the selected sites. Several models of strong balancing selection are used as examples, and the effects on the pattern of neutral polymorphism in the chromosomal region are discussed. We focus in particular on the statistical power to detect balancing selection when it is present.     

Positive Selection in Prion Protein

The prion diseases, such as Creutzfeldt-Jakob disease of humans and bovine spongiform encephalopathy, involve the aberrant metabolism and accumulation of prion protein PrP. There are three contradictory hypotheses about evolution of prion protein gene PRNP. Population genetic studies have proposed that PRNP could be under balancing selection, strong purifying selection, or mainly positive selection. We made use of the maximum likelihood tests for detection of positive selection at the amino acid level and present availability of PRNP coding sequences to contribute to these disagreements. Positive selection could occur at amino acids residing in active sites, and at amino acids involved in protein-protein interactions. Thus we tested a hypothesis that positive selection at the amino acid level in PrP might have taken place in human and related species from the superordinal group Euarchonta, as well as in bovine and related species from the superordinal clade Laurasiatheria. Our study and the present experimental evidences indicate that positive selection at the amino acid level might have taken place in the PrP signal sequences and conformationally plastic PrP regions, as well as at the protein X binding sites. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Prof. Vera Gamulin passed away.     

Evolution of <Emphasis Type="Italic">cd59</Emphasis> gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calculated. Results of synonymous and nonsynonymous substitution rates revealed that cd59 experienced negative selection in mammals overall. Four sites experiencing positive selection were found by using “site-specific” model in PAML software. These sites were distributed on the molecular surface, of which 2 sites located in the key functional domain. Furthermore, “branch-site-specific” model detected 1 positive site in cd59a and cd59b lineages which underwent accelerated evolution caused by positive selection after gene duplication in mouse.     

Footprints of intragenic recombination at HLA loci

 To evaluate the effect of balancing selection and intragenic recombination (or gene conversion) at six individual HLA loci, synonymous nucleotide diversity in different exon groups is examined within (π_w) and between (π_b) allelic lineages that may be defined by either serological or DNA sequence differences. Both π values are high in exons which encode for the peptide binding region (PBR) and tend to decrease in other exons. The value of π_w is significantly smaller than that of π_b in any exon of any locus. However, even π_w is much greater than nucleotide diversity at non-HLA loci. These observations provide additional strong evidence for the operation of balancing selection in PBR-encoding exons and its indirect effects on polymorphism at linked neighboring regions. It appears that allelic lineages have generally evolved in isolation but the linkage relationships within and between exons are incomplete throughout the long evolutionary history. To quantify intragenic recombination and account for the large discrepancy between the HLA and non-HLA diversity, a population genetics model is analyzed with special reference to the evolution of modern humans. The analysis suggests that the recombination rate between two sites 1000 base pairs apart is about 10^–5 per generation and that the effective size of human populations (equivalent roughly to the number of breeding individuals in a randomly mating population) has dropped from 10⁵ to 10⁴ in most of the Quaternary. One possibility for this reduction is discussed. Received: 11 August 1997 / Revised: 8 October 1997     

Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts

Background

Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.

Results

Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.

Conclusion

Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1369-8) contains supplementary material, which is available to authorized users.     

Characterization of expressed class II MHC sequences in the banner-tailed kangaroo rat (<Emphasis Type="Italic">Dipodomys spectabilis</Emphasis>) reveals multiple <Emphasis Type="Italic">DRB</Emphasis> loci

Genes of the major histocompatibility complex (MHC) are exceptionally polymorphic due to the combined effects of natural and sexual selection. Most research in wild populations has focused on the second exon of a single class II locus (DRB), but complete gene sequences can provide an illuminating backdrop for studies of intragenic selection, recombination, and organization. To this end, we characterized class II loci in the banner-tailed kangaroo rat (Dipodomys spectabilis). Seven DRB-like sequences (provisionally named MhcDisp-DRB*01 through *07) were isolated from spleen cDNA and most likely comprise ≥5 loci; this multiformity is quite unlike the situation in muroid rodents such as Mus, Rattus, and Peromyscus. In silico translation revealed the presence of important structural residues for glycosylation sites, salt bonds, and CD4+ T-cell recognition. Amino-acid distances varied widely among the seven sequences (2–34%). Nuclear DNA sequences from the Disp-DRB*07 locus (∼10 kb) revealed a conventional exon/intron structure as well as a number of microsatellites and short interspersed nuclear elements (B4, Alu, and IDL-Geo subfamilies). Rates of nucleotide substitution at Disp-DRB*07 are similar in both exons and introns (π = 0.015 and 0.012, respectively), which suggests relaxed selection and may indicate that this locus is an expressed pseudogene. Finally, we performed BLASTn searches against Dipodomys ordii genomic sequences (unassembled reads) and find 90–97% nucleotide similarity between the two kangaroo rat species. Collectively, these data suggest that class II diversity in heteromyid rodents is based on polylocism and departs from the muroid architecture. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession numbers EU817477–EU817485.