An STS-based radiation hybrid map of the human genome

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y-linked scaffolds using existing data. A total of 11.5 Mb Y-linked sequences were clustered into 33 scaffolds, and 187 protein-coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole-genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post-Neolithic period.

Construction and integration of three de novo Japanese human genome assemblies toward a population-specific reference

The complete human genome sequence is used as a reference for next-generation sequencing analyses. However, some ethnic ancestries are under-represented in the reference genome (e.g., GRCh37) due to its bias toward European and African ancestries. Here, we perform de novo assembly of three Japanese male genomes using > 100× Pacific Biosciences long reads and Bionano Genomics optical maps per sample. We integrate the genomes using the major allele for consensus and anchor the scaffolds using genetic and radiation hybrid maps to reconstruct each chromosome. The resulting genome sequence, JG1, is contiguous, accurate, and carries the Japanese major allele at most loci. We adopt JG1 as the reference for confirmatory exome re-analyses of seven rare-disease Japanese families and find that re-analysis using JG1 reduces total candidate variant calls versus GRCh37 while retaining disease-causing variants. These results suggest that integrating multiple genomes from a single population can aid genome analyses of that population.

Reduced Folate Carrier: an Entry Receptor for a Novel Feline Leukemia Virus Variant

Feline leukemia virus (FeLV) is horizontally transmitted among cats and causes a variety of hematopoietic disorders. Five subgroups of FeLV, A to D and T, each with distinct receptor usages, have been described. Recently, we identified a new FeLV Env (TG35-2) gene from a pseudotyped virus that does not belong to any known subgroup. FeLV-A is the primary virus from which other subgroups have emerged via mutation or recombination of the subgroup A env gene. Retrovirus entry into cells is mediated by the interaction of envelope protein (Env) with specific cell surface receptors. Here, phenotypic screening of a human/hamster radiation hybrid panel identified SLC19A1, a feline reduced folate carrier (RFC) and potential receptor for TG35-2-phenotypic virus. RFC is a multipass transmembrane protein. Feline and human RFC cDNAs conferred susceptibility to TG35-2-pseudotyped virus when introduced into nonpermissive cells but did not render these cells permissive to other FeLV subgroups or feline endogenous retrovirus. Moreover, human cells with genomic deletion of RFC were nonpermissive for TG35-2-pseudotyped virus infection, but the introduction of feline and human cDNAs rendered them permissive. Mutation analysis of FeLV Env demonstrated that amino acid substitutions within variable region A altered the specificity of the Env-receptor interaction. We isolated and reconstructed the full-length infectious TG35-2-phenotypic provirus from a naturally FeLV-infected cat, from which the FeLV Env (TG35-2) gene was previously isolated, and compared the replication of the virus in hematopoietic cell lines with that of FeLV-A 61E by measuring the viral RNA copy numbers. These results provide a tool for further investigation of FeLV infectious disease.IMPORTANCE Feline leukemia virus (FeLV) is a member of the genus Gammaretrovirus, which causes malignant diseases in cats. The most prevalent FeLV among cats is FeLV subgroup A (FeLV-A), and specific binding of FeLV-A Env to its viral receptor, thiamine transporter feTHTR1, is the first step of infection. In infected cats, novel variants of FeLV with altered receptor specificity for viral entry have emerged by mutation or recombination of the env gene. A novel FeLV variant arose from a subtle mutation of FeLV-A Env, which altered the specific interaction of the virus with its receptor. RFC, a folate transporter, is a potential receptor for the novel FeLV variant. The perturbation of specific retrovirus-receptor interactions under selective pressure by the host results in the emergence of novel viruses.

Dissecting Plant Chromosomes by the Use of Ionizing Radiation

Radiation treatment of genomes is used to generate chromosome breaks for numerous applications. This protocol describes the preparation of seeds and the determination of the optimal level of irradiation dosage for the creation of a radiation hybrid (RH) population. These RH lines can be used to generate high-resolution physical maps for the assembly of sequenced genomes as well as the fine mapping of genes. This procedure can also be used for mutation breeding and forward/reverse genetics.

Genetic screening reveals a link between Wnt signaling and antitubulin drugs

The antitubulin drugs, paclitaxel (PX) and colchicine (COL), inhibit cell growth and are therapeutically valuable. PX stabilizes microtubules, while COL promotes their depolymerization. But, the drug concentrations that alter tubulin polymerization are hundreds of times higher than their clinically useful levels. To map genetic targets for drug action at single-gene resolution, we used a human radiation hybrid panel. We identified loci that affected cell survival in the presence of five compounds of medical relevance. For PX and COL, the zinc and ring finger 3 (ZNRF3) gene dominated the genetic landscape at therapeutic concentrations. ZNRF3 encodes an R-spondin regulated receptor that inhibits Wingless/Int (Wnt) signaling. Overexpression of the ZNRF3 gene shielded cells from antitubulin drug action, while small interfering RNA knockdowns resulted in sensitization. Further a potent pharmacological inhibitor of Wnt signaling, Wnt-C59, protected cells from PX and COL. Our results suggest that the antitubulin drugs perturb microtubule dynamics, thereby influencing Wnt signaling.

Radiation Hybrid Map of Barley Chromosome 3H

Assembly of the barley (Hordeum vulgare L.) genome is complicated by its large size (5.1 Gb) and proportion of repetitive elements (84%). This process is facilitated by high resolution maps for aligning bacterial artificial chromosome (BAC) contigs along chromosomes. Available genetic maps, however, do not provide accurate information on the physical position of a large portion of the genome located in recombination-poor regions. Radiation hybrid (RH) mapping is an alternative approach, which is based on radiation-induced deletions along the length of chromosomes. In this study, the first RH map for barley chromosome 3H was developed. In total, 373 in vivo RH lines were generated by irradiating wheat (Triticum aestivum L.)-barley chromosome 3H addition lines and crossing them to a normal wheat cultivar. Each RH informative line (containing deletions) had, on average, three deletions. The induced deletion size varied from 36.58 Kb to 576.00 Mb, with an average length of 52.42 Mb. This initial chromosome 3H radiation hybrid (3H-RH) map had a 9.53× higher resolution than an analogous genetic map, reaching a maximum of >262.40× resolution in regions around the centromere. The final RH map was 3066.1 cR in length, with a 0.76 Mb resolution. It was estimated that the map resolution can be improved to an average of 30.34 Kb by saturating the 3H-RH map with molecular markers. The generated RH panel enabled alignment of BAC and sequenced contigs as small as 1.50 Kb in size. The high resolution and the coverage of poor-recombination regions make RH maps an ideal resource for barley genome assembly, as well as other genetic studies.

Reliable Radiation Hybrid Maps: An Efficient Scalable Clustering-Based Approach

The process of mapping markers from radiation hybrid mapping (RHM) experiments is equivalent to the traveling salesman problem and, thereby, has combinatorial complexity. As an additional problem, experiments typically result in some unreliable markers that reduce the overall quality of the map. We propose a clustering approach for addressing both problems efficiently by eliminating unreliable markers without the need for mapping the complete set of markers. Traditional approaches for eliminating markers use resampling of the full data set, which has an even higher computational complexity than the original mapping problem. In contrast, the proposed approach uses a divide-and-conquer strategy to construct framework maps based on clusters that exclude unreliable markers. Clusters are ordered using parallel processing and are then combined to form the complete map. We present three algorithms that explore the trade-off between the number of markers included in the map and placement accuracy. Using an RHM data set of the human genome, we compare the framework maps from our proposed approaches with published physical maps and with the results of using the Carthagene tool. Overall, our approaches have a very low computational complexity and produce solid framework maps with good chromosome coverage and high agreement with the physical map marker order.

An update of the goat genome assembly using dense radiation hybrid maps allows detailed analysis of evolutionary rearrangements in Bovidae

Background

The domestic goat (Capra hircus), an important livestock species, belongs to a clade of Ruminantia, Bovidae, together with cattle, buffalo and sheep. The history of genome evolution and chromosomal rearrangements on a small scale in ruminants remain speculative. Recently completed goat genome sequence was released but is still in a draft stage. The draft sequence used a variety of assembly packages, as well as a radiation hybrid (RH) map of chromosome 1 as part of its validation.

Results

Using an improved RH mapping pipeline, whole-genome dense maps of 45,953 SNP markers were constructed with statistical confidence measures and the saturated maps provided a fine map resolution of approximate 65 kb. Linking RH maps to the goat sequences showed that the assemblies of scaffolds/super-scaffolds were globally accurate. However, we observed certain flaws linked to the process of anchoring chromosome using conserved synteny with cattle. Chromosome assignments, long-range order, and orientation of the scaffolds were reassessed in an updated genome sequence version. We also present new results exploiting the updated goat genome sequence to understand genomic rearrangements and chromosome evolution between mammals during species radiations. The sequence architecture of rearrangement sites between the goat and cattle genomes presented abundant segmental duplication on regions of goat chromosome 9 and 14, as well as new insertions in homologous cattle genome regions. This complex interplay between duplicated sequences and Robertsonian translocations highlights the rearrangement mechanism of centromeric nonallelic homologous recombination (NAHR) in mammals. We observed that species-specific shifts in ANKRD26 gene duplication are coincident with breakpoint reuse in divergent lineages and this gene family may play a role in chromosome stabilization in chromosome evolution.

Conclusions

We generated dense maps of the complete whole goat genome. The chromosomal maps allowed us to anchor and orientate assembled genome scaffolds along the chromosomes, annotate chromosome rearrangements and thereby get a better understanding of the genome evolution of ruminants and other mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-625) contains supplementary material, which is available to authorized users.

Construction of a radiation hybrid panel and the first yellowtail (Seriola quinqueradiata) radiation hybrid map using a nanofluidic dynamic array

Background

Yellowtail (Seriola quinqueradiata) are an economically important species in Japan. However, there are currently no methods for captive breeding and early rearing for yellowtail. Thus, the commercial cultivation of this species is reliant upon the capture of wild immature fish. Given this, there is a need to develop captive breeding techniques to reduce pressure on wild stocks and facilitate the sustainable development of yellowtail aquaculture. We constructed a whole genome radiation hybrid (RH) panel for yellowtail gene mapping and developed a framework physical map using a nanofluidic dynamic array to use SNPs (single nucleotide polymorphisms) in ESTs (expressed sequence tags) for the DNA-assisted breeding of yellowtail.

Results

Clonal RH cell lines were obtained after ionizing radiation; specifically, 78, 64, 129, 55, 42, and 53 clones were isolated after treatment with 3,000, 4,000, 5,000, 6,000, 8,000, or 10,000 rads, respectively. A total of 421 hybrid cell lines were obtained by fusion with mouse B78 cells. Ninety-four microsatellite markers used in the genetic linkage map were genotyped using the 421 hybrid cell lines. Based upon marker retention and genome coverage, we selected 93 hybrid cell lines to form an RH panel. Importantly, we performed the first genotyping of yellowtail markers in an RH panel using a nanofluidic dynamic array (Fluidigm, CA, USA). Then, 580 markers containing ESTs and SNPs were mapped in the first yellowtail RH map.

Conclusions

We successfully developed a yellowtail RH panel to facilitate the localization of markers. Using this, a framework RH map was constructed with 580 markers. This high-density physical map will serve as a useful tool for the identification of genes related to important breeding traits using genetic structural information, such as conserved synteny. Moreover, in a comparison of 30 sequences in the RH group 1 (SQ1), yellowtail appeared to be evolutionarily closer to medaka and the green-spotted pufferfish than to zebrafish. We suggest that synteny analysis may be potentially useful as a tool to investigate chromosomal evolution by comparison with model fish.

Endosperm tolerance of paternal aneuploidy allows radiation hybrid mapping of the wheat D-genome and a measure of γ ray-induced chromosome breaks

Physical mapping and genome sequencing are underway for the ≈17 Gb wheat genome. Physical mapping methods independent of meiotic recombination, such as radiation hybrid (RH) mapping, will aid precise anchoring of BAC contigs in the large regions of suppressed recombination in Triticeae genomes. Reports of endosperm development following pollination with irradiated pollen at dosages that cause embryo abortion prompted us to investigate endosperm as a potential source of RH mapping germplasm. Here, we report a novel approach to construct RH based physical maps of all seven D-genome chromosomes of the hexaploid wheat ‘Chinese Spring’, simultaneously. An 81-member subset of endosperm samples derived from 20-Gy irradiated pollen was genotyped for deletions, and 737 markers were mapped on seven D-genome chromosomes. Analysis of well-defined regions of six chromosomes suggested a map resolution of ∼830 kb could be achieved; this estimate was validated with assays of markers from a sequenced contig. We estimate that the panel contains ∼6,000 deletion bins for D-genome chromosomes and will require ∼18,000 markers for high resolution mapping. Map-based deletion estimates revealed a majority of 1–20 Mb interstitial deletions suggesting mutagenic repair of double-strand breaks in pollen provides a useful resource for RH mapping and map based cloning studies.

Physical mapping resources for large plant genomes: radiation hybrids for wheat D-genome progenitor Aegilops tauschii

BACKGROUND

Development of a high quality reference sequence is a daunting task in crops like wheat with large (~17Gb), highly repetitive (>80%) and polyploid genome. To achieve complete sequence assembly of such genomes, development of a high quality physical map is a necessary first step. However, due to the lack of recombination in certain regions of the chromosomes, genetic mapping, which uses recombination frequency to map marker loci, alone is not sufficient to develop high quality marker scaffolds for a sequence ready physical map. Radiation hybrid (RH) mapping, which uses radiation induced chromosomal breaks, has proven to be a successful approach for developing marker scaffolds for sequence assembly in animal systems. Here, the development and characterization of a RH panel for the mapping of D-genome of wheat progenitor Aegilops tauschii is reported.

RESULTS

Radiation dosages of 350 and 450 Gy were optimized for seed irradiation of a synthetic hexaploid (AABBDD) wheat with the D-genome of Ae. tauschii accession AL8/78. The surviving plants after irradiation were crossed to durum wheat (AABB), to produce pentaploid RH1s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH1 plants was obtained, of which 592 plants were generated from the mature RH1 seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH1 seeds. This panel showed a homogenous marker loss (2.1%) after screening with SSR markers uniformly covering all the D-genome chromosomes. Different marker systems mostly detected different lines with deletions. Using markers covering known distances, the mapping resolution of this RH panel was estimated to be <140kb. Analysis of only 16 RH lines carrying deletions on chromosome 2D resulted in a physical map with cM/cR ratio of 1:5.2 and 15 distinct bins. Additionally, with this small set of lines, almost all the tested ESTs could be mapped. A set of 399 most informative RH lines with an average deletion frequency of ~10% were identified for developing high density marker scaffolds of the D-genome.

CONCLUSIONS

The RH panel reported here is the first developed for any wild ancestor of a major cultivated plant species. The results provided insight into various aspects of RH mapping in plants, including the genetically effective cell number for wheat (for the first time) and the potential implementation of this technique in other plant species. This RH panel will be an invaluable resource for mapping gene based markers, developing a complete marker scaffold for the whole genome sequence assembly, fine mapping of markers and functional characterization of genes and gene networks present on the D-genome.

Physical mapping resources for large plant genomes: radiation hybrids for wheat D-genome progenitor Aegilops tauschii

BACKGROUND

Development of a high quality reference sequence is a daunting task in crops like wheat with large (~17Gb), highly repetitive (>80%) and polyploid genome. To achieve complete sequence assembly of such genomes, development of a high quality physical map is a necessary first step. However, due to the lack of recombination in certain regions of the chromosomes, genetic mapping, which uses recombination frequency to map marker loci, alone is not sufficient to develop high quality marker scaffolds for a sequence ready physical map. Radiation hybrid (RH) mapping, which uses radiation induced chromosomal breaks, has proven to be a successful approach for developing marker scaffolds for sequence assembly in animal systems. Here, the development and characterization of a RH panel for the mapping of D-genome of wheat progenitor Aegilops tauschii is reported.

RESULTS

Radiation dosages of 350 and 450 Gy were optimized for seed irradiation of a synthetic hexaploid (AABBDD) wheat with the D-genome of Ae. tauschii accession AL8/78. The surviving plants after irradiation were crossed to durum wheat (AABB), to produce pentaploid RH1s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH1 plants was obtained, of which 592 plants were generated from the mature RH1 seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH1 seeds. This panel showed a homogenous marker loss (2.1%) after screening with SSR markers uniformly covering all the D-genome chromosomes. Different marker systems mostly detected different lines with deletions. Using markers covering known distances, the mapping resolution of this RH panel was estimated to be <140kb. Analysis of only 16 RH lines carrying deletions on chromosome 2D resulted in a physical map with cM/cR ratio of 1:5.2 and 15 distinct bins. Additionally, with this small set of lines, almost all the tested ESTs could be mapped. A set of 399 most informative RH lines with an average deletion frequency of ~10% were identified for developing high density marker scaffolds of the D-genome.

CONCLUSIONS

The RH panel reported here is the first developed for any wild ancestor of a major cultivated plant species. The results provided insight into various aspects of RH mapping in plants, including the genetically effective cell number for wheat (for the first time) and the potential implementation of this technique in other plant species. This RH panel will be an invaluable resource for mapping gene based markers, developing a complete marker scaffold for the whole genome sequence assembly, fine mapping of markers and functional characterization of genes and gene networks present on the D-genome.

Loss of TIMP-3 promotes tumor invasion via elevated IL-6 production and predicts poor survival and relapse in HPV-infected non-small cell lung cancer

Human papillomavirus (HPV) 16/18 E6 oncoprotein is expressed in lung tumors and is associated with p53 inactivation. The tissue inhibitor of metalloproteinase 3 (TIMP-3) is essential for limiting inflammation; therefore, we expected that TIMP-3 loss might induce chronic inflammation, thereby promoting tumor malignancy as well as poor survival and relapse in patients with HPV-infected non-small cell lung cancer. In this study, the loss of TIMP-3 by loss of heterozygosity and/or promoter hypermethylation was more frequent in HPV16/18 E6-positive tumors than in E6-negative tumors. To explore the possible underlying mechanism, E6-negative TL4 and CL1-0 cells were transfected with an E6 cDNA plasmid. A marked decrease in TIMP-3 expression was caused by promoter hypermethylation via increased DNA (cytosine-5-)-methyltransferase 1 (DNMT1) expression. Mechanistic studies indicated that TIMP-3 loss promoted interleukin-6 (IL-6) production, which led to cell invasion and anchorage-independent growth on soft agar plates. Kaplan-Meier and Cox regression models showed that patients with low-TIMP-3/high-IL-6 tumors had shorter overall survival and relapse-free survival periods when compared with patients with high-TIMP-3/low-IL-6 tumors. In summary, loss of TIMP-3 may increase IL-6 production via the tumor necrosis factor α/nuclear factor κB axis, thereby promoting tumor malignancy and subsequent relapse and poor survival in patients with HPV-infected non-small cell lung cancer.

INCONSISTENCIES BETWEEN MAPS OF HUMAN CHROMOSOME 22 CORRELATE WITH INCREASED FREQUENCY OF DISEASE-RELATED LOCI

The relationships between genetic or radiation hybrid (RH) and sequence maps of chromosome 22 have been reconsidered based on the sequence map. Integrated maps have been constructed by retaining only common markers between genetic or RH maps and the sequence map. Local inversions of markers have been detected. Ratios between either genetic or RH distances and sequence-based distances have been calculated for each map interval. Hot zones for recombination or radiation breakage have been delineated by merging together intervals displaying high distance ratios and located close to each other for sequence-constrained RH maps, and for female and male genetic maps. A statistically significant positive correlation was found between the distribution of disease-related genes and the hot zones for recombination or radiation breakage on both female genetic and Stanford-G3 RH maps. This observation indicates that investigation of chromosomal regions displaying inconsistencies between RH or genetic linkage and sequence-based maps can accelerate the initial phase of identification of disease-associated genes.

Effects of genome-wide copy number variation on expression in mammalian cells

Background

There is only a limited understanding of the relation between copy number and expression for mammalian genes. We fine mapped cis and trans regulatory loci due to copy number change for essentially all genes using a human-hamster radiation hybrid (RH) panel. These loci are called copy number expression quantitative trait loci (ceQTLs).

Results

Unexpected findings from a previous study of a mouse-hamster RH panel were replicated. These findings included decreased expression as a result of increased copy number for 30% of genes and an attenuated relationship between expression and copy number on the X chromosome suggesting an Xist independent form of dosage compensation. In a separate glioblastoma dataset, we found conservation of genes in which dosage was negatively correlated with gene expression. These genes were enriched in signaling and receptor activities. The observation of attenuated X-linked gene expression in response to increased gene number was also replicated in the glioblastoma dataset. Of 523 gene deserts of size > 600 kb in the human RH panel, 325 contained trans ceQTLs with -log₁₀ P > 4.1. Recently discovered genes, ultra conserved regions, noncoding RNAs and microRNAs explained only a small fraction of the results, suggesting a substantial portion of gene deserts harbor as yet unidentified functional elements.

Conclusion

Radiation hybrids are a useful tool for high resolution mapping of cis and trans loci capable of affecting gene expression due to copy number change. Analysis of two independent radiation hybrid panels show agreement in their findings and may serve as a discovery source for novel regulatory loci in noncoding regions of the genome.

A genome-wide map of human genetic interactions inferred from radiation hybrid genotypes

Using radiation hybrid genotyping data, 99% of all possible gene pairs across the mammalian genome were tested for interactions based on co-retention frequencies higher (attraction) or lower (repulsion) than chance. Gene interaction networks constructed from six independent data sets overlapped strongly. Combining the data sets resulted in a network of more than seven million interactions, almost all attractive. This network overlapped with protein-protein interaction networks on multiple measures and also confirmed the relationship between essentiality and centrality. In contrast to other biological networks, the radiation hybrid network did not show a scale-free distribution of connectivity but was Gaussian-like, suggesting a closer approach to saturation. The radiation hybrid (RH) network constitutes a platform for understanding the systems biology of the mammalian cell.

Methods for the analysis of large gene deletions and their application in some monogenic disorders

Számos monogénesen öröklődő kórképben a betegséget okozó gén teljes vagy részleges deletiója, illetve kópiaszámának megváltozása patogenetikai tényezőként jön számításba. A direkt DNS-szekvenálás nem alkalmas a gén nagy deletiójának, illetve kópiaszám-változásának kimutatására. Az összefoglalóban a szerzők áttekintik a nagy géndeletio vizsgálómódszereit, és két, monogénesen öröklődő betegségben végzett saját vizsgálataik példáján keresztül bemutatják a módszerek gyakorlati alkalmazásának lehetőségeit. Vázolják a géndeletio-vizsgálat hagyományos (kromoszóma-sávtechnika, Southern-blot, fluoreszcens in situ hibridizáció) és polimeráz láncreakcióra alapozott módszereit (denaturáló nagy felbontóképességű folyadékkromatográfia, kvantitatív valós idejű polimeráz láncreakció, mikroszatellitamarker-analízis, multiplex amplifikálhatópróba-hibridizáció, multiplex ligatióspróba-analízis), valamint a technikai és informatikai haladás legújabb vívmányait (komparatív genomhibridizálás, „array” analízis). Saját vizsgálataikban von Hippel–Lindau-szindrómában szenvedő betegekben kvantitatív valós idejű polimeráz láncreakció és multiplex ligatióspróba-amplifikálás alkalmazásával bemutatják a VHL, illetve congenitalis adrenalis hyperplasiás betegekben a CYP21A2 géndeletio-vizsgálat eredményeit és ezek klinikai jelentőségét.

A high-resolution radiation hybrid map of rhesus macaque chromosome 5 identifies rearrangements in the genome assembly

A 10,000-rad radiation hybrid (RH) cell panel of the rhesus macaque was generated to construct a comprehensive RH map of chromosome 5. The map represents 218 markers typed in 185 RH clones. The 4846-cR map has an average marker spacing of 798 kb. Alignments of the RH map to macaque and human genome sequences confirm a large inversion and reveal a previously unreported telomeric inversion. The macaque genome sequence indicates small translocations from the ancestral homolog of macaque chromosome 5 to macaque chromosomes 1 and 6. The RH map suggests that these are probably assembly artifacts. Unlike the genome sequence, the RH mapping data indicate the conservation of synteny between macaque chromosome 5 and human chromosome 4. This study shows that the 10,000-rad panel is appropriate for the generation of a high-resolution whole-genome RH map suitable for the verification of the rhesus genome assembly.

Improved assay-dependent searching of nucleic acid sequence databases

Nucleic acid-based biochemical assays are crucial to modern biology. Key applications, such as detection of bacterial, viral and fungal pathogens, require detailed knowledge of assay sensitivity and specificity to obtain reliable results. Improved methods to predict assay performance are needed for exploiting the exponentially growing amount of DNA sequence data and for reducing the experimental effort required to develop robust detection assays. Toward this goal, we present an algorithm for the calculation of sequence similarity based on DNA thermodynamics. In our approach, search queries consist of one to three oligonucleotide sequences representing either a hybridization probe, a pair of Padlock probes or a pair of PCR primers with an optional TaqMan™ probe (i.e. in silico or ‘virtual’ PCR). Matches are reported if the query and target satisfy both the thermodynamics of the assay (binding at a specified hybridization temperature and/or change in free energy) and the relevant biological constraints (assay sequences binding to the correct target duplex strands in the required orientations). The sensitivity and specificity of our method is evaluated by comparing predicted to known sequence tagged sites in the human genome. Free energy is shown to be a more sensitive and specific match criterion than hybridization temperature.

Radiation hybrid (RH) and HAPPY mapping in plants

Radiation hybrid (RH) and HAPPY mapping are two technologies used in animal systems that have attracted the attention of the plant genetics community because they bridge the resolution gap between meiotic and BAC-based physical mapping that would facilitate the analysis of plant species lacking substantial genomics resources. Research has shown that the essence of these approaches can be applied and that a variety of strategies can be used to produce mapping panels. Mapping panels composed of live plants, protoplast fusion cultures, and sub-genomic DNA samples have been described. The resolution achievable by RH mapping panels involving live-plant derivatives of a monosomic maize (Zea mays) chromosome 9 addition in allohexaploid oat (Avena sativa), a monosomic chromosome 1D addition in allotetraploid durum wheat (Triticum turgidum), and interspecific hybrids between two tetraploid cotton species (G. hirsutum and G. barbadense), has been estimated to range from 0.6 to 6 Mb. On the other hand, a more comprehensive evaluation of one panel from durum wheat suggests that a higher mapping resolution (approximately 200 kb) is possible. In cases involving RH mapping panels based on barley (Hordeum vulgare)-tobacco (Nicotiana tabacum) protoplast fusions or a HAPPY mapping panel based on genomic DNA from Arabidopsis thaliana, the potential mapping resolution appears to be higher (50 to 200 kb). Despite these encouraging results, the application of either RH or HAPPY mapping in plants is still in the experimental phase and additional work is clearly needed before these methods are more routinely utilized.

Identification of the myelin protein plasmolipin as the cell entry receptor for Mus caroli endogenous retrovirus

The Asian wild mouse species Mus caroli harbors an endogenous retrovirus (McERV) that is closely related to but distinct from the endogenous retrovirus family defined by the Mus dunni endogenous virus and the Mus musculus endogenous retrovirus. McERV could infect some cell types from humans, dogs, and rats, but not all, and did not infect any mouse cell line tested. Because of its interesting host range and proposed ancestral relationship to primate retroviruses and because none of the entry receptors for this family of retroviruses had been identified, we began a search for the McERV receptor. We determined the chromosomal location of the receptor gene in the human genome by phenotypic screening of the G3 human-hamster radiation hybrid cell line panel and confirmed the localization by assaying for receptor activity conferred by bacterial artificial chromosome (BAC) clones spanning the region. We next localized the gene more precisely in one positive BAC by assaying for receptor activity following BAC digestion with several restriction enzymes that cleaved different sets of genes, and we confirmed that the final candidate gene, plasmolipin (PLLP; TM4SF11), is the novel receptor by showing that the expression of the human PLLP cDNA renders hamster and mouse cells susceptible to McERV infection. PLLP functions as a voltage-dependent potassium ion channel and is expressed primarily in kidney and brain, helping to explain the limited range of cell types that McERV can infect. Interestingly, mouse PLLP also functioned well as a receptor for McERV but was simply not expressed in the mouse cell types that we originally tested.

Zebrafish ae2.2 encodes a second slc4a2 anion exchanger

The genome of zebrafish (Danio rerio) encodes two unlinked genes equally closely related to the SLC4A2/AE2 anion exchanger genes of mammals. One of these is the recently reported zebrafish ae2 gene (Shmukler BE, Kurschat CE, Ackermann GE, Jiang L, Zhou Y, Barut B, Stuart-Tilley AK, Zhao J, Zon LI, Drummond IA, Vandorpe DH, Paw BH, Alper SL. Am J Physiol Renal Physiol Renal Physiol 289: F835-F849, 2005), now called ae2.1. We now report the structural and functional characterization of Ae2.2, the product of the second zebrafish Ae2 gene, ae2.2. The ae2.2 gene of zebrafish linkage group 24 encodes a polypeptide of 1,232 aa in length, sharing 70% amino acid identity with zebrafish Ae2.1 and 67% identity with mouse AE2a. Zebrafish Ae2.2 expressed in Xenopus oocytes encodes a 135-kDa polypeptide that mediates bidirectional, DIDS-sensitive Cl(-)/Cl(-) exchange and Cl(-)/HCO3(-) exchange. Ae2.2-mediated Cl(-)/Cl(-) exchange is cation independent, voltage insensitive, and electroneutral. Acute regulation of anion exchange mediated by Ae2.2 includes activation by NH4+ and independent inhibition by acidic intracellular pH and by acidic extracellular pH. In situ hybridization reveals low-level expression of Ae2.2 mRNA in zebrafish embryo, most notably in posterior tectum, eye, pharynx, epidermal cells, and axial vascular structures, without notable expression in the Ae2.1-expressing pronephric duct. Knockdown of Ae2.2 mRNA, of Ae2.1 mRNA, or of both with nontoxic or minimally toxic levels of N-morpholino oligomers produced no grossly detectable morphological phenotype, and preserved normal structure of the head and the pronephric duct at 24 h postfertilization.

Radiation hybrid maps of Medaka chromosomes LG 12, 17, and 22

The Medaka is an excellent genetic system for studies of vertebrate development and disease and environmental and evolutionary biology studies. To facilitate the mapping of markers or the cloning of affected genes in Medaka mutants identified by forward-genetic screens, we have established a panel of whole-genome radiation hybrids (RHs) and RH maps for three Medaka chromosomes. RH mapping is useful, since markers to be mapped need not be polymorphic and one can establish the order of markers that are difficult to resolve by genetic mapping owing to low genetic recombination rates. RHs were generated by fusing the irradiated donor, OLF-136 Medaka cell line, with the host B78 mouse melanoma cells. Of 290 initial RH clones, we selected 93 on the basis of high retention of fragments of the Medaka genome to establish a panel that allows genotyping in the 96-well format. RH maps for linkage groups 12, 17, and 22 were generated using 159 markers. The average retention for the three chromosomes was 19% and the average break point frequency was ∼33 kb/cR. We estimate the potential resolution of the RH panel to be ∼186 kb, which is high enough for integrating RH data with bacterial artificial chromosome clones. Thus, this first RH panel will be a useful tool for mapping mutated genes in Medaka.

Zebrafish orthologs of human muscular dystrophy genes

Background

Human muscular dystrophies are a heterogeneous group of genetic disorders which cause decreased muscle strength and often result in premature death. There is no known cure for muscular dystrophy, nor have all causative genes been identified. Recent work in the small vertebrate zebrafish Danio rerio suggests that mutation or misregulation of zebrafish dystrophy orthologs can also cause muscular degeneration phenotypes in fish. To aid in the identification of new causative genes, this study identifies and maps zebrafish orthologs for all known human muscular dystrophy genes.

Results

Zebrafish sequence databases were queried for transcripts orthologous to human dystrophy-causing genes, identifying transcripts for 28 out of 29 genes of interest. In addition, the genomic locations of all 29 genes have been found, allowing rapid candidate gene discovery during genetic mapping of zebrafish dystrophy mutants. 19 genes show conservation of syntenic relationships with humans and at least two genes appear to be duplicated in zebrafish. Significant sequence coverage on one or more BAC clone(s) was also identified for 24 of the genes to provide better local sequence information and easy updating of genomic locations as the zebrafish genome assembly continues to evolve.

Conclusion

This resource supports zebrafish as a dystrophy model, suggesting maintenance of all known dystrophy-associated genes in the zebrafish genome. Coupled with the ability to conduct genetic screens and small molecule screens, zebrafish are thus an attractive model organism for isolating new dystrophy-causing genes/pathways and for use in high-throughput therapeutic discovery.

Radiation hybrid map of the porcine genome comprising 2035 EST loci

The IMpRH(7000-rad) radiation hybrid panel was used to map 2035 expressed sequence tags (ESTs) at a minimum LOD score of 4.0. A total of 134 linkage groups covers 57,192 cR or 78% of the predicted size of the porcine and 71% of the human genome, respectively. Approximately 81% (1649) of the porcine ESTs were annotated against the NCBI nonredundant database; 1422 mapped in silico to a location in build 35.1 of the human genome sequence (HGS) and 1185 to a gene and location in build 35.1 HGS. The map revealed 40 major breaks in synteny (1.00e (-25 )and lower) with the human genome, 37 of which fall within a single chromosome. At this improved level of resolution and coverage, porcine chromosomes (SSC) 2, 5, 6, 7, 12, and 14 remain "gene-rich" and homologous to human chromosomes (HSA) 17, 19, and 22, while SSC 1, 8, 11, and X have been confirmed to correspond to the "gene-deserts" on HSA 18, 4, 13, and X.

High-resolution radiation hybrid map of wheat chromosome 1D

Physical mapping methods that do not rely on meiotic recombination are necessary for complex polyploid genomes such as wheat (Triticum aestivum L.). This need is due to the uneven distribution of recombination and significant variation in genetic to physical distance ratios. One method that has proven valuable in a number of nonplant and plant systems is radiation hybrid (RH) mapping. This work presents, for the first time, a high-resolution radiation hybrid map of wheat chromosome 1D (D genome) in a tetraploid durum wheat (T. turgidum L., AB genomes) background. An RH panel of 87 lines was used to map 378 molecular markers, which detected 2312 chromosome breaks. The total map distance ranged from approximately 3,341 cR(35,000) for five major linkage groups to 11,773 cR(35,000) for a comprehensive map. The mapping resolution was estimated to be approximately 199 kb/break and provided the starting point for BAC contig alignment. To date, this is the highest resolution that has been obtained by plant RH mapping and serves as a first step for the development of RH resources in wheat.

Development of a deer mouse whole-genome radiation hybrid panel and comparative mapping of Mus chromosome 11 loci

A 5000-rad whole-genome radiation hybrid cell panel (BW5000) was developed for mapping the deer mouse (Peromyscus maniculatus bairdii) genome. The panel consists of 103 cell lines and has an estimated marker retention frequency of 63.9% (range, 28%-88%) based on PCR typing of 30 Type I (coding gene) and 25 Type II (microsatellite) markers. Using the composite Mus map, Type I markers were selected from six Mus chromosomes, 22 of which are on Mus Chr 11. Fifteen of the Mus Chr 11 markers were simultaneously mapped on an interspecific (P. maniculatus x P. polionotus) backcross panel to test the utility of the radiation hybrid panel, create a framework map, and help establish gene order. The radiation hybrids have effectively detected linkage in the deer mouse genome between markers as far apart as 6.7 cM and resolved markers that are, in the Mus genome, as close as 0.2 Mb. Combined results from both panels have indicated a high degree of gene order conservation of the telomeric 64 cM of Mus Chr 11 in the deer mouse genome. The remaining centromeric portion also shows gene order conservation with the deer mouse but as a separate linkage group. This indicates a translocation of that portion of Mus Chr 11 in P. maniculatus and is consistent with rearrangement breakpoints observed between Mus and other mammalian genomes, including rat and human. Furthermore, this separate linkage group is likely to reside in a chromosomal region of inversion polymorphism between P. maniculatus and P. polionotus.

A 12,000-rad porcine radiation hybrid (IMNpRH2) panel refines the conserved synteny between SSC12 and HSA17

Reverse or bidirectional Zoo-FISH suggests that synteny between porcine chromosome 12 (SSC12) and human chromosome 17 (HSA17) is completely conserved. The construction of a high-resolution radiation hybrid (RH) map for SSC12 provides a unique opportunity to determine whether chromosomal synteny is reflected at the molecular level by comparative gene mapping of SSC12 and HSA17. We report an initial, high-resolution RH map of SSC12 on the 12,000-rad IMNpRH2 panel using CarthaGene software. This map contains a total of 320 markers, including 20 microsatellites and 300 ESTs/genes, covering approximately 4836.9 cR12,000. The markers were ordered in 16 linkage groups at LOD 6.0 using framework markers previously mapped on the IMpRH7000-rad SSC12 and porcine genetic maps. Ten linkage groups ordered more than 10 markers, with the largest containing 101 STSs. The resolution of the current RH map is approximately 15.3 kb/cR on SSC12, a significant improvement over the second-generation EST SSC12 RH7000-rad map of 103 ESTs and 15 framework markers covering approximately 2287.2 cR7000. Compared to HSA17, six distinct segments were identified, revealing macro-rearrangements within the apparently complete synteny between SSC12 and HSA17. Further analysis of the order of 245 genes (ESTs) on HSA17 and SSC12 also revealed several micro-rearrangements within a synteny segment. A high-resolution SSC12 RH12,000-rad map will be useful in fine-mapping QTL and as a scaffold for sequencing this chromosome.

bloodthirsty, an RBCC/TRIM gene required for erythropoiesis in zebrafish

The Antarctic icefishes (family Channichthyidae, suborder Notothenioidei) constitute the only vertebrate taxon that fails to produce red blood cells. These fishes can be paired with closely related, but erythrocyte-producing, notothenioids to discover erythropoietic genes via representational difference analysis. Using a B30.2-domain-encoding DNA probe so derived from the hematopoietic kidney (pronephros) of a red-blooded Antarctic rockcod, Notothenia coriiceps, we discovered a related, novel gene, bloodthirsty (bty), that encoded a 547-residue protein that contains sequential RING finger, B Box, coiled-coil, and B30.2 domains. bty mRNA was expressed by the pronephric kidney of N. coriiceps at a steady-state level 10-fold greater than that found in the kidney of the icefish Chaenocephalus aceratus. To test the function of bty, we cloned the orthologous zebrafish gene from a kidney cDNA library. Whole-mount in situ hybridization of zebrafish embryos showed that bty mRNA was present throughout development and, after the mid-blastula transition, was expressed in the head and in or near the site of primitive erythropoiesis in the tail just prior to red cell production. One- to four-cell embryos injected with two distinct antisense morpholino oligonucleotides (MOs) targeted to the 5'-end of the bty mRNA failed to develop red cells, whereas embryos injected with 4- and 5-bp mismatch control MOs produced wild-type quantities of erythrocytes. The morphant phenotype was rescued by co-injection of synthetic bty mRNA containing an artificial 5'-untranslated region (UTR) with the antisense MO that bound the 5'-UTR of the wild-type bty transcript. Furthermore, the expression of genes that mark terminal erythroid differentiation was greatly reduced in the antisense-MO-treated embryos. We conclude that bty is likely to play a role in differentiation of the committed red cell progenitor.

Zebrafish slc4a2/ae2 anion exchanger: cDNA cloning, mapping, functional characterization, and localization

Although the zebrafish has been used increasingly for the study of pronephric kidney development, studies of renal ion transporters and channels of the zebrafish remain few. We report the cDNA cloning and characterization of the AE2 anion exchanger ortholog from zebrafish kidney, slc4a2/ae2. The ae2 gene in linkage group 2 encodes a polypeptide of 1,228 aa exhibiting 64% aa identity with mouse AE2a. The exon-intron boundaries of the zebrafish ae2 gene are nearly identical to those of the rodent and human genes. Whole-mount in situ hybridization detects ae2 mRNA in prospective midbrain as early as the five-somite stage, then later in the pronephric primordia and the forming pronephric duct, where it persists through 72 h postfertilization (hpf). Zebrafish Ae2 expressed in Xenopus laevis oocytes mediates Na(+)-independent, electroneutral (36)Cl(-)/Cl(-) exchange moderately sensitive to inhibition by DIDS, is inhibited by acidic intracellular pH and by acidic extracellular pH, but activated by (acidifying) ammonium and by hypertonicity. Zebrafish Ae2 also mediates Cl(-)/HCO(3)(-) exchange in X. laevis oocytes and accumulates in or near the plasma membrane in transfected HEK-293 cells. In 24-48 hpf zebrafish embryos, the predominant but not exclusive localization of Ae2 polypeptide is the apical membrane of pronephric duct epithelial cells. Thus Ae2 resembles its mammalian orthologs in function, mechanism, and acute regulation but differs in its preferentially apical expression in kidney. These results will inform tests of the role of Ae2 in zebrafish kidney development and function.

Radiation hybrid mapping of the species cytoplasm-specific (scsae) gene in wheat

Radiation hybrid (RH) mapping is based on radiation-induced chromosome breakage and analysis of chromosome segment retention or loss using molecular markers. In durum wheat (Triticum turgidum L., AABB), an alloplasmic durum line [(lo) durum] has been identified with chromosome 1D of T. aestivum L. (AABBDD) carrying the species cytoplasm-specific (scsae) gene. The chromosome 1D of this line segregates as a whole without recombination, precluding the use of conventional genome mapping. A radiation hybrid mapping population was developed from a hemizygous (lo) scsae--line using 35 krad gamma rays. The analysis of 87 individuals of this population with 39 molecular markers mapped on chromosome 1D revealed 88 radiation-induced breaks in this chromosome. This number of chromosome 1D breaks is eight times higher than the number of previously identified breaks and should result in a 10-fold increase in mapping resolution compared to what was previously possible. The analysis of molecular marker retention in our radiation hybrid mapping panel allowed the localization of scsae and 8 linked markers on the long arm of chromosome 1D. This constitutes the first report of using RH mapping to localize a gene in wheat and illustrates that this approach is feasible in a species with a large complex genome.

The feline major histocompatibility complex is rearranged by an inversion with a breakpoint in the distal class I region

In order to determine the genomic organization of the major histocompatibility complex (MHC) of the domestic cat (Felis catus), DNA probes for 61 markers were designed from human MHC reference sequences and used to construct feline MHC BAC contig map spanning ARE1 in the class II region to the olfactory receptor complex in the extended class I region. Selected BAC clones were then used to identify feline-specific probes for the three regions of the mammalian MHC (class II-class III-class I) for radiation hybrid mapping and fluorescent in situ hybridization to refine the organization of the domestic cat MHC. The results not only confirmed that the p-arm of domestic cat B2 is inverted relative to human Chromosome 6, but also demonstrated that one inversion breakpoint localized to the distal segment of the MHC class I between TRIM39 and TRIM26. The inversion thus disjoined the approximately 2.85 Mb of MHC containing class II-class III-class I (proximal region) from the approximately 0.50 Mb of MHC class I/extended class I region, such that TRIM39 is adjacent to the Chromosome B2 centromere and TRIM26 is adjacent to the B2 telomere in the domestic cat.

Characterization and localization of a human serine racemase

D-serine is present in the mammalian central nervous system, where it acts as one of the co-activators of N-methyl-D aspartate receptors. Synthesis of D-serine is catalyzed by the serine racemase enzyme. The current studies report on the isolation of a cDNA encoding a human serine racemase (SRR) from the human neuronal like cell line, NT2N. The SRR gene was localized on chromosome 17q13. The full-length cDNA has 1020 nucleotides which encode for a protein of 340 amino acids. The human protein shares 89% sequence identity with the mouse serine racemase. Human embryonic kidney 293 cells transiently transfected with this SRR gene were able to produce d-serine, indicating that the sequence encodes for an active enzyme. In Northern blot analysis the SRR mRNA was expressed in human brain, heart, skeletal muscle, kidney and liver tissues. Different splice forms of SRR were present in the peripheral tissues. Transcripts of at least three different sizes were present in heart and kidney, while in Western blot analysis multiple bands of different sizes were observed. Immunohistochemical studies, using a polyclonal anti-human serine racemase antibody, revealed a peripheral expression of serine racemase protein in human cardiac myocytes and convoluted tubules of the kidney. Experiments in non-human primate brain demonstrated the localization of SRR in amygdala nuclei, cortex, thalamus and hippocampus. Co-localization studies in the hippocampus demonstrated the exclusive expression of serine racemase in glial cells. The cloning of a functional human serine racemase and its expression in central nervous system of primates support a role for D-serine in neuronal activity. Furthermore, its presence in human periphery such as in heart and kidney suggest a potential biological role for D-serine in the regulation of N-methyl-D-aspartate (NMDA) receptor activity in these peripheral organs as well.

A high-resolution radiation hybrid map of chicken chromosome 5 and comparison with human chromosomes

BACKGROUND

The resolution of radiation hybrid (RH) maps is intermediate between that of the genetic and BAC (Bacterial Artificial Chromosome) contig maps. Moreover, once framework RH maps of a genome have been constructed, a quick location of markers by simple PCR on the RH panel is possible. The chicken ChickRH6 panel recently produced was used here to construct a high resolution RH map of chicken GGA5. To confirm the validity of the map and to provide valuable comparative mapping information, both markers from the genetic map and a high number of ESTs (Expressed Sequence Tags) were used. Finally, this RH map was used for testing the accuracy of the chicken genome assembly for chromosome 5.

RESULTS

A total of 169 markers (21 microsatellites and 148 ESTs) were typed on the ChickRH6 RH panel, of which 134 were assigned to GGA5. The final map is composed of 73 framework markers extending over a 1315.6 cR distance. The remaining 61 markers were placed alongside the framework markers within confidence intervals.

CONCLUSION

The high resolution framework map obtained in this study has markers covering the entire chicken chromosome 5 and reveals the existence of a high number of rearrangements when compared to the human genome. Only two discrepancies were observed in relation to the sequence assembly recently reported for this chromosome.

A high-resolution radiation hybrid map of chicken chromosome 5 and comparison with human chromosomes

Background

The resolution of radiation hybrid (RH) maps is intermediate between that of the genetic and BAC (Bacterial Artificial Chromosome) contig maps. Moreover, once framework RH maps of a genome have been constructed, a quick location of markers by simple PCR on the RH panel is possible. The chicken ChickRH6 panel recently produced was used here to construct a high resolution RH map of chicken GGA5. To confirm the validity of the map and to provide valuable comparative mapping information, both markers from the genetic map and a high number of ESTs (Expressed Sequence Tags) were used. Finally, this RH map was used for testing the accuracy of the chicken genome assembly for chromosome 5.

Results

A total of 169 markers (21 microsatellites and 148 ESTs) were typed on the ChickRH6 RH panel, of which 134 were assigned to GGA5. The final map is composed of 73 framework markers extending over a 1315.6 cR distance. The remaining 61 markers were placed alongside the framework markers within confidence intervals.

Conclusion

The high resolution framework map obtained in this study has markers covering the entire chicken chromosome 5 and reveals the existence of a high number of rearrangements when compared to the human genome. Only two discrepancies were observed in relation to the sequence assembly recently reported for this chromosome.

A web server for performing electronic PCR

'Electronic PCR' (e-PCR) refers to a computational procedure that is used to search DNA sequences for sequence tagged sites (STSs), each of which is defined by a pair of primer sequences and an expected PCR product size. To gain speed, our implementation extracts short 'words' from the 3' end of each primer and stores them in a sorted hash table that can be accessed efficiently during the search. One recent improvement is the use of overlapping discontinuous words to allow matches to be found despite the presence of a mismatch. Moreover, it is possible to allow gaps in the alignment between the primer and the sequence. The effect of these changes is to improve sensitivity without significantly affecting specificity. The new software provides a search mode using a query STS against a sequence database to augment the previously available mode using a query sequence against an STS database. Finally, e-PCR may now be used through a web service, with search results linked to other web resources such as the UniSTS database and the MapViewer genome browser. The e-PCR web server may be found at www.ncbi.nlm.nih.gov/sutils/e-pcr.

The use of DNA polymorphisms in genetic mapping

The introduction of molecular markers has revolutionized genetics. The range of polymorphisms that are available is increasing and the advent of large-scale cDNA and genomic sequencing is a source of an ever-increasing set of available markers. The ease with which any particular marker type can be applied to an experimental system depends, to some extent, on the amount of genomic information available for that system. However, comparative genomics is enabling a wider range of marker technology to be applied to relatively information-poor systems. The types of markers that are available include restriction fragment length polymorphisms, amplified fragment length polymorphisms, ransom amplified polymorphic DNAs, simple sequence repeats, single nucleotide polymorphisms and small insertions/deletions. The types of questions that can be addressed with these molecular markers include the generation of genetic and physical maps for the identification of interesting loci, the development of marker-based gene tags, map-based cloning of agronomically important genes, synteny mapping, marker-assisted selection and quantitative trait analysis. The continued development of technology including new high throughput methods, for example those being applied to single nucleotide polymorphisms, will change the ease with which current questions can be answered as well as enable new analyses that are presently impossible to undertake.

The first-generation whole-genome radiation hybrid map in the horse identifies conserved segments in human and mouse genomes

A first-generation radiation hybrid (RH) map of the equine (Equus caballus) genome was assembled using 92 horse x hamster hybrid cell lines and 730 equine markers. The map is the first comprehensive framework map of the horse that (1) incorporates type I as well as type II markers, (2) integrates synteny, cytogenetic, and meiotic maps into a consensus map, and (3) provides the most detailed genome-wide information to date on the organization and comparative status of the equine genome. The 730 loci (258 type I and 472 type II) included in the final map are clustered in 101 RH groups distributed over all equine autosomes and the X chromosome. The overall marker retention frequency in the panel is approximately 21%, and the possibility of adding any new marker to the map is approximately 90%. On average, the mapped markers are distributed every 19 cR (4 Mb) of the equine genome--a significant improvement in resolution over previous maps. With 69 new FISH assignments, a total of 253 cytogenetically mapped loci physically anchor the RH map to various chromosomal segments. Synteny assignments of 39 gene loci complemented the RH mapping of 27 genes. The results added 12 new loci to the horse gene map. Lastly, comparison of the assembly of 447 equine genes (256 linearly ordered RH-mapped and additional 191 FISH-mapped) with the location of draft sequences of their human and mouse orthologs provides the most extensive horse-human and horse-mouse comparative map to date. We expect that the foundation established through this map will significantly facilitate rapid targeted expansion of the horse gene map and consequently, mapping and positional cloning of genes governing traits significant to the equine industry.

The comprehensive mouse radiation hybrid map densely cross-referenced to the recombination map: a tool to support the sequence assemblies

We have developed a unique comprehensive mouse radiation hybrid (RH) map of nearly 23,000 markers integrating data from three international genome centers and over 400 independent laboratories. We have cross-referenced this map to the 0.5-cM resolution recombination-based Jackson Laboratory (TJL) backcross panel map, building a complete set of RH framework chromosome maps based on a high density of known-ordered anchor markers. We have systematically typed markers to improve coverage and resolve discrepancies, and have reanalyzed data sets as needed. The cross-linking of the RH and recombination maps has resulted in a highly accurate genome-wide map with consistent marker order. We have compared these linked framework maps to the Ensemble mouse genome sequence assembly, and show that they are a useful medium resolution tool for both validating sequence assembly and elucidating chromosome biology.

Linkage relationships of genes coding for alpha2-macroglobulin, C3 and C4 in the zebrafish: implications for the evolution of the complement and Mhc systems

The alpha2-macroglobulin (A2M) and the complement components C3 and C4 are related proteins derived from a common ancestor. Theoretically, this derivation could have occurred either by tandem duplications of their encoding genes or by polyploidization involving chromosomal segments, a chromosome or the whole genome. In tetrapods the A2M-, C3- and C4-encoding genes are generally each located on a different chromosome. This observation has been interpreted as supporting their origin by polyploidization. We identified and mapped (with the help of a radiation hybrid panel of cell lines) the A2M, C3 and C4 loci in the zebrafish, Danio rerio. Each of the three types of loci is present in the zebrafish in multiple copies, but all of the identified copies of a given type map to the same region in linkage groups 1 (C3) and 15 (A2M, C4). The A2M and C4 loci are mapped in the same region not linked to any of the class I or class II major histocompatibility complex (Mhc) loci. These observations are interpreted as supporting the origin of the A2M family of genes by tandem duplications, followed by the dispersal of the copies to different chromosomes. It is also argued that the association of C4 with the class I/II loci in tetrapods is accidental and without functional significance.

A candidate gene study in low HDL-cholesterol families provides evidence for the involvement of the APOA2 gene and the APOA1C3A4 gene cluster

In patients with premature coronary heart disease, the most common lipoprotein abnormality is high-density lipoprotein (HDL) deficiency. To assess the genetic background of the low HDL-cholesterol trait, we performed a candidate gene study in 25 families with low HDL, collected from the genetically isolated population of Finland. We studied 21 genes encoding essential proteins involved in the HDL metabolism by genotyping intragenic and flanking markers for these genes. We found suggestive evidence for linkage in two candidate regions: Marker D1S2844, in the apolipoprotein A-II (APOA2) region, yielded a LOD score of 2.14 and marker D11S939 flanking the apolipoprotein A-I/C-III/A-IV gene cluster (APOA1C3A4) produced a LOD score of 1.69. Interestingly, we identified potential shared haplotypes in these two regions in a subset of low HDL families. These families also contributed to the obtained positive LOD scores, whereas the rest of the families produced negative LOD scores. None of the remaining candidate regions provided any evidence for linkage. Since only a limited number of loci were tested in this candidate gene study, these LOD scores suggest significant involvement of the APOA2 gene and the APOA1C3A4 gene cluster, or loci in their immediate vicinity, in the pathogenesis of low HDL.

Developmental expression of functional cyclooxygenases in zebrafish

Study of the cyclooxygenases (COXs) has been limited by the role of COX-2 in murine reproduction and renal organogenesis. We sought to characterize COX expression and function in zebrafish (z). Full-length cDNAs of zCOX-1 and zCOX-2 were cloned and assigned to conserved regions of chromosomes 5 and 2, respectively. The deduced proteins are 67% homologous with their human orthologs. Prostaglandin (PG) E(2) is the predominant zCOX product detected by mass spectrometry. Pharmacological inhibitors demonstrate selectivity when directed against heterologously expressed zCOX isoforms. Zebrafish thrombocyte aggregation ex vivo and hemostasis in vivo are sensitive to inhibition of zCOX-1, but not zCOX-2. Both zCOXs were widely expressed during development, and knockdown of zCOX-1 causes growth arrest during early embryogenesis. zCOX-1 is widely evident in the embryonic vasculature, whereas zCOX-2 exhibits a more restricted pattern of expression. Both zCOX isoforms are genetically and functionally homologous to their mammalian orthologs. The zebrafish affords a tractable model system for the study of COX biology and development.

Molecular characterization of hiwi, a human member of the piwi gene family whose overexpression is correlated to seminomas

The piwi family genes are highly conserved during evolution and play essential roles in stem cell self-renewal, gametogenesis, and RNA interference in diverse organisms ranging from Drosophila melanogaster and C. elegans to Arabidopsis. Here we report the molecular characterization of hiwi, a human member of the piwi gene family. hiwi maps to the long arm of chromosome 12, band 12q24.33, a genomic region that displays genetic linkage to the development of testicular germ cell tumors of adolescents and adults (TGCTs), i.e., seminomas and nonseminomas. In addition, gain of this chromosomal region has been found in some TGCTs. hiwi encodes a 3.6 kb mRNA that is expressed abundantly in the adult testis. It encodes a highly basic 861-amino-acid protein that shares significant homology throughout its entire length with other members of the PIWI family proteins in Drosophila, C. elegans and mammals. In normal human testes, hiwi is specifically expressed in germline cells, with its expression detectable in spermatocytes and round spermatids during spermatogenesis. No expresssion was observed in testicular tumors of somatic origin, such as Sertoli cell and Leydig cell tumors. Enhanced expression was found in 12 out of 19 sampled testicular seminomas-tumors originating from embryonic germ cells with retention of germ cell phenotype. In contrast, no enhanced expression was detected in 10 nonseminomas-testicular tumors that originate from the same precursor cells as seminomas yet have lost their germ cell characteristics. Finally, no enhanced expression was detected in four spermatocytic seminomas-testicular tumors that most likely originate from germ cells capable of partial meiosis. Thus, hiwi is specifically expressed in both normal and malignant spermatogenic cells in a maturation stage-dependent pattern, in which it might function in germ cell proliferation.

Systematic evaluation of map quality: human chromosome 22

Marker positions on nine genetic linkage, radiation hybrid, and integrated maps of human chromosome 22 were compared with their corresponding positions in the completed DNA sequence. The proportion of markers whose map position is <250 kb from their respective sequence positions ranges from 100% to 35%. Several discordant markers were identified, as well as four regions that show common inconsistencies across multiple maps. These shared discordant regions surround duplicated DNA segments and may indicate mapping or assembly errors due to sequence homology. Recombination-rate distributions along the chromosome were also evaluated, with male and female meioses showing significantly different patterns of recombination, including an 8-Mb male recombination desert. The distributions of radiation-induced chromosome breakage for the GB4 and the G3 radiation hybrid panels were also evaluated. Both panels show fluctuations in breakage intensity, with different regions of significantly elevated rates of breakage. These results provide support for the common assumption that radiation-induced breaks are generally randomly distributed. The present studies detail the limitations of these important map resources and should prove useful for clarifying potential problems in the human maps and sequence assemblies, as well as for mapping and sequencing projects in and across other species.

Targeted construction of a high-resolution, integrated, comprehensive, and comparative map for a region specific to bovine chromosome 6 based on radiation hybrid mapping

To resolve a candidate chromosome region on the middle part of bovine chromosome 6 (BTA6) containing several different quantitative trait locus (QTL) intervals, we constructed a high-resolution, integrated, comprehensive, and comparative map using a 12,000-rad, whole-genome, cattle-hamster radiation hybrid (RH) panel. The RH map includes a total of 71 loci either selected from bovine and comparative maps or targeted directly from a microdissection library specific for the BTA6 region. All loci typed were placed in one linkage group at a lod score threshold of 4.0. The length of the comprehensive RH map, which is the first high-resolution RH map in cattle, spans 2568.8 cR(12,000). The order of markers obtained principally agrees with the order on published bovine genetic maps. Our RH map integrates markers as well as genes and ESTs available from several physical and genetic maps of BTA6 and the orthologous ovine chromosome 6, human chromosome 4, and mouse chromosomes 5/3. Comparative analysis confirms and refines current knowledge about conservation and rearrangements in corresponding chromosomal regions on BTA6. We identified and localized two new breakpoints for intrachromosomal rearrangements between human chromosome 4 and BTA6. This RH map is a powerful tool in all aspects of genetic, physical, transcript, and comparative mapping. Due to its links to the gene-dense maps of human and mouse, it can serve as a prerequisite to identify possible candidate genes for quantitative trait loci localized in the targeted BTA6 region.

Cloning, mapping and expression analysis of C15orf4, a novel human gene with homology to the yeast mitochondrial ribosomal protein Ym130 gene

We have identified C15orf4, a novel human gene showing homology to the yeast mitochondrial ribosomal protein YmL30. C15orf4 encodes a transcript of 1,006 nt with an ORF of 279 amino acids and a predicted protein size of 31.7 kDa. Expression of C15orf4 is enriched in testis. C15orf4 was positioned to chromosome 15q24 by radiation hybrid mapping. We have identified the C15orf4 mouse orthologue as well as homologues in other species.