Male mouse recombination maps for each autosome identified by chromosome painting

Linkage maps constructed from genetic analysis of gene order and crossover frequency provide few clues to the basis of genomewide distribution of meiotic recombination, such as chromosome structure, that influences meiotic recombination. To bridge this gap, we have generated the first cytological recombination map that identifies individual autosomes in the male mouse. We prepared meiotic chromosome (synaptonemal complex [SC]) spreads from 110 mouse spermatocytes, identified each autosome by multicolor fluorescence in situ hybridization of chromosome-specific DNA libraries, and mapped >2,000 sites of recombination along individual autosomes, using immunolocalization of MLH1, a mismatch repair protein that marks crossover sites. We show that SC length is strongly correlated with crossover frequency and distribution. Although the length of most SCs corresponds to that predicted from their mitotic chromosome length rank, several SCs are longer or shorter than expected, with corresponding increases and decreases in MLH1 frequency. Although all bivalents share certain general recombination features, such as few crossovers near the centromeres and a high rate of distal recombination, individual bivalents have unique patterns of crossover distribution along their length. In addition to SC length, other, as-yet-unidentified, factors influence crossover distribution leading to hot regions on individual chromosomes, with recombination frequencies as much as six times higher than average, as well as cold spots with no recombination. By reprobing the SC spreads with genetically mapped BACs, we demonstrate a robust strategy for integrating genetic linkage and physical contig maps with mitotic and meiotic chromosome structure.     

Using multiple sequence correlation analysis to characterize functionally important protein regions

Protein co-evolution under structural and functional constraints necessitates the preservation of important interactions. Identifying functionally important regions poses many obstacles in protein engineering efforts. In this paper, we present a bioinformatics-inspired approach (residue correlation analysis, RCA) for predicting functionally important domains from protein family sequence data. RCA is comprised of two major steps: (i) identifying pairs of residue positions that mutate in a coordinated manner, and (ii) using these results to identify protein regions that interact with an uncommonly high number of other residues. We hypothesize that strongly correlated pairs result not only from contacting pairs, but also from residues that participate in conformational changes involved during catalysis or important interactions necessary for retaining functionality. The results show that highly mobile loops that assist in ligand association/dissociation tend to exhibit high correlation. RCA results exhibit good agreement with the findings of experimental and molecular dynamics studies for the three protein families that are analyzed: (i) DHFR (dihydrofolate reductase), (ii) cyclophilin, and (iii) formyl-transferase. Specifically, the specificity (percentage of correct predictions) in all three cases is substantially higher than those obtained by entropic measures or contacting residue pairs. In addition, we use our approach in a predictive fashion to identify important regions of a transmembrane amino acid transporter protein for which there is limited structural and functional information available.     

Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches

Based on integrative computational hybrid approaches that combined statistical coupling analysis (SCA), molecular dynamics (MD), and normal mode analysis (NMA), evolutionarily coupled residues involved in functionally relevant motion in the adenylate kinase protein family were identified. The hybrids identified four top-ranking site pairs that belong to a conserved hydrogen bond network that is involved in the enzyme's flexibility. A second group of top-ranking site pairs was identified in critical regions for functional dynamics, such as those related to enzymatic turnover. The high consistency of the results obtained by SCA with NMA (SCA.NMA) and by SCA.MD hybrid analyses suggests that suitable replacement of the matrix of cross-correlation analysis of atomic fluctuations (derived by using NMA) with those based on MD contributes to the identification of such sites by means of a fast computational calculation. The analysis presented here strongly supports the hypothesis that evolutionary forces, such as coevolution at the sequence level, have promoted functional dynamic properties of the adenylate kinase protein family. Finally, these hybrid approaches can be used to identify, at the residue level, protein motion coordination patterns not previously observed, such as in hinge regions.     

Development of chimeric laccases by directed evolution

DNA recombination methods are useful tools to generate diversity in directed evolution protein engineering studies. We have designed an array of chimeric laccases with high‐redox potential by in vitro and in vivo DNA recombination of two fungal laccases (from Pycnoporus cinnabarinus and PM1 basidiomycete), which were previously tailored by laboratory evolution for functional expression in Saccharomyces cerevisiae. The laccase fusion genes (including the evolved α‐factor prepro‐leaders for secretion in yeast) were subjected to a round of family shuffling to construct chimeric libraries and the best laccase hybrids were identified in dual high‐throughput screening (HTS) assays. Using this approach, we identified chimeras with up to six crossover events in the whole sequence, and we obtained active hybrid laccases with combined characteristics in terms of pH activity and thermostability. Biotechnol. Bioeng. 2012; 109: 2978–2986. © 2012 Wiley Periodicals, Inc.     

Total degradation of pentachloroethane by an engineered Alcaligenes strain expressing a modified camphor monooxygenase and a hybrid dioxygenase

We engineered biphenyl-degrading Alcaligenes sp. strain KF711 for total degradation of pentachloroethane (PCA), which expresses a modified camphor monooxygenase and a hybrid dioxygenase consisting of TodC1 (a large subunit of toluene dioxygenase of Pseudomonas putida F1) and BphA2-BphA3-pbhA4 (a small subunit, ferredoxin and ferredoxin reductase of biphenyl dioxygenase, respectively, in strain KF707). Modified camphor monooxygenase genes (camCAB) were supplied as a plasmid and the todC1 gene was integrated within the chromosomal bph gene cluster by a single crossover recombination. The resultant strain KF711S-3cam dechlorinated PCA to trichloroethene by the action of the modified camphor monooxygenase under anaerobic conditions. The same strain subsequently degraded trichloroethene formed oxidatively by the action of the Tol-Bph hybrid dioxygenase under aerobic conditions. Thus sequential anaerobic and aerobic treatments of the KF711S-3cam resting cells resulted in efficient and total degradation of PCA.     

Conserved patterns of chromosome pairing and recombination in Brassica napus crosses.

The patterns of chromosome pairing and recombination in two contrasting Brassica napus F1 hybrids were deduced. One hybrid was from a winter oilseed rape (WOSR) x spring oilseed rape cross, the other from a resynthesized B. napus x WOSR cross. Segregation at 211 equivalent loci assayed in the population derived from each hybrid produced two collinear genetic maps. Alignment of the maps indicated that B. napus chromosomes behaved reproducibly as 19 homologous pairs and that the 19 distinct chromosomes of B. napus each recombined with unique chromosomes from the interspecific hybrid between Brassica rapa and Brassica oleracea. This result indicated that the genomes of the diploid progenitors of amphidiploid B. napus have remained essentially unaltered since the formation of the species and that the progenitor genomes were similar to those of modern-day B. rapa and B. oleracea. The frequency and distribution of crossovers were almost indistinguishable in the two populations, suggesting that the recombination machinery of B. napus could cope easily with different degrees of genetic divergence between homologous chromosomes. Efficient recombination in wide crosses will facilitate the introgression of novel alleles into oilseed rape from B. rapa and B. oleracea (via resynthesized B. napus) and reduce linkage drag.     

Prediction of protein residue contacts with a PDB-derived likelihood matrix

Proteins with similar folds often display common patterns of residue variability. A widely discussed question is how these patterns can be identified and deconvoluted to predict protein structure. In this respect, correlated mutation analysis (CMA) has shown considerable promise. CMA compares multiple members of a protein family and detects residues that remain constant or mutate in tandem. Often this behavior points to structural or functional interdependence between residues. CMA has been used to predict pairs of amino acids that are distant in the primary sequence but likely to form close contacts in the native three-dimensional structure. Until now these methods have used evolutionary or biophysical models to score the fit between residues. We wished to test whether empirical methods, derived from known protein structures, would provide useful predictive power for CMA. We analyzed 672 known protein structures, derived contact likelihood scores for all possible amino acid pairs, and used these scores to predict contacts. We then tested the method on 118 different protein families for which structures have been solved to atomic resolution. The mean performance was almost seven times better than random prediction. Used in concert with secondary structure prediction, the new CMA method could supply restraints for predicting still undetermined structures.     

Analysis of shuffled gene libraries

In vitro recombination of homologous genes (family shuffling) has been proposed as an effective search strategy for laboratory evolution of genes and proteins. Few data are available, however, on the composition of shuffled gene libraries, from which one could assess the efficiency of recombination and optimize protocols. Here, probe hybridization is used in a macroarray format to analyze chimeric DNA libraries created by DNA shuffling. Characterization of hundreds of shuffled genes encoding dioxygenases has elucidated important biases in the shuffling reaction. As expected, crossovers are favored in regions of high sequence identity. A sequence-based model of homologous recombination that captures this observed bias was formulated using the experimental results. The chimeric genes were found to show biases in the incorporation of sequences from certain parents, even before selection. Statistically different patterns of parental incorporation in genes expressing functional proteins can help to identify key sequence-function relationships.     

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.     

Relationships between rust resistance genes at the M locus in flax

Genes at the M locus in flax ( Linum usitatissimum ) that confer resistance to flax rust ( Melampsora lini ) occur in complex haplotypes containing up to 15 related genes or gene fragments. We have cloned two additional functional resistance genes at this locus, M1 and M3 , by transposon tagging and candidate gene approaches, and investigated the genetic relationships between four genes ( M , M1 , M3 and M4 ) by recombination analysis. M1 and M3 , like M , are members of the nucleotide binding site, leucine-rich repeat (NBS-LRR) family. Comparisons of the predicted M1 and M3 amino acid sequences with M and L6 reveal that: (i) M1 contains four additional LRRs, probably as a result of an unequal crossover event between duplicated regions; (ii) M1 shares large segments of exact identity with M and M3, indicative of intragenic recombination events; and (iii) a large number of amino acid differences are scattered throughout the M, M1 and M3 proteins. Recombination analysis (here and in previous studies) has revealed that M readily recombines with M1 , M3 and M4 , whereas these three genes fail to recombine despite large family sizes (>5800) in two test-cross families, suggesting that they may occupy allelic positions in the gene cluster. Several restriction fragment length polymorphism markers within or near the M locus were mapped with respect to seven crossover events between M and M1 . The results of this and previous studies provide evidence of structural differences between: (i) homoeologous loci in the different genomes of flax; (ii) different haplotypes at the M locus; (iii) different resistance genes in the M group; and (iv) the flanking regions downstream of M locus resistance genes.     

Comparative analysis of sequence covariation methods to mine evolutionary hubs: Examples from selected GPCR families

Covariation between positions in a multiple sequence alignment may reflect structural, functional, and/or phylogenetic constraints and can be analyzed by a wide variety of methods. We explored several of these methods for their ability to identify covarying positions related to the divergence of a protein family at different hierarchical levels. Specifically, we compared seven methods on a model system composed of three nested sets of G‐protein‐coupled receptors (GPCRs) in which a divergence event occurred. The covariation methods analyzed were based on: χ² test, mutual information, substitution matrices, and perturbation methods. We first analyzed the dependence of the covariation scores on residue conservation (measured by sequence entropy), and then we analyzed the networking structure of the top pairs. Two methods out of seven—OMES (Observed minus Expected Squared) and ELSC (Explicit Likelihood of Subset Covariation)—favored pairs with intermediate entropy and a networking structure with a central residue involved in several high‐scoring pairs. This networking structure was observed for the three sequence sets. In each case, the central residue corresponded to a residue known to be crucial for the evolution of the GPCR family and the subfamily specificity. These central residues can be viewed as evolutionary hubs, in relation with an epistasis‐based mechanism of functional divergence within a protein family. Proteins 2014; 82:2141–2156. © 2014 Wiley Periodicals, Inc.     

Hybrid genes between HLA-A2 and HLA-A3 constructed by in vivo recombination allow mapping of HLA-A2 and HLA-A3 polymorphic antigenic determinants

HLA-A2 and -A3 genes have been modified in their third exon (second domain) by using in vivo recombination. In this method Escherichia coli are transfected with a plasmid which contains two highly homologous sequences (e.g., the third exons of HLA-A2 and -A3) and has been linearized by cleavage between these two sequences. Circularization takes place in the bacteria by homologous recombination leading to hybrid A2-A3 sequences. The analysis by DNA sequencing of a number of such recombinants shows that they indeed occur by homologous recombination (no insertions or deletions) and that the probability of crossing over decreases as the distance from the free end of DNA in the homologous region increases. No double recombinants were observed. These hybrid exons were reinserted into either HLA-A2 or HLA-A3 genes, thus generating a panel of functional hybrid genes containing one or several HLA-A2 specific substitutions in an HLA-A3 background or vice versa. These genes were expressed by transfection into murine P815-high transfection efficiency recipient cells. Serologic analysis leads to the conclusion that expression of polymorphic antigenic determinants specific for HLA-A2 (detected with M58, A2A28M1, and CR11.351 mAb) is linked to the presence of threonine residue (amino acid (AA) 142) and/or histidine residue (AA 145) and valine residue (AA 152). The expression of specific HLA-A3 polymorphic determinants (recognized by GAP-A3 mAb) is correlated with the existence of a asparagine residue (AA 127) and a aspartic residue (AA 161). But aspartic residue 161 contributes with glutamic acid residue 152 in the formation of the A3 epitope recognized by the anti-A3 mAb X1.23.2.     

Overlapping alternative donor splice sites in the human genome

Over 50% of donor splice sites in the human genome have a potential alternative donor site at a distance of three to six nucleotides. Conservation of these potential sites is determined by the consensus requirements and by its exonic or intronic location. Several hundred pairs of overlapping sites are confirmed to be alternatively spliced as both sites in a pair are supported by a protein, by a full-length mRNA, or by expressed sequence tags (ESTs) from at least two independent clone libraries. Overlapping sites may clash with consensus requirements. Pairs with a site shift of four nucleotides are the most abundant, despite the frameshift in the protein-coding region that they introduce. The site usage in pairs is usually uneven, and the major site is more frequently conserved in other mammalian genomes. Overlapping alternative donor sites and acceptor sites may have different functional roles: alternative splicing of overlapping acceptor sites leads mainly to microvariations in protein sequences; whereas alternative donor sites often lead to frameshifts and thus either yield major differences in the protein sequence and structure, or generate nonsense-mediated decay-inducing mRNA isoforms likely involved in regulated unproductive splicing pathways.     

SIDEpro: a novel machine learning approach for the fast and accurate prediction of side-chain conformations

Accurate protein side-chain conformation prediction is crucial for protein modeling and existing methods for the task are widely used; however, faster and more accurate methods are still required. Here we present a new machine learning approach to the problem where an energy function for each rotamer in a structure is computed additively over pairs of contacting atoms. A family of 156 neural networks indexed by amino acid and contacting atom types is used to compute these rotamer energies as a function of atomic contact distances. Although direct energy targets are not available for training, the neural networks can still be optimized by converting the energies to probabilities and optimizing these probabilities using Markov Chain Monte Carlo methods. The resulting predictor SIDEpro makes predictions by initially setting the rotamer probabilities for each residue from a backbone-dependent rotamer library, then iteratively updating these probabilities using the trained neural networks. After convergences of the probabilities, the side-chains are set to the highest probability rotamer. Finally, a post processing clash reduction step is applied to the models. SIDEpro represents a significant improvement in speed and a modest, but statistically significant, improvement in accuracy when compared with the state-of-the-art for rapid side-chain prediction method SCWRL4 on the following datasets: (1) 379 protein test set of SCWRL4; (2) 94 proteins from CASP9; (3) a set of seven large protein-only complexes; and (4) a ribosome with and without the RNA. Using the SCWRL4 test set, SIDEpro's accuracy (χ(1) 86.14%, χ(1+2) 74.15%) is slightly better than SCWRL4-FRM (χ(1) 85.43%, χ(1+2) 73.47%) and it is 7.0 times faster. On the same test set SIDEpro is clearly more accurate than SCWRL4-rigid rotamer model (RRM) (χ(1) 84.15%, χ(1+2) 71.24%) and 2.4 times faster. Evaluation on the additional test sets yield similar accuracy results with SIDEpro being slightly more accurate than SCWRL4-flexible rotamer model (FRM) and clearly more accurate than SCWRL4-RRM; however, the gap in CPU time is much more significant when the methods are applied to large protein complexes. SIDEpro is part of the SCRATCH suite of predictors and available from: http://scratch.proteomics.ics.uci.edu/.     

Random Field Model Reveals Structure of the Protein Recombinational Landscape

We are interested in how intragenic recombination contributes to the evolution of proteins and how this mechanism complements and enhances the diversity generated by random mutation. Experiments have revealed that proteins are highly tolerant to recombination with homologous sequences (mutation by recombination is conservative); more surprisingly, they have also shown that homologous sequence fragments make largely additive contributions to biophysical properties such as stability. Here, we develop a random field model to describe the statistical features of the subset of protein space accessible by recombination, which we refer to as the recombinational landscape. This model shows quantitative agreement with experimental results compiled from eight libraries of proteins that were generated by recombining gene fragments from homologous proteins. The model reveals a recombinational landscape that is highly enriched in functional sequences, with properties dominated by a large-scale additive structure. It also quantifies the relative contributions of parent sequence identity, crossover locations, and protein fold to the tolerance of proteins to recombination. Intragenic recombination explores a unique subset of sequence space that promotes rapid molecular diversification and functional adaptation.     

Development and Linkage Analysis of 104 New Microsatellite Markers for Bay Scallop (Argopecten irradians)

For genetic analysis and linkage mapping of bay scallop (Argopecten irradians), a set of 120 novel simple sequence repeat markers were developed from microsatellite-enriched libraries and expressed sequence tags. An inter-subspecies hybrid bay scallop family (CC5) of 46 progeny was analyzed as the reference population to confirm polymorphism and test the segregation patterns of these loci. A total of 104 microsatellite markers were polymorphic in the reference family, among which 36 in female, 28 in male, and 40 in both parents, respectively. Linkage analysis allowed mapping these markers to 15 linkage groups, which is close to the haploid chromosome number of bay scallop (n = 16). Analysis of the 40 markers segregating in both parents showed a higher recombination rate in the female parent, with the average of female-to-male recombination ratio of 1.09:1 between linked pairs of markers. When null alleles were considered, there were 17 loci showing segregation distortion at the 5% significance level using the chi-square test. The microsatellite markers developed in this study provide a useful resource for future linkage mapping and quantitative loci analysis in A. irradians.     

The beta recombinase of plasmid pSM19035 binds to two adjacent sites, making different contacts at each of them.

The beta recombinase from plasmid pSM19035 catalyzes intramolecular site-specific recombination between two directly or inversely oriented six sites in the presence of a chromatin-associated protein (Hbsu, HU or HMG-1). The six site is a DNA segment containing two binding sites (I and II) for beta protein dimers. We show that beta recombinase binds sequentially to both sites, having a different affinity for each one. Hydroxyl radical footprints show a different protection pattern at each site. Positions critical for beta protein binding have been identified by methylation interference and missing nucleoside assays. The results indicate that the protein recognizes each site in a different way. Comparison of the beta protein recombination site with that of DNA resolvases and DNA invertases of the Tn3 family, to which it belongs, shows that these sequences can be divided into two regions. One corresponds to the crossover point and is similar for all recombinases of the family. The other region differs in the different subfamilies and seems to have an architectural role in aligning the crossover sites at the synaptic complex. The different ways to assemble this complex could explain why each system leads to a particular recombination event: DNA resolution (resolvases), inversion (invertases) or both (beta recombinase).     

Genetic and physical mapping of two centromere-proximal regions of chromosome IV in Aspergillus nidulans

Chromosome IV is the smallest chromosome of Aspergillus nidulans. The centromere-proximal portion of the chromosome was mapped physically using overlapping clones of a cosmid genomic library. Two contiguous segments of a physical map, based on restriction mapping of cosmid clones, were generated, together covering more than 0.4 Mb DNA. A reverse genetic mapping approach was used to establish a correlation between physical and genetic maps; i.e., marker genes were integrated into physically mapped segments and subsequently mapped by mitotic and meiotic recombination. The resulting data, together with additional classical genetic mapping, lead to a substantial revision of the genetic map of the chromosome, including the position of the centromere. Comparison of physical and genetic maps indicates that meiotic recombination is low in subcentromeric DNA, its frequency being reduced from 1 crossover per 0.8 Mb to approximately 1 crossover per 5 Mb per meiosis. The portion of the chromosome containing the functional centromere was not mapped because repeat-rich regions hindered further chromosome walking. The size of the missing segment was estimated to be between 70 and 400 kb.     

The ssDNA-binding protein MEIOB acts as a dosage-sensitive regulator of meiotic recombination

Meiotic recombination enables reciprocal exchange of genetic information between parental chromosomes and is essential for fertility. MEIOB, a meiosis-specific ssDNA-binding protein, regulates early meiotic recombination. Here we report that the human infertility-associated missense mutation (N64I) in MEIOB causes protein degradation and reduced crossover formation in mouse testes. Although the MEIOB N64I substitution is associated with human infertility, the point mutant mice are fertile despite meiotic defects. Meiob mutagenesis identifies serine 67 as a critical residue for MEIOB. Biochemically, these two mutations (N64I and S67 deletion) cause self-aggregation of MEIOB and sharply reduced protein half-life. Molecular genetic analyses of both point mutants reveal an important role for MEIOB in crossover formation in late meiotic recombination. Furthermore, we find that the MEIOB protein levels directly correlate with the severity of meiotic defects. Our results demonstrate that MEIOB regulates meiotic recombination in a dosage-dependent manner.     

Protein building blocks preserved by recombination

Borrowing concepts from the schema theory of genetic algorithms, we have developed a computational algorithm to identify the fragments of proteins, or schemas, that can be recombined without disturbing the integrity of the three-dimensional structure. When recombination leaves these schemas undisturbed, the hybrid proteins are more likely to be folded and functional. Crossovers found by screening libraries of several randomly shuffled proteins for functional hybrids strongly correlate with those predicted by this approach. Experimental results from the construction of hybrids of two beta-lactamases that share 40% amino acid identity demonstrate a threshold in the amount of schema disruption that the hybrid protein can tolerate. To the extent that introns function to promote recombination within proteins, natural selection would serve to bias their locations to schema boundaries.