Genome-wide analysis of adaptive molecular evolution in the carnivorous plant Utricularia gibba

The genome of the bladderwort Utricularia gibba provides an unparalleled opportunity to uncover the adaptive landscape of an aquatic carnivorous plant with unique phenotypic features such as absence of roots, development of water-filled suction bladders, and a highly ramified branching pattern. Despite its tiny size, the U. gibba genome accommodates approximately as many genes as other plant genomes. To examine the relationship between the compactness of its genome and gene turnover, we compared the U. gibba genome with that of four other eudicot species, defining a total of 17,324 gene families (orthogroups). These families were further classified as either 1) lineage-specific expanded/contracted or 2) stable in size. The U. gibba-expanded families are generically related to three main phenotypic features: 1) trap physiology, 2) key plant morphogenetic/developmental pathways, and 3) response to environmental stimuli, including adaptations to life in aquatic environments. Further scans for signatures of protein functional specialization permitted identification of seven candidate genes with amino acid changes putatively fixed by positive Darwinian selection in the U. gibba lineage. The Arabidopsis orthologs of these genes (AXR, UMAMIT41, IGS, TAR2, SOL1, DEG9, and DEG10) are involved in diverse plant biological functions potentially relevant for U. gibba phenotypic diversification, including 1) auxin metabolism and signal transduction, 2) flowering induction and floral meristem transition, 3) root development, and 4) peptidases. Taken together, our results suggest numerous candidate genes and gene families as interesting targets for further experimental confirmation of their functional and adaptive roles in the U. gibba's unique lifestyle and highly specialized body plan.

HUMAN GENOME VARIATIONS AND EVOLUTION WITH A FOCUS ON THE ANALYSIS OF TRANSPOSABLE ELEMENTS

Genome sequence varies in numerous ways among individuals although the gross architecture is fixed for all humans. Retrotransposons create one of the most abundant structural variants in the human genome and are divided in many families, with certain members in some families, e.g., L1, Alu, SVA, and HERV-K, remaining active for transposition. Along with other types of genomic variants, retrotransponson-derived variants contribute to the whole spectrum of genome variants in humans. With the advancement of sequencing techniques, many human genomes are being sequenced at the individual level, fueling the comparative research on these variants among individuals. In this thesis, the evolution and functional impact of structural variations is examined primarily focusing on retrotransposons in the context of human evolution. The thesis comprises of three different studies on the topics that are presented in three data chapters. First, the recent evolution of all human specific AluYb members, representing the second most active subfamily of Alus, was tracked to identify their source/master copy using a novel approach. All human-specific AluYb elements from the reference genome were extracted, aligned with one another to construct clusters of similar copies and each cluster was analyzed to generate the evolutionary relationship between the members of the cluster. The approach resulted in identification of one major driver copy of all human specific Yb8 and the source copy of the Yb9 lineage. Three new subfamilies within the AluYb family – Yb8a1, Yb10 and Yb11 were also identified, with Yb11 being the youngest and most polymorphic. Second, an attempt to construct a relation between transposable elements (TEs) and tandem repeats (TRs) was made at a genome-wide scale for the first time. Upon sequence comparison, positional cross-checking and other relevant analyses, it was observed that over 20% of all TRs are derived from TEs. This result established the first connection between these two types of repetitive elements, and extends our appreciation for the impact of TEs on genomes. Furthermore, only 6% of these TE-derived TRs follow the already postulated initiation and expansion mechanisms, suggesting that the others are likely to follow a yet-unidentified mechanism. Third, by taking a combination of multiple computational approaches involving all types of genetic variations published so far including transposable elements, the first whole genome sequence of the most recent common ancestor of all modern human populations that diverged into different populations around 125,000-100,000 years ago was constructed. The study shows that the current reference genome sequence is 8.89 million base pairs larger than our common ancestor’s genome, contributed by a whole spectrum of genetic mechanisms. The use of this ancestral reference genome to facilitate the analysis of personal genomes was demonstrated using an example genome and more insightful recent evolutionary analyses involving the Neanderthal genome. The three data chapters presented in this thesis conclude that the tandem repeats and transposable elements are not two entirely distinctly isolated elements as over 20% TRs are actually derived from TEs. Certain subfamilies of TEs themselves are still evolving with the generation of newer subfamilies. The evolutionary analyses of all TEs along with other genomic variants helped to construct the genome sequence of the most recent common ancestor to all modern human populations which provides a better alternative to human reference genome and can be a useful resource for the study of personal genomics, population genetics, human and primate evolution.

Post-genomics approaches towards monitoring changes within the microbial ecology of the gut

The human gut microbiota, comprising many hundreds of different microbial species, has closely co-evolved with its human host over the millennia. Diet has been a major driver of this co-evolution, in particular dietary non-digestible carbohydrates. This dietary fraction reaches the colon and becomes available for microbial fermentation, and it is in the colon that the great diversity of gut microorganisms resides. For the vast majority of our evolutionary history humans followed hunter-gatherer life-styles and consumed diets with many times more non-digestible carbohydrates, fiber and whole plant polyphenol rich foods than typical Western style diets today.

Bacterial pathogen evolution: breaking news

The immense social and economic impact of bacterial pathogens, from drug-resistant infections in hospitals to the devastation of agricultural resources, has resulted in major investment to understand the causes and conse- quences of pathogen evolution. Recent genome se- quencing projects have provided insight into the evolution of bacterial genome structures; revealing the impact of mobile DNA on genome restructuring and pathogenicity. Sequencing of multiple genomes of relat- ed strains has enabled the delineation of pathogen evo- lution and facilitated the tracking of bacterial pathogens globally. Other recent theoretical and empirical studies have shown that pathogen evolution is significantly influenced by ecological factors, such as the distribution of hosts within the environment and the effects of co- infection. We suggest that the time is ripe for experi- mentalists to use genomics in conjunction with evolu- tionary ecology experiments to further understanding of how bacterial pathogens evolve.

Comparative genomics suggests that an ancestral polyploidy event leads to enhanced root nodule symbiosis in the Papilionoideae

Root nodule symbiosis (RNS) is one of the most efficient biological systems for nitrogen fixation and it occurs in 90% of genera in the Papilionoideae, the largest subfamily of legumes. Most papilionoid species show evidence of a polyploidy event occurred approximately 58 million years ago. Although polyploidy is considered to be an important evolutionary force in plants, the role of this papilionoid polyploidy event, especially its association with RNS, is not understood. In this study, we explored this role using an integrated comparative genomic approach and conducted gene expression comparisons and gene ontology enrichment analyses. The results show the following: (1) approximately a quarter of the papilionoid-polyploidy-derived duplicate genes are retained; (2) there is a striking divergence in the level of expression of gene duplicate pairs derived from the polyploidy event; and (3) the retained duplicates are frequently involved in the processes crucial for RNS establishment, such as symbiotic signalling, nodule organogenesis, rhizobial infection and nutrient exchange and transport. Thus, we conclude that the papilionoid polyploidy event might have further refined RNS and induced a more robust and enhanced symbiotic system. This conclusion partly explains the widespread occurrence of the Papilionoideae.

Divergent evolution of the activity and regulation of the glutamate decarboxylase systems in Listeria monocytogenes EGD-e and 10403S : roles in virulence and acid tolerance

The glutamate decarboxylase (GAD) system has been shown to be important for the survival of Listeria monocytogenes in low pH environments. The bacterium can use this faculty to maintain pH homeostasis under acidic conditions. The accepted model for the GAD system proposes that the antiport of glutamate into the bacterial cell in exchange for γ-aminobutyric acid (GABA) is coupled to an intracellular decarboxylation reaction of glutamate into GABA that consumes protons and therefore facilitates pH homeostasis. Most strains of L. monocytogenes possess three decarboxylase genes (gadD1, D2 & D3) and two antiporter genes (gadT1 & gadT2). Here, we confirm that the gadD3 encodes a glutamate decarboxylase dedicated to the intracellular GAD system (GADi), which produces GABA from cytoplasmic glutamate in the absence of antiport activity. We also compare the functionality of the GAD system between two commonly studied reference strains, EGD-e and 10403S with differences in terms of acid resistance. Through functional genomics we show that EGD-e is unable to export GABA and relies exclusively in the GADi system, which is driven primarily by GadD3 in this strain. In contrast 10403S relies upon GadD2 to maintain both an intracellular and extracellular GAD system (GADi/GADe). Through experiments with a murinised variant of EGD-e (EGDm) in mice, we found that the GAD system plays a significant role in the overall virulence of this strain. Double mutants lacking either gadD1D3 or gadD2D3 of the GAD system displayed reduced acid tolerance and were significantly affected in their ability to cause infection following oral inoculation. Since EGDm exploits GADi but not GADe the results indicate that the GADi system makes a contribution to virulence within the mouse. Furthermore, we also provide evidence that there might be a separate line of evolution in the GAD system between two commonly used reference strains.

NAD Biosynthesis Evolution in Bacteria: Lateral Gene Transfer of Kynurenine Pathway in Xanthomonadales and Flavobacteriales

The biosynthesis of quinolinate, the de novo precursor of nicotinamide adenine dinucleotide (NAD), may be performed by two distinct pathways, namely, the bacterial aspartate (aspartate-to-quinolinate) and the eukaryotic kynurenine (tryptophan-to-quinolinate). Even though the separation into eukaryotic and bacterial routes is long established, recent genomic surveys have challenged this view, because certain bacterial species also carry the genes for the kynurenine pathway. In this work, both quinolinate biosynthetic pathways were investigated in the Bacteria clade and with special attention to Xanthomonadales and Bacteroidetes, from an evolutionary viewpoint. Genomic screening has revealed that a small number of bacterial species possess some of the genes for the kynurenine pathway, which is complete in the genus Xanthomonas and in the order Flavobacteriales, where the aspartate pathway is absent. The opposite pattern (presence of the aspartate pathway and absence of the kynurenine pathway) in close relatives (Xylella ssp. and the order Bacteroidales, respectively) points to the idea of a recent acquisition of the kynurenine pathway through lateral gene transfer in these bacterial groups. In fact, sequence similarity comparison and phylogenetic reconstruction both suggest that at least part of the genes of the kynurenine pathway in Xanthomonas and Flavobacteriales is shared by eukaryotes. These results reinforce the idea of the role that lateral gene transfer plays in the configuration of bacterial genomes, thereby providing alternative metabolic pathways, even with the replacement of primary and essential cell functions, as exemplified by NAD biosynthesis.

GRASSIUS: A Platform for Comparative Regulatory Genomics across the Grasses

Transcription factors (TFs) are major players in gene regulatory networks and interactions between TFs and their target genes furnish spatiotemporal patterns of gene expression. Establishing the architecture of regulatory networks requires gathering information on TFs, their targets in the genome, and the corresponding binding sites. We have developed GRASSIUS (Grass Regulatory Information Services) as a knowledge-based Web resource that integrates information on TFs and gene promoters across the grasses. In its initial implementation, GRASSIUS consists of two separate, yet linked, databases. GrassTFDB holds information on TFs from maize (Zea mays), sorghum (Sorghum bicolor), sugarcane (Saccharum spp.), and rice (Oryza sativa). TFs are classified into families and phylogenetic relationships begin to uncover orthologous relationships among the participating species. This database also provides a centralized clearinghouse for TF synonyms in the grasses. GrassTFDB is linked to the grass TFome collection, which provides clones in recombination-based vectors corresponding to full-length open reading frames for a growing number of grass TFs. GrassPROMDB contains promoter and cis-regulatory element information for those grass species and genes for which enough data are available. The integration of GrassTFDB and GrassPROMDB will be accomplished through GrassRegNet as a first step in representing the architecture of grass regulatory networks. GRASSIUS can be accessed from www.grassius.org.

Comparative chromosome mapping of repetitive sequences. Implications for genomic evolution in the fish, Hoplias malabaricus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Vertical inheritance and bursts of transposition have shaped the evolution of the BS non-LTR retrotransposon in Drosophila

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Chromosomal evolution of neotropical cichlids: The role of repetitive DNA sequences in the organization and structure of karyotype

Cichlids are important in the aquaculture and ornamental fish trade and are considered models for evolutionary biology. However, most studies of cichlids have investigated African species, and the South American cichlids remain poorly characterized. Studies in neotropical regions have focused almost exclusively on classical cytogenetic approaches without investigating physical chromosomal mapping of specific sequences. The aim of the present study is to investigate the genomic organization of species belonging to different tribes of the subfamily Cichlinae (Cichla monoculus, Astronotus ocellatus, Geophagus proximus, Acaronia nassa, Bujurquina peregrinabunda, Hoplarchus psittacus, Hypselecara coryphaenoides, Hypselecara temporalis, Caquetaia spectabilis, Uaru amphiacanthoides, Pterophyllum leopoldi, Pterophyllum scalare, and Symphysodon discus) and reexamine the karyotypic evolutionary patterns proposed for this group. Variations in some cytogenetic markers were observed, although no trends were found in terms of the increase, decrease, or maintenance of the basal diploid chromosome number 2n = 48 in the tribes. Several species were observed to have 18S rDNA genetic duplications, as well as multiple rDNA loci. In most of the taxa analyzed, the 5S rDNA was located in the interstitial region of a pair of homologous chromosomes, although variations from this pattern were observed. Interstitial telomere sites were also observed and appear to be involved in chromosomal rearrangement events and the accumulation of repeat-rich satellite DNA sequences. Our data demonstrated the karyotypic diversity that exists among neotropical cichlids, suggesting that most of this diversity is due to the repetitive sequences present in heterochromatic regions and that repeat sequences have greatly influenced the karyotypic evolution of these fishes. © 2012 Springer Science+Business Media B.V.

Evolution and genomic organization of muscle microRNAs in fish genomes

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Comparative Genomics of Early-Diverging Brucella Strains Reveals a Novel Lipopolysaccharide Biosynthesis Pathway

Brucella species are Gram-negative bacteria that infect mammals. Recently, two unusual strains (Brucella inopinata BO1(T) and B. inopinata-like BO2) have been isolated from human patients, and their similarity to some atypical brucellae isolated from Australian native rodent species was noted. Here we present a phylogenomic analysis of the draft genome sequences of BO1(T) and BO2 and of the Australian rodent strains 83-13 and NF2653 that shows that they form two groups well separated from the other sequenced Brucella spp. Several important differences were noted. Both BO1(T) and BO2 did not agglutinate significantly when live or inactivated cells were exposed to monospecific A and Mantisera against O-side chain sugars composed of N-formyl-perosamine. While BO1(T) maintained the genes required to synthesize a typical Brucella O-antigen, BO2 lacked many of these genes but still produced a smooth LPS (lipopolysaccharide). Most missing genes were found in the wbk region involved in O-antigen synthesis in classic smooth Brucella spp. In their place, BO2 carries four genes that other bacteria use for making a rhamnose-based O-antigen. Electrophoretic, immunoblot, and chemical analyses showed that BO2 carries an antigenically different O-antigen made of repeating hexose-rich oligosaccharide units that made the LPS water-soluble, which contrasts with the homopolymeric O-antigen of other smooth brucellae that have a phenol-soluble LPS. The results demonstrate the existence of a group of early-diverging brucellae with traits that depart significantly from those of the Brucella species described thus far. IMPORTANCE This report examines differences between genomes from four new Brucella strains and those from the classic Brucella spp. Our results show that the four new strains are outliers with respect to the previously known Brucella strains and yet are part of the genus, forming two new clades. The analysis revealed important information about the evolution and survival mechanisms of Brucella species, helping reshape our knowledge of this important zoonotic pathogen. One discovery of special importance is that one of the strains, BO2, produces an O-antigen distinct from any that has been seen in any other Brucella isolates to date.