Morphology and molecular phylogeny of a marine interstitial tetraflagellate with putative endosymbionts: <Emphasis Type="Italic">Auranticordis quadriverberis</Emphasis> n. gen. et sp. (Cercozoa)

Background  

Comparative morphological studies and environmental sequencing surveys indicate that marine benthic environments contain a diverse assortment of microorganisms that are just beginning to be explored and characterized. The most conspicuous predatory flagellates in these habitats range from about 20–150 μm in size and fall into three major groups of eukaryotes that are very distantly related to one another: dinoflagellates, euglenids and cercozoans. The Cercozoa is a diverse group of amoeboflagellates that cluster together in molecular phylogenies inferred mainly from ribosomal gene sequences. These molecular phylogenetic studies have demonstrated that several enigmatic taxa, previously treated as Eukaryota insertae sedis, fall within the Cercozoa, and suggest that the actual diversity of this group is largely unknown. Improved knowledge of cercozoan diversity is expected to help resolve major branches in the tree of eukaryotes and demonstrate important cellular innovations for understanding eukaryote evolution.     

Microbial communities in low permeability,high pH uranium mine tailings: characterization and potential effects

Aims

To describe the diversity and metabolic potential of microbial communities in uranium mine tailings characterized by high pH, high metal concentration and low permeability.

Methods and Results

To assess microbial diversity and their potential to influence the geochemistry of uranium mine tailings using aerobic and anaerobic culture‐based methods, in conjunction with next generation sequencing and clone library sequencing targeting two universal bacterial markers (the 16S rRNA and cpn60 genes). Growth assays revealed that 69% of the 59 distinct culturable isolates evaluated were multiple‐metal resistant, with 15% exhibiting dual‐metal hypertolerance. There was a moderately positive correlation coefficient (R = 0·43, P < 0·05) between multiple‐metal resistance of the isolates and their enzyme expression profile. Of the isolates tested, 17 reduced amorphous iron, 22 reduced molybdate and seven oxidized arsenite. Based on next generation sequencing, tailings depth was shown to influence bacterial community composition, with the difference in the microbial diversity of the upper (0–20 m) and middle (20–40 m) tailings zones being highly significant (P < 0·01) from the lower zone (40–60 m) and the difference in diversity of the upper and middle tailings zone being significant (P < 0·05). Phylotypes closely related to well‐known sulfate‐reducing and iron‐reducing bacteria were identified with low abundance, yet relatively high diversity.

Conclusions

The presence of a population of metabolically‐diverse, metal‐resistant micro‐organisms within the tailings environment, along with their demonstrated capacity for transforming metal elements, suggests that these organisms have the potential to influence the long‐term geochemistry of the tailings.

Significance and Impact of the study

This study is the first investigation of the diversity and functional potential of micro‐organisms present in low permeability, high pH uranium mine tailings.     

Distribution of Extensive nifH Gene Diversity Across Physical Soil Microenvironments

The diversity of nitrogen-fixing bacteria is well described for aquatic environments; however, terrestrial analyses remain mostly biased to rhizobial plant–microbe associations. We maximized the level of resolution for this study through the use of nucleotide sequence information extracted from a series of soil microenvironments, ranging from macroaggregates at 2000 μm to the clay fraction at <75 μm in diameter. In addition, we attempted to create an overview of the distribution of terrestrial nitrogen fixers across such microenvironments by combining culture-independent techniques with a suite of natural soil environments from uniquely different origins. Soil diazotroph diversity was analyzed phylogenetically for 600 terrestrial nifH sequences from 12 midsized clone libraries based on microenvironments of three separate soils across a global scale. Statistical analyses of nifH gene clone libraries were used to estimate coverage, establish degrees of sequence overlap, and compare cluster distributions. These analyses revealed an extensive diversity in a tropical (19 phylotypes) and an arctic soil (17 phylotypes), and moderate diversity in a temperate soil (11 phylotypes). Within each soil, comparisons across aggregate size fractions delineated nifH gene cluster shifts within populations and degrees of sequence overlap that ranged from significantly different (arctic, tropical) to significantly similar (temperate). We suggest that this is due to population separation across aggregates of different size classes, which results from differences in the temporal stability of aggregates as niches for microbial communities. This study not only provides new knowledge of the arrangement of diazotrophic communities at the soil microscale, but it also contributes to the underrepresented knowledge of soil nifH sequences in the public databases.     

Mitochondrial DNA and Y-chromosomal diversity in ancient populations of domestic sheep (Ovis aries) in Finland: comparison with contemporary sheep breeds

Background

Several molecular and population genetic studies have focused on the native sheep breeds of Finland. In this work, we investigated their ancestral sheep populations from Iron Age, Medieval and Post-Medieval periods by sequencing a partial mitochondrial DNA D-loop and the 5’-promoter region of the SRY gene. We compared the maternal (mitochondrial DNA haplotypes) and paternal (SNP oY1) genetic diversity of ancient sheep in Finland with modern domestic sheep populations in Europe and Asia to study temporal changes in genetic variation and affinities between ancient and modern populations.

Results

A 523-bp mitochondrial DNA sequence was successfully amplified for 26 of 36 sheep ancient samples i.e. five, seven and 14 samples representative of Iron Age, Medieval and Post-Medieval sheep, respectively. Genetic diversity was analyzed within the cohorts. This ancient dataset was compared with present-day data consisting of 94 animals from 10 contemporary European breeds and with GenBank DNA sequence data to carry out a haplotype sharing analysis. Among the 18 ancient mitochondrial DNA haplotypes identified, 14 were present in the modern breeds. Ancient haplotypes were assigned to the highly divergent ovine haplogroups A and B, haplogroup B being the major lineage within the cohorts. Only two haplotypes were detected in the Iron Age samples, while the genetic diversity of the Medieval and Post-Medieval cohorts was higher. For three of the ancient DNA samples, Y-chromosome SRY gene sequences were amplified indicating that they originated from rams. The SRY gene of these three ancient ram samples contained SNP G-oY1, which is frequent in modern north-European sheep breeds.

Conclusions

Our study did not reveal any sign of major population replacement of native sheep in Finland since the Iron Age. Variations in the availability of archaeological remains may explain differences in genetic diversity estimates and patterns within the cohorts rather than demographic events that occurred in the past. Our ancient DNA results fit well with the genetic context of domestic sheep as determined by analyses of modern north-European sheep breeds.     

Integrating metagenomic and amplicon databases to resolve the phylogenetic and ecological diversity of the Chlamydiae

In the era of metagenomics and amplicon sequencing, comprehensive analyses of available sequence data remain a challenge. Here we describe an approach exploiting metagenomic and amplicon data sets from public databases to elucidate phylogenetic diversity of defined microbial taxa. We investigated the phylum Chlamydiae whose known members are obligate intracellular bacteria that represent important pathogens of humans and animals, as well as symbionts of protists. Despite their medical relevance, our knowledge about chlamydial diversity is still scarce. Most of the nine known families are represented by only a few isolates, while previous clone library-based surveys suggested the existence of yet uncharacterized members of this phylum. Here we identified more than 22 000 high quality, non-redundant chlamydial 16S rRNA gene sequences in diverse databases, as well as 1900 putative chlamydial protein-encoding genes. Even when applying the most conservative approach, clustering of chlamydial 16S rRNA gene sequences into operational taxonomic units revealed an unexpectedly high species, genus and family-level diversity within the Chlamydiae, including 181 putative families. These in silico findings were verified experimentally in one Antarctic sample, which contained a high diversity of novel Chlamydiae. In our analysis, the Rhabdochlamydiaceae, whose known members infect arthropods, represents the most diverse and species-rich chlamydial family, followed by the protist-associated Parachlamydiaceae, and a putative new family (PCF8) with unknown host specificity. Available information on the origin of metagenomic samples indicated that marine environments contain the majority of the newly discovered chlamydial lineages, highlighting this environment as an important chlamydial reservoir.     

Unveiling trophic functions of uncultured protist taxa by incubation experiments in the brackish baltic sea

Background

Our knowledge of the phylogeny and diversity of aquatic protists is rapidly increasing due to molecular surveys and next-generation sequencing approaches. This has led to a considerable discrepancy between the taxa known from cultures and those known from environmental 18S rRNA gene sequences. Hence, it is generally difficult to assign ecological functions to new taxa detected by culture-independent molecular approaches.

Methodology/Principal Findings

A combination of unamended dark incubations and 18S rRNA sequencing was chosen to link molecular diversity data of uncultured protists with heterotrophic, presumably bacterivorous, growth. The incubations, conducted with Baltic Sea brackish water, resulted in a consistent shift from a protistan community dominated by phototrophs to one in which heterotrophs predominated. This was determined on the basis of cell abundance and 18S rRNA sequences derived from fingerprint analysis and clone libraries. The bulk of enriched phylotypes after incubation were related to hitherto uncultured marine taxa within chrysophytes, ochrophytes, choanoflagellates, cercozoans, and picobiliphytes, mostly represented in recently established or here defined environmental clades. Their growth in the dark, together with coinciding results from studies with a similar objective, provides evidence that these uncultured taxa represent heterotrophic or mixotrophic species.

Conclusions/Significance

These findings shed some light into the trophic role of diverse uncultured protists especially within functionally heterogeneous groups (e.g., chrysophytes, ochrophytes) and groups that appear to be puzzling with regard to their nutrition (picobiliphytes). Additionally, our results indicate that the heterotrophic flagellate community in the southwestern Baltic Sea is dominated by species of marine origin. The combination of unamended incubations with molecular diversity analysis is thus confirmed as a promising approach to explore the trophic mode of environmentally relevant protist taxa for which only sequence data are currently available.     

Development of a low bias method for characterizing viral populations using next generation sequencing technology

Background

With an estimated 38 million people worldwide currently infected with human immunodeficiency virus (HIV), and an additional 4.1 million people becoming infected each year, it is important to understand how this virus mutates and develops resistance in order to design successful therapies.

Methodology/Principal Findings

We report a novel experimental method for amplifying full-length HIV genomes without the use of sequence-specific primers for high throughput DNA sequencing, followed by assembly of full length viral genome sequences from the resulting large dataset. Illumina was chosen for sequencing due to its ability to provide greater coverage of the HIV genome compared to prior methods, allowing for more comprehensive characterization of the heterogeneity present in the HIV samples analyzed. Our novel amplification method in combination with Illumina sequencing was used to analyze two HIV populations: a homogenous HIV population based on the canonical NL4-3 strain and a heterogeneous viral population obtained from a HIV patient''s infected T cells. In addition, the resulting sequence was analyzed using a new computational approach to obtain a consensus sequence and several metrics of diversity.

Significance

This study demonstrates how a lower bias amplification method in combination with next generation DNA sequencing provides in-depth, complete coverage of the HIV genome, enabling a stronger characterization of the quasispecies present in a clinically relevant HIV population as well as future study of how HIV mutates in response to a selective pressure.     

Deep Sequencing of the Trypanosoma cruzi GP63 Surface Proteases Reveals Diversity and Diversifying Selection among Chronic and Congenital Chagas Disease Patients

Background

Chagas disease results from infection with the diploid protozoan parasite Trypanosoma cruzi. T. cruzi is highly genetically diverse, and multiclonal infections in individual hosts are common, but little studied. In this study, we explore T. cruzi infection multiclonality in the context of age, sex and clinical profile among a cohort of chronic patients, as well as paired congenital cases from Cochabamba, Bolivia and Goias, Brazil using amplicon deep sequencing technology.

Methodology/ Principal Findings

A 450bp fragment of the trypomastigote TcGP63I surface protease gene was amplified and sequenced across 70 chronic and 22 congenital cases on the Illumina MiSeq platform. In addition, a second, mitochondrial target—ND5—was sequenced across the same cohort of cases. Several million reads were generated, and sequencing read depths were normalized within patient cohorts (Goias chronic, n = 43, Goias congenital n = 2, Bolivia chronic, n = 27; Bolivia congenital, n = 20), Among chronic cases, analyses of variance indicated no clear correlation between intra-host sequence diversity and age, sex or symptoms, while principal coordinate analyses showed no clustering by symptoms between patients. Between congenital pairs, we found evidence for the transmission of multiple sequence types from mother to infant, as well as widespread instances of novel genotypes in infants. Finally, non-synonymous to synonymous (dn:ds) nucleotide substitution ratios among sequences of TcGP63Ia and TcGP63Ib subfamilies within each cohort provided powerful evidence of strong diversifying selection at this locus.

Conclusions/Significance

Our results shed light on the diversity of parasite DTUs within each patient, as well as the extent to which parasite strains pass between mother and foetus in congenital cases. Although we were unable to find any evidence that parasite diversity accumulates with age in our study cohorts, putative diversifying selection within members of the TcGP63I gene family suggests a link between genetic diversity within this gene family and survival in the mammalian host.     

Cryptic diversity in the pen shell Atrina pectinata (Bivalvia: Pinnidae): high divergence and hybridization revealed by molecular and morphological data

Cryptic species have been increasingly revealed in the marine realm through an analytical approach incorporating multiple lines of evidence (e.g., mtDNA, nuclear genes and morphology). Illustrations of cryptic taxa improve our understanding of species diversity and evolutionary histories within marine animals. The pen shell Atrina pectinata is known to exhibit extensive morphological variations that may harbour cryptic diversity. In this study, we investigated A. pectinata populations along the coast of China and one from Japan to explore possible cryptic diversity and hybridization using a combination of mitochondrial (cytochrome c oxidase subunit I, mtCOI) and nuclear (ribosomal internal transcribed spacer, nrITS) genes as well as morphology. Phylogenetic analyses of mtCOI ‘DNA barcoding gene’ sequences resolved six divergent lineages with intralineage divergences between 0.4% and 0.8%. Interlineage sequence differences ranged from 4.3% to 22.0%, suggesting that six candidate cryptic species are present. The nrITS gene revealed five deep lineages with Kimura 2‐parameter distances of 3.7–30.3%. The five nuclear lineages generally corresponded to mtCOI lineages 1–4 and (5 + 6), suggestive of five distinct evolutionary lineages. Multiple nrITS sequences of significant variance were found within an individual, clearly implying recent hybridization events between/among the evolutionary lineages, which contributed to cytonuclear discordance. Morphologically, five morphotypes matched the five genetic lineages, although the intermediates may well blur the boundaries of different morphotypes. This study demonstrates the importance of combining multiple lines of evidence to explore species cryptic diversity and past evolutionary histories.     

Massively parallel tag sequencing reveals the complexity of anaerobic marine protistan communities

Background  

Recent advances in sequencing strategies make possible unprecedented depth and scale of sampling for molecular detection of microbial diversity. Two major paradigm-shifting discoveries include the detection of bacterial diversity that is one to two orders of magnitude greater than previous estimates, and the discovery of an exciting 'rare biosphere' of molecular signatures ('species') of poorly understood ecological significance. We applied a high-throughput parallel tag sequencing (454 sequencing) protocol adopted for eukaryotes to investigate protistan community complexity in two contrasting anoxic marine ecosystems (Framvaren Fjord, Norway; Cariaco deep-sea basin, Venezuela). Both sampling sites have previously been scrutinized for protistan diversity by traditional clone library construction and Sanger sequencing. By comparing these clone library data with 454 amplicon library data, we assess the efficiency of high-throughput tag sequencing strategies. We here present a novel, highly conservative bioinformatic analysis pipeline for the processing of large tag sequence data sets.     

Fast-evolving noncoding sequences in the human genome

Background

Gene regulation is considered one of the driving forces of evolution. Although protein-coding DNA sequences and RNA genes have been subject to recent evolutionary events in the human lineage, it has been hypothesized that the large phenotypic divergence between humans and chimpanzees has been driven mainly by changes in gene regulation rather than altered protein-coding gene sequences. Comparative analysis of vertebrate genomes has revealed an abundance of evolutionarily conserved but noncoding sequences. These conserved noncoding (CNC) sequences may well harbor critical regulatory variants that have driven recent human evolution.

Results

Here we identify 1,356 CNC sequences that appear to have undergone dramatic human-specific changes in selective pressures, at least 15% of which have substitution rates significantly above that expected under neutrality. The 1,356 'accelerated CNC' (ANC) sequences are enriched in recent segmental duplications, suggesting a recent change in selective constraint following duplication. In addition, single nucleotide polymorphisms within ANC sequences have a significant excess of high frequency derived alleles and high F _ST values relative to controls, indicating that acceleration and positive selection are recent in human populations. Finally, a significant number of single nucleotide polymorphisms within ANC sequences are associated with changes in gene expression. The probability of variation in an ANC sequence being associated with a gene expression phenotype is fivefold higher than variation in a control CNC sequence.

Conclusion

Our analysis suggests that ANC sequences have until very recently played a role in human evolution, potentially through lineage-specific changes in gene regulation.     

Introduction of a novel 18S rDNA gene arrangement along with distinct ITS region in the saline water microalga Dunaliella

Background

Sundarban is the world's largest coastal sediment comprising of mangrove forest which covers about one million hectares in the south-eastern parts of India and southern parts of Bangladesh. The microbial diversity in this sediment is largely unknown till date. In the present study an attempt has been made to understand the microbial diversity in this sediment using a cultivation-independent molecular approach.

Results

Two 16 S rRNA gene libraries were constructed and partial sequencing of the selected clones was carried out to identify bacterial strains present in the sediment. Phylogenetic analysis of partially sequenced 16 S rRNA gene sequences revealed the diversity of bacterial strains in the Sundarban sediment. At least 8 different bacterial phyla were detected. The major divisions of detected bacterial phyla were Proteobacteria (alpha, beta, gamma, and delta), Flexibacteria (CFB group), Actinobacteria, Acidobacteria, Chloroflexi, Firmicutes, Planctomycetes and Gammatimonadates.

Conclusion

The gammaproteobacteria were found to be the most abundant bacterial group in Sundarban sediment. Many clones showed similarity with previously reported bacterial lineages recovered from various marine sediments. The present study indicates a probable hydrocarbon and oil contamination in this sediment. In the present study, a number of clones were identified that have shown similarity with bacterial clones or isolates responsible for the maintenance of the S-cycle in the saline environment.     

Comparison of homoeolocus organisation in paired BAC clones from white clover (<Emphasis Type="Italic">Trifolium repens</Emphasis> L.) and microcolinearity with model legume species

Background  

White clover (Trifolium repens L.) is an outbreeding allotetraploid species and an important forage legume in temperate grassland agriculture. Comparison of sub-genome architecture and study of nucleotide sequence diversity within allopolyploids provides insight into evolutionary divergence mechanisms, and is also necessary for the development of whole-genome sequencing strategies. This study aimed to evaluate the degree of divergence between the O and P' sub-genomes of white clover through sequencing of BAC clones containing paired homoeoloci. The microsyntenic relationships between the genomes of white clover and the model legumes Lotus japonicus and Medicago truncatula as well as Arabidopsis thaliana were also characterised.     

Analysis of the sequence variations in the Mhc DRB1-like gene of the endangered Humboldt penguin (Spheniscus humboldti)

The Major Histocompatibility Complex (Mhc) genomic region of many vertebrates is known to contain at least one highly polymorphic class II gene that is homologous in sequence to one or other of the human Mhc DRB1 class II genes. The diversity of the avian Mhc class II gene sequences have been extensively studied in chickens, quails, and some songbirds, but have been largely ignored in the oceanic birds, including the flightless penguins. We have previously reported that several penguin species have a high degree of polymorphism on exon 2 of the Mhc class II DRB1-like gene. In this study, we present for the first time the complete nucleotide sequences of exon 2, intron 2, and exon 3 of the DRB1-like gene of 20 Humboldt penguins, a species that is presently vulnerable to the dangers of extinction. The Humboldt DRB1-like nucleotide and amino acid sequences reveal at least eight unique alleles. Phylogenetic analysis of all the available avian DRB-like sequences showed that, of five penguin species and nine other bird species, the sequences of the Humboldt penguins grouped most closely to the Little penguin and the mallard, respectively. The present analysis confirms that the sequence variations of the Mhc class II gene, DRB1, are useful for discriminating among individuals within the same penguin population as well those within different penguin population groups and species.The nucleotide sequence and amino acid sequence data reported in this paper have been submitted to the DDBJ database and have been assigned the accession numbers AB088371–AB088374, AB089199, AB154393–AB154399, and AB162144.     

Phylogenetic diversity of transition and anoxic zone bacterial communities within a near‐shore anoxic basin: Nitinat Lake

At two stations surveyed in Nitinat Lake, a ～200‐m‐deep anoxic tidal fjord, sulfide was detected as close as 15 m from the surface. Biological characterization, determined from small subunit ribosomal RNA gene sequencing, of the chemocline and anaerobic zone revealed many sequences related to sulfur‐oxidizing bacteria, suggesting that sulfur cycling is a dominant process. γ‐ and ε‐Proteobacteria related to thiotrophic symbionts, as well as Chlorobium sp., dominated the transition zone. These are expected to play a role in dark and phototrophic CO₂ fixation, respectively. ε‐Proteobacteria phylotype abundance increased with depth, eventually comprising 69–97% of all sequences recovered from the anoxic zone. The vast majority (74%) of these phylotypes were affiliated with a novel Acrobacter sp. group (NITEP5). Quantification of NITEP5 revealed that up to 2.8 × 10⁵ cells ml^?1 were present in the anoxic zone. Surprisingly, although sequences related to known sulfate‐reducing bacteria were recovered from the transition zone, quantification of the dsr gene and ³⁵SO₄^2? uptake tests suggest that sulfate‐reduction within the water column is negligible. Overall, sequence diversity between different vertical zones was high, although the spatial segregation of γ‐Proteobacteria, Chlorobi, and ε‐Proteobacteria did not appear to vary significantly between seasons.     

Removing Noise From Pyrosequenced Amplicons

Background  

In many environmental genomics applications a homologous region of DNA from a diverse sample is first amplified by PCR and then sequenced. The next generation sequencing technology, 454 pyrosequencing, has allowed much larger read numbers from PCR amplicons than ever before. This has revolutionised the study of microbial diversity as it is now possible to sequence a substantial fraction of the 16S rRNA genes in a community. However, there is a growing realisation that because of the large read numbers and the lack of consensus sequences it is vital to distinguish noise from true sequence diversity in this data. Otherwise this leads to inflated estimates of the number of types or operational taxonomic units (OTUs) present. Three sources of error are important: sequencing error, PCR single base substitutions and PCR chimeras. We present AmpliconNoise, a development of the PyroNoise algorithm that is capable of separately removing 454 sequencing errors and PCR single base errors. We also introduce a novel chimera removal program, Perseus, that exploits the sequence abundances associated with pyrosequencing data. We use data sets where samples of known diversity have been amplified and sequenced to quantify the effect of each of the sources of error on OTU inflation and to validate these algorithms.     

Phylogenetic relationships between the chlorophyll a/b binding protein (CAB) multigene family: An intra- and interspecies study

Summary The genome ofGlycine max (L.) Merr. cv. Dare contains a chlorophyll a/b binding (Cab) protein gene family consisting of 10 genes. The primary structures of two linkedCab genes (Cab 4 andCab 5) were determined. A comparison of the nucleic acid and predicted amino acid sequences ofCab 4 andCab 5 revealed a high degree of similarity (96% and 98%, respectively). Phylogenetic inferences drawn from sequence comparisons between previously characterized soybeanCab 1, 2, and 3 andCab 4 and 5 suggested that soybeanCab 3 was an evolutionarily distant member within this family. We further investigated the molecular evolution of theCab gene family by comparing nucleotide sequences from 25 differentCab genes representing diverse phylogenetic taxa including moncot and dicot species. Phylogenetic inferences from these data support existing morphological phylogenies in that all species within one family clustered together. These data suggested that the Solanaceae were more evolutionarily distant from the monocots than the Fabaceae and Brassicaceae. In addition, these data supported the theory thatCab Type I and II genes originated prior to divergence of the monocots and dicots.     

RiTE database: a resource database for genus-wide rice genomics and evolutionary biology

Background

Comparative evolutionary analysis of whole genomes requires not only accurate annotation of gene space, but also proper annotation of the repetitive fraction which is often the largest component of most if not all genomes larger than 50 kb in size.

Results

Here we present the Rice TE database (RiTE-db) - a genus-wide collection of transposable elements and repeated sequences across 11 diploid species of the genus Oryza and the closely-related out-group Leersia perrieri. The database consists of more than 170,000 entries divided into three main types: (i) a classified and curated set of publicly-available repeated sequences, (ii) a set of consensus assemblies of highly-repetitive sequences obtained from genome sequencing surveys of 12 species; and (iii) a set of full-length TEs, identified and extracted from 12 whole genome assemblies.

Conclusions

This is the first report of a repeat dataset that spans the majority of repeat variability within an entire genus, and one that includes complete elements as well as unassembled repeats. The database allows sequence browsing, downloading, and similarity searches. Because of the strategy adopted, the RiTE-db opens a new path to unprecedented direct comparative studies that span the entire nuclear repeat content of 15 million years of Oryza diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1762-3) contains supplementary material, which is available to authorized users.