Innovation in Nucleotide-Binding Oligomerization-Like Receptor and Toll-Like Receptor Sensing Drives the Major Histocompatibility Complex-II Free Atlantic Cod Immune System.

The absence of MHC class II antigen presentation and multiple pathogen recognition receptors in the Atlantic cod has not impaired its immune response however how underlying mechanisms have adapted remains largely unknown. In this study, ex vivo cod macrophages were challenged with various bacterial and viral microbe-associated molecular patterns (MAMP) to identify major response pathways. Cytosolic MAMP-PRR pathways based upon the NOD-like receptors (NLRs) and RIG-I-like receptors (RLRs) were identified as the critical response pathways. Our analyses suggest that internalization of exogenous ligands through scavenger receptors drives both pathways activating transcription factors like NF-kB (Nuclear factor-kappa B) and interferon regulatory factors (IRFs). Further, ligand-dependent differential expression of a unique TLR25 isoform and multiple NLR paralogues suggests (sub)neofunctionalization toward specific immune defensive strategies. Our results further demonstrate that the unique immune system of the Atlantic cod provides an unprecedented opportunity to explore the evolutionary history of PRR-based signaling in vertebrate immunity.

The Grayling Genome Reveals Selection on Gene Expression Regulation after Whole-Genome Duplication.

Whole-genome duplication (WGD) has been a major evolutionary driver of increased genomic complexity in vertebrates. One such event occurred in the salmonid family ∼80 Ma (Ss4R) giving rise to a plethora of structural and regulatory duplicate-driven divergence, making salmonids an exemplary system to investigate the evolutionary consequences of WGD. Here, we present a draft genome assembly of European grayling (Thymallus thymallus) and use this in a comparative framework to study evolution of gene regulation following WGD. Among the Ss4R duplicates identified in European grayling and Atlantic salmon (Salmo salar), one-third reflect nonneutral tissue expression evolution, with strong purifying selection, maintained over ∼50 Myr. Of these, the majority reflect conserved tissue regulation under strong selective constraints related to brain and neural-related functions, as well as higher-order protein-protein interactions. A small subset of the duplicates have evolved tissue regulatory expression divergence in a common ancestor, which have been subsequently conserved in both lineages, suggestive of adaptive divergence following WGD. These candidates for adaptive tissue expression divergence have elevated rates of protein coding- and promoter-sequence evolution and are enriched for immune- and lipid metabolism ontology terms. Lastly, lineage-specific duplicate divergence points toward underlying differences in adaptive pressures on expression regulation in the nonanadromous grayling versus the anadromous Atlantic salmon. Our findings enhance our understanding of the role of WGD in genome evolution and highlight cases of regulatory divergence of Ss4R duplicates, possibly related to a niche shift in early salmonid evolution.

From gene trees to a dated allopolyploid network: insights from the angiosperm genus Viola (Violaceae).

Allopolyploidization accounts for a significant fraction of speciation events in many eukaryotic lineages. However, existing phylogenetic and dating methods require tree-like topologies and are unable to handle the network-like phylogenetic relationships of lineages containing allopolyploids. No explicit framework has so far been established for evaluating competing network topologies, and few attempts have been made to date phylogenetic networks. We used a four-step approach to generate a dated polyploid species network for the cosmopolitan angiosperm genus Viola L. (Violaceae Batch.). The genus contains ca 600 species and both recent (neo-) and more ancient (meso-) polyploid lineages distributed over 16 sections. First, we obtained DNA sequences of three low-copy nuclear genes and one chloroplast region, from 42 species representing all 16 sections. Second, we obtained fossil-calibrated chronograms for each nuclear gene marker. Third, we determined the most parsimonious multilabeled genome tree and its corresponding network, resolved at the section (not the species) level. Reconstructing the "correct" network for a set of polyploids depends on recovering all homoeologs, i.e., all subgenomes, in these polyploids. Assuming the presence of Viola subgenome lineages that were not detected by the nuclear gene phylogenies ("ghost subgenome lineages") significantly reduced the number of inferred polyploidization events. We identified the most parsimonious network topology from a set of five competing scenarios differing in the interpretation of homoeolog extinctions and lineage sorting, based on (i) fewest possible ghost subgenome lineages, (ii) fewest possible polyploidization events, and (iii) least possible deviation from expected ploidy as inferred from available chromosome counts of the involved polyploid taxa. Finally, we estimated the homoploid and polyploid speciation times of the most parsimonious network. Homoploid speciation times were estimated by coalescent analysis of gene tree node ages. Polyploid speciation times were estimated by comparing branch lengths and speciation rates of lineages with and without ploidy shifts. Our analyses recognize Viola as an old genus (crown age 31 Ma) whose evolutionary history has been profoundly affected by allopolyploidy. Between 16 and 21 allopolyploidizations are necessary to explain the diversification of the 16 major lineages (sections) of Viola, suggesting that allopolyploidy has accounted for a high percentage-between 67% and 88%-of the speciation events at this level. The theoretical and methodological approaches presented here for (i) constructing networks and (ii) dating speciation events within a network, have general applicability for phylogenetic studies of groups where allopolyploidization has occurred. They make explicit use of a hitherto underexplored source of ploidy information from chromosome counts to help resolve phylogenetic cases where incomplete sequence data hampers network inference. Importantly, the coalescent-based method used herein circumvents the assumption of tree-like evolution required by most techniques for dating speciation events.

Ancient hybridizations among the ancestral genomes of bread wheat.

The allohexaploid bread wheat genome consists of three closely related subgenomes (A, B, and D), but a clear understanding of their phylogenetic history has been lacking. We used genome assemblies of bread wheat and five diploid relatives to analyze genome-wide samples of gene trees, as well as to estimate evolutionary relatedness and divergence times. We show that the A and B genomes diverged from a common ancestor ~7 million years ago and that these genomes gave rise to the D genome through homoploid hybrid speciation 1 to 2 million years later. Our findings imply that the present-day bread wheat genome is a product of multiple rounds of hybrid speciation (homoploid and polyploid) and lay the foundation for a new framework for understanding the wheat genome as a multilevel phylogenetic mosaic.

Alterations of microbiota in urine from women with interstitial cystitis.

BACKGROUND: Interstitial Cystitis (IC) is a chronic inflammatory condition of the bladder with unknown etiology. The aim of this study was to characterize the microbial community present in the urine from IC female patients by 454 high throughput sequencing of the 16S variable regions V1V2 and V6. The taxonomical composition, richness and diversity of the IC microbiota were determined and compared to the microbial profile of asymptomatic healthy female (HF) urine. RESULTS: The composition and distribution of bacterial sequences differed between the urine microbiota of IC patients and HFs. Reduced sequence richness and diversity were found in IC patient urine, and a significant difference in the community structure of IC urine in relation to HF urine was observed. More than 90% of the IC sequence reads were identified as belonging to the bacterial genus Lactobacillus, a marked increase compared to 60% in HF urine. CONCLUSION: The 16S rDNA sequence data demonstrates a shift in the composition of the bacterial community in IC urine. The reduced microbial diversity and richness is accompanied by a higher abundance of the bacterial genus Lactobacillus, compared to HF urine. This study demonstrates that high throughput sequencing analysis of urine microbiota in IC patients is a powerful tool towards a better understanding of this enigmatic disease.

Inferring species networks from gene trees in high-polyploid North American and Hawaiian violets (Viola, Violaceae).

The phylogenies of allopolyploids take the shape of networks and cannot be adequately represented as bifurcating trees. Especially for high polyploids (i.e., organisms with more than six sets of nuclear chromosomes), the signatures of gene homoeolog loss, deep coalescence, and polyploidy may become confounded, with the result that gene trees may be congruent with more than one species network. Herein, we obtained the most parsimonious species network by objective comparison of competing scenarios involving polyploidization and homoeolog loss in a high-polyploid lineage of violets (Viola, Violaceae) mostly or entirely restricted to North America, Central America, or Hawaii. We amplified homoeologs of the low-copy nuclear gene, glucose-6-phosphate isomerase (GPI), by single-molecule polymerase chain reaction (PCR) and the chloroplast trnL-F region by conventional PCR for 51 species and subspecies. Topological incongruence among GPI homoeolog subclades, owing to deep coalescence and two instances of putative loss (or lack of detection) of homoeologs, were reconciled by applying the maximum tree topology for each subclade. The most parsimonious species network and the fossil-based calibration of the homoeolog tree favored monophyly of the high polyploids, which has resulted from allodecaploidization 9-14 Ma, involving sympatric ancestors from the extant Viola sections Chamaemelanium (diploid), Plagiostigma (paleotetraploid), and Viola (paleotetraploid). Although two of the high-polyploid lineages (Boreali-Americanae, Pedatae) remained decaploid, recurrent polyploidization with tetraploids of section Plagiostigma within the last 5 Ma has resulted in two 14-ploid lineages (Mexicanae, Nosphinium) and one 18-ploid lineage (Langsdorffianae). This implies a more complex phylogenetic and biogeographic origin of the Hawaiian violets (Nosphinium) than that previously inferred from rDNA data and illustrates the necessity of considering polyploidy in phylogenetic and biogeographic reconstruction.

A metagenomic study of methanotrophic microorganisms in Coal Oil Point seep sediments.

BACKGROUND: Methane oxidizing prokaryotes in marine sediments are believed to function as a methane filter reducing the oceanic contribution to the global methane emission. In the anoxic parts of the sediments, oxidation of methane is accomplished by anaerobic methanotrophic archaea (ANME) living in syntrophy with sulphate reducing bacteria. This anaerobic oxidation of methane is assumed to be a coupling of reversed methanogenesis and dissimilatory sulphate reduction. Where oxygen is available aerobic methanotrophs take part in methane oxidation. In this study, we used metagenomics to characterize the taxonomic and metabolic potential for methane oxidation at the Tonya seep in the Coal Oil Point area, California. Two metagenomes from different sediment depth horizons (0-4 cm and 10-15 cm below sea floor) were sequenced by 454 technology. The metagenomes were analysed to characterize the distribution of aerobic and anaerobic methanotrophic taxa at the two sediment depths. To gain insight into the metabolic potential the metagenomes were searched for marker genes associated with methane oxidation. RESULTS: Blast searches followed by taxonomic binning in MEGAN revealed aerobic methanotrophs of the genus Methylococcus to be overrepresented in the 0-4 cm metagenome compared to the 10-15 cm metagenome. In the 10-15 cm metagenome, ANME of the ANME-1 clade, were identified as the most abundant methanotrophic taxon with 8.6% of the reads. Searches for particulate methane monooxygenase (pmoA) and methyl-coenzyme M reductase (mcrA), marker genes for aerobic and anaerobic oxidation of methane respectively, identified pmoA in the 0-4 cm metagenome as Methylococcaceae related. The mcrA reads from the 10-15 cm horizon were all classified as originating from the ANME-1 clade. CONCLUSIONS: Most of the taxa detected were present in both metagenomes and differences in community structure and corresponding metabolic potential between the two samples were mainly due to abundance differences. The results suggests that the Tonya Seep sediment is a robust methane filter, where taxa presently dominating this process could be replaced by less abundant methanotrophic taxa in case of changed environmental conditions.

Evolution of plant RNA polymerase IV/V genes: evidence of subneofunctionalization of duplicated NRPD2/NRPE2-like paralogs in Viola (Violaceae).

BACKGROUND: DNA-dependent RNA polymerase IV and V (Pol IV and V) are multi-subunit enzymes occurring in plants. The origin of Pol V, specific to angiosperms, from Pol IV, which is present in all land plants, is linked to the duplication of the gene encoding the largest subunit and the subsequent subneofunctionalization of the two paralogs (NRPD1 and NRPE1). Additional duplication of the second-largest subunit, NRPD2/NRPE2, has happened independently in at least some eudicot lineages, but its paralogs are often subject to concerted evolution and gene death and little is known about their evolution nor their affinity with Pol IV and Pol V. RESULTS: We sequenced a approximately 1500 bp NRPD2/E2-like fragment from 18 Viola species, mostly paleopolyploids, and 6 non-Viola Violaceae species. Incongruence between the NRPD2/E2-like gene phylogeny and species phylogeny indicates a first duplication of NRPD2 relatively basally in Violaceae, with subsequent sorting of paralogs in the descendants, followed by a second duplication in the common ancestor of Viola and Allexis. In Viola, the mutation pattern suggested (sub-) neofunctionalization of the two NRPD2/E2-like paralogs, NRPD2/E2-a and NRPD2/E2-b. The dN/dS ratios indicated that a 54 bp region exerted strong positive selection for both paralogs immediately following duplication. This 54 bp region encodes a domain that is involved in the binding of the Nrpd2 subunit with other Pol IV/V subunits, and may be important for correct recognition of subunits specific to Pol IV and Pol V. Across all Viola taxa 73 NRPD2/E2-like sequences were obtained, of which 23 (32%) were putative pseudogenes - all occurring in polyploids. The NRPD2 duplication was conserved in all lineages except the diploid MELVIO clade, in which NRPD2/E2-b was lost, and its allopolyploid derivates from hybridization with the CHAM clade, section Viola and section Melanium, in which NRPD2/E2-a occurred in multiple copies while NRPD2/E2-b paralogs were either absent or pseudogenized. CONCLUSIONS: Following the relatively recent split of Pol IV and Pol V, our data indicate that these two multi-subunit enzymes are still in the process of specialization and each acquiring fully subfunctionalized copies of their subunit genes. Even after specialization, the NRPD2/E2-like paralogs are prone to pseudogenization and gene conversion and NRPD2 and NRPE2 copy number is a highly dynamic process modulated by allopolyploidy and gene death.

Genome signatures of colon carcinoma cell lines.

In cancer biology, cell lines are often used instead of primary tumors because of their widespread availability and close reflection of the in vivo state. Cancer is a genetic disease, commonly caused by small- and large-scale DNA rearrangements. Therefore, it is essential to know the genomic profiles of tumor cell lines to enable their correct and efficient use as experimental tools. Here, we present a comprehensive study of the genomic profiles of 20 colon cancer cell lines combining conventional karyotyping (G-banding), comparative genomic hybridization (CGH), and multicolor fluorescence in situ hybridization (M-FISH). Major differences between the microsatellite instability (MSI) and chromosome instability (CIN) cell lines are shown; the CIN cell lines exhibited complex karyotypes involving many chromosomes (mean: 8.5 copy number changes), whereas the MSI cell lines showed considerably fewer aberrations (mean: 2.6). The 3 techniques complement each other to provide a detailed picture of the numerical and structural chromosomal changes that characterize cancer cells. Therefore, 7 of the cell lines (Colo320, EB, Fri, IS2, IS3, SW480, and V9P) are here completely karyotyped for the first time and, among these, 5 have not previously been cytogenetically described. By hierarchical cluster analysis, we show that the cell lines are representative models for primary carcinomas at the genome level. We also present the genomic profiles of an experimental model for tumor progression, including 3 cell lines (IS1, IS2, and IS3) established from a primary carcinoma, its corresponding liver- and peritoneal metastasis from the same patient. To address the question of clonality, we compared the genome of 3 common cell lines grown in 2 laboratories. Finally, we compared all our results with previously published CGH data and karyotypes of colorectal cell lines. In conclusion, the large variation in genetic complexity of the cell lines highlights the importance of a comprehensive reference of genomic profiles for investigators engaged in functional studies using these research tools.

Genome profiling of breast cancer cells selected against in vitro shows copy number changes.

About 20% of breast carcinomas show no clonal chromosome abnormalities when analyzed after short-term culturing. An interesting question is whether this subset of breast carcinomas really is karyotypically normal or if selection for normal cells occurred in vitro. To address this issue, 26 breast carcinomas that had shown no cytogenetic changes by chromosome banding analysis were examined by comparative genomic hybridization (CGH), a technique that does not require culturing or tumor metaphase cells. All but one case showed copy number changes by CGH (median, four). A comparison of these findings with those of a karyotypically abnormal series analyzed using the same CGH protocol found that the cytogenetically "normal" cases were typically genetically less complex (median, four and eight, respectively; P = 0.0058). Although largely the same alterations were found in both series, some differences with respect to the frequencies of specific imbalances were seen. Gains of 3p and 6q and losses of 10q, 14q, and 17p more often were found in the cytogenetically abnormal series than in the normal tumors. We conclude that in most instances cells found to be normal by chromosome banding analysis after short-term culture do not belong to the tumor parenchyma. Furthermore, when we compared the distribution of the number of imbalances detected by CGH in the total data set according to the mitotic index in vivo (scored from 1 to 3), the median values were three, seven, and 18, respectively (P < 0.001). These data indicate not only that karyotypically normal breast carcinomas may represent a genetically simpler subgroup that grows poorly in vitro but also that this subset of tumors already has a slow growth rate in vivo.