Long-Read Sequencing Genome Assembly of Ceratocystis fimbriata Enables Development of Molecular Diagnostics for Sweetpotato Black Rot.

Ceratocystis fimbriata is a destructive fungal pathogen of sweetpotato (Ipomoea batatas) that leads to losses at all stages of sweetpotato production. Accurate detection of C. fimbriata would allow for more efficient deployment of management tactics in sweetpotato production. To develop a diagnostic assay, a hybrid genome assembly of C. fimbriata isolate AS236 was generated. The resulting 31.7-MB assembly was near-chromosome level, with 18 contigs, 6,481 predicted genes, and a BUSCO completion score of 98.4% when compared with the fungus-specific lineage database. Additional Illumina DNA reads from C. manginecans, C. platani, and a second C. fimbriata isolate (C1421) were then mapped to the assembled genome using BOWTIE2 and counted using HTSeq, which identified 148 genes present only within C. fimbriata as molecular diagnostic candidates; 6 single-copy and 35 highly multi-copy (>40 BLAST hits), as determined through a self-BLAST-P alignment. Primers for PCR were designed in the 200-bp flanking region of the first exon for each candidate, and the candidates were validated against a diverse DNA panel containing Ceratocystis species, sweetpotato pathogens, and plants. After validation, two diagnostic candidates amplified only C. fimbriata DNA and were considered to be highly specific to the species. These genetic markers will serve as valuable diagnostic tools with multiple applications including the detection of C. fimbriata in seed, soil, and wash water in sweetpotato production.

Development, Validation, and Utility of Species-Specific Diagnostic Markers for Detection of Peronospora belbahrii.

Peronospora belbahrii is an oomycete and the cause of basil downy mildew, one of the most destructive diseases affecting basil production worldwide. Disease management is challenging due to wind-dispersed sporangia and contaminated seed; therefore, identifying P. belbahrii in seed lots before sale or planting or in the field before symptoms develop could allow for timely deployment of disease management strategies. In this study, a draft genome assembly and next-generation sequencing reads for P. belbahrii, as well as publicly available DNA-seq and RNA-seq reads of several other downy mildew pathogens, were incorporated into a bioinformatics pipeline to predict P. belbahrii-specific diagnostic markers. The specificity of each candidate marker was validated against a diverse DNA collection of P. belbahrii, host tissue, and related oomycetes using PCR. Two species-specific markers were identified and used as templates to develop a highly sensitive probe-based real-time quantitative PCR (qPCR) assay that could detect P. belbahrii in leaf tissue and seed samples. Both markers were capable of reliably detecting as low as 500 fg/µl of P. belbahrii genomic DNA and as few as 10 sporangia. The qPCR assay was then validated with seed samples collected from a basil cultivar experiment. In total, 48 seed samples were collected and tested; P. belbahrii was detected in samples of all cultivars at estimated concentrations of 600 fg/µl up to 250 pg/µl and at as few as 10 sporangia up to >1,000 sporangia. The markers and assays are valuable for diagnostics and identifying P. belbahrii-contaminated seed lots to mitigate the effects of future basil downy mildew epidemics.

Genome Sequencing and Transcriptome Analysis of the Hop Downy Mildew Pathogen Pseudoperonospora humuli Reveal Species-Specific Genes for Molecular Detection.

Pseudoperonospora humuli is an obligate oomycete pathogen of hop (Humulus lupulus) that causes downy mildew, an important disease in most production regions in the Northern Hemisphere. The pathogen can cause a systemic infection in hop, overwinter in the root system, and infect propagation material. Substantial yield loss may occur owing to P. humuli infection of strobiles (seed cones), shoots, and cone-bearing branches. Fungicide application and cultural practices are the primary methods to manage hop downy mildew. However, effective, sustainable, and cost-effective management of downy mildew can be improved by developing early detection systems to inform on disease risk and timely fungicide application. However, no species-specific diagnostic assays or genomic resources are available for P. humuli. The genome of the P. humuli OR502AA isolate was partially sequenced using Illumina technology and assembled with ABySS. The assembly had a minimum scaffold length of 500 bp and an N50 (median scaffold length of the assembled genome) of 19.2 kbp. A total number of 18,656 genes were identified using MAKER standard gene predictions. Additionally, transcriptome assemblies were generated using RNA-seq and Trinity for seven additional P. humuli isolates. Bioinformatics analyses of next generation sequencing reads of P. humuli and P. cubensis (a closely related sister species) identified 242 candidate species-specific P. humuli genes that could be used as diagnostic molecular markers. These candidate genes were validated using polymerase chain reaction against a diverse collection of isolates from P. humuli, P. cubensis, and other oomycetes. Overall, four diagnostic markers were found to be uniquely present in P. humuli. These candidate markers identified through comparative genomics can be used for pathogen diagnostics in propagation material, such as rhizomes and vegetative cuttings, or adapted for biosurveillance of airborne sporangia, an important source of inoculum in hop downy mildew epidemics.

Serum NGAL can act as an early renal safety biomarker during long-term nucleos(t)ide analogue antiviral therapy in chronic hepatitis B.

Tubular renal toxicity is a side-effect of long-term therapy with nucleos(t)ide analogue(s) (NA) in chronic hepatitis B (CHB). There are no established surrogate markers in plasma of early NA-related toxicity. Neutrophil gelatinase-associated lipocalin (NGAL) is a protein produced by tubular cells following renal damage. We aimed therefore to retrospectively compare conventional renal markers (estimated glomerular filtration rates (eGFR) and urinary protein/creatinine ratio uPCR) with a sensitive biomarker (NGAL) in CHB patients on long-term NA therapy and assess the ability of new markers to predict NA-related renal toxicity (new onset of nonalbumin proteinuria). A total of 192 naïve CHB patients (median age 41 years, 78% males, 25% HBeAg+, 35% cirrhosis) were NA treated for at least 5 years (median 8.34 years, range 5.54-11.1 years). The eGFR and uPCR were compared at baseline and last clinical visit with serum NGAL concentrations measured by ELISA at same time-points and assessed according to the presence/absence of nonalbumin proteinuria at last visit. While baseline and last visit eGFR were similar (median:78 vs 84 mL/min), serum NGAL concentrations increased during therapy (median:9.4 vs 16.4 ng/mL, P < .05). The proportion of patients with proteinuria (uPCR > 15) increased between baseline and last visit (4.6% vs 21.4%, P < .05), with 30 (16%) patients having de novo nonalbumin proteinuria at last visit. High baseline NGAL concentrations were exclusive to patients with de novo nonalbumin proteinuria (median:31.7 vs 7.8 ng/mL, P < .01) and baseline NGAL levels >25 mg/mL were predictive of nonalbumin proteinuria at last visit (AUROC = 0.813). In conclusion, serum NGAL can act as a surrogate marker of early renal injury (de novo nonalbumin proteinuria) in CHB on long-term NA therapy.

Effects of testosterone on the transgender male voice.

There is sparse prospective data on the effects of testosterone therapy on the voices of transgender men (also referred to as trans men or female-to-male transsexuals). Our aim was to investigate the timing and degree of voice deepening over 12 months among an ethnically diverse sample of transgender men. This was a prospective 12-month study at an academic outpatient endocrinology clinic and speech and voice center. The participants were seven transgender men naïve to testosterone therapy. All patients received two voice assessments at baseline and one assessment at 3, 6, 9, and 12 months while on intramuscular testosterone esters. Serum testosterone and estradiol concentrations were measured at baseline and every 3 months. All seven transgender men reached a cisgender male mean fundamental frequency (MF0) within 6 months of testosterone therapy and four continued to experience a decrease after 6 months. The mean decrease in MF0 after 12 months of testosterone therapy was 6.4 semitones (49 Hz). Several patterns emerged regarding the extent and timing. For example, some participants showed no decrease in MF0 within the first 3 months of testosterone therapy, whereas others showed their greatest decrease in MF0. We concluded that transgender men who start testosterone therapy display different patterns of voice lowering. Clinicians should counsel transgender men that they may or may not experience voice lowering within the first 3 months of testosterone therapy and that the majority of voice deepening will occur within 6-9 months.

Detection of the NS3 Q80K polymorphism by Sanger and deep sequencing in hepatitis C virus genotype 1a strains in the UK.

The Q80K polymorphism in the hepatitis C virus (HCV) NS3 enzyme reduces susceptibility to simeprevir and other novel protease inhibitors. The aims of this study were to determine the prevalence of Q80K in treatment-naïve HCV-1a carriers in the North West region (NW) and South East region (SE) of England, investigate the occurrence of Q80K as a minority variant, and characterize viral phylogeny. Plasma samples from subjects who were naïve to anti-HCV therapy were subjected to conventional (Sanger) and deep (Illumina-Miseq, 1% interpretative cut-off) sequencing of NS3. Q80K occurred in 44 of 238 subjects (18.5%, 95% CI 13.6-23.4%), including 19 of 70 (27.1%) in the NW and 25 of 168 (14.9%) in the SE (p 0.0425), with no difference in HCV RNA load or human immunodeficiency virus (HIV) status. Q80K frequencies in reads of samples subjected to Illumina sequencing were >40% in all cases. Among subjects with Q80K, five of 44 (11.4%) showed one additional major resistance-associated mutation in NS3, detected at frequencies of >10% (V36L and V55A) or <10% (V36M). Phylogenetic analyses identified the two recognized HCV-1a lineages with (clade I) and without (clade II) Q80K. Overall, 148 of 238 (62.2%) sequences occurred within regional or inter-regional clusters, each comprising 3-20 sequences. There was no unique clustering of English sequences relative to strains from continental Europe and North America. In conclusion, Q80K was found at a high prevalence among treatment-naïve HCV-1a carriers in England, and was reliably detected by conventional sequencing, with no increased detection by deep sequencing. English sequences were highly interspersed with sequences from elsewhere in Europe (clade II) and North America (clade I), and their phylogeny was consistent with multiple introductions from different areas.

Intersex in the clam Scrobicularia plana (Da Costa): Widespread occurrence in English Channel estuaries and surrounding areas.

Estuarine clams Scrobicularia plana were sampled from 108 intertidal locations around the English Channel and adjacent areas. Although S. plana is believed to be a strict gonochorist, 58% of the populations sampled included intersexed individuals (described as male clams exhibiting ovotestis). Over the entire region, on average, 8.6% of male clams exhibited intersex, although proportions of affected males ranged from 0% to 53% depending on location. The severity of intersex was assessed using a simple classification scale, with the majority of individuals showing low levels of impact. Sex ratios were significantly skewed at some sites. There were no significant relationships between incidence or severity of intersex; or with size or parasitism of individual clams. Intersex in S. plana is a useful tool to assess endocrine disruptive effects in estuaries, although mechanisms of impact and causative agents remain uncertain.

Mechanism of mda-5 Inhibition by paramyxovirus V proteins.

The RNA helicases encoded by melanoma differentiation-associated gene 5 (mda-5) and retinoic acid-inducible gene I (RIG-I) detect foreign cytoplasmic RNA molecules generated during the course of a virus infection, and their activation leads to induction of type I interferon synthesis. Paramyxoviruses limit the amount of interferon produced by infected cells through the action of their V protein, which binds to and inhibits mda-5. Here we show that activation of both mda-5 and RIG-I by double-stranded RNA (dsRNA) leads to the formation of homo-oligomers through self-association of the helicase domains. We identify a region within the helicase domain of mda-5 that is targeted by all paramyxovirus V proteins and demonstrate that they inhibit activation of mda-5 by blocking dsRNA binding and consequent self-association. In addition to this commonly targeted domain, some paramyxovirus V proteins target additional regions of mda-5. In contrast, V proteins cannot bind to RIG-I and consequently have no effect on the ability of RIG-I to bind dsRNA or to form oligomers.

Mapuera virus, a rubulavirus that inhibits interferon signalling in a wide variety of mammalian cells without degrading STATs.

Mapuera virus (MPRV) is a paramyxovirus that was originally isolated from bats, but its host range remains unknown. It was classified as a member of the genus Rubulavirus on the basis of structural and genetic features. Like other rubulaviruses it encodes a V protein (MPRV/V) that functions as an interferon (IFN) antagonist. Here we show that MPRV/V differs from the IFN antagonists of other rubulaviruses in that it does not induce the proteasomal degradation of STAT proteins, key factors in the IFN signalling cascade. Rather, MPRV/V prevents the nuclear translocation of STATs in response to IFN stimulation and inhibits the formation of the transcription factor complex ISGF3. We also show that MPRV/V blocks IFN signalling in cells from diverse mammalian species and discuss the IFN response as a barrier to cross-species infections.

Molecular characterization of Escherichia coli O157:H7 hide contamination routes: feedlot to harvest.

This study was conducted to identify the origin of Escherichia coli O157:H7 contamination on steer hides at the time of harvest. Samples were collected from the feedlot, transport trailers, and packing plant holding pens and from the colons and hides of feedlot steers. A total of 50 hide samples were positive for E. coli O157:H7 in two geographical locations: the Midwest (25 positive hides) and Southwest (25 positive hides). Hide samples were screened, and the presence of E. coli O157: H7 was confirmed. E. coli O157:H7 isolates were fingerprinted by pulsed-field gel electrophoresis and subjected to multiplex PCR procedures for amplification of E. coli O157:H7 genes stx1, stx2, eaeA, fliC, rfbEO157, and hlyA. Feedlot water trough, pen floor, feed bunk, loading chute, truck trailer side wall and floor, packing plant holding pen floor and side rail, and packing plant cattle drinking water samples were positive for E. coli O157:H7. Pulsed-field gel electrophoresis banding patterns were analyzed after classifying isolates according to the marker genes present and according to packing plant. In this study, hide samples positive for E. coli O157:H7 were traced to other E. coli O157:H7-positive hide, colon, feedlot pen floor fecal, packing plant holding pen drinking water, and transport trailer side wall samples. Links were found between packing plant side rails, feedlot loading chutes, and feedlot pens and between truck trailer, different feedlots, and colons of multiple cattle. This study is the first in which genotypic matches have been made between E. coli O157:H7 isolates obtained from transport trailer side walls and those from cattle hide samples within the packing plant.

FINDMOL: automated identification of macromolecules in electron-density maps.

Automating the determination of novel macromolecular structures via X-ray crystallographic methods involves building a model into an electron-density map. Unfortunately, the conventional crystallographic asymmetric unit volumes are usually not well matched to the biological molecular units. In most cases, the facets of the asymmetric unit cut the molecules into a number of disconnected fragments, rendering interpretation by the crystallographer significantly more difficult. The FINDMOL algorithm is designed to quickly parse the arrangement of trace points (pseudo-atoms) derived from a skeletonized electron-density map without requiring higher level prior information such as sequence information or number of molecules in the asymmetric unit. The algorithm was tested with a variety of density-modified maps computed with medium- to low-resolution data. Typically, the resulting volume resembles the biological unit. In the remaining cases the number of disconnected fragments is very small. In all examples, secondary-structural elements such as alpha-helices or beta-sheets are easily identifiable in the defragmented arrangement. FINDMOL can greatly assist a crystallographer during manual model building or in cases where automatic model building can only build partial models owing to limitations of the data such as low resolution and/or poor phases.

Simian virus 5 V protein acts as an adaptor, linking DDB1 to STAT2, to facilitate the ubiquitination of STAT1.

The V protein of simian virus 5 (SV5) facilitates the ubiquitination and subsequent proteasome-mediated degradation of STAT1. Here we show, by visualizing direct protein-protein interactions and by using the yeast two-hybrid system, that while the SV5 V protein fails to bind to STAT1 directly, it binds directly and independently to both DDB1 and STAT2, two cellular proteins known to be essential for SV5-mediated degradation of STAT1. We also demonstrate that STAT1 and STAT2 interact independently of SV5 V and show that SV5 V protein acts as an adaptor molecule linking DDB1 to STAT2/STAT1 heterodimers, which in the presence of additional accessory cellular proteins, including Cullin 4a, can ubiquitinate STAT1. Additionally, we show that the avidity of STAT2 for V is relatively weak but is significantly enhanced by the presence of both STAT1 and DDB1, i.e., the complex of STAT1, STAT2, DDB1, and SV5 V is more stable than a complex of STAT2 and V. From these studies we propose a dynamic model in which SV5 V acts as a bridge, bringing together a DDB1/Cullin 4a-containing ubiquitin ligase complex and STAT1/STAT2 heterodimers, which leads to the degradation of STAT1. The loss of STAT1 results in a decrease in affinity of binding of STAT2 for V such that STAT2 either dissociates from V or is displaced from V by STAT1/STAT2 complexes, thereby ensuring the cycling of the DDB1 and SV5 V containing E3 complex for continued rounds of STAT1 ubiquitination and degradation.

The V proteins of paramyxoviruses bind the IFN-inducible RNA helicase, mda-5, and inhibit its activation of the IFN-beta promoter.

Most paramyxoviruses circumvent the IFN response by blocking IFN signaling and limiting the production of IFN by virus-infected cells. Here we report that the highly conserved cysteine-rich C-terminal domain of the V proteins of a wide variety of paramyxoviruses binds melanoma differentiation-associated gene 5 (mda-5) product. mda-5 is an IFN-inducible host cell DExD/H box helicase that contains a caspase recruitment domain at its N terminus. Overexpression of mda-5 stimulated the basal activity of the IFN-beta promoter in reporter gene assays and significantly enhanced the activation of the IFN-beta promoter by intracellular dsRNA. Both these activities were repressed by coexpression of the V proteins of simian virus 5, human parainfluenza virus 2, mumps virus, Sendai virus, and Hendra virus. Similar results to the reporter assays were obtained by measuring IFN production. Inhibition of mda-5 by RNA interference or by dominant interfering forms of mda-5 significantly inhibited the activation of the IFN-beta promoter by dsRNA. It thus appears that mda-5 plays a central role in an intracellular signal transduction pathway that can lead to the activation of the IFN-beta promoter, and that the V proteins of paramyxoviruses interact with mda-5 to block its activity.

A molecular cytogenetic map of sorghum chromosome 1. Fluorescence in situ hybridization analysis with mapped bacterial artificial chromosomes.

We used structural genomic resources for Sorghum bicolor (L.) Moench to target and develop multiple molecular cytogenetic probes that would provide extensive coverage for a specific chromosome of sorghum. Bacterial artificial chromosome (BAC) clones containing molecular markers mapped across sorghum linkage group A were labeled as probes for fluorescence in situ hybridization (FISH). Signals from single-, dual-, and multiprobe BAC-FISH to spreads of mitotic chromosomes and pachytene bivalents were associated with the largest sorghum chromosome, which bears the nucleolus organizing region (NOR). The order of individual BAC-FISH loci along the chromosome was fully concordant to that of marker loci along the linkage map. In addition, the order of several tightly linked molecular markers was clarified by FISH analysis. The FISH results indicate that markers from the linkage map positions 0.0-81.8 cM reside in the short arm of chromosome 1 whereas markers from 81.8-242.9 cM are located in the long arm of chromosome 1. The centromere and NOR were located in a large heterochromatic region that spans approximately 60% of chromosome 1. In contrast, this region represents only 0.7% of the total genetic map distance of this chromosome. Variation in recombination frequency among euchromatic chromosomal regions also was apparent. The integrated data underscore the value of cytological data, because minor errors and uncertainties in linkage maps can involve huge physical regions. The successful development of multiprobe FISH cocktails suggests that it is feasible to develop chromosome-specific "paints" from genomic resources rather than flow sorting or microdissection and that when applied to pachytene chromatin, such cocktails provide an especially powerful framework for mapping. Such a molecular cytogenetic infrastructure would be inherently cross-linked with other genomic tools and thereby establish a cytogenomics system with extensive utility in development and application of genomic resources, cloning, transgene localization, development of plant "chromonomics," germplasm introgression, and marker-assisted breeding. In combination with previously reported work, the results indicate that a sorghum cytogenomics system would be partially applicable to other gramineous genera.

Development and mapping of AFLP markers linked to the sorghum fertility restorer gene rf4.

The restoration of male fertility in the sorghum IS1112 C (A3) male-sterile cytoplasm is through a two-gene gametophytic system involving complementary action of the restoring alleles Rf3 and Rf4. To develop markers suitable for mapping rf4, AFLP technology was applied to bulks of sterile and fertile individuals from a segregating BC(3)F(1) population. Three AFLP markers linked to rf4were identified and subsequently converted to STS/CAPS markers, two of which are co-dominant. Based on a population of 378 BC(1)F(1) individuals, two STS/CAPS markers, LW7 and LW8, mapped to within 5.31 and 3.18 cM, respectively, of rf4, while an STS marker, LW9, was positioned 0.79 cM on the flanking side of rf4. Markers LW8 and LW9 were used to screen sorghum BAC libraries to identify the genomic region encoding rf4. A series of BAC clones shown to represent a genomic region of linkage group E were identified by the rf4-linked markers. A contig of BAC clones flanking the LW9 marker represent seed clones on linkage group E, from which fine mapping of the rf4 locus and chromosome walking can be initiated.

Verification of different glial fibrillary acidic protein (GFAP) analyses as accurate detectors of central nervous system tissue in advanced meat recovery (AMR) products.

A glial fibrillary acidic protein (GFAP) fluorescent enzyme linked immunosorbant assay (ELISA) was compared with an ELISA test kit for GFAP to determine the level of central nervous system (CNS) tissue in advanced meat recovery (AMR) products. The test kit results were highly correlated (r=0.975) with the fluorescent ELISA. Meat cuts and AMR were analyzed on site at 14 meat plants utilizing the test kits. In seven of the plants all AMR samples had less than 1 ng GFAP. Seven of the plants had greater than 1 ng GFAP in AMR samples. Development of proper process controls to eliminate inclusion of spinal cord in AMR materials should bring all values to less than 1 ng GFAP, a level slightly above background.

Mapping genes on an integrated sorghum genetic and physical map using cDNA selection technology.

Sorghum is an important target of plant genomics. This cereal has unusual tolerance to adverse environments, a small genome (750 Mbp) relative to most other grasses, a diverse germplasm, and utility for comparative genomics with rice, maize and other grasses. In this study, a modified cDNA selection protocol was developed to aid the discovery and mapping of genes across an integrated genetic and physical map of the sorghum genome. BAC DNA from the sorghum genome map was isolated and covalently bound in arrayed tubes for efficient liquid handling. Amplifiable cDNA sequence tags were isolated by hybridization to individual sorghum BACs, cloned and sequenced. Analysis of a fully sequenced sorghum BAC indicated that about 80% of known or predicted genes were detected in the sequence tags, including multiple tags from different regions of individual genes. Data from cDNA selection using the fully sequenced BAC indicate that the occurrence of mislocated cDNA tags is very low. Analysis of 35 BACs (5.25 Mb) from sorghum linkage group B revealed (and therefore mapped) two sorghum genes and 58 sorghum ESTs. Additionally, 31 cDNA tags that had significant homologies to genes from other species were also isolated. The modified cDNA selection procedure described here will be useful for genome-wide gene discovery and EST mapping in sorghum, and for comparative genomics of sorghum, rice, maize and other grasses.

The detection of central nervous system tissue on beef carcasses and in comminuted beef.

We report the development and validation of a fluorescent enzyme-linked immunosorbent assay (ELISA) for glial fibrillary acidic protein (GFAP), which can be used as a rapid and sensitive method to detect CNS tissue in meat products. The fluorometric assay is sensitive to 0.2 ng GFAP and has an intra-assay coefficient of variation (CV) of 2.0% and an interassay CV of 14.1%. Bovine spinal cord and brain demonstrate dose-response curves that are parallel to GFAP standards, whereas peripheral sciatic nerve and cervical ganglia also cross-react at high tissue levels. The use of another central nervous system marker, syntaxin 1-B, was not effective for neural tissue detection. Less than 1.0 ng GFAP per mg tissue was found on most beef subprimals and advanced meat recovery (AMR) product. Occasional samples contained higher levels of GFAP, probably because of contamination by the carcass-splitting saw, incomplete removal of the spinal cord, or a chance sampling of a major nerve. Further reduction of CNS content was facilitated by removal of the cervical vertebrae and the spinal canal prior to processing beef chuck bones through AMR equipment. The presence of GFAP was very low (0.037 ng/mg) in beef patties collected from major processors throughout the USA. The presence of normal sausage ingredients or heating the product to 80 degrees C for 60 min did not affect the detection of GFAP. Heating the product to 115 degrees C for 100 min eliminated the detectability of GFAP.

A high-throughput AFLP-based method for constructing integrated genetic and physical maps: progress toward a sorghum genome map.

Sorghum is an important target for plant genomic mapping because of its adaptation to harsh environments, diverse germplasm collection, and value for comparing the genomes of grass species such as corn and rice. The construction of an integrated genetic and physical map of the sorghum genome (750 Mbp) is a primary goal of our sorghum genome project. To help accomplish this task, we have developed a new high-throughput PCR-based method for building BAC contigs and locating BAC clones on the sorghum genetic map. This task involved pooling 24,576 sorghum BAC clones ( approximately 4x genome equivalents) in six different matrices to create 184 pools of BAC DNA. DNA fragments from each pool were amplified using amplified fragment length polymorphism (AFLP) technology, resolved on a LI-COR dual-dye DNA sequencing system, and analyzed using Bionumerics software. On average, each set of AFLP primers amplified 28 single-copy DNA markers that were useful for identifying overlapping BAC clones. Data from 32 different AFLP primer combinations identified approximately 2400 BACs and ordered approximately 700 BAC contigs. Analysis of a sorghum RIL mapping population using the same primer pairs located approximately 200 of the BAC contigs on the sorghum genetic map. Restriction endonuclease fingerprinting of the entire collection of sorghum BAC clones was applied to test and extend the contigs constructed using this PCR-based methodology. Analysis of the fingerprint data allowed for the identification of 3366 contigs each containing an average of 5 BACs. BACs in approximately 65% of the contigs aligned by AFLP analysis had sufficient overlap to be confirmed by DNA fingerprint analysis. In addition, 30% of the overlapping BACs aligned by AFLP analysis provided information for merging contigs and singletons that could not be joined using fingerprint data alone. Thus, the combination of fingerprinting and AFLP-based contig assembly and mapping provides a reliable, high-throughput method for building an integrated genetic and physical map of the sorghum genome.