Genomic analysis and marker development for the Tsn1 locus in wheat using bin-mapped ESTs and flanking BAC contigs

The wheat Tsn1 gene confers sensitivity to the host-selective toxin Ptr ToxA produced by the tan spot fungus (Pyrenophora tritici-repentis). The long-term goal of this research is to isolate Tsn1 using a positional cloning approach. Here, we evaluated 54 ESTs (expressed sequence tags) physically mapped to deletion bin 5BL 0.75–0.76, which is a gene-rich region containing Tsn1. Twenty-three EST loci were mapped as either PCR-based single-stranded conformational polymorphism or RFLP markers in a low-resolution wheat population. The genetic map corresponding to the 5BL 0.75–0.76 deletion bin spans 18.5 cM and contains 37 markers for a density of 2 markers/cM. The EST-based genetic map will be useful for tagging other genes, establishing colinearity with rice, and anchoring sequence ready BAC contigs of the 5BL 0.75–0.76 deletion bin. High-resolution mapping showed that EST-derived markers together with previously developed AFLP-derived markers delineated Tsn1 to a 0.8 cM interval. Flanking markers were used to screen the Langdon durum BAC library and contigs of 205 and 228 kb flanking Tsn1 were assembled, sequenced, and anchored to the genetic map. Recombination frequency averaged 760 kb/cM across the 228 kb contig, but no recombination was observed across the 205 kb contig resulting in an expected recombination frequency of more than 10 Mb/cM. Therefore, chromosome walking within the Tsn1 region may be difficult. However, the sequenced BACs allowed the identification of one microsatellite in each contig for which markers were developed and shown to be highly suitable for marker-assisted selection of Tsn1.     

High-resolution comparative physical mapping of mouse Chromosome 10 in the region of homology with human Chromosome 21

Comparative mapping of human and mouse chromosomes can be used to predict locations of homologous loci between the species, provides the substrate to examine the process of chromosomal evolution, and facilitates the continuing development of mouse genetic models for human disorders. A YAC contig of the region of mouse Chromosome (Chr) 10 (MMU10) that demonstrates conserved linkage with the distal portion of human Chr 21 (HSA21) has been constructed. The contig contains all known genes mapped in both species, defines the proximal region of homology between MMU10 and HSA22, and contains the evolutionary junction between HSA21 and HSA22 on MMU10. It consists of 23 YACs and 2 PACs, and covers 3.2 Mb of MMU10. The average marker density for this region is 1 marker/69 kb. Nine of 22 expressed sequences are mapped here for the first time in mouse, and two are newly characterized expressed sequences. The contig also contains 12 simple sequence repeats (SSRs) and 16 YAC and PAC endclone markers. YAC fragmentation analysis was used to create a physical map for the proximal 2.2 Mb of the contig. Cloning of the corresponding region of HSA21 has proven difficult, and the mouse contig includes segments absent from previously described sequence ready maps of HSA21. Received: 22 July 1998 / Accepted: 13 November 1998     

Sequence-ready physical map of the mouse Chromosome 16 region with conserved synteny to the human Velocardiofacial syndrome region on 22q11.2

Proximal mouse Chromosome (Chr) 16 shows conserved synteny with human Chrs 16, 8, 22, and 3. The mouse Chr 16/human Chr 22 conserved synteny region includes the DiGeorge/Velocardiofacial syndrome region of human Chr 22q11.2. A physical map of the entire mouse Chr 16/human Chr 22 region of conserved synteny has been constructed to provide a substrate for gene discovery, genomic sequencing, and animal model development. A YAC contig was constructed that extends ca. 5.4 Mb from a region of conserved synteny with human Chr 8 at Prkdc through the region conserved with human Chr 3 at DVL3. Sixty-one markers including 37 genes are mapped with average marker spacing of 90 kb. Physical distance was determined across the 2.6-Mb region from D16Mit74 to Hira with YAC fragmentation. The central region from D16Jhu28 to Igl-C1 was converted into BAC and PAC clones, further refining the physical map and providing sequence-ready template. The gene content and borders of three blocks of conserved linkage between human Chr 22q11.2 mouse Chr 16 are refined. Received: 4 November 1998 / Accepted: 21 December 1998     

Genomic organization and mapping of the human and mouse neuronal β2-nicotinic acetylcholine receptor genes

As a first step in determining whether there are polymorphisms in the nicotinic acetylcholine receptor (nAChR) genes that are associated with nicotine addiction, we isolated genomic clones of the β2-nAChR genes from human and mouse BAC libraries. Although cDNA sequences were available for the human gene, only the promoter sequence had been reported for the mouse gene. We determined the genomic structures by sequencing 12 kb of the human gene and over 7 kb of the mouse gene. While the sizes of exons in the mouse and human genes are the same, the introns differ in size. Both promoters have a high GC content (60–80%) proximal to the AUG and share a neural-restrictive silencer element (NRSE), but overall sequence identity is only 72%. Using a 6-Mb YAC contig of Chr 1, we mapped the human β2-nAChR gene, CHRNB2, to 1q21.3 with the order of markers cen, FLG, IVL, LOR, CHRNB2, tel. The mouse gene, Acrb2, had previously been mapped to Chr 3 in a region orthologous to human Chr 1. We refined mapping of the mouse gene and other markers on a radiation hybrid panel of Chr 3 and found the order cen, Acrb2, Lor, Iv1, Flg, tel. Our results indicate that this cluster of markers on human Chr 1 is inverted with respect to its orientation on the chromosome compared with markers in the orthologous region of mouse Chr 3. Received: 26 January 1999 / Accepted: 10 May 1999     

Characterization of a cluster comprising ∼100 odorant receptor genes in mouse

With ∼1000 genes, the odorant receptor (OR) gene repertoire is the largest gene family in the mouse genome. Here we have established a 129/Sv BAC contig for mouse OR gene cluster 7 (Olfr7) on Chromosome (Chr) 9. The assembled ∼2-Mb contig consists of 75 BACs and may contain as many as 100 OR genes, or ∼10% of the mouse repertoire. Facilitated by the lack of introns in the coding region, we have determined the nucleotide sequence of 37 full-length, 2 partial, and 3 pseudo coding regions. These 42 OR genes and 3 additional OR genes previously mapped to the mouse Olfr7 cluster can be organized into 13 classes based on OR probe cross-hybridizations with 129/Sv mouse genomic DNA. OR genes belonging to the same class tend to be located next to each other within the cluster. Comparison of published full-length mouse and rat OR coding sequences with those identified here shows that the Olfr7 OR genes are highly related to each other, clustering on two major branches of an unrooted phylogenetic tree. Eight ORs contain an unusual NXC sequon at the amino-terminal extracellular domain that may represent a novel N-linked glycosylation site. The BAC contig presented here provides the substrate for sequencing of the cluster. Received: 27 June 2000 / Accepted: 17 August 2000     

A BAC contig map of the Ly49 gene cluster in 129 mice reveals extensive differences in gene content relative to C57BL/6 mice

The murine Ly49 gene family is functionally analogous to the human killer cell Ig-like receptor (KIR) family of class I major histocompatibility complex (MHC) receptors. The number of KIR genes varies dramatically between individuals; however, the organization of the Ly49 genes has only been determined for the C57BL/6 (B6) mouse. The organization of the 129 Ly49 loci was determined from a BAC contig map by PCR and Southern blot analysis. In addition to the 10 Ly49 genes known from previous studies of the 129/J strain, 8 new genes were localized to the 129 Ly49 cluster. A gene order of Ly49q(1), e, (v, q(2)), e/c(2), l/r, s, t, e/c(1), r, u, u/i, i(1), g, p/d, (i(2), p), and o was determined. The 129 Ly49 gene cluster is predicted to span approximately 600 kb. These results indicate that Ly49 gene numbers can be significantly different between inbred mouse strains, analogous to the haplotype differences observed in the human KIR genes.     

A sequence-ready map of the human chromosome 1q telomere

A 260-kb half-YAC clone derived from human chromosome 1q was mapped at high resolution using cosmid subclone fingerprint analysis and was integrated with overlapping clones from the telomeric end of a separately derived 1q44 BAC contig to create a sequence-ready map extending to the molecular telomere of 1q. Analysis of 100 kb of sample sequences from across the 260-kb region encompassed by the half-YAC revealed the presence of EST sequence matches corresponding to 12 separate Unigene clusters and to 12 separate unclustered EST sequences. Low-copy subtelomeric repeats typical of many human telomere regions are present within the distal-most 30 kb of 1q. The previously isolated and radiation hybrid-mapped markers Bda84F03, 1QTEL019, and WI11861 localized at distances approximately 32, 88, and 99 kb, respectively, from the 1q terminus. This sequence-ready map permits high-resolution integration of genetic maps with the DNA sequences directly adjacent to the tip of human chromosome 1q and will enable telomeric closure of the human chromosome 1q DNA reference sequence by connecting the molecular 1q telomere to an internal BAC contig.     

Physical mapping across an avirulence locus of Phytophthora infestans using a highly representative, large-insert bacterial artificial chromosome library

The oomycete plant pathogen Phytophthora infestans is the causal agent of late blight, one of the most devastating diseases of potato worldwide. As part of efforts to clone avirulence (Avr) genes and pathogenicity factors from P. infestans, we have constructed a bacterial artificial chromosome (BAC) library from an isolate containing six Avr genes. The BAC library comprises clones with an average insert size of 98 kb and represents an estimated 10 genome equivalents. A three-dimensional pooling strategy was developed to screen the BAC library for amplified fragment length polymorphism (AFLP) markers, as this type of marker has been extensively used in construction of a P. infestans genetic map. Multiple positive clones were identified for each AFLP marker tested. The pools were used to construct a contig of 11 BAC clones in a region of the P. infestans genome containing a cluster of three avirulence genes. The BAC contig is predicted to encompass the Avr11 locus but mapping of the BAC ends will be required to determine if the Avr3 and Avr10 loci are also present in the BAC contig. These results are an important step towards the positional cloning of avirulence genes from P. infestans, and the BAC library represents a valuable resource for largescale studies of oomycete genome organisation and gene content.     

Identification of 129S6/SvEvTac-specific polymorphisms on mouse chromosome 11

Polymorphisms such as single-nucleotide polymorphisms (SNPs) and insertions/deletions (Indels) can be associated with phenotypic traits and be used as markers for disease diagnosis. Identification of these genetic variations within laboratory mice is crucial to improve our understanding of the genetic background of the mice used for research. As part of a positional cloning project, we sequenced six genes (Mettl16, Evi2a, Psmd11, Cct6d, Rffl, and Ap2b1) within a 6.8-Mb domain of mmu chr 11 in the C57BL/6J and 129S6/SvEvTac inbred strains. Although 129S6/SvEvTac is widely used in the mouse community, there is very little current (or projected future) sequence information available for this strain. We identified 6 Indels and 21 novel SNPs and confirmed genotype information for 114 additional SNPs in these 6 genes. Mettl16 and Ap2b1 contained the largest numbers of variants between the C57BL/6J and 129S6/SvEvTac strains. In addition, we found five new SNPs between 129S6/SvEvTac and 129S1/SvImJ within the Ap2b1 locus. Although we did not detect differences between C57BL/6J and 129S6/SvEvTac within Evi2a, this locus contains a relatively high SNP density compared with the surrounding sequence. Our study highlights the genetic differences among three inbred mouse strains (C57BL/6J, 129S6/SvEvTac, and 129S1/SvImJ) and provides valuable sequence information that can be used to track alleles in genomics-based studies.     

Sequence and characterization of the Ig heavy chain constant and partial variable region of the mouse strain 129S1

Although the entire mouse genome has been sequenced, there remain challenges concerning the elucidation of particular complex and polymorphic genomic loci. In the murine Igh locus, different haplotypes exist in different inbred mouse strains. For example, the Igh(b) haplotype sequence of the Mouse Genome Project strain C57BL/6 differs considerably from the Igh(a) haplotype of BALB/c, which has been widely used in the analyses of Ab responses. We have sequenced and annotated the 3' half of the Igh(a) locus of 129S1/SvImJ, covering the C(H) region and approximately half of the V(H) region. This sequence comprises 128 V(H) genes, of which 49 are judged to be functional. The comparison of the Igh(a) sequence with the homologous Igh(b) region from C57BL/6 revealed two major expansions in the germline repertoire of Igh(a). In addition, we found smaller haplotype-specific differences like the duplication of five V(H) genes in the Igh(a) locus. We generated a V(H) allele table by comparing the individual V(H) genes of both haplotypes. Surprisingly, the number and position of D(H) genes in the 129S1 strain differs not only from the sequence of C57BL/6 but also from the map published for BALB/c. Taken together, the contiguous genomic sequence of the 3' part of the Igh(a) locus allows a detailed view of the recent evolution of this highly dynamic locus in the mouse.     

An integrated physical, genetic and cytogenetic map around the sunn locus of Medicago truncatula.

The sunn mutation of Medicago truncatula is a single-gene mutation that confers a novel supernodulation phenotype in response to inoculation with Sinorhizobium meliloti. We took advantage of the publicly available codominant PCR markers, the high-density genetic map, and a linked cytogenetic map to define the physical and genetic region containing sunn. We determined that sunn is located at the bottom of linkage group 4, where a fine-structure genetic map was used to place the locus within a approximately 400-kb contig of bacterial artificial chromosome (BAC) clones. Genetic analyses of the sunn contig, as well as of a second, closely linked BAC contig designated NUM1, indicate that the physical to genetic distance within this chromosome region is in the range of 1000 -1100 kb.cM-1. The ratio of genetic to cytogenetic distance determined across the entire region is 0.3 cM.microm(-1). These estimates are in good agreement with the empirically determined value of approximately 300 kb.microm(-1) measured for the NUM1 contig. The assignment of sunn to a defined physical interval should provide a basis for sequencing and ultimately cloning the responsible gene.     

Generation of a 5.5-Mb BAC/PAC contig of pig chromosome 6q1.2 and its integration with existing RH, genetic and comparative maps

We generated a sequence-ready BAC/PAC contig spanning approximately 5.5 Mb on porcine chromosome 6q1.2, which represents a very gene-rich genome region. STS content mapping was used as the main strategy for the assembly of the contig and a total of 6 microsatellite markers, 53 gene-related STS and 116 STS corresponding to BAC and PAC end sequences were analyzed. The contig comprises 316 BAC and PAC clones covering the region between the genes GPI and LIPE. The correct contig assembly was verified by RH-mapping of STS markers and comparative mapping of BAC/PAC end sequences using BLAST searches. The use of microsatellite primer pairs allowed the integration of the physical maps with the genetic map of this region. Comparative mapping of the porcine BAC/PAC contig with respect to the gene-rich region on the human chromosome 19q13.1 map revealed a completely conserved gene order of this segment, however, physical distances differ somewhat between HSA19q13.1 and SSC6q1.2. Three major differences in DNA content between human and pig are found in two large intergenic regions and in one region of a clustered gene family, respectively. While there is a complete conservation of gene order between pig and human, the comparative analysis with respect to the rodent species mouse and rat shows one breakpoint where a genome segment is inverted.     

A bacterial artificial chromosome library for the Chinese alligator (Alligator sinensis)

Chinese alligator (Alligator sinensis) is a rare and endangered species endemic to China. To better understand genetic details of the Chinese alligator genomic structure, a highly redundant bacterial artificial chromosome (BAC) library was constructed. This library consists of 216,238 clones with an average insert size of about 90kb, indicating that the library contains 6.8-fold genome equivalents. Subsequently, we constructed a 516kb contig map for the Chinese alligator olfactory receptor (OR) genes, which spans nine BAC clones, and subjected the BACs to full sequencing. The sequence analysis revealed that this contig contained 16 OR functional genes and meanwhile demonstrated that the nine BACs, which constituted the contig, overlapped correctly, proving the usability of this genome library. As a result, this BAC library could provide a useful platform for physical mapping, genome sequencing or complex analysis of targeted genomic regions for this rare species.     

Chromosome Walking in the Petunia Inflata Self-Incompatibility (S-) Locus and Gene Identification in an 881-kb Contig Containing S 2 -RNase

Self-incompatibility (SI) in the Solanaceae, Rosaceae and Scrophulariaceae is controlled by the polymorphic S locus, which contains two separate genes encoding pollen and pistil determinants in SI interactions. The S-RNase gene encodes the pistil determinant, whereas the pollen determinant gene, named the pollen S gene, has not yet been identified. Here, we set out to construct an integrated genetic and physical map of the S locus of Petunia inflata and identify any additional genes located at this locus. We first conducted chromosome walking at the S2 locus using BAC clones that contained either S2-RNase or one of the nine markers tightly linked to the S locus. Ten separate contigs were constructed, which collectively spanned 4.4 Mb. To identify additional genes located at the S2 locus, a 328-kb region (part of an 881-kb BAC contig) containing S2-RNase was completely sequenced. Approximately 76% of the region contained repetitive sequences, including transposon-like sequences. Other than S2-RNase, an F-box gene, named PiSLF2 (S2-allele of P. inflata S-locus F-box gene), was the only predicted gene whose deduced amino acid sequence was similar to the sequences of known proteins in the database. Two different cDNA selection methods were used to identify additional genes in the 881-kb contig; 11 groups of cDNA clones were identified in addition to those for S2-RNase and PiSLF2. RT-PCR analysis of expression profiles and PCR analysis of BAC clones and genomic DNA confirmed that seven of these 11 newly identified genes were located in the 881-kb contig.     

Assessment of the gene content of the chromosomal regions flanking bovine DGAT1

As a first step towards verifying the candidate status of DGAT1 as the causal gene for milk fat percentage in cattle, we constructed a bovine BAC contig spanning 576 kb of the chromosomal region containing DGAT1. High content of NotI sites (21 within the contig) indicated that the region is gene-rich. Twenty-three genes neighboring DGAT1 were mapped, including two bovine cDNA sequences that have no orthologous sequences within the NCBI sequence databases. On average, 2015 bp for each of the 23 neighboring genes were sequenced and entered into EMBL. Likewise, 10 new STS markers were established by BAC-end sequencing. Within the genes and STS markers, 55 polymorphisms were discovered. These will form the basis of future linkage disequilibrium studies to test whether any genes neighboring DGAT1 are associated with variation in milk fat percentage, thereby testing the candidate status of DGAT1.     

Assignment of 280 swine genomic inserts including 31 microsatellites from BAC clones to the swine RH map (IMpRH map)

A precise genetic map containing anonymous markers and genes is indispensable for the efficient selection of candidate gene(s) responsible for quantitative trait loci (QTL) traits. For this purpose, a first version of a radiation hybrid cell (RH) map has been constructed by using the INRA-University of Minnesota RH panel for 757 markers (IMpRH) (Hawken et al. 1999, Mamm. Genome 10: 824–830). In this study, 280 swine genomic fragments in BAC clones were assigned to the IMpRH map; 255 BAC clones were successfully linked to first-generation linkage groups (LOD > 4.8). The remaining 25 clones could not be mapped, because their lod-scores to the closest markers in the first generation map were less than 4.8. In addition, 16 BAC clones, mapped to swine Chromosome (Chr) 1 by IMpRH mapping, were subjected to isolation of microsatellites (MSs). Thirty-one MSs were isolated from 15 BAC clones, and 24 of 31 (77%) MSs derived from 14 clones were found to be polymorphic. We also mapped both termini of 12 BAC clones to the IMpRH map, in order to measure resolution of the IMpRH map; the resolution was found to range from 8 kb/centiRay to more than 126 kb/centiRay depending on the region. Received: 21 June 2001 / Accepted: 28 September 2001     

Recent developments from the Leishmania genome project

A first generation cosmid contig map of the Leishmania major Friedlin genome has been constructed, and genomic sequencing is well underway. Chromosome 1 (Chr1) and Chr3 have been completely sequenced, and Chr4 is virtually complete. Sequencing of several other chromosomes is in progress and the complete genome sequence may be available as soon as 2003. More than 600 completely sequenced new genes have been identified, representing approximately 8% of the total gene complement (approximately 8,600 genes) of Leishmania. Notably, a large proportion (approximately 69%) of the genes remain unclassified, with 40% of these being potentially Leishmania- (or kinetoplastid-) specific. Most interestingly, the genes are organized into large (>100-300 kb) polycistronic clusters of adjacent genes on the same DNA strand. Chr1 contains two such clusters organized in a 'divergent' manner, whereas Chr3 contains two 'convergent' clusters, with a single 'divergent' gene at one telomere, with the two large clusters separated by a tRNA gene. Statistical analyses of Chr1 show that the 'divergent junction' region between the two polycistronic gene clusters may be a candidate for an origin of DNA replication.