Electronic PCR: bridging the gap between genome mapping and genome sequencing

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat

The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

Development of genic SSR marker resources from RNA-seq data in Camellia japonica and their application in the genus Camellia

Camellia is a genus of flowering plants in the family Theaceae, and several species in this genus have economic importance. Although a great deal of molecular makers has been developed for molecular assisted breeding in genus Camellia in the past decade, the number of simple sequence repeats (SSRs) publicly available for plants in this genus is insufficient. In this study, a total of 28,854 potential SSRs were identified with a frequency of 4.63 kb. A total of 172 primer pairs were synthesized and preliminarily screened in 10 C. japonica accessions, and of these primer pairs, 111 were found to be polymorphic. Fifty-one polymorphic SSR markers were randomly selected to perform further analysis of the genetic relationships of 89 accessions across the genus Camellia. Cluster analysis revealed major clusters corresponding to those based on taxonomic classification and geographic origin. Furthermore, all the genotypes of C. japonica separated and consistently grouped well in the genetic structure analysis. The results of the present study provide high-quality SSR resources for molecular genetic breeding studies in camellia plants.

In silico and in vitro evaluation of primers for molecular differentiation of Leishmania species

Leishmaniasis is a zoonotic disease caused by over 20 species of protozoan parasites of the genus Leishmania. Infection is commonly spread by sandflies and produces a wide spectrum of clinical signs and symptoms. Therefore, from an epidemiological and therapeutic standpoint, it is important to detect and differentiate Leishmania spp. The objective of this study was to combinate in silico and in vitro strategies to evaluate the analytical specificity of primers previously described in the literature. According to electronic PCR (e-PCR) analysis, 23 out of 141 pairs of primers selected through literature search matched their previously reported analytical specificity. In vitro evaluation of nine of these primer pairs by quantitative PCR (qPCR) confirmed the analytical specificity of five of them at the level of Leishmania spp., L. mexicana complex or Leishmania and Viannia subgenera. Based on these findings, the combination of e-PCR and qPCR is suggested to be a valuable approach to maximize the specificity of new primer pairs for the laboratory diagnosis of infections with Leishmania spp.

Identification of reference genes for real-time polymerase chain reaction gene expression studies in Nile rats fed Water-Soluble Palm Fruit Extract

The Nile rat (Arvicanthis niloticus) is a novel diurnal carbohydrate-sensitive rodent useful for studies on type 2 diabetes mellitus (T2DM) and the metabolic syndrome. Hepatic responses to T2DM and any interventions thereof can be evaluated via transcriptomic gene expression analysis. However, the study of gene expression via real-time reverse transcription quantitative polymerase chain reaction (RT-qPCR) requires identification of stably expressed reference genes for accurate normalisation. This study describes the evaluation and identification of stable reference genes in the livers from Control Nile rats as well as those supplemented with Water-Soluble Palm Fruit Extract, which has been previously shown to attenuate T2DM in this animal model. Seven genes identified as having stable expression in RNA-Sequencing transcriptome analysis were chosen for verification using real-time RT-qPCR. Six commonly used reference genes from previous literature and two genes from a previous microarray gene expression study in Nile rats were also evaluated. The expression data of these 15 candidate reference genes were analysed using the RefFinder software which incorporated analyses performed by various algorithms. The Hpd, Pnpla6 and Vpp2 genes were identified as the most stable across the 36 samples tested. Their applicability was demonstrated through the normalisation of the gene expression profiles of two target genes, Cela1 and Lepr. In conclusion, three novel reference genes which can be used for robust normalisation of real-time RT-qPCR data were identified, thereby facilitating future hepatic gene expression studies in the Nile rat.

MultiplexSSR: A pipeline for developing multiplex SSR-PCR assays from resequencing data

Next-generation sequencing has greatly promoted the investigation of single nucleotide polymorphisms, while studies of simple sequence repeats are sharply decreasing. However, simple sequence repeats still present some advantages in conservation genetics. In this study, an end-to-end pipeline referred to as MultiplexSSR was established to develop multiplex PCR assays in batches with highly polymorphic simple sequence repeats for capillary platforms from resequencing data. The distribution of single sequence repeats in the genome, the error profiles of genotypes and allelotypes, and the increase in the allele length range depending on the number of individuals were investigated. A total of 98% of single sequence repeats presented lengths of less than 100 bp. The error rate of the genotyping and allelotyping of dimeric patterns was ten times higher than those for other patterns. The error rate of allelotyping was less than that of genotyping. The allele length range reached approximate saturation with 10 individuals. This pipeline uses allele numbers to select highly polymorphic loci, masks loci with variation, and applies in silico PCR to improve primer specificity. The application of the developed multiplex SSR-PCR assays validated the pipeline's robustness, showing higher polymorphism and stability for the developed simple sequence repeats and a lower cost for genotyping and providing low-depth resequencing data from less than a dozen individuals for the development of markers. This pipeline fills the gap between next-generation sequencing and multiplex SSR-PCR.

Genome-Wide Mapping of Loci Associated With Resistance to Clubroot in Brassica napus ssp. napobrassica (Rutabaga) Accessions From Nordic Countries

Rutabaga [Brassica napus ssp. napobrassica (L.) Hanelt] is reported to be an excellent source of clubroot (Plasmodiophora brassicae) resistance genes. In this study, 124 rutabaga accessions from the Nordic countries (Norway, Sweden, Finland, Denmark, and Iceland) were evaluated for their reaction to five single-spore isolates representing P. brassicae pathotypes 2F, 3H, 5I, 6M, and 8N and 12 field isolates representing pathotypes 2B, 3A, 3O, 5C, 5G, 5K, 5L, 5X (two isolates, L-G2 and L-G3), 8E, 8J, and 8P. The accessions were also genotyped using a 15K Brassica SNP array and 60 PCR-based primers linked to previously identified clubroot resistance genes. Six thousand eight hundred sixty-one SNP markers were retained after filtering with TASSEL 5.0, and used to evaluate four general linear models (GLM) and four mixed linear models (MLM). The PCA + K and Q + K MLM models gave the minimal deviance of the observed from the expected distribution in quantile-quantile plots, and hence were used for SNP-clubroot association analyses. In addition, 108 alleles derived from the PCR-based markers and the phenotypic data were analyzed with the PCA + K model. Forty-five SNPs and four PCR-based markers were identified to be associated strongly with resistance to isolates representing 13 pathotypes (2F, 3H, 5I, 6M, 8N, 2B, 3A, 3O, 5C, 5G, 5K, 5L, and 8P). These markers revealed the top and bottom segments of rutabaga chromosome A03 and the middle segment of chromosome A08 as genomic hotspots associated with resistance to the different P. brassicae pathotypes.

Genome-wide identification of microsatellite markers from cultivated peanut (Arachis hypogaea L.)

Background

Microsatellites, or simple sequence repeats (SSRs), represent important DNA variations that are widely distributed across the entire plant genome and can be used to develop SSR markers, which can then be used to conduct genetic analyses and molecular breeding. Cultivated peanut (A. hypogaea L.), an important oil crop worldwide, is an allotetraploid (AABB, 2n = 4× = 40) plant species. Because of its complex genome, genomic marker development has been very challenging. However, sequencing of cultivated peanut genome allowed us to develop genomic markers and construct a high-density physical map.

Results

A total of 8,329,496 SSRs were identified, including 3,772,653, 4,414,961, and 141,882 SSRs that were distributed in subgenome A, B, and nine scaffolds, respectively. Based on the flanking sequences of the identified SSRs, a total of 973,984 newly developed SSR markers were developed in subgenome A (462,267), B (489,394), and nine scaffolds (22,323), with an average density of 392.45 markers per Mb. In silico PCR evaluation showed that an average of 88.32% of the SSR markers generated only one in silico-specific product in two tetraploid A. hypogaea varieties, Tifrunner and Shitouqi. A total of 39,599 common SSR markers were identified among the two A. hypogaea varieties and two progenitors, A. duranensis and A. ipaensis. Additionally, an amplification effectiveness of 44.15% was observed by real PCR validation. Moreover, a total of 1276 public SSR loci were integrated with the newly developed SSR markers. Finally, a previously known leaf spot quantitative trait locus (QTL), qLLS_T13_A05_7, was determined to be in a 1.448-Mb region on chromosome A05. In this region, a total of 819 newly developed SSR markers were located and 108 candidate genes were detected.

Conclusions

The availability of these newly developed and public SSR markers both provide a large number of molecular markers that could potentially be used to enhance the process of trait genetic analyses and improve molecular breeding strategies for cultivated peanut.

Genome-wide analysis of SSR and ILP markers in trees: diversity profiling, alternate distribution, and applications in duplication

Molecular markers are efficient tools for breeding and genetic studies. However, despite their ecological and economic importance, their development and application have long been hampered. In this study, we identified 524,170 simple sequence repeat (SSR), 267,636 intron length polymorphism (ILP), and 11,872 potential intron polymorphism (PIP) markers from 16 tree species based on recently available genome sequences. Larger motifs, including hexamers and heptamers, accounted for most of the seven different types of SSR loci. Within these loci, A/T bases comprised a significantly larger proportion of sequence than G/C. SSR and ILP markers exhibited an alternative distribution pattern. Most SSRs were monomorphic markers, and the proportions of polymorphic markers were positively correlated with genome size. By verifying with all 16 tree species, 54 SSR, 418 ILP, and four PIP universal markers were obtained, and their efficiency was examined by PCR. A combination of five SSR and six ILP markers were used for the phylogenetic analysis of 30 willow samples, revealing a positive correlation between genetic diversity and geographic distance. We also found that SSRs can be used as tools for duplication analysis. Our findings provide important foundations for the development of breeding and genetic studies in tree species.

Genomic survey sequencing for development and validation of single-locus SSR markers in peanut (Arachis hypogaea L.)

BACKGROUND

Single-locus markers have many advantages compared with multi-locus markers in genetic and breeding studies because their alleles can be assigned to particular genomic loci in diversity analyses. However, there is little research on single-locus SSR markers in peanut. Through the de novo assembly of DNA sequencing reads of A. hypogaea, we developed single-locus SSR markers in a genomic survey for better application in genetic and breeding studies of peanut.

RESULTS

In this study, DNA libraries with four different insert sizes were used for sequencing with 150 bp paired-end reads. Approximately 237 gigabases of clean data containing 1,675,631,984 reads were obtained after filtering. These reads were assembled into 2,102,446 contigs with an N50 length of 1,782 bp, and the contigs were further assembled into 1,176,527 scaffolds with an N50 of 3,920 bp. The total length of the assembled scaffold sequences was 2.0 Gbp, and 134,652 single-locus SSRs were identified from 375,180 SSRs. Among these developed single-locus SSRs, trinucleotide motifs were the most abundant, followed by tetra-, di-, mono-, penta- and hexanucleotide motifs. The most common motif repeats for the various types of single-locus SSRs have a tendency to be A/T rich. A total of 1,790 developed in silico single-locus SSR markers were chosen and used in PCR experiments to confirm amplification patterns. Of them, 1,637 markers that produced single amplicons in twelve inbred lines were considered putative single-locus markers, and 290 (17.7 %) showed polymorphisms. A further F2 population study showed that the segregation ratios of the 97 developed SSR markers, which showed polymorphisms between the parents, were consistent with the Mendelian inheritance law for single loci (1:2:1). Finally, 89 markers were assigned to an A. hypogaea linkage map. A subset of 100 single-locus SSR markers was shown to be highly stable and universal in a collection of 96 peanut accessions. A neighbor-joining tree of this natural population showed that genotypes have obviously correlation with botanical varieties.

CONCLUSIONS

We have shown that the detection of single-locus SSR markers from a de novo genomic assembly of a combination of different-insert-size libraries is highly efficient. This is the first report of the development of genome-wide single-locus markers for A. hypogaea, and the markers developed in this study will be useful for gene tagging, sequence scaffold assignment, linkage map construction, diversity analysis, variety identification and association mapping in peanut.

Genome sequencing and comparative genomics of the broad host-range pathogen Rhizoctonia solani AG8

Rhizoctonia solani is a soil-borne basidiomycete fungus with a necrotrophic lifestyle which is classified into fourteen reproductively incompatible anastomosis groups (AGs). One of these, AG8, is a devastating pathogen causing bare patch of cereals, brassicas and legumes. R. solani is a multinucleate heterokaryon containing significant heterozygosity within a single cell. This complexity posed significant challenges for the assembly of its genome. We present a high quality genome assembly of R. solani AG8 and a manually curated set of 13,964 genes supported by RNA-seq. The AG8 genome assembly used novel methods to produce a haploid representation of its heterokaryotic state. The whole-genomes of AG8, the rice pathogen AG1-IA and the potato pathogen AG3 were observed to be syntenic and co-linear. Genes and functions putatively relevant to pathogenicity were highlighted by comparing AG8 to known pathogenicity genes, orthology databases spanning 197 phytopathogenic taxa and AG1-IA. We also observed SNP-level "hypermutation" of CpG dinucleotides to TpG between AG8 nuclei, with similarities to repeat-induced point mutation (RIP). Interestingly, gene-coding regions were widely affected along with repetitive DNA, which has not been previously observed for RIP in mononuclear fungi of the Pezizomycotina. The rate of heterozygous SNP mutations within this single isolate of AG8 was observed to be higher than SNP mutation rates observed across populations of most fungal species compared. Comparative analyses were combined to predict biological processes relevant to AG8 and 308 proteins with effector-like characteristics, forming a valuable resource for further study of this pathosystem. Predicted effector-like proteins had elevated levels of non-synonymous point mutations relative to synonymous mutations (dN/dS), suggesting that they may be under diversifying selection pressures. In addition, the distant relationship to sequenced necrotrophs of the Ascomycota suggests the R. solani genome sequence may prove to be a useful resource in future comparative analysis of plant pathogens.

Supermodels: sorghum and maize provide mutual insight into the genetics of flowering time

Nested association mapping (NAM) offers power to dissect complex, quantitative traits. This study made use of a recently developed sorghum backcross (BC)-NAM population to dissect the genetic architecture of flowering time in sorghum; to compare the QTL identified with other genomic regions identified in previous sorghum and maize flowering time studies and to highlight the implications of our findings for plant breeding. A subset of the sorghum BC-NAM population consisting of over 1,300 individuals from 24 families was evaluated for flowering time across multiple environments. Two QTL analysis methodologies were used to identify 40 QTLs with predominately small, additive effects on flowering time; 24 of these co-located with previously identified QTL for flowering time in sorghum and 16 were novel in sorghum. Significant synteny was also detected with the QTL for flowering time detected in a comparable NAM resource recently developed for maize (Zea mays) by Buckler et al. (Science 325:714-718, 2009). The use of the sorghum BC-NAM population allowed us to catalogue allelic variants at a maximal number of QTL and understand their contribution to the flowering time phenotype and distribution across diverse germplasm. The successful demonstration of the power of the sorghum BC-NAM population is exemplified not only by correspondence of QTL previously identified in sorghum, but also by correspondence of QTL in different taxa, specifically maize in this case. The unification across taxa of the candidate genes influencing complex traits, such as flowering time can further facilitate the detailed dissection of the genetic control and causal genes.

Expedited batch processing and analysis of transposon insertions

Background

With advances in sequencing technology, greater and greater amounts of eukaryotic genome data are becoming available. Often, large portions of these genomes consist of transposable elements, frequently accounting for 50% or more in vertebrates. Each transposable element family may have thousands or tens of thousands of individual copies within a given genome, and therefore it can take an exorbitant amount of time and effort to process data in a meaningful fashion.

Findings

In order to combat this problem, we developed a set of bioinformatics techniques and programs to streamline the analysis. This includes a unique Perl script which automates the process of taking BLAST, Repeatmasker and similar data to extract and manipulate the hit sequences from the genome. This script, called Process_hits uses an object-oriented methodology to compile all hit locations from a given file for processing, organize this data into useable categories, and output it in multiple formats.

Conclusions

The program proved capable of handling large amounts of transposon data in an efficient fashion. It is equipped with a number of useful sub-functions, each of which is contained within its own sub-module to allow for greater expandability and as a foundation for future program design.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

Using the NCBI Map Viewer to browse genomic sequence data

This unit includes a basic protocol with an introduction to the Map Viewer, describing how to perform a simple text-based search of genome annotations to view the genomic context of a gene, navigate along a chromosome, zoom in and out, and change the displayed maps to hide and show information. It also describes some of NCBI's sequence-analysis tools, which are provided as links from the Map Viewer. The alternate protocols describe different ways to query the genome sequence, and also illustrate additional features of the Map Viewer. Alternate Protocol 1 shows how to perform and interpret the results of a BLAST search against the human genome. Alternate Protocol 2 demonstrates how to retrieve a list of all genes between two STS markers. Finally, Alternate Protocol 3 shows how to find all annotated members of a gene family.

Location of major effect genes in sorghum (Sorghum bicolor (L.) Moench)

Major effect genes are often used for germplasm identification, for diversity analyses and as selection targets in breeding. To date, only a few morphological characters have been mapped as major effect genes across a range of genetic linkage maps based on different types of molecular markers in sorghum (Sorghum bicolor (L.) Moench). This study aims to integrate all available previously mapped major effect genes onto a complete genome map, linked to the whole genome sequence, allowing sorghum breeders and researchers to link this information to QTL studies and to be aware of the consequences of selection for major genes. This provides new opportunities for breeders to take advantage of readily scorable morphological traits and to develop more effective breeding strategies. We also provide examples of the impact of selection for major effect genes on quantitative traits in sorghum. The concepts described in this paper have particular application to breeding programmes in developing countries where molecular markers are expensive or impossible to access.

Using the NCBI map viewer to browse genomic sequence data

This unit includes a Basic Protocol with an introduction to the Map Viewer, describing how to perform a simple text-based search of genome annotations to view the genomic context of a gene, navigate along a chromosome, zoom in and out, and change the displayed maps to hide and show information. It also describes some of NCBI's sequence-analysis tools, which are provided as links from the Map Viewer. The Alternate Protocols describe different ways to query the genome sequence, and also illustrate additional features of the Map Viewer. Alternate Protocol 1 shows how to perform and interpret the results of a BLAST search against the human genome. Alternate Protocol 2 demonstrates how to retrieve a list of all genes between two STS markers. Finally, Alternate Protocol 3 shows how to find all annotated members of a gene family.

QTL global meta-analysis: are trait determining genes clustered?

Background

A key open question in biology is if genes are physically clustered with respect to their known functions or phenotypic effects. This is of particular interest for Quantitative Trait Loci (QTL) where a QTL region could contain a number of genes that contribute to the trait being measured.

Results

We observed a significant increase in gene density within QTL regions compared to non-QTL regions and/or the entire bovine genome. By grouping QTL from the Bovine QTL Viewer database into 8 categories of non-redundant regions, we have been able to analyze gene density and gene function distribution, based on Gene Ontology (GO) with relation to their location within QTL regions, outside of QTL regions and across the entire bovine genome. We identified a number of GO terms that were significantly over represented within particular QTL categories. Furthermore, select GO terms expected to be associated with the QTL category based on common biological knowledge have also proved to be significantly over represented in QTL regions.

Conclusion

Our analysis provides evidence of over represented GO terms in QTL regions. This increased GO term density indicates possible clustering of gene functions within QTL regions of the bovine genome. Genes with similar functions may be grouped in specific locales and could be contributing to QTL traits. Moreover, we have identified over-represented GO terminology that from a biological standpoint, makes sense with respect to QTL category type.

Understanding the development of human bladder cancer by using a whole-organ genomic mapping strategy

The search for the genomic sequences involved in human cancers can be greatly facilitated by maps of genomic imbalances identifying the involved chromosomal regions, particularly those that participate in the development of occult preneoplastic conditions that progress to clinically aggressive invasive cancer. The integration of such regions with human genome sequence variation may provide valuable clues about their overall structure and gene content. By extension, such knowledge may help us understand the underlying genetic components involved in the initiation and progression of these cancers. We describe the development of a genome-wide map of human bladder cancer that tracks its progression from in situ precursor conditions to invasive disease. Testing for allelic losses using a genome-wide panel of 787 microsatellite markers was performed on multiple DNA samples, extracted from the entire mucosal surface of the bladder and corresponding to normal urothelium, in situ preneoplastic lesions, and invasive carcinoma. Using this approach, we matched the clonal allelic losses in distinct chromosomal regions to specific phases of bladder neoplasia and produced a detailed genetic map of bladder cancer development. These analyses revealed three major waves of genetic changes associated with growth advantages of successive clones and reflecting a stepwise conversion of normal urothelial cells into cancer cells. The genetic changes map to six regions at 3q22-q24, 5q22-q31, 9q21-q22, 10q26, 13q14, and 17p13, which may represent critical hits driving the development of bladder cancer. Finally, we performed high-resolution mapping using single nucleotide polymorphism markers within one region on chromosome 13q14, containing the model tumor suppressor gene RB1, and defined a minimal deleted region associated with clonal expansion of in situ neoplasia. These analyses provided new insights on the involvement of several non-coding sequences mapping to the region and identified novel target genes, termed forerunner (FR) genes, involved in early phases of cancer development.

Using the NCBI Map Viewer to browse genomic sequence data

This unit includes an introduction to the Map Viewer, which describes how to perform a simple text-based search of genome annotations to view the genomic context of a gene, navigate along a chromosome, zoom in and out, and change the displayed maps to hide and show information. It also describes some of NCBI's sequence-analysis tools, which are provided as links from the Map Viewer. The Alternate Protocols describe different ways to query the genome sequence, and also illustrate additional features of the Map Viewer. Alternate Protocol 1 shows how to perform and interpret the results of a BLAST search against the human genome. Alternate Protocol 2 demonstrates how to retrieve a list of all genes between two STS markers. Finally, Alternate Protocol 3 shows how to find all annotated members of a gene family.

Using the Ensembl genome server to browse genomic sequence data

The Ensembl genome Web browser (http://www.ensembl.org) provides a comprehensive source of automatic annotation of the human genome sequence (as well as other species of biomedical interest), with confirmed gene predictions that have been integrated with external data sources. This unit describes how to use the Ensembl browser, how to find your gene or protein of interest and get information and external links about them, and how to use the comparative genomic data.

14th International HLA and Immunogenetics Workshop: report on mapping microsatellite markers in the major histocompatibility complex region

This paper describes the use of the program e-pcr to localize 687 known major histocompatability complex (MHC) microsatellite primer pairs to their sequence positions in several genomic assemblies across the MHC region. The sequences used were the Sequences of Sanger Institute's MHC Haplotype Project: COX, PGF, QBL, as well as the Celera, and Reference (PGF across extended MHC) sequences from the NCBI genomic build 36. More than 95% (664/687) of the markers mapped unambiguously to the Reference assembly sequence. All primer pairs used in this analysis, and those were previously unknown to UniSTS, have now been assigned permanent public UniSTS identifiers. Mapping and descriptive data for each primer pair are available at the publicly accessible dbMHC microsatellite resource: http://www.ncbi.nlm.nih.gov/projects/mhc/xslcgi.fcgi?cmd=mssearch.

Identification and fine mapping of AvrPi15, a novel avirulence gene of Magnaporthe grisea

Avirulence of Magnaporthe grisea isolate CHL346 on rice cultivar GA25 was studied with 242 ascospore progenies derived from the cross CHL346 x CHL42. Segregation analysis of the avirulence in the progeny population was in agreement with the existence of a single avirulence (Avr) gene, designated as AvrPi15. For mapping the Avr gene, we developed a total of 121 microsatellite DNA markers [simple sequence repeat (SSR)], which evenly distributed in the whole-genome of M. grisea through bioinformatics analysis (BIA) using the publicly available sequence. Linkage analysis of the AvrPi15 gene with these SSR markers showed that six markers on chromosome 6, MS6-1, MS6-2, MS6-3, MS6-7, MS6-8 and MS6-10, were linked to the AvrPi15 locus. To further define the chromosomal location of the AvrPi15 locus, two additional markers, MS6-17 and STS6-6, which were developed based on the sequences of telomeric region 11 (TEL11), were subjected to linkage analysis. The results showed that MS6-17 and STS6-6 were associated with the locus by 3.3 and 0.8 cM, respectively. To finely map the Avr gene, two additional candidate avirulence gene (CAG) markers, CAG6-1 and CAG6-2, were developed based on the gene annotation of the sequence of TEL 11. Linkage analysis of the Avr gene with these two markers revealed that both of them completely cosegregated with the AvrPi15 locus. Finally, this locus was physically mapped into approximately 7.2-kb interval of the TEL11 by BIA using these sequence-ready markers. This is the key step toward positional cloning of the AvrPi15 gene.

Genetic analysis and fine mapping of a rice brown planthopper (Nilaparvata lugens Stål) resistance gene bph19(t)

Genetic analysis and fine mapping of a resistance gene against brown planthopper (BPH) biotype 2 in rice was performed using two F(2) populations derived from two crosses between a resistant indica cultivar (cv.), AS20-1, and two susceptible japonica cvs., Aichi Asahi and Lijiangxintuanheigu. Insect resistance was evaluated using F(1) plants and the two F(2) populations. The results showed that a single recessive gene, tentatively designated as bph19(t), conditioned the resistance in AS20-1. A linkage analysis, mainly employing microsatellite markers, was carried out in the two F(2) populations through bulked segregant analysis and recessive class analysis (RCA), in combination with bioinformatics analysis (BIA). The resistance gene locus bph19(t) was finely mapped to a region of about 1.0 cM on the short arm of chromosome 3, flanked by markers RM6308 and RM3134, where one known marker RM1022, and four new markers, b1, b2, b3 and b4, developed in the present study were co-segregating with the locus. To physically map this locus, the bph19(t)-linked markers were landed on bacterial artificial chromosome or P1 artificial chromosome clones of the reference cv., Nipponbare, released by the International Rice Genome Sequencing Project. Sequence information of these clones was used to construct a physical map of the bph19(t) locus, in silico, by BIA. The bph19(t) locus was physically defined to an interval of about 60 kb. The detailed genetic and physical maps of the bph19(t) locus will facilitate marker-assisted gene pyramiding and cloning.

Systematic recovery and analysis of full-ORF human cDNA clones

The Mammalian Gene Collection (MGC) consortium (http://mgc.nci.nih.gov) seeks to establish publicly available collections of full-ORF cDNAs for several organisms of significance to biomedical research, including human. To date over 15,200 human cDNA clones containing full-length open reading frames (ORFs) have been identified via systematic expressed sequence tag (EST) analysis of a diverse set of cDNA libraries; however, further systematic EST analysis is no longer an efficient method for identifying new cDNAs. As part of our involvement in the MGC program, we have developed a scalable method for targeted recovery of cDNA clones to facilitate recovery of genes absent from the MGC collection. First, cDNA is synthesized from various RNAs, followed by polymerase chain reaction (PCR) amplification of transcripts in 96-well plates using gene-specific primer pairs flanking the ORFs. Amplicons are cloned into a sequencing vector, and full-length sequences are obtained. Sequences are processed and assembled using Phred and Phrap, and analyzed using Consed and a number of bioinformatics methods we have developed. Sequences are compared with the Reference Sequence (RefSeq) database, and validation of sequence discrepancies is attempted using other sequence databases including dbEST and dbSNP. Clones with identical sequence to RefSeq or containing only validated changes will become part of the MGC human gene collection. Clones containing novel splice variants or polymorphisms have also been identified. Our approach to clone recovery, applied at large scale, has the potential to recover many and possibly most of the genes absent from the MGC collection.

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.

Identification of a purine-rich intronic enhancer element in the mouse eosinophil-associated ribonuclease 2 (mEar 2) gene

The Mus musculus eosinophil-associated ribonuclease (mEar) gene cluster includes multiple distinct coding sequences that are highly divergent orthologs of the human eosinophil ribonucleases, eosinophil-derived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3). We present a transcriptional analysis of the gene encoding mEar 2, the only member of this cluster with a well-defined expression profile. In this work, we demonstrate that the presence of non-coding exon 1 and the intron in tandem with a 361-bp 5' promoter of mEar 2 results in enhanced reporter gene expression, as much as 6-to 10-fold over the activity observed with the 5' promoter alone. We have identified a conserved purine-rich element in the intron of the mEar 2 gene that is necessary for maximum transcription and that interacts specifically with NFAT-binding proteins in nuclear extracts derived from the mouse LA4 epithelial cell line. Similar intronic enhancers have been described as regulating transcription of the human EDN gene, suggesting an overall conservation of an important regulatory strategy.

1-Mb resolution array-based comparative genomic hybridization using a BAC clone set optimized for cancer gene analysis

Array-based comparative genomic hybridization (aCGH) is a recently developed tool for genome-wide determination of DNA copy number alterations. This technology has tremendous potential for disease-gene discovery in cancer and developmental disorders as well as numerous other applications. However, widespread utilization of a CGH has been limited by the lack of well characterized, high-resolution clone sets optimized for consistent performance in aCGH assays and specifically designed analytic software. We have assembled a set of approximately 4100 publicly available human bacterial artificial chromosome (BAC) clones evenly spaced at approximately 1-Mb resolution across the genome, which includes direct coverage of approximately 400 known cancer genes. This aCGH-optimized clone set was compiled from five existing sets, experimentally refined, and supplemented for higher resolution and enhancing mapping capabilities. This clone set is associated with a public online resource containing detailed clone mapping data, protocols for the construction and use of arrays, and a suite of analytical software tools designed specifically for aCGH analysis. These resources should greatly facilitate the use of aCGH in gene discovery.

Replication timing of the human genome

We have developed a directly quantitative method utilizing genomic clone DNA microarrays to assess the replication timing of sequences during the S phase of the cell cycle. The genomic resolution of the replication timing measurements is limited only by the genomic clone size and density. We demonstrate the power of this approach by constructing a genome-wide map of replication timing in human lymphoblastoid cells using an array with clones spaced at 1 Mb intervals and a high-resolution replication timing map of 22q with an array utilizing overlapping sequencing tile path clones. We show a positive correlation, both genome-wide and at a high resolution, between replication timing and a range of genome parameters including GC content, gene density and transcriptional activity.

A functional genomic analysis of cell morphology using RNA interference

Background

The diversity of metazoan cell shapes is influenced by the dynamic cytoskeletal network. With the advent of RNA-interference (RNAi) technology, it is now possible to screen systematically for genes controlling specific cell-biological processes, including those required to generate distinct morphologies.

Results

We adapted existing RNAi technology in Drosophila cell culture for use in high-throughput screens to enable a comprehensive genetic dissection of cell morphogenesis. To identify genes responsible for the characteristic shape of two morphologically distinct cell lines, we performed RNAi screens in each line with a set of double-stranded RNAs (dsRNAs) targeting 994 predicted cell shape regulators. Using automated fluorescence microscopy to visualize actin filaments, microtubules and DNA, we detected morphological phenotypes for 160 genes, one-third of which have not been previously characterized in vivo. Genes with similar phenotypes corresponded to known components of pathways controlling cytoskeletal organization and cell shape, leading us to propose similar functions for previously uncharacterized genes. Furthermore, we were able to uncover genes acting within a specific pathway using a co-RNAi screen to identify dsRNA suppressors of a cell shape change induced by Pten dsRNA.

Conclusions

Using RNAi, we identified genes that influence cytoskeletal organization and morphology in two distinct cell types. Some genes exhibited similar RNAi phenotypes in both cell types, while others appeared to have cell-type-specific functions, in part reflecting the different mechanisms used to generate a round or a flat cell morphology.

Locating sequence on FPC maps and selecting a minimal tiling path

This study discusses three software tools, the first two aid in integrating sequence with an FPC physical map and the third automatically selects a minimal tiling path given genomic draft sequence and BAC end sequences. The first tool, FSD (FPC Simulated Digest), takes a sequenced clone and adds it back to the map based on a fingerprint generated by an in silico digest of the clone. This allows verification of sequenced clone positions and the integration of sequenced clones that were not originally part of the FPC map. The second tool, BSS (Blast Some Sequence), takes a query sequence and positions it on the map based on sequence associated with the clones in the map. BSS has multiple uses as follows: (1) When the query is a file of marker sequences, they can be added as electronic markers. (2) When the query is draft sequence, the results of BSS can be used to close gaps in a sequenced clone or the physical map. (3) When the query is a sequenced clone and the target is BAC end sequences, one may select the next clone for sequencing using both sequence comparison results and map location. (4) When the query is whole-genome draft sequence and the target is BAC end sequences, the results can be used to select many clones for a minimal tiling path at once. The third tool, pickMTP, automates the majority of this last usage of BSS. Results are presented using the rice FPC map, BAC end sequences, and whole-genome shotgun from Syngenta.

Amplification of marine methanotrophic enrichment DNA with 16S rDNA PCR primers for type II alpha proteobacteria methanotrophs

Type II alpha proteobacteria methanotrophs are capable of a wide range of cometabolic transformations of chlorinated solvents and polycyclic aromatic hydrocarbons (PAHs), and this activity has been exploited in many terrestrial bioremediation systems. However, at present, all known obligately marine methanotrophic isolates are Type I gamma proteobacteria which do not have this activity to the extent of Type II methanotrophs. In previous work in our laboratory, determining the presence of Type II alpha proteobacteria methanotrophs in marine enrichment cultures that co-metabolized PAHs required a more sensitive assay. 16S rDNA PCR primers were designed based on oligonucleotide probes for serine pathway methanotrophs and serine pathway methylotrophs with an approximate amplification fragment size of 870 base pairs. Comparison of the primers using double primer BLAST searches in established nucleotide databases showed potential amplification with all Methylocystis and Methylosinus spp., as well as potential amplification with Methylocella palustrus. DNA from Methylosinus trichosporium OB3b, a Type II methanotroph, amplified with the primers with a fragment size of approximately 850 base pairs, whereas DNA extracted from Methylomonas methanica, a Type I methanotroph, did not. The primers were used to amplify DNA extracted from two marine methanotrophic enrichment cultures: a low nitrogen/low copper enrichment to select for Type II methanotrophs and a high nitrogen/high copper enrichment to select for Type I methanotrophs. Although DNA from both cultures amplified with the PCR primers, amplification was stronger in cultures that were specifically enriched for Type II methanotrophs, suggesting the presence of higher numbers of Type II methanotrophs. These results provide further evidence for the existence of Type II marine methanotrophs, suggesting the possibility of exploiting cometabolic activity in marine systems.

DNannotator: Annotation software tool kit for regional genomic sequences

Sequence annotation is essential for genomics-based research. Investigators of a specific genomic region who have developed abundant local discoveries such as genes and genetic markers, or have collected annotations from multiple resources, can be overwhelmed by the difficulty in creating local annotation and the complexity of integrating all the annotations. Presenting such integrated data in a form suitable for data mining and high-throughput experimental design is even more daunting. DNannotator, a web application, was designed to perform batch annotation on a sizeable genomic region. It takes annotation source data, such as SNPs, genes, primers, and so on, prepared by the end-user and/or a specified target of genomic DNA, and performs de novo annotation. DNannotator can also robustly migrate existing annotations in GenBank format from one sequence to another. Annotation results are provided in GenBank format and in tab-delimited text, which can be imported and managed in a database or spreadsheet and combined with existing annotation as desired. Graphic viewers, such as Genome Browser or Artemis, can display the annotation results. Reference data (reports on the process) facilitating the user's evaluation of annotation quality are optionally provided. DNannotator can be accessed at http://sky.bsd.uchicago.edu/DNannotator.htm.

The microcell mediated transfer of human chromosome 8 into highly metastatic rat liver cancer cell line C5F

AIM: Our previous research on the surgical samples of primary liver cancer with CGH showed that the loss of human chromosome 8p had correlation with the metastatic phenotype of liver cancer. In order to seek the functional evidence that there could be a metastatsis suppressor gene (s) for liver cancer on human chromosome 8, we tried to transfer normal human chromosome 8 into rat liver cancer cell line C5F, which had high metastatic potential to lung.

METHODS: Human chromosome 8 randomly marked with neo gene was introduced into C5F cell line by MMCT and positive microcell hybrids were screened by double selections of G418 and HAT. Single cell isolation cloning was applied to clone microcell hybrids. Finally, STS-PCR and WCP-FISH were used to confirm the introduction.

RESULTS: Microcell hybrids resistant to HAT and G418 were obtained and 15 clones were obtained by single-cell isolation cloning. STS-PCR and WCP-FISH proved that human chromosome 8 had been successfully introduced into rat liver cancer cell line C5F. STS-PCR detected a random loss in the chromosome introduced and WCP-FISH found a consistent recombination of the introduced human chromosome with the rat chromosome.

CONCLUSION: The successful introduction of human chromosome 8 into highly metastatic rat liver cancer cell line builds the basis for seeking functional evidence of a metastasis suppressor gene for liver cancer harboring on human chromosome 8 and its subsequent cloning.

A transcript map encompassing a susceptibility locus for bipolar affective disorder on chromosome 4q35

Bipolar affective disorder is one of the most common mental illnesses with a population prevalence of approximately 1%. The disorder is genetically complex, with an increasing number of loci being implicated through genetic linkage studies. However, the specific genetic variations and molecules involved in bipolar susceptibility and pathogenesis are yet to be identified. Genetic linkage analysis has identified a bipolar disorder susceptibility locus on chromosome 4q35, and the interval harbouring this susceptibility gene has been narrowed to a size that is amenable to positional cloning. We have used the resources of the Human Genome Project (HGP) and Celera Genomics to identify overlapping sequenced BAC clones and sequence contigs that represent the region implicated by linkage analysis. A combination of bioinformatic tools and laboratory techniques have been applied to annotate this DNA sequence data and establish a comprehensive transcript map that spans approximately 5.5 Mb. This map encompasses the chromosome 4q35 bipolar susceptibility locus, which localises to a "most probable" candidate interval of approximately 2.3 Mb, within a more conservative candidate interval of approximately 5 Mb. Localised within this map are 11 characterised genes and eight novel genes of unknown function, which together provide a collection of candidate transcripts that may be investigated for association with bipolar disorder. Overall, this region was shown to be very gene-poor, with a high incidence of pseudogenes, and redundant and novel repetitive elements. Our analysis of the interval has demonstrated a significant difference in the extent to which the current HGP and Celera sequence data sets represent this region.

CRIM1 is involved in endothelial cell capillary formation in vitro and is expressed in blood vessels in vivo

In endothelial cells that form capillary-like structures in vitro a variety of genes is upregulated as we have demonstrated previously. In addition to well known genes, we also identified genes never described in endothelial cells before. Here, we report the further characterization of one selected gene called cysteine-rich motor neuron 1 (CRIM1). CRIM1 is strongly upregulated in endothelial cells during tube formation and is expressed by a variety of adherent growing cell lines whereas cell lines grown in suspension do not express CRIM1. By using antisense technology we were able to inhibit CRIM1 expression and demonstrate impaired formation of capillary-like structures in vitro in transfected endothelial cells. Furthermore, we show that CRIM1 is a glycosylated type I transmembrane protein, that accumulates at sites of close cell-to-cell contact upon stimulation. Finally, we found CRIM1 protein to be expressed by endothelial cells of the inner lining of blood vessels in vivo. Taken together our results imply a possible role of CRIM1 in capillary formation and maintainance during angiogenesis.

Accessing the human genome

The majority of the sequence for the human genome is now available. Regardless of the researcher's area of interest, it is quite likely that they will want to use some aspect of this data. This unit helps researchers achieve that goal. It presents the gene models available at NCBI, the UCSC Genome Browser, and Ensembl. It reviews the features and options available from the three web sites to query, display, and download the data. In addition, the unit illustrates how to query each of the databases in order to identify information, such as the genomic location of a novel cDNA, a BAC clone that contains a particular gene, and homologous human genes to a particular protein sequence from a different organism.