Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

Genomic characterization of closely related species in the Rumoiensis clade infers ecogenomic signatures to non-marine environments

Members of the family Vibrionaceae are generally found in marine and brackish environments, playing important roles in nutrient cycling. The Rumoiensis clade is an unconventional group in the genus Vibrio, currently comprising six species from different origins including two species isolated from non-marine environments. In this study, we performed comparative genome analysis of all six species in the clade using their complete genome sequences. We found that two non-marine species, Vibrio casei and Vibrio gangliei, lacked the genes responsible for algal polysaccharide degradation, while a number of glycoside hydrolase genes were enriched in these two species. Expansion of insertion sequences was observed in V. casei and Vibrio rumoiensis, which suggests ongoing genomic changes associated with niche adaptations. The genes responsible for the metabolism of glucosylglycerate, a compound known to play a role as compatible solutes under nitrogen limitation, were conserved across the clade. These characteristics, along with genes encoding species-specific functions, may reflect the habit expansion which has led to the current distribution of Rumoiensis clade species. Genome analysis of all species in a single clade give us valuable insights into the genomic background of the Rumoiensis clade species and emphasize the genomic diversity and versatility of Vibrionaceae.     

Analysis of mitochondrial and chloroplast genomes in two volvocine algae: Eudorina elegans and Eudorina cylindrica (Volvocaceae,Chlorophyta)

The colonial volvocine algae span the full range of organizational complexity, from four-celled species to multicellular species, and this group of algae is often used for the study of evolution. In recent years, many organelle genomes have been sequenced using the application of next generation sequencing technology; however, only a few organelle genomes have been reported in colonial volvocine algae. In this study, we determined the organelle genomes of Eudorina elegans and Eudorina cylindrica and analysed the organelle genome size, structure and gene content between these volvocine species. This provided useful information to help us understand the composition of colonial volvocine organelle genomes. Based on the chloroplast genome protein-coding genes, we conducted a phylogenomics analysis of the volvocine algae. The result revealed an unexpected phylogenetic relationship, namely, E. elegans is more closely related to Pleodorina starrii than to E. cylindrica. The substitution rate of volvocine algae was then calculated based on organelle genome protein-coding genes; our analysis suggested the possibility that the two Eudorina species may be under similar evolutionary pressure. Lastly, the synteny analysis of the mitochondrial genome showed that gene arrangements and contents are highly conserved in the family Volvocaceae, and the synteny analysis of the chloroplast genome indicated that the genus Eudorina may have experienced genomic changes.     

Horizontal gene transfer of a genetic island encoding a type III secretion system distributed in Vibrio cholerae

Twelve Vibrio cholerae isolates with genes for a type III secretion system (T3SS) were detected among 110 environmental and 14 clinical isolates. T3SS‐related genes were distributed among the various serogroups and pulsed‐field gel electrophoresis of NotI‐digested genomes showed genetic diversity in these strains. However, the restriction fragment length polymorphism profiles of the T3SS‐related genes had similar patterns. Additionally, naturally competent T3SS‐negative V. cholerae incorporated the ca. 47 kb gene cluster of T3SS, which had been integrated into a site on the chromosome by recombination. Therefore, it is suggested that horizontal gene transfer of T3SS‐related genes occurs among V. cholerae in natural ecosystems.     

Organization of the variable region of the immunoglobulin heavy-chain gene locus of the rat

We have mapped and annotated the variable region of the immunoglobulin heavy (IGH) gene locus of the Brown Norway (BN) rat (assembly V3.4; Rat Genomic Sequence Consortium). In addition to known variable region genes, we found 12 novel previously unidentified functional IGHV genes and 1 novel functional IGHD gene. In total, the variable region of the rat IGH locus is composed of at least 353 unique IGHV genes, 21 IGHD genes, and 5 IGHJ genes, of which 131, 14, and 4 are potentially functional genes, respectively. Of all species studied so far, the rat seems to have the highest number of functional IGHV genes in the genome. Rat IGHV genes can be classified into 13 IGHV families based on nucleotide sequence identity. The variable region of the BN rat spans a total length of approximately 4.9 Mb and is organized in a typical translocon organization. Like the mouse, members of the various IGHV gene families are more or less grouped together on the genome, albeit some members of IGHV gene families are found intermingled with each other. In the rat, the largest IGHV gene families are IGHV1, IGHV2, and IGHV5. The overall conclusion is that the genomic organization of the variable region of the rat IGH locus is strikingly similar to that of the mouse, illustrating the close evolutionary relationship between these two species.     

Analysis of the Codon Usage Pattern in the Vibrio Cholerae Genome

Abstract

The codon usage in the Vibrio cholerae genome is analyzed in this paper. Although there are much more genes on the chromosome 1 than on chromosome 2, the codon usage patterns of genes on the two chromosomes are quite similar, indicating that the two chromosomes may have coexisted in the same cell for a very long history. Unlike the base frequency pattern observed in other genomes, the G+C content at the third codon position of the V. cholerae genome varies in a rather small interval. The most notable feature of codon usage of V. cholerae genome is that there is a fraction of genes show significant bias in base choice at the second codon position. The 2006 known genes can be classified into two clusters according to the base frequencies at this position. The smaller cluster contains 227 genes, most of which code for proteins involved in transport and binding functions. The encoding products of these genes have significant bias in amino acids composition as compared with other genes. The codon usage patterns for the 1836 function unknown ORFs are also analyzed, which is useful to study their functions.     

Sugar Composition of Lipopolysaccharides of Family Vibrionaceae

A comparative study was carried out on the sugar composition of lipopolysaccharides (LPS) isolated from representative strains of members of the family Vibrionaceae including all of the constituting genera, i.e., Vibrio, Aeromonas, Photobacterium, Plesiomonas, and Lucibacterium. More than 100 strains were examined. It was found that, with the exception of Vibrio parahaemolyticus 06, 2-keto-3-deoxyoctonate (KDO), known generally as a component sugar in the core region of usual gram-negative bacterial LPS, is virtually absent from LPS of the Vibrionaceae strains so far examined. Furthermore, mannose was also lacking in LPS of Vibrionaceae strains with the exception of only one strain, A. anaerogenes (ATCC 15467). Instead, some KDO-like substances were found in LPS from Vibrio (“Beneckea”) nereida (ATCC 25917) and Plesiomonas shigelloides including the type strain (ATCC 14029), the same as those found in LPS from V. parahaemolyticus O7 and O12, and three strains of V. alginolyticus. These substances were strongly positive in the periodate-thiobarbituric acid test, yielding a color with maximum absorption at 549 nm. The spectra were identical to that of KDO, whereas substances differed from KDO at least in behavior in high-voltage paper electrophoresis and thin-layer chromatography. A particularly interesting feature from the chemotaxonomical point of view was found in the sugar composition of LPS isolated from V. cholerae. Fructose was present exclusively in LPS of V. cholerae (both O1 and non-O1 groups and classical and eltor biotypes) with the exception of one strain of Photobacterium phosphoreum (NCMB 844). In addition, a pair of rarely occurring amino sugars, perosamine and quinovosamine, was found in LPS from O1 group V. cholerae regardless of either the biotype (classical or eltor) or the serotype (Inaba or Ogawa), whereas this pair was not present in non-O1 group V. cholerae (the so-called NAG vibrios). This feature was confirmed with LPS from more than 30 additional strains of O1 group V. cholerae isolated from patients. The virtual absence of KDO in LPS of the family Vibrionaceae was demonstrated for the first time in this study. These results are compatible with the interpretation that the absence of KDO in LPS can be used as one of the taxonomical characteristics of Vibrionaceae in addition to (G+C) content, DNA (or RNA) homology, and numerical analysis data.     

Using Illumina next generation sequencing technologies to sequence multigene families in de novo species

The advent of Next Generation Sequencing Technology (NGST) has revolutionized molecular biology research, allowing for rapid gene/genome sequencing from a multitude of diverse species. As high throughput sequencing becomes more accessible, more efficient workflows must be developed to deal with the amounts of data produced and better assemble the genomes of de novo lineages. We combine traditional laboratory methods with Illumina NGST to amplify and sequence the largest mammalian multigene family, the Olfactory Receptor gene family, for species with and without a reference genome. We develop novel assembly methods to annotate and filter these data, which can be utilized for any gene family or any species. We find no significant difference between the ratio of genes within their respective gene families of our data compared with available genomic data. Using simulated data we explore the limitations of short‐read sequence data and our assembly in recovering this gene family. We highlight the benefits and shortcomings of these methods. Compared with data generated from traditional polymerase chain reaction, cloning and Sanger sequencing methodologies, sequence data generated using our pipeline increases yield and sequencing efficiency without reducing the number of unique genes amplified. A cloning step is not required, therefore shortening data generation time. The novel downstream methodologies and workflows described provide a tool to be utilized by many fields of biology, to access and analyze the vast quantities of data generated. By combining laboratory and in silico methods, we provide a means of extracting genomic information for multigene families without complete genome sequencing.     

Evolutionary genomics of the Fox genes: Origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.     

OrthoParaMap: Distinguishing orthologs from paralogs by integrating comparative genome data and gene phylogenies

Background  

In eukaryotic genomes, most genes are members of gene families. When comparing genes from two species, therefore, most genes in one species will be homologous to multiple genes in the second. This often makes it difficult to distinguish orthologs (separated through speciation) from paralogs (separated by other types of gene duplication). Combining phylogenetic relationships and genomic position in both genomes helps to distinguish between these scenarios. This kind of comparison can also help to describe how gene families have evolved within a single genome that has undergone polyploidy or other large-scale duplications, as in the case of Arabidopsis thaliana – and probably most plant genomes.     

A naturally occurring point mutation in the 13-mer R repeat affects the <Emphasis Type="Italic">oriC</Emphasis> function of the large chromosome of <Emphasis Type="Italic">Vibrio cholerae</Emphasis> O1 classical biotype

The genome of Vibrio cholerae consists of two circular chromosomes of different sizes. Here, a comparative analysis of the replication origins of the large chromosomes (oriCI_VC) of classical and El Tor biotypes of the pathogen is reported. Extensive nucleotide sequence analyses revealed that the oriCI_VC region has six DnaA boxes instead of the five found in Escherichia coli oriC. The additional DnaA box, designated R_V, was unique in V. cholerae as well as in other members of the family Vibrionaceae. However, R_V was not found to be essential for the autonomous replication function of the 307-bp oriCI_VC minimal region. In contrast to El Tor and the recently evolved V. cholerae O139 strains, the oriCI_VC region of the classical biotype showed only a single base transition (TG) in a highly conserved AT-rich 13-mer R repeat region. From the minichromosome copy number and its transformational efficiency analyses, it appears that the single base substitution in the oriCI_VC of the classical biotype has a significant effect on its replication initiation.     

Competence and natural transformation in vibrios

Natural transformation is a major mechanism of horizontal gene transfer in bacteria. By incorporating exogenous DNA elements into chromosomes, bacteria are able to acquire new traits that can enhance their fitness in different environments. Within the past decade, numerous studies have revealed that natural transformation is prevalent among members of the Vibrionaceae, including the pathogen Vibrio cholerae. Four environmental factors: (i) nutrient limitation, (ii) availability of extracellular nucleosides, (iii) high cell density and (iv) the presence of chitin, promote genetic competence and natural transformation in Vibrio cholerae by co‐ordinating expression of the regulators CRP, CytR, HapR and TfoX respectively. Studies of other Vibrionaceae members highlight the general importance of natural transformation within this bacterial family.     

Chromosome-level de novo genome assembly of two conifer-parasitic wasps,Megastigmus duclouxiana and Megastigmus sabinae,reveals genomic imprints of adaptation to hosts

Conifers make up about one third of global forests but are threatened by seed parasitoid wasp species. Many of these wasps belong to the genus Megastigmus, yet little is known about their genomic background. In this study, we provide chromosome-level genome assemblies for two oligophagous conifer parasitoid species of Megastigmus, which represent the first two chromosome-level genomes of the genus. The assembled genomes of Megastigmus duclouxiana and M. sabinae are 878.48 Mb (scaffold N50 of 215.60 Mb) and 812.98 Mb (scaffold N50 of 139.16 Mb), respectively, which are larger than the genome size of most hymenopterans due to the expansion of transposable elements. Expanded gene families highlight the difference in sensory-related genes between the two species, reflecting the difference in their hosts. We further found that these two species have fewer family members but more single-gene duplications than polyphagous congeners in the gene families of ATP-binding cassette transporter (ABC), cytochrome P450 (P450) and olfactory receptors (OR). These findings shed light on the pattern of adaptation to a narrow spectrum of hosts in oligophagous parasitoids. Our findings suggest potential drivers underlying genome evolution and parasitism adaptation, and provide valuable resources for understanding the ecology, genetics and evolution of Megastigmus, as well as for the research and biological control of global conifer forest pests.     

Worldwide Occurrence of Integrative Conjugative Element Encoding Multidrug Resistance Determinants in Epidemic Vibrio cholerae O1

In the last decades, there has been an increase of cholera epidemics caused by multidrug resistant strains. Particularly, the integrative and conjugative element (ICE) seems to play a major role in the emergence of multidrug resistant Vibrio cholerae. This study fully characterized, by whole genome sequencing, new ICEs carried by multidrug resistant V. cholerae O1 strains from Nigeria (2010) (ICEVchNig1) and Nepal (1994) (ICEVchNep1). The gene content and gene order of these two ICEs are the same, and identical to ICEVchInd5, ICEVchBan5 and ICEVchHai1 previously identified in multidrug resistant V. cholerae O1. This ICE is characterized by dfrA1, sul2, strAB and floR antimicrobial resistance genes, and by unique gene content in HS4 and HS5 ICE regions. Screening for ICEs, in publicly available V. cholerae genomes, revealed the occurrence and widespread distribution of this ICE among V. cholerae O1. Metagenomic analysis found segments of this ICE in marine environments far from the direct influence of the cholera epidemic. Therefore, this study revealed the epidemiology of a spatio-temporal prevalent ICE in V. cholerae O1. Its occurrence and dispersion in V. cholerae O1 strains from different continents throughout more than two decades can be indicative of its role in the fitness of the current pandemic lineage.     

Genomic and Phenotypic Diversity of Coastal Vibrio cholerae Strains Is Linked to Environmental Factors

Studies of Vibrio cholerae diversity have focused primarily on pathogenic isolates of the O1 and O139 serotypes. However, autochthonous environmental isolates of this species routinely display more extensive genetic diversity than the primarily clonal pathogenic strains. In this study, genomic and metabolic profiles of 41 non-O1/O139 environmental isolates from central California coastal waters and four clinical strains are used to characterize the core genome and metabolome of V. cholerae. Comparative genome hybridization using microarrays constructed from the fully sequenced V. cholerae O1 El Tor N16961 genome identified 2,787 core genes that approximated the projected species core genome within 1.6%. Core genes are almost universally present in strains with widely different niches, suggesting that these genes are essential for persistence in diverse aquatic environments. In contrast, the dispensable genes and phenotypic traits identified in this study should provide increased fitness for certain niche environments. Environmental parameters, measured in situ during sample collection, are correlated to the presence of specific dispensable genes and metabolic capabilities, including utilization of mannose, sialic acid, citrate, and chitosan oligosaccharides. These results identify gene content and metabolic pathways that are likely selected for in certain coastal environments and may influence V. cholerae population structure in aquatic environments.     

Cyclitol metabolism is a central feature of Burkholderia leaf symbionts

The symbioses between plants of the Rubiaceae and Primulaceae families with Burkholderia bacteria represent unique and intimate plant–bacterial relationships. Many of these interactions have been identified through PCR-dependent typing methods, but there is little information available about their functional and ecological roles. We assembled 17 new endophyte genomes representing endophytes from 13 plant species, including those of two previously unknown associations. Genomes of leaf endophytes belonging to Burkholderia s.l. show extensive signs of genome reduction, albeit to varying degrees. Except for one endophyte, none of the bacterial symbionts could be isolated on standard microbiological media. Despite their taxonomic diversity, all endophyte genomes contained gene clusters linked to the production of specialized metabolites, including genes linked to cyclitol sugar analog metabolism and in one instance non-ribosomal peptide synthesis. These genes and gene clusters are unique within Burkholderia s.l. and are likely horizontally acquired. We propose that the acquisition of secondary metabolite gene clusters through horizontal gene transfer is a prerequisite for the evolution of a stable association between these endophytes and their hosts.     

Limited functional conservation of a global regulator among related bacterial genera: Lrp in <Emphasis Type="Italic">Escherichia</Emphasis>, <Emphasis Type="Italic">Proteus</Emphasis> and <Emphasis Type="Italic">Vibrio</Emphasis>

Background  

Bacterial genome sequences are being determined rapidly, but few species are physiologically well characterized. Predicting regulation from genome sequences usually involves extrapolation from better-studied bacteria, using the hypothesis that a conserved regulator, conserved target gene, and predicted regulator-binding site in the target promoter imply conserved regulation between the two species. However many compared organisms are ecologically and physiologically diverse, and the limits of extrapolation have not been well tested. In E. coli K-12 the leucine-responsive regulatory protein (Lrp) affects expression of ~400 genes. Proteus mirabilis and Vibrio cholerae have highly-conserved lrp orthologs (98% and 92% identity to E. coli lrp). The functional equivalence of Lrp from these related species was assessed.