Comprehensive characterization and RNA-Seq profiling of the HD-Zip transcription factor family in soybean (Glycine max) during dehydration and salt stress

BACKGROUND

The homeodomain leucine zipper (HD-Zip) transcription factor family is one of the largest plant specific superfamilies, and includes genes with roles in modulation of plant growth and response to environmental stresses. Many HD-Zip genes are characterized in Arabidopsis (Arabidopsis thaliana), and members of the family are being investigated for abiotic stress responses in rice (Oryza sativa), maize (Zea mays), poplar (Populus trichocarpa) and cucumber (Cucmis sativus). Findings in these species suggest HD-Zip genes as high priority candidates for crop improvement.

RESULTS

In this study we have identified members of the HD-Zip gene family in soybean cv. 'Williams 82', and characterized their expression under dehydration and salt stress. Homology searches with BLASTP and Hidden Markov Model guided sequence alignments identified 101 HD-Zip genes in the soybean genome. Phylogeny reconstruction coupled with domain and gene structure analyses using soybean, Arabidopsis, rice, grape (Vitis vinifera), and Medicago truncatula homologues enabled placement of these sequences into four previously described subfamilies. Of the 101 HD-Zip genes identified in soybean, 88 exist as whole-genome duplication-derived gene pairs, indicating high retention of these genes following polyploidy in Glycine ~13 Mya. The HD-Zip genes exhibit ubiquitous expression patterns across 24 conditions that include 17 tissues of soybean. An RNA-Seq experiment performed to study differential gene expression at 0, 1, 6 and 12 hr soybean roots under dehydration and salt stress identified 20 differentially expressed (DE) genes. Several of these DE genes are orthologs of genes previously reported to play a role under abiotic stress, implying conservation of HD-Zip gene functions across species. Screening of HD-Zip promoters identified transcription factor binding sites that are overrepresented in the DE genes under both dehydration and salt stress, providing further support for the role of HD-Zip genes in abiotic stress responses.

CONCLUSIONS

We provide a thorough description of soybean HD-Zip genes, and identify potential candidates with probable roles in dehydration and salt stress. Expression profiles generated for all soybean genes, under dehydration and salt stress, at four time points, will serve as an important resource for the soybean research community, and will aid in understanding plant responses to abiotic stress.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Candidate genes and genetic architecture of symbiotic and agronomic traits revealed by whole-genome, sequence-based association genetics in Medicago truncatula

Genome-wide association study (GWAS) has revolutionized the search for the genetic basis of complex traits. To date, GWAS have generally relied on relatively sparse sampling of nucleotide diversity, which is likely to bias results by preferentially sampling high-frequency SNPs not in complete linkage disequilibrium (LD) with causative SNPs. To avoid these limitations we conducted GWAS with >6 million SNPs identified by sequencing the genomes of 226 accessions of the model legume Medicago truncatula. We used these data to identify candidate genes and the genetic architecture underlying phenotypic variation in plant height, trichome density, flowering time, and nodulation. The characteristics of candidate SNPs differed among traits, with candidates for flowering time and trichome density in distinct clusters of high linkage disequilibrium (LD) and the minor allele frequencies (MAF) of candidates underlying variation in flowering time and height significantly greater than MAF of candidates underlying variation in other traits. Candidate SNPs tagged several characterized genes including nodulation related genes SERK2, MtnodGRP3, MtMMPL1, NFP, CaML3, MtnodGRP3A and flowering time gene MtFD as well as uncharacterized genes that become candidates for further molecular characterization. By comparing sequence-based candidates to candidates identified by in silico 250K SNP arrays, we provide an empirical example of how reliance on even high-density reduced representation genomic makers can bias GWAS results. Depending on the trait, only 30–70% of the top 20 in silico array candidates were within 1 kb of sequence-based candidates. Moreover, the sequence-based candidates tagged by array candidates were heavily biased towards common variants; these comparisons underscore the need for caution when interpreting results from GWAS conducted with sparsely covered genomes.

Comparative genomics of the core and accessory genomes of 48 Sinorhizobium strains comprising five genospecies

BACKGROUND

The sinorhizobia are amongst the most well studied members of nitrogen-fixing root nodule bacteria and contribute substantial amounts of fixed nitrogen to the biosphere. While the alfalfa symbiont Sinorhizobium meliloti RM 1021 was one of the first rhizobial strains to be completely sequenced, little information is available about the genomes of this large and diverse species group.

RESULTS

Here we report the draft assembly and annotation of 48 strains of Sinorhizobium comprising five genospecies. While S. meliloti and S. medicae are taxonomically related, they displayed different nodulation patterns on diverse Medicago host plants, and have differences in gene content, including those involved in conjugation and organic sulfur utilization. Genes involved in Nod factor and polysaccharide biosynthesis, denitrification and type III, IV, and VI secretion systems also vary within and between species. Symbiotic phenotyping and mutational analyses indicated that some type IV secretion genes are symbiosis-related and involved in nitrogen fixation efficiency. Moreover, there is a correlation between the presence of type IV secretion systems, heme biosynthesis and microaerobic denitrification genes, and symbiotic efficiency.

CONCLUSIONS

Our results suggest that each Sinorhizobium strain uses a slightly different strategy to obtain maximum compatibility with a host plant. This large genome data set provides useful information to better understand the functional features of five Sinorhizobium species, especially compatibility in legume-Sinorhizobium interactions. The diversity of genes present in the accessory genomes of members of this genus indicates that each bacterium has adopted slightly different strategies to interact with diverse plant genera and soil environments.

Phylogenetic signal variation in the genomes of Medicago (Fabaceae)

Genome-scale data offer the opportunity to clarify phylogenetic relationships that are difficult to resolve with few loci, but they can also identify genomic regions with evolutionary history distinct from that of the species history. We collected whole-genome sequence data from 29 taxa in the legume genus Medicago, then aligned these sequences to the Medicago truncatula reference genome to confidently identify 87 596 variable homologous sites. We used this data set to estimate phylogenetic relationships among Medicago species, to investigate the number of sites needed to provide robust phylogenetic estimates and to identify specific genomic regions supporting topologies in conflict with the genome-wide phylogeny. Our full genomic data set resolves relationships within the genus that were previously intractable. Subsampling the data reveals considerable variation in phylogenetic signal and power in smaller subsets of the data. Even when sampling 5000 sites, no random sample of the data supports a topology identical to that of the genome-wide phylogeny. Phylogenetic relationships estimated from 500-site sliding windows revealed genome regions supporting several alternative species relationships among recently diverged taxa, consistent with the expected effects of deep coalescence or introgression in the recent history of Medicago.

Prevalence of single nucleotide polymorphism among 27 diverse alfalfa genotypes as assessed by transcriptome sequencing

BACKGROUND

Alfalfa, a perennial, outcrossing species, is a widely planted forage legume producing highly nutritious biomass. Currently, improvement of cultivated alfalfa mainly relies on recurrent phenotypic selection. Marker assisted breeding strategies can enhance alfalfa improvement efforts, particularly if many genome-wide markers are available. Transcriptome sequencing enables efficient high-throughput discovery of single nucleotide polymorphism (SNP) markers for a complex polyploid species.

RESULT

The transcriptomes of 27 alfalfa genotypes, including elite breeding genotypes, parents of mapping populations, and unimproved wild genotypes, were sequenced using an Illumina Genome Analyzer IIx. De novo assembly of quality-filtered 72-bp reads generated 25,183 contigs with a total length of 26.8 Mbp and an average length of 1,065 bp, with an average read depth of 55.9-fold for each genotype. Overall, 21,954 (87.2%) of the 25,183 contigs represented 14,878 unique protein accessions. Gene ontology (GO) analysis suggested that a broad diversity of genes was represented in the resulting sequences. The realignment of individual reads to the contigs enabled the detection of 872,384 SNPs and 31,760 InDels. High resolution melting (HRM) analysis was used to validate 91% of 192 putative SNPs identified by sequencing. Both allelic variants at about 95% of SNP sites identified among five wild, unimproved genotypes are still present in cultivated alfalfa, and all four US breeding programs also contain a high proportion of these SNPs. Thus, little evidence exists among this dataset for loss of significant DNA sequence diversity from either domestication or breeding of alfalfa. Structure analysis indicated that individuals from the subspecies falcata, the diploid subspecies caerulea, and the tetraploid subspecies sativa (cultivated tetraploid alfalfa) were clearly separated.

CONCLUSION

We used transcriptome sequencing to discover large numbers of SNPs segregating in elite breeding populations of alfalfa. Little loss of SNP diversity was evident between unimproved and elite alfalfa germplasm. The EST and SNP markers generated from this study are publicly available at the Legume Information System ( http://medsa.comparative-legumes.org/) and can contribute to future alfalfa research and breeding applications.

Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici

The oomycete vegetable pathogen Phytophthora capsici has shown remarkable adaptation to fungicides and new hosts. Like other members of this destructive genus, P. capsici has an explosive epidemiology, rapidly producing massive numbers of asexual spores on infected hosts. In addition, P. capsici can remain dormant for years as sexually recombined oospores, making it difficult to produce crops at infested sites, and allowing outcrossing populations to maintain significant genetic variation. Genome sequencing, development of a high-density genetic map, and integrative genomic or genetic characterization of P. capsici field isolates and intercross progeny revealed significant mitotic loss of heterozygosity (LOH) in diverse isolates. LOH was detected in clonally propagated field isolates and sexual progeny, cumulatively affecting >30% of the genome. LOH altered genotypes for more than 11,000 single-nucleotide variant sites and showed a strong association with changes in mating type and pathogenicity. Overall, it appears that LOH may provide a rapid mechanism for fixing alleles and may be an important component of adaptability for P. capsici.

A comprehensive transcriptome assembly of Pigeonpea (Cajanus cajan L.) using sanger and second-generation sequencing platforms

A comprehensive transcriptome assembly for pigeonpea has been developed by analyzing 128.9 million short Illumina GA IIx single end reads, 2.19 million single end FLX/454 reads, and 18 353 Sanger expressed sequenced tags from more than 16 genotypes. The resultant transcriptome assembly, referred to as CcTA v2, comprised 21 434 transcript assembly contigs (TACs) with an N50 of 1510 bp, the largest one being ~8 kb. Of the 21 434 TACs, 16 622 (77.5%) could be mapped on to the soybean genome build 1.0.9 under fairly stringent alignment parameters. Based on knowledge of intron junctions, 10 009 primer pairs were designed from 5033 TACs for amplifying intron spanning regions (ISRs). By using in silico mapping of BAC-end-derived SSR loci of pigeonpea on the soybean genome as a reference, putative mapping positions at the chromosome level were predicted for 6284 ISR markers, covering all 11 pigeonpea chromosomes. A subset of 128 ISR markers were analyzed on a set of eight genotypes. While 116 markers were validated, 70 markers showed one to three alleles, with an average of 0.16 polymorphism information content (PIC) value. In summary, the CcTA v2 transcript assembly and ISR markers will serve as a useful resource to accelerate genetic research and breeding applications in pigeonpea.

Population genomics of the facultatively mutualistic bacteria Sinorhizobium meliloti and S. medicae

The symbiosis between rhizobial bacteria and legume plants has served as a model for investigating the genetics of nitrogen fixation and the evolution of facultative mutualism. We used deep sequence coverage (>100×) to characterize genomic diversity at the nucleotide level among 12 Sinorhizobium medicae and 32 S. meliloti strains. Although these species are closely related and share host plants, based on the ratio of shared polymorphisms to fixed differences we found that horizontal gene transfer (HGT) between these species was confined almost exclusively to plasmid genes. Three multi-genic regions that show the strongest evidence of HGT harbor genes directly involved in establishing or maintaining the mutualism with host plants. In both species, nucleotide diversity is 1.5–2.5 times greater on the plasmids than chromosomes. Interestingly, nucleotide diversity in S. meliloti but not S. medicae is highly structured along the chromosome – with mean diversity (θ_π) on one half of the chromosome five times greater than mean diversity on the other half. Based on the ratio of plasmid to chromosome diversity, this appears to be due to severely reduced diversity on the chromosome half with less diversity, which is consistent with extensive hitchhiking along with a selective sweep. Frequency-spectrum based tests identified 82 genes with a signature of adaptive evolution in one species or another but none of the genes were identified in both species. Based upon available functional information, several genes identified as targets of selection are likely to alter the symbiosis with the host plant, making them attractive targets for further functional characterization.

Facultative mutualisms are relationships between two species that can live independently, but derive benefits when living together with their mutualistic partners. The facultative mutualism between rhizobial bacteria and legume plants contributes approximately half of all biologically fixed nitrogen, an essential plant nutrient, and is an important source of nitrogen to both natural and agricultural ecosystems. We resequenced the genomes of 44 strains of two closely related species of the genus Sinorhizobium that form facultative mutualisms with the model legme Medicago truncatula. These data provide one of the most complete examinations of genomic diversity segregating within microbial species that are not causative agents of human illness. Our analyses reveal that horizontal gene transfer, a common source of new genes in microbial species, disproportionately affects genes with direct roles in the rhizobia-plant symbiosis. Analyses of nucleotide diversity segregating within each species suggests that strong selection, along with genetic hitchhiking has sharply reduced diversity along an entire chromosome half in S. meliloti. Despite the two species' ecological similarity, we did not find evidence for selection acting on the same genetic targets. In addition to providing insight into the evolutionary history of rhizobial, this study shows the feasibility and potential power of applying population genomic analyses to microbial species.

Genome-enabled insights into legume biology

Legumes are the third-largest family of angiosperms, the second-most-important crop family, and a key source of biological nitrogen in agriculture. Recently, the genome sequences of Glycine max (soybean), Medicago truncatula, and Lotus japonicus were substantially completed. Comparisons among legume genomes reveal a key role for duplication, especially a whole-genome duplication event approximately 58 Mya that is shared by most agriculturally important legumes. A second and more recent genome duplication occurred only in the lineage leading to soybean. Outcomes of genome duplication, including gene fractionation and sub- and neofunctionalization, have played key roles in shaping legume genomes and in the evolution of legume-specific traits. Analysis of legume genome sequences also enables the discovery of legume-specific gene families and provides a framework for genome-wide association mapping that will target phenotypes of special importance in legumes. Translating genomic resources from sequenced species to less studied but still important "orphan" legumes will enhance prospects for world food production.

The Medicago genome provides insight into the evolution of rhizobial symbioses

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers

Pigeonpea is an important legume food crop grown primarily by smallholder farmers in many semi-arid tropical regions of the world. We used the Illumina next-generation sequencing platform to generate 237.2 Gb of sequence, which along with Sanger-based bacterial artificial chromosome end sequences and a genetic map, we assembled into scaffolds representing 72.7% (605.78 Mb) of the 833.07 Mb pigeonpea genome. Genome analysis predicted 48,680 genes for pigeonpea and also showed the potential role that certain gene families, for example, drought tolerance-related genes, have played throughout the domestication of pigeonpea and the evolution of its ancestors. Although we found a few segmental duplication events, we did not observe the recent genome-wide duplication events observed in soybean. This reference genome sequence will facilitate the identification of the genetic basis of agronomically important traits, and accelerate the development of improved pigeonpea varieties that could improve food security in many developing countries.

The expansion of the PRAME gene family in Eutheria

The PRAME gene family belongs to the group of cancer/testis genes whose expression is restricted primarily to the testis and a variety of cancers. The expansion of this gene family as a result of gene duplication has been observed in primates and rodents. We analyzed the PRAME gene family in Eutheria and discovered a novel Y-linked PRAME gene family in bovine, PRAMEY, which underwent amplification after a lineage-specific, autosome-to-Y transposition. Phylogenetic analyses revealed two major evolutionary clades. Clade I containing the amplified PRAMEYs and the unamplified autosomal homologs in cattle and other eutherians is under stronger functional constraints; whereas, Clade II containing the amplified autosomal PRAMEs is under positive selection. Deep-sequencing analysis indicated that eight of the identified 16 PRAMEY loci are active transcriptionally. Compared to the bovine autosomal PRAME that is expressed predominantly in testis, the PRAMEY gene family is expressed exclusively in testis and is up-regulated during testicular maturation. Furthermore, the sense RNA of PRAMEY is expressed specifically whereas the antisense RNA is expressed predominantly in spermatids. This study revealed that the expansion of the PRAME family occurred in both autosomes and sex chromosomes in a lineage-dependent manner. Differential selection forces have shaped the evolution and function of the PRAME family. The positive selection observed on the autosomal PRAMEs (Clade II) may result in their functional diversification in immunity and reproduction. Conversely, selective constraints have operated on the expanded PRAMEYs to preserve their essential function in spermatogenesis.

A comparative study on PVC based sensors in determination of

The proposed work described the synthesis of neutral ionophores; 4,8-diaza-3,3,10,10-tetramethyl-1,2-dithiacyclodecane (S1), 5,8-diaza-3,3,10,10-tetramethyl-1,2 dithiacyclodecane-N,N''-diacetic acid (S2) and N,N''-bis(2,2-dimethyl-2-mercaptoethyl)ethylenediamine-N,N-diacetic acid (S3) used for comparative analysis in determination of molybdenum(v) as PVC-based membrane sensors. The best result was obtained with sensor no. 7 based on S3 ionophore with membrane composition (w/w, mg%); 5.0(S) : 30.0(PVC) : 5.0(KTpClPB) : 60.0(o-NPOE) showing a working range of 2.3 × 10-7-1.0 × 10-2 M with a detection limit of 1.2 × 10-7 M and a toleration of non-aqueous media up to 15% with slope of 11.7 mV/decade of activity. The sensor no. 7 can also be used to assess the Mo(v) concentration in different natural samples (water bodies, soils, root nodules and urine samples) with comparative analysis of ETAAS and spectrophotometric methods. The proposed sensor no. 7 can be used 2.5 months without any distortion in results after that leaching of ionophore was observed from the membrane phase in aqueous solutions of Mo(v). The proposed sensor has shown a good dynamic response time of 11 s.

ZNF280BY and ZNF280AY: autosome derived Y-chromosome gene families in Bovidae

Background

Recent progress in exploring the Y-chromosome gene content in humans, mice and cats have suggested that "autosome-to-Y" transposition of the male fertility genes is a recurrent theme during the mammalian Y-chromosome evolution. These transpositions are lineage-dependent. The purpose of this study is to investigate the lineage-specific Y-chromosome genes in bovid.

Results

We took a direct testis cDNA selection strategy and discovered two novel gene families, ZNF280BY and ZNF280AY, on the bovine (Bos taurus) Y-chromosome (BTAY), which originated from the transposition of a gene block on the bovine chromosome 17 (BTA17) and subsequently amplified. Approximately 130 active ZNF280BY loci (and ~240 pseudogenes) and ~130 pseudogenized ZNF280AY copies are present over the majority of the male-specific region (MSY). Phylogenetic analysis indicated that both gene families fit with the "birth-and-death" model of evolution. The active ZNF280BY loci share high sequence similarity and comprise three major genomic structures, resulted from insertions/deletions (indels). Assembly of a 1.2 Mb BTAY sequence in the MSY ampliconic region demonstrated that ZNF280BY and ZNF280AY, together with HSFY and TSPY families, constitute the major elements within the repeat units. The ZNF280BY gene family was found to express in different developmental stages of testis with sense RNA detected in all cell types of the seminiferous tubules while the antisense RNA detected only in the spermatids. Deep sequencing of the selected cDNAs revealed that different loci of ZNF280BY were differentially expressed up to 60-fold. Interestingly, different copies of the ZNF280AY pseudogenes were also found to differentially express up to 10-fold. However, expression level of the ZNF280AY pseudogenes was almost 6-fold lower than that of the ZNF280BY genes. ZNF280BY and ZNF280AY gene families are present in bovid, but absent in other mammalian lineages.

Conclusions

ZNF280BY and ZNF280AY are lineage-specific, multi-copy Y-gene families specific to Bovidae, and are derived from the transposition of an autosomal gene block. The temporal and spatial expression patterns of ZNF280BYs in testis suggest a role in spermatogenesis. This study offers insights into the genomic organization of the bovine MSY and gene regulation in spermatogenesis, and provides a model for studying evolution of multi-copy gene families in mammals.

The Sorghum bicolor genome and the diversification of grasses

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

Methylation-sensitive linking libraries enhance gene-enriched sequencing of complex genomes and map DNA methylation domains

Background

Many plant genomes are resistant to whole-genome assembly due to an abundance of repetitive sequence, leading to the development of gene-rich sequencing techniques. Two such techniques are hypomethylated partial restriction (HMPR) and methylation spanning linker libraries (MSLL). These libraries differ from other gene-rich datasets in having larger insert sizes, and the MSLL clones are designed to provide reads localized to "epigenetic boundaries" where methylation begins or ends.

Results

A large-scale study in maize generated 40,299 HMPR sequences and 80,723 MSLL sequences, including MSLL clones exceeding 100 kb. The paired end reads of MSLL and HMPR clones were shown to be effective in linking existing gene-rich sequences into scaffolds. In addition, it was shown that the MSLL clones can be used for anchoring these scaffolds to a BAC-based physical map. The MSLL end reads effectively identified epigenetic boundaries, as indicated by their preferential alignment to regions upstream and downstream from annotated genes. The ability to precisely map long stretches of fully methylated DNA sequence is a unique outcome of MSLL analysis, and was also shown to provide evidence for errors in gene identification. MSLL clones were observed to be significantly more repeat-rich in their interiors than in their end reads, confirming the correlation between methylation and retroelement content. Both MSLL and HMPR reads were found to be substantially gene-enriched, with the SalI MSLL libraries being the most highly enriched (31% align to an EST contig), while the HMPR clones exhibited exceptional depletion of repetitive DNA (to ~11%). These two techniques were compared with other gene-enrichment methods, and shown to be complementary.

Conclusion

MSLL technology provides an unparalleled approach for mapping the epigenetic status of repetitive blocks and for identifying sequences mis-identified as genes. Although the types and natures of epigenetic boundaries are barely understood at this time, MSLL technology flags both approximate boundaries and methylated genes that deserve additional investigation. MSLL and HMPR sequences provide a valuable resource for maize genome annotation, and are a uniquely valuable complement to any plant genome sequencing project. In order to make these results fully accessible to the community, a web display was developed that shows the alignment of MSLL, HMPR, and other gene-rich sequences to the BACs; this display is continually updated with the latest ESTs and BAC sequences.

Genetic analysis of opaque2 modifier loci in quality protein maize

Quality protein maize (QPM) was created by selecting genetic modifiers that convert the starchy endosperm of an opaque2 (o2) mutant to a hard, vitreous phenotype. Genetic analysis has shown that there are multiple, unlinked o2 modifiers (Opm), but their identity and mode of action are unknown. Using two independently developed QPM lines, we mapped several major Opm QTLs to chromosomes 1, 7 and 9. A microarray hybridization performed with RNA obtained from true breeding o2 progeny with vitreous and opaque kernel phenotypes identified a small group of differentially expressed genes, some of which map at or near the Opm QTLs. Several of the genes are associated with ethylene and ABA signaling and suggest a potential linkage of o2 endosperm modification with programmed cell death.

Physical and genetic structure of the maize genome reflects its complex evolutionary history

Maize (Zea mays L.) is one of the most important cereal crops and a model for the study of genetics, evolution, and domestication. To better understand maize genome organization and to build a framework for genome sequencing, we constructed a sequence-ready fingerprinted contig-based physical map that covers 93.5% of the genome, of which 86.1% is aligned to the genetic map. The fingerprinted contig map contains 25,908 genic markers that enabled us to align nearly 73% of the anchored maize genome to the rice genome. The distribution pattern of expressed sequence tags correlates to that of recombination. In collinear regions, 1 kb in rice corresponds to an average of 3.2 kb in maize, yet maize has a 6-fold genome size expansion. This can be explained by the fact that most rice regions correspond to two regions in maize as a result of its recent polyploid origin. Inversions account for the majority of chromosome structural variations during subsequent maize diploidization. We also find clear evidence of ancient genome duplication predating the divergence of the progenitors of maize and rice. Reconstructing the paleoethnobotany of the maize genome indicates that the progenitors of modern maize contained ten chromosomes.

As a cash crop and a model biological system, maize is of great public interest. To facilitate maize molecular breeding and its basic biology research, we built a high-resolution physical map with two different fingerprinting methods on the same set of bacterial artificial chromosome clones. The physical map was integrated to a high-density genetic map and further serves as a framework for the maize genome-sequencing project. Comparative genomics showed that the euchromatic regions between rice and maize are very conserved. Physically we delimited these conserved regions and thus detected many genome rearrangements. We defined extensively the duplication blocks within the maize genome. These blocks allowed us to reconstruct the chromosomes of the maize progenitor. We detected that maize genome has experienced two rounds of genome duplications, an ancient one before maize–rice divergence and a recent one after tetraploidization.

Analysis of EST sequences suggests recent origin of allotetraploid colonial and creeping bentgrasses

Advances in plant genomics have permitted the analysis of several members of the grass family, including the major domesticated species, and provided new insights into the evolution of the major crops on earth. Two members, colonial bentgrass (Agrostis capillaris L.) and creeping bentgrass (A. stolonifera L.) have only recently been domesticated and provide an interesting case of polyploidy and comparison to crops that have undergone human selection for thousands of years. As an initial step of characterizing these genomes, we have sampled roughly 10% of their gene content, thereby also serving as a starting point for the construction of their physical and genetic maps. Sampling mRNA from plants subjected to environmental stress showed a remarkable increase in transcription of transposable elements. Both colonial and creeping bentgrass are allotetraploids and are considered to have one genome in common, designated the A2 genome. Analysis of conserved genes present among the ESTs suggests the colonial and creeping bentgrass A2 genomes diverged from a common ancestor approximately 2.2 million years ago (MYA), thereby providing an enhanced evolutionary zoom in respect to the origin of maize, which formed 4.8 MYA, and tetraploid wheat, which formed only 0.5 MYA and is the progenitor of domesticated hexaploid wheat.

Uneven chromosome contraction and expansion in the maize genome

Maize (Zea mays or corn), both a major food source and an important cytogenetic model, evolved from a tetraploid that arose about 4.8 million years ago (Mya). As a result, maize has extensive duplicated regions within its genome. We have sequenced the two copies of one such region, generating 7.8 Mb of sequence spanning 17.4 cM of the short arm of chromosome 1 and 6.6 Mb (25.6 cM) from the long arm of chromosome 9. Rice, which did not undergo a similar whole genome duplication event, has only one orthologous region (4.9 Mb) on the short arm of chromosome 3, and can be used as reference for the maize homoeologous regions. Alignment of the three regions allowed identification of syntenic blocks, and indicated that the maize regions have undergone differential contraction in genic and intergenic regions and expansion by the insertion of retrotransposable elements. Approximately 9% of the predicted genes in each duplicated region are completely missing in the rice genome, and almost 20% have moved to other genomic locations. Predicted genes within these regions tend to be larger in maize than in rice, primarily because of the presence of predicted genes in maize with larger introns. Interestingly, the general gene methylation patterns in the maize homoeologous regions do not appear to have changed with contraction or expansion of their chromosomes. In addition, no differences in methylation of single genes and tandemly repeated gene copies have been detected. These results, therefore, provide new insights into the diploidization of polyploid species.

Structure and architecture of the maize genome

Maize (Zea mays or corn) plays many varied and important roles in society. It is not only an important experimental model plant, but also a major livestock feed crop and a significant source of industrial products such as sweeteners and ethanol. In this study we report the systematic analysis of contiguous sequences of the maize genome. We selected 100 random regions averaging 144 kb in size, representing about 0.6% of the genome, and generated a high-quality dataset for sequence analysis. This sampling contains 330 annotated genes, 91% of which are supported by expressed sequence tag data from maize and other cereal species. Genes averaged 4 kb in size with five exons, although the largest was over 59 kb with 31 exons. Gene density varied over a wide range from 0.5 to 10.7 genes per 100 kb and genes did not appear to cluster significantly. The total repetitive element content we observed (66%) was slightly higher than previous whole-genome estimates (58%-63%) and consisted almost exclusively of retroelements. The vast majority of genes can be aligned to at least one sequence read derived from gene-enrichment procedures, but only about 30% are fully covered. Our results indicate that much of the increase in genome size of maize relative to rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) is attributable to an increase in number of both repetitive elements and genes.

Whole-genome validation of high-information-content fingerprinting

Fluorescent-based high-information-content fingerprinting (HICF) techniques have recently been developed for physical mapping. These techniques make use of automated capillary DNA sequencing instruments to enable both high-resolution and high-throughput fingerprinting. In this article, we report the construction of a whole-genome HICF FPC map for maize (Zea mays subsp. mays cv B73), using a variant of HICF in which a type IIS restriction enzyme is used to generate the fluorescently labeled fragments. The HICF maize map was constructed from the same three maize bacterial artificial chromosome libraries as previously used for the whole-genome agarose FPC map, providing a unique opportunity for direct comparison of the agarose and HICF methods; as a result, it was found that HICF has substantially greater sensitivity in forming contigs. An improved assembly procedure is also described that uses automatic end-merging of contigs to reduce the effects of contamination and repetitive bands. Several new features in FPC v7.2 are presented, including shared-memory multiprocessing, which allows dramatically faster assemblies, and automatic end-merging, which permits more accurate assemblies. It is further shown that sequenced clones may be digested in silico and located accurately on the HICF assembly, despite size deviations that prevent the precise prediction of experimental fingerprints. Finally, repetitive bands are isolated, and their effect on the assembly is studied.

Characterization of the maize endosperm transcriptome and its comparison to the rice genome

The cereal endosperm is a major organ of the seed and an important component of the world's food supply. To understand the development and physiology of the endosperm of cereal seeds, we focused on the identification of genes expressed at various times during maize endosperm development. We constructed several cDNA libraries to identify full-length clones and subjected them to a twofold enrichment. A total of 23,348 high-quality sequence-reads from 5'- and 3'-ends of cDNAs were generated and assembled into a unigene set representing 5326 genes with paired sequence-reads. Additional sequencing yielded a total of 3160 (59%) completely sequenced, full-length cDNAs. From 5326 unigenes, 4139 (78%) can be aligned with 5367 predicted rice genes and by taking only the "best hit" be mapped to 3108 positions on the rice genome. The 22% unigenes not present in rice indicate a rapid change of gene content between rice and maize in only 50 million years. Differences in rice and maize gene numbers also suggest that maize has lost a large number of duplicated genes following tetraploidization. The larger number of gene copies in rice suggests that as many as 30% of its genes arose from gene amplification, which would extrapolate to a significant proportion of the estimated 44,027 candidate genes of its entire genome. Functional classification of the maize endosperm unigene set indicated that more than a fourth of the novel functionally assignable genes found in this study are involved in carbohydrate metabolism, consistent with its role as a storage organ.

Sequence composition and genome organization of maize

Zea mays L. ssp. mays, or corn, one of the most important crops and a model for plant genetics, has a genome approximately 80% the size of the human genome. To gain global insight into the organization of its genome, we have sequenced the ends of large insert clones, yielding a cumulative length of one-eighth of the genome with a DNA sequence read every 6.2 kb, thereby describing a large percentage of the genes and transposable elements of maize in an unbiased approach. Based on the accumulative 307 Mb of sequence, repeat sequences occupy 58% and genic regions occupy 7.5%. A conservative estimate predicts approximately 59,000 genes, which is higher than in any other organism sequenced so far. Because the sequences are derived from bacterial artificial chromosome clones, which are ordered in overlapping bins, tagged genes are also ordered along continuous chromosomal segments. Based on this positional information, roughly one-third of the genes appear to consist of tandemly arrayed gene families. Although the ancestor of maize arose by tetraploidization, fewer than half of the genes appear to be present in two orthologous copies, indicating that the maize genome has undergone significant gene loss since the duplication event.

Sequence analysis of the long arm of rice chromosome 11 for rice-wheat synteny

The DNA sequence of 106 BAC/PAC clones in the minimum tiling path (MTP) of the long arm of rice chromosome 11, between map positions 57.3 and 116.2 cM, has been assembled to phase 2 or PLN level. This region has been sequenced to 10x redundancy by the Indian Initiative for Rice Genome Sequencing (IIRGS) and is now publicly available in GenBank. The region, excluding overlaps, has been predicted to contain 2,932 genes using different software. A gene-by-gene BLASTN search of the NCBI wheat EST database of over 420,000 cDNA sequences revealed that 1,143 of the predicted rice genes (38.9%) have significant homology to wheat ESTs (bit score >/= 100). Further BLASTN search of these 1,143 rice genes with the GrainGenes database of sequence contigs containing bin-mapped wheat ESTs allowed 113 of the genes to be placed in bins located on wheat chromosomes of different homoeologous groups. The largest number of genes, about one-third, mapped to the homoeologous group 4 chromosomes of wheat, suggesting a common evolutionary origin. The remaining genes were located on wheat chromosomes of different groups with significantly higher numbers for groups 3 and 5. Location of bin-mapped wheat contigs to chromosomes of all the seven homoeologous groups can be ascribed to movement of genes (transpositions) or chromosome segments (translocations) within rice or the hexaploid wheat genomes. Alternatively, it could be due to ancient duplications in the common ancestral genome of wheat and rice followed by selective elimination of genes in the wheat and rice genomes. While there exists definite conservation of gene sequences and the ancestral chromosomal identity between rice and wheat, there is no obvious conservation of the gene order at this level of resolution. Lack of extensive colinearity between rice and wheat genomes suggests that there have been many insertions, deletions, duplications and translocations that make the synteny comparisons much more complicated than earlier thought. However, enhanced resolution of comparative sequence analysis may reveal smaller conserved regions of colinearity, which will facilitate selection of markers for saturation mapping and sequencing of the gene-rich regions of the wheat genome.

Mutants of Arabidopsis as tools to understand the regulation of phenylpropanoid pathway and UVB protection mechanisms

Plants accumulate certain phenylpropanoid compounds in the vacuoles of their epidermal and subepidermal cell layers thereby protecting the underlying tissue against UVB-induced damage. However, a number of mutants of Arabidopsis thaliana are known that fail to synthesize these protective pigments, thereby allowing harmful UVB radiation to penetrate into their dermal layers. Study of several of these nonlethal mutants, defective in various aspects of flavonoid and lignin biosynthesis, has led to a better understanding of the coordinate regulation and expression of important genes as well as of mechanisms involved in plant defense against UVB radiation. The characteristics of the various phenylpropanoid mutants of Arabidopsis, viz. flavonoid mutants (banyuls [ban]; increased chalcone synthase expression 1 [icx1]; transparent testa [tt] and ultraviolet sensitive [uvs]) and hydroxycinnamic acid ester mutants (ferulic acid hydroxylase 1 [fah1] and sinapoylglucose accumulator 1 [sng1]) are discussed in detail. We have briefly touched upon, wherever relevant, the unique aspects in other plant species too.

Involvement of amino acids 361 to 364 of human topoisomerase I in camptothecin resistance and enzyme catalysis

Camptothecins are antineoplastic drugs that specifically target the enzyme DNA topoisomerase I. Prior work has identified a human topoisomerase I mutation, F361S, that confers resistance to camptothecin. We now demonstrate that substitutions in the 361-364 region can alter DNA cleavage/ligation by the enzyme. The defective catalysis exhibited by certain mutants likely relates to an impaired interaction with DNA, since these enzymes are more sensitive to the inhibitory effects of DNA binding ligands. Moreover, studies with peptides and fusion proteins suggest that the 361-364 region may bind DNA directly. The finding that the 361-364 region is involved in both enzyme catalysis and camptothecin resistance suggests that this region is part of the active site of human topoisomerase I and that camptothecin may interact with the enzyme at this site.

Identification of a nucleolin binding site in human topoisomerase I

DNA topoisomerase I (topo I) is involved in the regulation of DNA supercoiling, gene transcription, and rDNA recombination. However, little is known about interactions between topo I and other nuclear proteins. We used affinity chromatography with a topo I fusion protein to screen U-937 leukemic cell extracts and have identified nucleolin as a topo I-binding protein. Coimmunoprecipitation and other studies demonstrate that the interaction between topo I and nucleolin is direct. Furthermore, deletion analyses have identified the 166-210-amino acid region of topo I as sufficient for the interaction with nucleolin. Since nucleolin has been implicated in nuclear transport and in a variety of transcriptional processes, the interaction with topo I may relate to the cellular localization of topo I or to the known role of this topoisomerase in transcription.