Origin of domesticated water chestnuts (Trapa bispinosa Roxb.) and genetic variation in wild water chestnuts.

The water chestnut Trapa bispinosa Roxb. has been domesticated in China and has been reported as the only domesticated species of this genus. To understand the origin of T. bispinosa and its evolution pathway, we compared the genetic similarity and seed morphology of domesticated water chestnut T. bispinosa with three wild species T. natans, T. incisa, and T. japonica along with archeological seed samples from the Tianluoshan site (approximately 7000-6300 cal BP) in China. The largest seed size was observed only in the domesticated species, whereas other wild species showed smaller size including T. natans L. genetically close to the domesticated type, and T. incisa was the smallest in size. The volumes of the seed capsule and endosperm were measured using X ray CT scans, showing the ratios of total volumes between T. bispinosa and wild species ranged from 4.2 to 4.5. The ratios of endosperm volume ranged from 3.3 to 3.7. Both measurements showed domesticated species have larger seed volume. Genome size was indirectly estimated by flow cytometry. Domesticated species with larger seed size was estimated as diploid, as were the wild species except for tetraploid species T. japonica. Domesticated species clearly showed the largest edible organs, but it was not a result of ploidy level changes. Maternal lineages traced using complete whole chloroplast sequences, suggested that T. natans is the closest to T. bispinosa, both of which are close to T. japonica. The result was confirmed by PCR genotyping with chloroplast insertion/deletion (cpINDEL) markers developed in the study. T. incisa showed distinct plastid types within the species, and T. japonica showed a unique plastid genotype. Our study concludes the largest volumes for the edible endosperm have been accomplished through nearly 6000 years of artificial selection, but the domestication did not involve ploidy level changes.

Metagenome-assembled genomes reveal greatly expanded taxonomic and functional diversification of the abundant marine Roseobacter RCA cluster.

BACKGROUND: The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster. RESULTS: The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44-3.12 Mbp) and higher G + C contents (46.36-53.70%) than the four Pseudoplanktomarina species (2.26-2.72 Mbp, 42.22-43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species. CONCLUSIONS: Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches. Video Abstract.

Genome organization and genomics in Chlamydia: whole genome sequencing increases understanding of chlamydial virulence, evolution, and phylogeny.

The genus Chlamydia contains important obligate intracellular bacterial pathogens to humans and animals, including C. trachomatis and C. pneumoniae. Since 1998, when the first Chlamydia genome was published, our understanding of how these microbes interact, evolved and adapted to different intracellular host environments has been transformed due to the expansion of chlamydial genomes. This review explores the current state of knowledge in Chlamydia genomics and how whole genome sequencing has revolutionised our understanding of Chlamydia virulence, evolution, and phylogeny over the past two and a half decades. This review will also highlight developments in multi-omics and other approaches that have complemented whole genome sequencing to advance knowledge of Chlamydia pathogenesis and future directions for chlamydial genomics.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization.

BACKGROUND: Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization. RESULTS: Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events. CONCLUSION: These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Escherichia coli Chromosome Copy Number Measurement Using Flow Cytometry Analysis.

Flow cytometry is a high-throughput technique that analyzes individual particles as they pass through a laser beam. These particles can be individual cells and by detecting cell-scattered light their number and relative size can be measured as they pass through the beam. Labeling of molecules, usually via a fluorescent reporter, allows the amount of these molecules per cell to be quantified. DNA content can be estimated using this approach and here we describe how flow cytometry can be used to assess the DNA content of Escherichia coli cells.

The impact of flanking sequence features on DNA CpG methylation.

Epigenetics and DNA methylation play a pivotal role in many processes of the cell and we often observe that an aberrant methylation pattern characterizes pathologies. In this work we investigate the role that the flanking sequences of CGs play in the methylation process in human. We built four different CG datasets: methylated, unmethylated, and two randomly extracted ones. We evaluated features associated to the flanking sequences of those CG sets, for different size around the CG, through five measures accounting for different aspects of sequence composition complexity and structure. The analysis performed through those measures revealed evident different behaviors between methylated and unmethylated probe sets. Major differences were observed for GC content and CG dinucleotide frequency in a window size of 300-400 bp and for CG self-attraction in 3K bp. It is remarkable as the effect of methylated CG lasts much more than expected far from the CG.

AGeNNT: annotation of enzyme families by means of refined neighborhood networks.

BACKGROUND: Large enzyme families may contain functionally diverse members that give rise to clusters in a sequence similarity network (SSN). In prokaryotes, the genome neighborhood of a gene-product is indicative of its function and thus, a genome neighborhood network (GNN) deduced for an SSN provides strong clues to the specific function of enzymes constituting the different clusters. The Enzyme Function Initiative ( http://enzymefunction.org/ ) offers services that compute SSNs and GNNs. RESULTS: We have implemented AGeNNT that utilizes these services, albeit with datasets purged with respect to unspecific protein functions and overrepresented species. AGeNNT generates refined GNNs (rGNNs) that consist of cluster-nodes representing the sequences under study and Pfam-nodes representing enzyme functions encoded in the respective neighborhoods. For cluster-nodes, AGeNNT summarizes the phylogenetic relationships of the contributing species and a statistic indicates how unique nodes and GNs are within this rGNN. Pfam-nodes are annotated with additional features like GO terms describing protein function. For edges, the coverage is given, which is the relative number of neighborhoods containing the considered enzyme function (Pfam-node). AGeNNT is available at https://github.com/kandlinf/agennt . CONCLUSIONS: An rGNN is easier to interpret than a conventional GNN, which commonly contains proteins without enzymatic function and overly specific neighborhoods due to phylogenetic bias. The implemented filter routines and the statistic allow the user to identify those neighborhoods that are most indicative of a specific metabolic capacity. Thus, AGeNNT facilitates to distinguish and annotate functionally different members of enzyme families.

Genome content analysis yields new insights into the relationship between the human malaria parasite Plasmodium falciparum and its anopheline vectors.

BACKGROUND: The persistent and growing gap between the availability of sequenced genomes and the ability to assign functions to sequenced genes led us to explore ways to maximize the information content of automated annotation for studies of anopheline mosquitos. Specifically, we use genome content analysis of a large number of previously sequenced anopheline mosquitos to follow the loss and gain of protein families over the evolutionary history of this group. The importance of this endeavor lies in the potential for comparative genomic studies between Anopheles and closely related non-vector species to reveal ancestral genome content dynamics involved in vector competence. In addition, comparisons within Anopheles could identify genome content changes responsible for variation in the vectorial capacity of this family of important parasite vectors. RESULTS: The competence and capacity of P. falciparum vectors do not appear to be phylogenetically constrained within the Anophelinae. Instead, using ancestral reconstruction methods, we suggest that a previously unexamined component of vector biology, anopheline nucleotide metabolism, may contribute to the unique status of anophelines as P. falciparum vectors. While the fitness effects of nucleotide co-option by P. falciparum parasites on their anopheline hosts are not yet known, our results suggest that anopheline genome content may be responding to selection pressure from P. falciparum. Whether this response is defensive, in an attempt to redress improper nucleotide balance resulting from P. falciparum infection, or perhaps symbiotic, resulting from an as-yet-unknown mutualism between anophelines and P. falciparum, is an open question that deserves further study. CONCLUSIONS: Clearly, there is a wealth of functional information to be gained from detailed manual genome annotation, yet the rapid increase in the number of available sequences means that most researchers will not have the time or resources to manually annotate all the sequence data they generate. We believe that efforts to maximize the amount of information obtained from automated annotation can help address the functional annotation deficit that most evolutionary biologists now face, and here demonstrate the value of such an approach.

A Potential New Human Pathogen Belonging to Helicobacter Genus, Identified in a Bloodstream Infection.

We isolated from aerobic and anaerobic blood culture bottles from a febrile patient, a Helicobacter-like Gram negative, rod-shaped bacterium that MALDI-TOF MS failed to identify. Blood agar cultures incubated in a microaerobic atmosphere revealed a motile Gram negative rod, which was oxidase, catalase, nitrate reductase, esterase, and alkaline phosphatase positive. It grew at 42°C with no detectable urease activity. Antimicrobial susceptibility testing showed that the organism was susceptible to beta-lactams, gentamicin, erythromycin, and tetracycline but resistant to ciprofloxacin. Electronic microscopy analysis revealed a 3 × 0.5 µm curved rod bacterium harboring two sheathed amphitrichous flagella. Whole genome sequencing revealed a genome 1,708,265 base-pairs long with a GC content of 37.80% and a total of 1,697 coding sequences. The genomic analyses using the nucleotide sequences of the 16S rRNA gene, hsp60 and gyrB genes, as well as the GyrA protein sequence, and the results of Average Nucleotide Identity and in silico DNA-DNA hybridization suggest evidence for a novel Helicobacter species close to Helicobacter equorum and belonging to the group of enterohepatic Helicobacter species. As soon as the particular peptide mass fingerprint of this pathogen is added to the spectral databases, MALDI-TOF MS technology will improve its identification from clinical specimens, especially in case of "sterile infection". We propose to associate the present strain with the Latin name of the place of isolation; Caesarodunum (Tours, France) and suggest "Helicobacter caesarodunensis" for further description of this new bacterium.

Escherichia coli Chromosome Copy Number Measurement Using Flow Cytometry Analysis.

Flow cytometry is a high-throughput technique that analyzes individual particles as they pass through a laser beam. These particles can be individual cells and by detecting cell-scattered light their number and relative size can be measured as they pass through the beam. Labeling of molecules, usually via a fluorescent reporter, allows the amount of these molecules per cell to be quantified. DNA content can be estimated using this approach and here we describe how flow cytometry can be used to assess the DNA content of Escherichia coli cells.

Comparative Genomics of a Parthenogenesis-Inducing Wolbachia Symbiont.

Wolbachia is an intracellular symbiont of invertebrates responsible for inducing a wide variety of phenotypes in its host. These host-Wolbachia relationships span the continuum from reproductive parasitism to obligate mutualism, and provide a unique system to study genomic changes associated with the evolution of symbiosis. We present the genome sequence from a parthenogenesis-inducing Wolbachia strain (wTpre) infecting the minute parasitoid wasp Trichogramma pretiosum The wTpre genome is the most complete parthenogenesis-inducing Wolbachia genome available to date. We used comparative genomics across 16 Wolbachia strains, representing five supergroups, to identify a core Wolbachia genome of 496 sets of orthologous genes. Only 14 of these sets are unique to Wolbachia when compared to other bacteria from the Rickettsiales. We show that the B supergroup of Wolbachia, of which wTpre is a member, contains a significantly higher number of ankyrin repeat-containing genes than other supergroups. In the wTpre genome, there is evidence for truncation of the protein coding sequences in 20% of ORFs, mostly as a result of frameshift mutations. The wTpre strain represents a conversion from cytoplasmic incompatibility to a parthenogenesis-inducing lifestyle, and is required for reproduction in the Trichogramma host it infects. We hypothesize that the large number of coding frame truncations has accompanied the change in reproductive mode of the wTpre strain.

Quantifying the Relative Importance of Phylogeny and Environmental Preferences As Drivers of Gene Content in Prokaryotic Microorganisms.

Two complementary forces shape microbial genomes: vertical inheritance of genes by phylogenetic descent, and acquisition of new genes related to adaptation to particular habitats and lifestyles. Quantification of the relative importance of each driving force proved difficult. We determined the contribution of each factor, and identified particular genes or biochemical/cellular processes linked to environmental preferences (i.e., propensity of a taxon to live in particular habitats). Three types of data were confronted: (i) complete genomes, which provide gene content of different taxa; (ii) phylogenetic information, via alignment of 16S rRNA sequences, which allowed determination of the distance between taxa, and (iii) distribution of species in environments via 16S rRNA sampling experiments, reflecting environmental preferences of different taxa. The combination of these three datasets made it possible to describe and quantify the relationships among them. We found that, although phylogenetic descent was responsible for shaping most genomes, a discernible part of the latter was correlated to environmental adaptations. Particular families of genes were identified as environmental markers, as supported by direct studies such as metagenomic sequencing. These genes are likely important for adaptation of bacteria to particular conditions or habitats, such as carbohydrate or glycan metabolism genes being linked to host-associated environments.

Mitochondrial genome of Opsaridium microlepis (cypriniformes: cyprinidae).

The complete mitochondrial genome of Opsaridium microlepis was determined by using llumina sequencing method. The genome is 168 25 bp in length, comprising 13 protein-coding genes, 20 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region (D-loop region). Most genes are encoded on the heavy strand (H-strand), except for nad6 and eight tRNA genes. Base composition of the H-strand are A (27.97%), C (26.34%), G (18.37%) and T (27.32%) with AT bias of 55.29%.

Pseudo-backcrossing design for rapidly pyramiding multiple traits into a preferential rice variety.

BACKGROUND: Pyramiding multiple genes into a desirable genetic background can take years to accomplish. In this paper, a pseudo-backcrossing scheme was designed to shorten the backcrossing cycle needed. PinK3, an aromatic and potentially high-yielding rice variety-although one that is intolerant to flash flooding (Sub) and susceptible to bacterial leaf blight (BB), leaf-neck blast (BL) and the brown planthopper (BPH)-was used as a genetic basis for significant improvements through gene pyramiding. RESULTS: Four resistance donors with five target genes (Sub1A-C, xa5, Xa21, TPS and SSIIa) and three QTLs (qBph3, qBL1 and qBL11) were backcrossed individually using markers into the pseudo-recurrent parent 'PinK3' via one cycle of backcrossing followed by two cycles of pseudo-backcrossing and three selfings with rigorous foreground marker-assisted selection. In total, 29 pseudo-backcross inbred lines (BILs) were developed. Genome composition was surveyed using 61 simple sequence repeats (SSRs), 35 of which were located on six carrier chromosomes, with the remainder located on six non-carrier chromosomes. The recurrent genome content (%RGC) and donor genome content (%DGC), which were based on the physical positions of BC1F2, ranged from 69.99 to 88.98% and 11.02 to 30.01%, respectively. For the pseudo-BC3F3BILs, the %RGC and %DGC ranged from 74.50 to 81.30% and 18.70 to 25.50%, respectively. These results indicated that without direct background selection, no further increases in %RGC were obtained during pseudo-backcrossing, whereas rigorous foreground marker-assisted selection tended to reduce linkage drag during pseudo-backcrossing. The evaluation of new traits in selected pseudo-BC3F3BILs indicated significant improvements in resistance to BB, BL, BPH and Sub compared with PinK3, as well as significant improvements in grain yield (21-68%) over the donors, although yield was 7-26% lower than in 'PinK3'. All pyramided lines were aromatic and exhibited improved starch profiles, rendering them suitable for industrial food applications. CONCLUSIONS: Results show that our new pyramiding platform, which is based on marker-assisted pseudo-backcrossing, can fix five target genes and three QTLs into a high-yielding pseudo-recurrent background within seven breeding cycles in four years. This multiple pseudo-backcrossing platform decreases the time required to generate new rice varieties exhibiting complex, durable resistance to biotic and abiotic stresses in backgrounds with desirable qualities.

Streptococcus pneumoniae Supragenome Hybridization Arrays for Profiling of Genetic Content and Gene Expression.

There is extensive genomic diversity among Streptococcus pneumoniae isolates. Approximately half of the comprehensive set of genes in the species (the supragenome or pangenome) is present in all the isolates (core set), and the remaining is unevenly distributed among strains (distributed set). The Streptococcus pneumoniae Supragenome Hybridization (SpSGH) array provides coverage for an extensive set of genes and polymorphisms encountered within this species, capturing this genomic diversity. Further, the capture is quantitative. In this manner, the SpSGH array allows for both genomic and transcriptomic analyses of diverse S. pneumoniae isolates on a single platform. In this unit, we present the SpSGH array, and describe in detail its design and implementation for both genomic and transcriptomic analyses. The methodology can be applied to construction and modification of SpSGH array platforms, as well to other bacterial species as long as multiple whole-genome sequences are available that collectively capture the vast majority of the species supragenome.

Emergence of a novel subpopulation of CC398 Staphylococcus aureus infecting animals is a serious hazard for humans.

Until recently, Staphylococcus aureus from clonal complex (CC)398 were mostly described as colonizing asymptomatic raised pigs and pig-farmers. Currently, the epidemiology of the CC398 lineage is becoming more complex. CC398 human-adapted isolates are increasingly being identified in bloodstream infections in humans living in animal-free environments. In addition, CC398 isolates are increasingly responsible for invasive infections in various animals. CC398 isolates that colonize asymptomatic pigs and the isolates that infect humans living in animal-free environments (human-adapted isolates) both lack several clinically important S. aureus-associated virulence factors but differ on the basis of their prophage content. Recent findings have provided insight into the influence of a φMR11-like helper prophage on the ability of CC398 isolates to infect humans. To assess the recent spread of the CC398 lineage to various animal species and to investigate the links between the φMR11-like prophage and the emergence of CC398 isolates infecting animals, we studied 277 isolates causing infections in unrelated animals. The prevalence of CC398 isolates increased significantly between 2007 and 2013 (p < 0.001); 31.8% of the animal isolates harbored the φMR11-like prophage. High-density DNA microarray experiments with 37 representative infected-animal isolates positive for φMR11-like DNA established that most infected-animal isolates carried many genetic elements related to antimicrobial resistance and virulence genes, and a φ3 prophage encoding immune-modulating proteins and associated with animal-to-human jumps. Our findings suggest recent clonal expansion and dissemination of a new subpopulation of CC398 isolates, responsible for invasive infections in various animals, with a considerable potential to colonize and infect humans, probably greater than that of human-adapted CC398 isolates, justifying active surveillance.

Historical constraints on vertebrate genome evolution.

Recent analyses indicated that genes with larger effect of knockout or mutation and with larger probability to revert to single copy after whole genome duplication are expressed earlier in development. Here, we further investigate whether tissue specificity of gene expression is constrained by the age of origin of the corresponding genes. We use 38 metazoan genomes and a comparative genomic application system to integrate inference of gene duplication with expression data from 17,503 human genes into a strictly phylogenetic framework. We show that the number of anatomical systems in which genes are expressed decreases steadily with decreased age of the genes' first appearance in the phylogeny: the oldest genes are expressed, on average, in twice as many anatomical systems than the genes gained recently in evolution. These results are robust to different sources of expression data, to different levels of the anatomical system hierarchy, and to the use of gene families rather than duplication events. Finally, we show that the rate of increase in gene tissue specificity correlates with the relative rate of increase in the maximum number of cell types in the corresponding taxa. Although subfunctionalization and increase in cell type number throughout evolution could constitute, respectively, the proximal and ultimate causes of this correlation, the two phenomena are intermingled. Our analyses identify a striking historical constraint in gene expression: the number of cell types in existence at the time of a gene appearance (through duplication or de novo origination) tends to determine its level of tissue specificity for tens or hundreds of millions of years.