Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

Background: Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource. Results: The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten. Conclusions: It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.

Hydrochemical Atlas of the Arctic Ocean

Introduction: Chemical composition of water determines its physical properties and character of processes proceeding in it: freezing temperature, volume of evaporation, density, color, transparency, filtration capacity, etc. Presence of chemical elements in water solution confers waters special physical properties exerting significant influence on their circulation, creates necessary conditions for development and inhabitance of flora and fauna, and imparts to the ocean waters some chemical features that radically differ them from the land waters (Alekin & Liakhin, 1984). Hydrochemical information helps to determine elements of water circulation, convection depth, makes it easier to distinguish water masses and gives additional knowledge of climatic variability of ocean conditions. Hydrochemical information is a necessary part of biological research. Water chemical composition can be the governing characteristics determining possibility and limits of use of marine objects, both stationary and moving in sea water. Subject of investigation of hydrochemistry is study of dynamics of chemical composition, i.e. processes of its formation and hydrochemical conditions of water bodies (Alekin & Liakhin 1984). The hydrochemical processes in the Arctic Ocean are the least known. Some information on these processes can be obtained in odd publications. A generalizing study of hydrochemical conditions in the Arctic Ocean based on expeditions conducted in the years 1948-1975 has been carried out by Rusanov et al. (1979). The "Atlas of the World Ocean: the Arctic Ocean" contains a special section "Hydrochemistry" (Gorshkov, 1980). Typical vertical profiles, transects and maps for different depths - 0, 100, 300, 500, 1000, 2000, 3000 m are given in this section for the following parameters: dissolved oxygen, phosphate, silicate, pH and alkaline-chlorine coefficient. The maps were constructed using the data of expeditions conducted in the years 1948-1975. The illustrations reflect main features of distribution of the hydrochemical elements for multi-year period and represent a static image of hydrochemical conditions. Distribution of the hydrochemical elements on the ocean surface is given for two seasons - winter and summer, for the other depths are given mean annual fields. Aim of the present Atlas is description of hydrochemical conditions in the Arctic Ocean on the basis of a greater body of hydrochemical information for the years 1948-2000 and using the up-to-date methods of analysis and electronic forms of presentation of hydrochemical information. The most wide-spread characteristics determined in water samples were used as hydrochemical indices. They are: dissolved oxygen, phosphate, silicate, pH, total alkalinity, nitrite and nitrate. An important characteristics of water salt composition - "salinity" has been considered in the Oceanographic Atlas of the Arctic Ocean (1997, 1998). Presentation of the hydrochemical characteristics in this Hydrochemical Atlas is wider if compared with that of the former Atlas (Gorshkov, 1980). Maps of climatic distribution of the hydrochemical elements were constructed for all the standard depths, and seasonal variability of the hydrochemical parameters is given not only for the surface, but also for the underlying standard depths up to 400 m and including. Statistical characteristics of the hydrochemical elements are given for the first time. Detailed accuracy estimates of initial data and map construction are also given in the Atlas. Calculated values of mean-root deviations, maximum and minimum values of the parameters demonstrate limits of their variability for the analyzed period of observations. Therefore, not only investigations of chemical statics are summarized in the Atlas, but also some elements of chemical dynamics are demonstrated. Digital arrays of the hydrochemical elements obtained in nodes of a regular grid are the new form of characteristics presentation in the Atlas. It should be mentioned that the same grid and the same boxes were used in the Atlas, as those that had been used by creation of the US-Russian climatic Oceanographic Atlas. It allows to combine hydrochemical and oceanographic information of these Atlases. The first block of the digital arrays contains climatic characteristics calculated using direct observational data. These climatic characteristics were not calculated in the regions without observations, and the information arrays for these regions have gaps. The other block of climatic information in a gridded form was obtained with the help of objective analysis of observational data. Procedure of the objective analysis allowed us to obtain climatic estimates of the hydrochemical characteristics for the whole water area of the Arctic Ocean including the regions not covered by observations. Data of the objective analysis can be widely used, in particular, in hydrobiological investigations and in modeling of hydrochemical conditions of the Arctic Ocean. Array of initial measurements is a separate block. It includes all the available materials of hydrochemical observations in the form, as they were presented in different sources. While keeping in mind that this array contains some amount of perverted information, the authors of the Atlas assumed it necessary to store this information in its primary form. Methods of data quality control can be developed in future in the process of hydrochemical information accumulation. It can be supposed that attitude can vary in future to the data that were rejected according to the procedure accepted in the Atlas. The hydrochemical Atlas of the Arctic Ocean is the first specialized and electronic generalization of hydrochemical observations in the Arctic Ocean and finishes the program of joint efforts of Russian and US specialists in preparation of a number of atlases for the Arctic. The published Oceanographic Atlas (1997, 1998), Atlas of Arctic Meteorology and Climate (2000), Ice Atlas of the Arctic Ocean prepared for publication and Hydrochemical Atlas of the Arctic Ocean represent a united series of fundamental generalizations of empirical knowledge of Arctic Ocean nature at climatic level. The Hydrochemical Atlas of the Arctic Ocean was elaborated in the result of joint efforts of the SRC of the RF AARI and IARC. Dr. Ye. Nikiforov was scientific supervisor of the Atlas, Dr. R. Colony was manager on behalf of the USA and Dr. L. Timokhov - on behalf of Russia.

Double strand breaks at the HIS2 recombination hot spot in Saccharomyces cerevisiae

Double strand breaks (DSBs) have been found at several meiotic recombination hot spots in Saccharomyces cerevisiae; more global studies have found that they occur at many places along several yeast chromosomes during meiosis. Indeed, the number of breaks found is consistent with the number of recombination events predicted from the genetic map. We have previously demonstrated that the HIS2 gene is a recombination hot spot, exhibiting a high frequency of gene conversion and associated crossing over. This paper shows that DSBs occur in meiosis at a site in the coding region and at a site downstream of the HIS2 gene and that the DSBs are dependent upon genes required for recombination. The frequency of DSBs at HIS2 increases when the gene conversion frequency is increased by alterations in the DNA around HIS2, and vice versa. A deletion that increases both DSBs and conversion can stimulate both when heterozygous; that is, it is semidominant and acts to stimulate DSBs in trans. These data are consistent with the view that homologous chromosomes associate with each other before the formation of the DSBs.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

A novel modifier gene for plasma von Willebrand factor level maps to distal mouse chromosome 11

Type 1 von Willebrand disease (VWD), characterized by reduced levels of plasma von Willebrand factor (VWF), is the most common inherited bleeding disorder in humans. Penetrance of VWD is incomplete, and expression of the bleeding phenotype is highly variable. In addition, plasma VWF levels vary widely among normal individuals. To identify genes that influence VWF level, we analyzed a genetic cross between RIIIS/J and CASA/Rk, two strains of mice that exhibit a 20-fold difference in plasma VWF level. DNA samples from F2 progeny demonstrating either extremely high or extremely low plasma VWF levels were pooled and genotyped for 41 markers spanning the autosomal genome. A novel locus accounting for 63% of the total variance in VWF level was mapped to distal mouse chromosome 11, which is distinct from the murine Vwf locus on chromosome 6. We designated this locus Mvwf for “modifier of VWF.” Additional genotyping of as many as 2407 meioses established a high resolution genetic map with gene order Cola1-Itg3a-Ngfr-Mvwf/Gip-Hoxb9-Hoxb1-Cbx·rs2-Cox5a-Gfap. The Mvwf candidate interval between Ngfr and Hoxb9 is ≈0.5 centimorgan (cM). These results demonstrate that a single dominant gene accounts for the low VWF phenotype of RIIIS/J mice in crosses with several other strains. The pattern of inheritance suggests a gain-of-function mutation in a unique component of VWF biosynthesis or processing. Characterization of the human homologue for Mvwf may have relevance for a subset of type 1 VWD cases and may define an important genetic factor modifying penetrance and expression of mutations at the VWF locus.

High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato

A high-resolution physical and genetic map of a major fruit weight quantitative trait locus (QTL), fw2.2, has been constructed for a region of tomato chromosome 2. Using an F2 nearly isogenic line mapping population (3472 individuals) derived from Lycopersicon esculentum (domesticated tomato) × Lycopersicon pennellii (wild tomato), fw2.2 has been placed near TG91 and TG167, which have an interval distance of 0.13 ± 0.03 centimorgan. The physical distance between TG91 and TG167 was estimated to be ≤ 150 kb by pulsed-field gel electrophoresis of tomato DNA. A physical contig composed of six yeast artificial chromosomes (YACs) and encompassing fw2.2 was isolated. No rearrangements or chimerisms were detected within the YAC contig based on restriction fragment length polymorphism analysis using YAC-end sequences and anchored molecular markers from the high-resolution map. Based on genetic recombination events, fw2.2 could be narrowed down to a region less than 150 kb between molecular markers TG91 and HSF24 and included within two YACs: YAC264 (210 kb) and YAC355 (300 kb). This marks the first time, to our knowledge, that a QTL has been mapped with such precision and delimited to a segment of cloned DNA. The fact that the phenotypic effect of the fw2.2 QTL can be mapped to a small interval suggests that the action of this QTL is likely due to a single gene. The development of the high-resolution genetic map, in combination with the physical YAC contig, suggests that the gene responsible for this QTL and other QTLs in plants can be isolated using a positional cloning strategy. The cloning of fw2.2 will likely lead to a better understanding of the molecular biology of fruit development and to the genetic engineering of fruit size characteristics.

Subgenome chromosome walking in wheat: A 450-kb physical contig in Triticum monococcum L. spans the Lr10 resistance locus in hexaploid wheat (Triticum aestivum L.)

For many agronomically important plant genes, only their position on a genetic map is known. In the absence of an efficient transposon tagging system, such genes have to be isolated by map-based cloning. In bread wheat Triticum aestivum, the genome is hexaploid, has a size of 1.6 × 1010 bp, and contains more than 80% of repetitive sequences. So far, this genome complexity has not allowed chromosome walking and positional cloning. Here, we demonstrate that chromosome walking using bacterial artificial chromosome (BAC) clones is possible in the diploid wheat Triticum monococcum (Am genome). BAC end sequences were mostly repetitive and could not be used for the first walking step. New probes corresponding to rare low-copy sequences were efficiently identified by low-pass DNA sequencing of the BACs. Two walking steps resulted in a physical contig of 450 kb on chromosome 1AmS. Genetic mapping of the probes derived from the BAC contig demonstrated perfect colinearity between the physical map of T. monococcum and the genetic map of bread wheat on chromosome 1AS. The contig genetically spans the Lr10 leaf rust disease resistance locus in bread wheat, with 0.13 centimorgans corresponding to 300 kb between the closest flanking markers. Comparison of the genetic to physical distances has shown large variations within 350 kb of the contig. The physical contig can now be used for the isolation of the orthologous regions in bread wheat. Thus, subgenome chromosome walking in wheat can produce large physical contigs and saturate genomic regions to support positional cloning.

Multiple pathways for SOS-induced mutagenesis in Escherichia coli: An overexpression of dinB/dinP results in strongly enhancing mutagenesis in the absence of any exogenous treatment to damage DNA

dinP is an Escherichia coli gene recently identified at 5.5 min of the genetic map, whose product shows a similarity in amino acid sequence to the E. coli UmuC protein involved in DNA damage-induced mutagenesis. In this paper we show that the gene is identical to dinB, an SOS gene previously localized near the lac locus at 8 min, the function of which was shown to be required for mutagenesis of nonirradiated λ phage infecting UV-preirradiated bacterial cells (termed λUTM for λ untargeted mutagenesis). A newly constructed dinP null mutant exhibited the same defect for λUTM as observed previously with a dinB::Mu mutant, and the defect was complemented by plasmids carrying dinP as the only intact bacterial gene. Furthermore, merely increasing the dinP gene expression, without UV irradiation or any other DNA-damaging treatment, resulted in a strong enhancement of mutagenesis in F′lac plasmids; at most, 800-fold increase in the G6-to-G5 change. The enhanced mutagenesis did not depend on recA, uvrA, or umuDC. Thus, our results establish that E. coli has at least two distinct pathways for SOS-induced mutagenesis: one dependent on umuDC and the other on dinB/P.

The optimal measure of allelic association

Allelic association between pairs of loci is derived in terms of the association probability ρ as a function of recombination θ, effective population size N, linear systematic pressure v, and time t, predicting both ρrt, the decrease of association from founders and ρct, the increase by genetic drift, with ρt = ρrt + ρct. These results conform to the Malecot equation, with time replaced by distance on the genetic map, or on the physical map if recombination in the region is uniform. Earlier evidence suggested that ρ is less sensitive to variations in marker allele frequencies than alternative metrics for which there is no probability theory. This robustness is confirmed for six alternatives in eight samples. In none of these 48 tests was the residual variance as small as for ρ. Overall, efficiency was less than 80% for all alternatives, and less than 30% for two of them. Efficiency of alternatives did not increase when information was estimated simultaneously. The swept radius within which substantial values of ρ are conserved lies between 385 and 893 kb, but deviation of parameters between measures is enormously significant. The large effort now being devoted to allelic association has little value unless the ρ metric with the strongest theoretical basis and least sensitivity to marker allele frequencies is used for mapping of marker association and localization of disease loci.

Unique chromosomal regions associated with virulence of an avian pathogenic Escherichia coli strain.

The avian pathogenic Escherichia coli strain (chi)7122 (serotype O78:K80:H9) causes airsacculitis and colisepticemia in chickens. To identify genes associated with avian disease, a genomic subtraction technique was performed between strain (chi)7122 and the E. coli K-12 strain (chi)289. The DNA isolated using this method was found only in strain (chi)7122 and was used to identify cosmid clones carrying unique DNA from a library of (chi)7122 that were then used to map the position of unique DNA on the E. coli chromosome. A total of 12 unique regions were found, 5 of which correspond to previously identified positions for unique DNA sequence in E. coli strains. To assess the role each unique region plays in virulence, mutants of (chi)7122 were constructed in which a segment of unique DNA was replaced with E. coli K-12 DNA by cotransduction of linked transposon insertions in DNA flanking the unique sequence. The resulting replacement mutants were assessed for inability to colonize the air sac and cause septicemia in 2-week-old white Leghorn chickens. Two mutants were found to be avirulent when injected into the right caudal air sac of 2-week-old chickens. One avirulent mutant, designated (chi)7145, carries a replacement of the rfb locus at 44 min, generating a rough phenotype. The second mutant is designated (chi)7146, and carries a replacement at position 0.0 min on the genetic map. Both mutants could be complemented to partial virulence by cosmids carrying sequences unique to (chi)7122.

Assignment of Rfp-Y to the chicken major histocompatibility complex/NOR microchromosome and evidence for high-frequency recombination associated with the nucleolar organizer region.

Rfp-Y is a second region in the genome of the chicken containing major histocompatibility complex (MHC) class I and II genes. Haplotypes of Rfp-Y assort independently from haplotypes of the B system, a region known to function as a MHC and to be located on chromosome 16 (a microchromosome) with the single nucleolar organizer region (NOR) in the chicken genome. Linkage mapping with reference populations failed to reveal the location of Rfp-Y, leaving Rfp-Y unlinked in a map containing >400 markers. A possible location of Rfp-Y became apparent in studies of chickens trisomic for chromosome 16 when it was noted that the intensity of restriction fragments associated with Rfp-Y increased with increasing copy number of chromosome 16. Further evidence that Rfp-Y might be located on chromosome 16 was obtained when individuals trisomic for chromosome 16 were found to transmit three Rfp-Y haplotypes. Finally, mapping of cosmid cluster III of the molecular map of chicken MHC genes (containing a MHC class II gene and two rRNA genes) to Rfp-Y validated the assignment of Rfp-Y to the MHC/NOR microchromosome. A genetic map can now be drawn for a portion of chicken chromosome 16 with Rfp-Y, encompassing two MHC class I and three MHC class II genes, separated from the B system by a region containing the NOR and exhibiting highly frequent recombination.

Mutator tRNAs are encoded by the Escherichia coli mutator genes mutA and mutC: a novel pathway for mutagenesis.

We have previously described the mutator alleles mutA and mutC, which map at 95 minutes and 42 minutes, respectively, on the Escherichia coli genetic map and which stimulate transversions; the A.T-->T.A and G.C-->T.A substitutions are the most prominent. In this study we show that both mutA and mutC result from changes in the anticodon in one of four copies of the same glycine tRNA, at either the glyV or the glyW locus. This change results in a tRNA that inserts glycine at aspartic acid codons. In view of previous studies of missense suppressor tRNAs, the mistranslation of aspartic acid codons is assumed to occur at approximately 1-2%. We postulate that the mutator tRNA effect is exerted by generating a mutator polymerase and suggest that the epsilon subunit of DNA polymerase, which provides a proofreading function, is the most likely target. The implications of these findings for the contribution of mistranslation to observed spontaneous mutation rates in wild-type strains, as well as other cellular phenomena such as aging, are discussed.

Mapping the mouse genome: current status and future prospects.

The mouse is the best model system for the study of mammalian genetics and physiology. Because of the feasibility and importance of studying genetic crosses, the mouse genetic map has received tremendous attention in recent years. It currently contains over 14,000 genetically mapped markers, including 700 mutant loci, 3500 genes, and 6500 simple sequence length polymorphisms (SSLPs). The mutant loci and genes allow insights and correlations concerning physiology and development. The SSLPs provide highly polymorphic anchor points that allow inheritance to be traced in any cross and provide a scaffold for assembling physical maps. Adequate physical mapping resources--notably large-insert yeast artificial chromosome (YAC) libraries--are available to support positional cloning projects based on the genetic map, but a comprehensive physical map is still a few years away. Large-scale sequencing efforts have not yet begun in mouse, but comparative sequence analysis between mouse and human is likely to provide tremendous information about gene structure and regulation.

Analysing diversity in sugarcane resistance gene analogues

As resistance genes have been shown to contain conserved motifs and cluster in many plant genomes, the identification of resistance gene analogues can be used as a strategy for both the discovery of DNA markers linked to disease resistance loci and the map-based cloning of disease resistance genes. Sugarcane suffers from many important diseases and an analysis of resistance gene analogues offers a means to identify DNA markers linked to resistance loci. However, sugarcane has the most complex genome of any crop plant and initially it is important to understand the extent of resistance gene analogue diversity in the sugarcane genome before genetic analysis. We review herein how more than 100 expressed sequence tags with homology to different resistance genes have been identified in sugarcane with many mapped as single-dose restriction fragment length polymorphism markers. Importantly, some of these resistance gene analogues have been shown to be linked to disease resistance genes or disease quantitative trait loci. In an attempt to more efficiently analyse additional resistance gene analogues in sugarcane, we report on experiments aimed at investigating the molecular diversity of several resistance gene analogue families using a modified form of a technique termed Ecotilling. Using Ecotilling, we were able to rapidly detect single nucleotide polymorphisms in fragments amplified by PCR from four different resistance gene analogue families, SoRP1D, SoPTO, SoXa21 and SoHs1pro-1. An analysis of a diverse set of sugarcane varieties, including modern sugarcane cultivars and several S. officinarum and S. spontaneum clones, indicated that all amplicons, apart from SoHs1pro-1, contained significant polymorphism within the gene region studied. However, a comparison among these sugarcane clones, including between the parents of two sugarcane mapping populations, indicated that most polymorphisms were multi-dose, not single-dose, preventing their genetic map location or association with disease susceptibility or resistance from being determined.

The identification and characterisation of alleles of sucrose phosphate synthase gene family III in sugarcane

Little is known about the extent of allelic diversity of genes in the complex polyploid, sugarcane. Using sucrose phosphate synthase (SPS) Gene (SPS) Family III as an example, we have amplified and sequenced a 400 nt region from this gene from two sugarcane lines that are parents of a mapping population. Ten single nucleotide polymorphisms (SNPs) were identified within the 400 nt region of which seven were present in both lines. In the elite commercial cultivar Q165(A), 10 sequence haplotypes were identified, with four haplotypes recovered at 9% or greater frequency. Based on SNP presence, two clusters of haplotypes were observed. In IJ76-514, a Saccharum officinarum accession, 8 haplotypes were identified with 4 haplotypes recovered at 13% or greater frequency. Again, two clusters of haplotypes were observed. The results suggest that there may be two SPS Gene Family III genes per genome in sugarcane, each with different numbers of different alleles. This suggestion is supported by sequencing results in an elite parental sorghum line, 403463-2-1, in which 4 haplotypes, corresponding to two broad types, were also identified. Primers were designed to the sugarcane SNPs and screened over bulked DNA from high and low Sucrose-containing progeny from a cross between Q165(A) and IJ76-514. The SNP frequency did not vary in the two bulked DNA samples, suggesting that these SNPs from this SPS gene family are not associated with variation in sucrose content. Using an ecotilling approach, two of the SPS Gene Family III haplotypes were mapped to two different linkage groups in homology group 1 in Q165(A). Both haplotypes mapped near QTLs for increased sucrose content but were not themselves associated with any sugar-related trait.