Genetic linkage map of the edible basidiomycete Pleurotus ostreatus

We have constructed a genetic linkage map of the edible basidiomycete Pleurotus ostreatus (var. Florida). The map is based on the segregation of 178 random amplified polymorphic DNA and 23 restriction fragment length polymorphism markers; four hydrophobin, two laccase, and two manganese peroxidase genes; both mating type loci; one isozyme locus (est1); the rRNA gene sequence; and a repetitive DNA sequence in a population of 80 sibling monokaryons. The map identifies 11 linkage groups corresponding to the chromosomes of P. ostreatus, and it has a total length of 1,000.7 centimorgans (cM) with an average of 35.1 kbp/cM. The map shows a high correlation (0.76) between physical and genetic chromosome sizes. The number of crossovers observed per chromosome per individual cell is 0.89. This map covers nearly the whole genome of P. ostreatus.     

The longest (A+T) and (G+C) blocks in the human and other genomes

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are imbedded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

Genetic Linkage Map of the Edible Basidiomycete Pleurotus ostreatus

We have constructed a genetic linkage map of the edible basidiomycete Pleurotus ostreatus (var. Florida). The map is based on the segregation of 178 random amplified polymorphic DNA and 23 restriction fragment length polymorphism markers; four hydrophobin, two laccase, and two manganese peroxidase genes; both mating type loci; one isozyme locus (est1); the rRNA gene sequence; and a repetitive DNA sequence in a population of 80 sibling monokaryons. The map identifies 11 linkage groups corresponding to the chromosomes of P. ostreatus, and it has a total length of 1,000.7 centimorgans (cM) with an average of 35.1 kbp/cM. The map shows a high correlation (0.76) between physical and genetic chromosome sizes. The number of crossovers observed per chromosome per individual cell is 0.89. This map covers nearly the whole genome of P. ostreatus.     

The Longest (A+T) and (G+C) Blocks in the Human and Other Genomes

Abstract

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16^th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are embeeded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

Efficient reconstruction of haplotype structure via perfect phylogeny

Each person's genome contains two copies of each chromosome, one inherited from the father and the other from the mother. A person's genotype specifies the pair of bases at each site, but does not specify which base occurs on which chromosome. The sequence of each chromosome separately is called a haplotype. The determination of the haplotypes within a population is essential for understanding genetic variation and the inheritance of complex diseases. The haplotype mapping project, a successor to the human genome project, seeks to determine the common haplotypes in the human population. Since experimental determination of a person's genotype is less expensive than determining its component haplotypes, algorithms are required for computing haplotypes from genotypes. Two observations aid in this process: first, the human genome contains short blocks within which only a few different haplotypes occur; second, as suggested by Gusfield, it is reasonable to assume that the haplotypes observed within a block have evolved according to a perfect phylogeny, in which at most one mutation event has occurred at any site, and no recombination occurred at the given region. We present a simple and efficient polynomial-time algorithm for inferring haplotypes from the genotypes of a set of individuals assuming a perfect phylogeny. Using a reduction to 2-SAT we extend this algorithm to handle constraints that apply when we have genotypes from both parents and child. We also present a hardness result for the problem of removing the minimum number of individuals from a population to ensure that the genotypes of the remaining individuals are consistent with a perfect phylogeny. Our algorithms have been tested on real data and give biologically meaningful results. Our webserver (http://www.cs.columbia.edu/compbio/hap/) is publicly available for predicting haplotypes from genotype data and partitioning genotype data into blocks.     

Predicting recombination frequency from map distance

Map distance is one of the key measures in genetics and indicates the expected number of crossovers between two loci. Map distance is estimated from the observed recombination frequency using mapping functions, the most widely used of those, Haldane and Kosambi, being developed at the time when the number of markers was low and unobserved crossovers had a substantial effect on the recombination fractions. In contemporary high-density marker data, the probability of multiple crossovers between adjacent loci is negligible and different mapping functions yield the same result, that is, the recombination frequency between adjacent loci is equal to the map distance in Morgans. However, high-density linkage maps contain an interpretation problem: the map distance over a long interval is additive and its association with recombination frequency is not defined. Here, we demonstrate with high-density linkage maps from humans and stickleback fishes that the inverses of Haldane’s and Kosambi’s mapping functions systematically underpredict recombination frequencies from map distance. To remedy this, we formulate a piecewise function that yields more accurate predictions of recombination frequency from map distance. Our results demonstrate that the association between map distance and recombination frequency is context-dependent and without a universal solution.Subject terms: Genetic linkage study, Evolutionary genetics     

Comparison of a high-density genetic linkage map to genome features in the model grass <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F₂ mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.     

Constructing a dense genetic linkage map and mapping QTL for the traits of flower development in Brassica carinata

Key message

An integrated dense genetic linkage map was constructed in a B. carinata population and used for comparative genome analysis and QTL identification for flowering time.

Abstract

An integrated dense linkage map of Brassica carinata (BBCC) was constructed in a doubled haploid population based on DArT-Seq^TM markers. A total of 4,031 markers corresponding to 1,366 unique loci were mapped including 639 bins, covering a genetic distance of 2,048 cM. We identified 136 blocks and islands conserved in Brassicaceae, which showed a feature of hexaploidisation representing the suggested ancestral crucifer karyotype. The B and C genome of B. carinata shared 85 % of commonly conserved blocks with the B genome of B. nigra/B. juncea and 80 % of commonly conserved blocks with the C genome of B. napus, and shown frequent structural rearrangements such as insertions and inversions. Up to 24 quantitative trait loci (QTL) for flowering and budding time were identified in the DH population. Of these QTL, one consistent QTL (qFT.B4-2) for flowering time was identified in all of the environments in the J block of the B4 linkage group, where a group of genes for flowering time were aligned in A. thaliana. Another major QTL for flowering time under a winter-cropped environment was detected in the E block of C6, where the BnFT-C6 gene was previously localised in B. napus. This high-density map would be useful not only to reveal the genetic variation in the species with QTL analysis and genome sequencing, but also for other applications such as marker-assisted selection and genomic selection, for the African mustard improvement.     

Synorth: exploring the evolution of synteny and long-range regulatory interactions in vertebrate genomes

Genomic regulatory blocks are chromosomal regions spanned by long clusters of highly conserved noncoding elements devoted to long-range regulation of developmental genes, often immobilizing other, unrelated genes into long-lasting syntenic arrangements. Synorth is a web resource for exploring and categorizing the syntenic relationships in genomic regulatory blocks across multiple genomes, tracing their evolutionary fate after teleost whole genome duplication at the level of genomic regulatory block loci, individual genes, and their phylogenetic context.     

Genetic analysis and genome mapping in Raphanus.

The first genetic map of the Raphanus genome was developed based on meiosis in a hybrid between Raphanus sativus (cultivated radish) and Raphanus raphanistrum (wild radish). This hybrid was used to produce a BC1 population of 54 individuals and an F2 population of 85 individuals. A total of 236 marker loci were assayed in these populations using a set of 144 informative Brassica RFLP probes previously used for genetic mapping in other crucifer species. The genetic maps derived from the BC1 and F2 populations were perfectly collinear and were integrated to produce a robust Raphanus map. Cytological observations demonstrated strict bivalent pairing in the R. sativus x R. raphanistrum hybrids. Productive pairing along the length of each chromosome was confirmed by the identification of nine extensive linkage groups and the lack of clustering of marker loci. Indeed, the distributions of both marker loci and crossovers was more random than those reported for other crop species. The genetic markers and the reference map of Raphanus will be of considerable value for future trait mapping and marker-assisted breeding in this crop, as well as in the intergenomic transfer of Raphanus genes into Brassica crops. The future benefits of comparative mapping with Arabidopsis and Brassica species are also discussed.     

A highly recombined,high‐density,eight‐founder wheat MAGIC map reveals extensive segregation distortion and genomic locations of introgression segments

Multiparent Advanced Generation Intercross (MAGIC) mapping populations offer unique opportunities and challenges for marker and QTL mapping in crop species. We have constructed the first eight‐parent MAGIC genetic map for wheat, comprising 18 601 SNP markers. We validated the accuracy of our map against the wheat genome sequence and found an improvement in accuracy compared to published genetic maps. Our map shows a notable increase in precision resulting from the three generations of intercrossing required to create the population. This is most pronounced in the pericentromeric regions of the chromosomes. Sixteen percent of mapped markers exhibited segregation distortion (SD) with many occurring in long (>20 cM) blocks. Some of the longest and most distorted blocks were collinear with noncentromeric high‐marker‐density regions of the genome, suggesting they were candidates for introgression fragments introduced into the bread wheat gene pool from other grass species. We investigated two of these linkage blocks in detail and found strong evidence that one on chromosome 4AL, showing SD against the founder Robigus, is an interspecific introgression fragment. The completed map is available from http://www.niab.com/pages/id/326/Resources .     

Recombinants between herpes simplex virus types 1 and 2: analyses of genome structures and expression of immediate early polypeptides.

Recombinants between temperature-sensitive mutants of herpes simplex virus types 1 (HSV-1) and 2 (HSV-2) were constructed. Using restriction endonucleases, we analyzed the genome composition of 17 intertypic recombinants and detected crossovers in every region of the genome. The virion DNA of one recombinant appeared to be largely "frozen" in two of the four possible genome arrangements of HSV. Knowledge of the genome structures of recombinants enabled us to physically map immediate early polypeptides. We present evidence that the immediate early polypeptide Vmw IE 110 of HSV-1 and its functionally equivalent polypeptide, Vmw IE 118, of HSV-2 may map in the repetitive sequences bounding the long unique region of HSV.     

Haplotype block structure is conserved across mammals

Genetic variation in genomes is organized in haplotype blocks, and species-specific block structure is defined by differential contribution of population history effects in combination with mutation and recombination events. Haplotype maps characterize the common patterns of linkage disequilibrium in populations and have important applications in the design and interpretation of genetic experiments. Although evolutionary processes are known to drive the selection of individual polymorphisms, their effect on haplotype block structure dynamics has not been shown. Here, we present a high-resolution haplotype map for a 5-megabase genomic region in the rat and compare it with the orthologous human and mouse segments. Although the size and fine structure of haplotype blocks are species dependent, there is a significant interspecies overlap in structure and a tendency for blocks to encompass complete genes. Extending these findings to the complete human genome using haplotype map phase I data reveals that linkage disequilibrium values are significantly higher for equally spaced positions in genic regions, including promoters, as compared to intergenic regions, indicating that a selective mechanism exists to maintain combinations of alleles within potentially interacting coding and regulatory regions. Although this characteristic may complicate the identification of causal polymorphisms underlying phenotypic traits, conservation of haplotype structure may be employed for the identification and characterization of functionally important genomic regions.     

Using Markers in Gene Introgression Breeding Programs

We investigate the use of markers to hasten the recovery of the recipient genome during an introgression breeding program. The effects of time and intensity of selection, population size, number and position of selected markers are studied for chromosomes either carrying or not carrying the introgressed gene. We show that marker assisted selection may lead to a gain in time of about two generations, an efficiency below previous theoretical predictions. Markers are most useful when their map position is known. In the early generations, it is shown that increasing the number of markers over three per non-carrier chromosome is not efficient, that the segment surrounding the introgressed gene is better controlled by rather distant markers unless high selection intensity can be applied, and that selection on this segment first can reduce the selection intensity available for selection on non-carrier chromosomes. These results are used to propose an optimal strategy for selection on the whole genome, making the most of available material and conditions (e.g., population size and fertility, genetic map).     

The breakdown of genomic ancestry blocks in hybrid lineages given a finite number of recombination sites

When a lineage originates from hybridization genomic blocks of contiguous ancestry from different ancestors are fragmented through genetic recombination. The resulting blocks are delineated by so called junctions, which accumulate with every generation that passes. Modeling the accumulation of ancestry block junctions can elucidate processes and timeframes of genomic admixture. Previous models have not addressed ancestry block dynamics for chromosomes that consist of a finite number of recombination sites. However, genomic data typically consist of informative markers that are interspersed with fragments for which no ancestry information is available. Hence, repeated recombination events may occur between markers, effectively removing existing junctions. Here, we present an analytical treatment of the dynamics of the mean number of junctions over time, taking into account the number of recombination sites per chromosome, population size, genetic map length, and the frequency of the ancestral species in the founding hybrid swarm. We describe the expected number of junctions using equidistant molecular markers and estimate the number of junctions using random markers. This extended theory of junctions thus reflects properties of empirical data and can serve to study the genomic patterns following admixture.     

Large-scale heterospecific segregation distortion in <Emphasis Type="Italic">Populus</Emphasis> revealed by a dense genetic map

We report the most complete genetic map to have been constructed for the genus Populus. This map includes 544 markers mapped onto 19 linkage groups, equivalent to the Populus chromosome number, with all markers displaying internally consistent linkage patterns. We estimate the genome length to be between 2,300 and 2,500 cM, based both on the observed number of crossovers in the maternal haplotypes, as well as the total observed map length. Genome coverage was estimated to be greater than 99.9% at 20 cM per marker. We did not detect obvious recombination repression in the maternal tree (a hybrid of Populus trichocarpa Hooker × P. deltoides Marsh.) compared to the paternal tree (pure P. deltoides). Finally, most markers exhibiting segregation distortion were derived from the donor parent in this backcross, and generally occurred in large contiguous blocks on two linkage groups. We hypothesize that divergent selection has occurred on chromosomal scales among the parental species used to create this pedigree, and explore the evolutionary implications of this observation. This genetic linkage map provides the most comprehensive view of the Populus genome reported to date and will prove invaluable for future inquiries into the structural and functional genomics, evolutionary biology, and genetic improvement of this ecologically important model species.Electronic Supplementary Material Supplementary material is available for this article at Communicated by D.B. Neale     

A high‐resolution genetic map of the cereal crown rot pathogen Fusarium pseudograminearum provides a near‐complete genome assembly

Fusarium pseudograminearum is an important pathogen of wheat and barley, particularly in semi‐arid environments. Previous genome assemblies for this organism were based entirely on short read data and are highly fragmented. In this work, a genetic map of F. pseudograminearum has been constructed for the first time based on a mapping population of 178 individuals. The genetic map, together with long read scaffolding of a short read‐based genome assembly, was used to give a near‐complete assembly of the four F. pseudograminearum chromosomes. Large regions of synteny between F. pseudograminearum and F. graminearum, the related pathogen that is the primary causal agent of cereal head blight disease, were previously proposed in the core conserved genome, but the construction of a genetic map to order and orient contigs is critical to the validation of synteny and the placing of species‐specific regions. Indeed, our comparative analyses of the genomes of these two related pathogens suggest that rearrangements in the F. pseudograminearum genome have occurred in the chromosome ends. One of these rearrangements includes the transposition of an entire gene cluster involved in the detoxification of the benzoxazolinone (BOA) class of plant phytoalexins. This work provides an important genomic and genetic resource for F. pseudograminearum, which is less well characterized than F. graminearum. In addition, this study provides new insights into a better understanding of the sexual reproduction process in F. pseudograminearum, which informs us of the potential of this pathogen to evolve.     

COMPAM :visualization of combining pairwise alignments for multiple genomes

COMPAM is a tool for visualizing relationships among multiple whole genomes by combining all pairwise genome alignments. It displays shared conserved regions (blocks) and where these blocks occur (edges) as block relation graphs which can be explored interactively. An unannotated genome, e.g. can then be explored using information from well-annotated genomes, COG-based genome annotation and genes. COMPAM can run either as a stand-alone application or through an applet that is provided as service to PLATCOM, a toolset for whole genome comparative analysis, where a wide variety of genomes can be easily selected. Features provided by COMPAM include the ability to export genome relationship information into file formats that can be used by other existing tools. AVAILABILITY: http://bio.informatics.indiana.edu/projects/compam/     

Haplotyping as perfect phylogeny: a direct approach.

A full haplotype map of the human genome will prove extremely valuable as it will be used in large-scale screens of populations to associate specific haplotypes with specific complex genetic-influenced diseases. A haplotype map project has been announced by NIH. The biological key to that project is the surprising fact that some human genomic DNA can be partitioned into long blocks where genetic recombination has been rare, leading to strikingly fewer distinct haplotypes in the population than previously expected (Helmuth, 2001; Daly et al., 2001; Stephens et al., 2001; Friss et al., 2001). In this paper we explore the algorithmic implications of the no-recombination in long blocks observation, for the problem of inferring haplotypes in populations. This assumption, together with the standard population-genetic assumption of infinite sites, motivates a model of haplotype evolution where the haplotypes in a population are assumed to evolve along a coalescent, which as a rooted tree is a perfect phylogeny. We consider the following algorithmic problem, called the perfect phylogeny haplotyping problem (PPH), which was introduced by Gusfield (2002) - given n genotypes of length m each, does there exist a set of at most 2n haplotypes such that each genotype is generated by a pair of haplotypes from this set, and such that this set can be derived on a perfect phylogeny? The approach taken by Gusfield (2002) to solve this problem reduces it to established, deep results and algorithms from matroid and graph theory. Although that reduction is quite simple and the resulting algorithm nearly optimal in speed, taken as a whole that approach is quite involved, and in particular, challenging to program. Moreover, anyone wishing to fully establish, by reading existing literature, the correctness of the entire algorithm would need to read several deep and difficult papers in graph and matroid theory. However, as stated by Gusfield (2002), many simplifications are possible and the list of "future work" in Gusfield (2002) began with the task of developing a simpler, more direct, yet still efficient algorithm. This paper accomplishes that goal, for both the rooted and unrooted PPH problems. It establishes a simple, easy-to-program, O(nm(2))-time algorithm that determines whether there is a PPH solution for input genotypes and produces a linear-space data structure to represent all of the solutions. The approach allows complete, self-contained proofs. In addition to algorithmic simplicity, the approach here makes the representation of all solutions more intuitive than in Gusfield (2002), and solves another goal from that paper, namely, to prove a nontrivial upper bound on the number of PPH solutions, showing that that number is vastly smaller than the number of haplotype solutions (each solution being a set of n pairs of haplotypes that can generate the genotypes) when the perfect phylogeny requirement is not imposed.     

A statistical framework for haplotype block inference

The existence of haplotype blocks transmitted from parents to offspring has been suggested recently. This has created an interest in the inference of the block structure and length. The motivation is that haplotype blocks that are characterized well will make it relatively easier to quickly map all the genes carrying human diseases. To study the inference of haplotype block systematically, we propose a statistical framework. In this framework, the optimal haplotype block partitioning is formulated as the problem of statistical model selection; missing data can be handled in a standard statistical way; population strata can be implemented; block structure inference/hypothesis testing can be performed; prior knowledge, if present, can be incorporated to perform a Bayesian inference. The algorithm is linear in the number of loci, instead of NP-hard for many such algorithms. We illustrate the applications of our method to both simulated and real data sets.