Gene Synteny Research Articles

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma.

Bees provide indispensable pollination services to both agricultural crops and wild plant populations, and several species of bees have become important models for the study of learning and memory, plant–insect interactions, and social behavior. Orchid bees (Apidae: Euglossini) are especially important to the fields of pollination ecology, evolution, and species conservation. Here we report the nuclear and mitochondrial genome sequences of the orchid bee Euglossa dilemma Bembé & Eltz. E. dilemma was selected because it is widely distributed, highly abundant, and it was recently naturalized in the southeastern United States. We provide a high-quality assembly of the 3.3 Gb genome, and an official gene set of 15,904 gene annotations. We find high conservation of gene synteny with the honey bee throughout 80 MY of divergence time. This genomic resource represents the first draft genome of the orchid bee genus Euglossa, and the first draft orchid bee mitochondrial genome, thus representing a valuable resource to the research community.

Plastid genome structure and phylogenomics of Nymphaeales: conserved gene order and new insights into relationships

The plastid genomes of early-diverging angiosperms were among the first land plant plastomes investigated. Despite their importance to understanding angiosperm evolution, no investigation has so far compared gene content or gene synteny of these plastid genomes with a focus on the Nymphaeales. Here, we report an evaluation and comparison of gene content, gene synteny and inverted repeat length for a set of 15 plastid genomes of early-diverging angiosperms. Seven plastid genomes of the Nymphaeales were newly sequenced for this investigation. We compare gene order and inverted repeat (IR) length across all genomes, review the gene annotations of previously published genomes, generate a multi-gene alignment of 77 plastid-encoded genes and reconstruct the phylogenetic relationships of the taxa under study. Our results show that gene content and synteny are highly conserved across early-diverging angiosperms: All species analyzed display complete gene synteny when accounting for expansions and contractions of the IRs. This conservation was initially obscured by ambiguous and potentially incorrect gene annotations in previously published genomes. We also report the presence of intact open reading frames across all taxa analyzed. The multi-gene phylogeny displays maximum support for the families Cabombaceae and Hydatellaceae, but no support for a clade of all Nymphaeaceae. It further indicates that the genus Victoria is embedded within Nymphaea. Plastid genomes of Trithuria were found to deviate by numerous substitutions and length changes in the IRs. Phylogenetic analyses further indicate that a previously published plastome named Nymphaea mexicana falls into a clade of N. odorata and should be re-evaluated.

Evolutionary divergence of the vertebrate TNFAIP8 gene family: Applying the spotted gar orthology bridge to understand ohnolog loss in teleosts.

Comparative functional genomic studies require the proper identification of gene orthologs to properly exploit animal biomedical research models. To identify gene orthologs, comprehensive, conserved gene synteny analyses are necessary to unwind gene histories that are convoluted by two rounds of early vertebrate genome duplication, and in the case of the teleosts, a third round, the teleost genome duplication (TGD). Recently, the genome of the spotted gar, a holostean outgroup to the teleosts that did not undergo this third genome duplication, was sequenced and applied as an orthology bridge to facilitate the identification of teleost orthologs to human genes and to enhance the power of teleosts as biomedical models. In this study, we apply the spotted gar orthology bridge to help describe the gene history of the vertebrate TNFAIP8 family. Members of the TNFAIP8 gene family have been linked to regulation of immune function and homeostasis and the development of multiple cancer types. Through a conserved gene synteny analysis, we identified zebrafish orthologs to human TNFAIP8L1 and TNFAIP8L3 genes and two co-orthologs to human TNFAIP8L2, but failed to identify an ortholog to human TNFAIP8. Through the application of the orthology bridge, we determined that teleost orthologs to human TNFAIP8 genes were likely lost in a genome inversion event after their divergence from their common ancestor with spotted gar. These findings demonstrate the value of this enhanced approach to gene history analysis and support the development of teleost models to study complex questions related to an array of biomedical issues, including immunity and cancer.

Novel resistance to Cydia pomonella granulovirus (CpGV) in codling moth shows autosomal and dominant inheritance and confers cross-resistance to different CpGV genome groups.

Commercial Cydia pomonella granulovirus (CpGV) products have been successfully applied to control codling moth (CM) in organic and integrated fruit production for more than 30 years. Since 2005, resistance against the widely used isolate CpGV-M has been reported from different countries in Europe. The inheritance of this so-called type I resistance is dominant and linked to the Z chromosome. Recently, a second form (type II) of CpGV resistance in CM was reported from a field population (NRW-WE) in Germany. Type II resistance confers reduced susceptibility not only to CpGV-M but to most known CpGV isolates and it does not follow the previously described Z-linked inheritance of type I resistance. To further analyze type II resistance, two CM strains, termed CpR5M and CpR5S, were generated from parental NRW-WE by repeated mass crosses and selection using the two isolates CpGV-M and CpGV-S, respectively. Both CpR5M and CpR5S were considered to be genetically homogeneous for the presence of the resistance allele(s). By crossing and backcrossing experiments with a susceptible CM strain, followed by resistance testing of the offspring, an autosomal dominant inheritance of resistance was elucidated. In addition, cross-resistance to CpGV-M and CpGV-S was detected in both strains, CpR5M and CpR5S. To test the hypothesis that the autosomal inheritance of type II resistance was caused by a large interchromosomal rearrangement involving the Z chromosome, making type I resistance appear to be autosomal in these strains; fluorescence in situ hybridization with bacterial artificial chromosome probes (BAC-FISH) was used to physically map the Z chromosomes of different CM strains. Conserved synteny of the Z-linked genes in CpR5M and other CM strains rejects this hypothesis and argues for a novel genetic and functional mode of resistance in CM populations with type II resistance.

Genomic characterization of two novel SAR11 isolates from the Red Sea, including the first strain of the SAR11 Ib clade.

The SAR11 clade (Pelagibacterales) is a diverse group that forms a monophyletic clade within the Alphaproteobacteria, and constitutes up to one third of all prokaryotic cells in the photic zone of most oceans. Pelagibacterales are very abundant in the warm and highly saline surface waters of the Red Sea, raising the question of adaptive traits of SAR11 populations in this water body and warmer oceans through the world. In this study, two pure cultures were successfully obtained from surface waters on the Red Sea: one isolate of subgroup Ia and one of the previously uncultured SAR11 Ib lineage. The novel genomes were very similar to each other and to genomes of isolates of SAR11 subgroup Ia (Ia pan-genome), both in terms of gene content and synteny. Among the genes that were not present in the Ia pan-genome, 108 (RS39, Ia) and 151 genes (RS40, Ib) were strain specific. Detailed analyses showed that only 51 (RS39, Ia) and 55 (RS40, Ib) of these strain-specific genes had not reported before on genome fragments of Pelagibacterales. Further analyses revealed the potential production of phosphonates by some SAR11 members and possible adaptations for oligotrophic life, including pentose sugar utilization and adhesion to marine particulate matter.

A Chromosome-Scale Assembly of the Bactrocera cucurbitae Genome Provides Insight to the Genetic Basis of white pupae.

Genetic sexing strains (GSS) used in sterile insect technique (SIT) programs are textbook examples of how classical Mendelian genetics can be directly implemented in the management of agricultural insect pests. Although the foundation of traditionally developed GSS are single locus, autosomal recessive traits, their genetic basis are largely unknown. With the advent of modern genomic techniques, the genetic basis of sexing traits in GSS can now be further investigated. This study is the first of its kind to integrate traditional genetic techniques with emerging genomics to characterize a GSS using the tephritid fruit fly pest Bactrocera cucurbitae as a model. These techniques include whole-genome sequencing, the development of a mapping population and linkage map, and quantitative trait analysis. The experiment designed to map the genetic sexing trait in B. cucurbitae, white pupae (wp), also enabled the generation of a chromosome-scale genome assembly by integrating the linkage map with the assembly. Quantitative trait loci analysis revealed SNP loci near position 42 MB on chromosome 3 to be tightly linked to wp. Gene annotation and synteny analysis show a near perfect relationship between chromosomes in B. cucurbitae and Muller elements A–E in Drosophila melanogaster. This chromosome-scale genome assembly is complete, has high contiguity, was generated using a minimal input DNA, and will be used to further characterize the genetic mechanisms underlying wp. Knowledge of the genetic basis of genetic sexing traits can be used to improve SIT in this species and expand it to other economically important Diptera.

Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species

BackgroundNucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.ResultsBased on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.ConclusionsThe structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.

A Novel Long Non-coding RNA, durga Modulates Dendrite Density and Expression of kalirin in Zebrafish.

Kalirin, a key player in axonal development, nerve growth and synaptic re-modeling, is implicated in many pathological conditions like schizophrenia and autism-spectrum disorders. Alternative promoters and splicing lead to functionally distinct isoforms, but the post-transcriptional regulation of Kalirin has not been studied. Here, we report a novel non-coding RNA, which we name durga, arising from the first exon of kalirin a (kalrna) in the antisense orientation in zebrafish. The kalrna and durga transcripts are barely detectable during early development, but steadily increase by 24 hours post-fertilization (hpf) as the brain develops. Over-expression of durga in the zebrafish embryo led to an increase in kalrna expression. The morphology of the neurons cultured from durga injected embryos had significantly fewer and shorter dendrites. Although durga has no apparent sequence homolog in mammals, based on gene synteny, we found a non-coding RNA arising from the 5′ end of the human Kalrn gene and expressed in the human neuronal cell line, SH-SY5Y. We propose that the zebrafish lncRNA durga maintains dendritic length and density through regulation of kalrna expression and this may have further implications in mammalian systems.

On the Molecular Evolution of Leptin, Leptin Receptor, and Endospanin.

Over a decade passed between Friedman’s discovery of the mammalian leptin gene (1) and its cloning in fish (2) and amphibians (3). Since 2005, the concept of gene synteny conservation (vs. gene sequence homology) was instrumental in identifying leptin genes in dozens of species, and we now have leptin genes from all major classes of vertebrates. This database of LEP (leptin), LEPR (leptin receptor), and LEPROT (endospanin) genes has allowed protein structure modeling, stoichiometry predictions, and even functional predictions of leptin function for most vertebrate classes. Here, we apply functional genomics to model hundreds of LEP, LEPR, and LEPROT proteins from both vertebrates and invertebrates. We identify conserved structural motifs in each of the three leptin signaling proteins and demonstrate Drosophila Dome protein’s conservation with vertebrate leptin receptors. We model endospanin structure for the first time and identify endospanin paralogs in invertebrate genomes. Finally, we argue that leptin is not an adipostat in fishes and discuss emerging knockout models in fishes.

Nine co-localized cytochrome P450 genes of the CYP2N, CYP2AD, and CYP2P gene families in the mangrove killifish Kryptolebias marmoratus genome: Identification and expression in response to B[α]P, BPA, OP, and NP

The CYP2 genes are the largest and most diverse cytochrome P450 (CYP) subfamily in vertebrates. We have identified nine co-localized CYP2 genes (∼55kb) in a new cluster in the genome of the highly resilient ecotoxicological fish model Kryptolebias marmoratus. Molecular characterization, temporal and tissue-specific expression pattern, and response to xenobiotics of these genes were examined. The CYP2 gene clusters were characterized and designated CYP2N22-23, CYP2AD12, and CYP2P16-20. Gene synteny analysis confirmed that the cluster in K. marmoratus is similar to that found in other teleost fishes, including zebrafish. A gene duplication event with diverged catalytic function was observed in CYP2AD12. Moreover, a high level of divergence in expression was observed among the co-localized genes. Phylogeny of the cluster suggested an orthologous relationship with similar genes in zebrafish and Japanese medaka. Gene expression analysis showed that CYP2P19 and CYP2N20 were consecutively expressed throughout embryonic development, whereas CYP2P18 was expressed in all adult tissues, suggesting that members of each CYP2 gene family have different physiological roles even though they are located in the same cluster. Among endocrine-disrupting chemicals (EDCs), benzo[α]pyrene (B[α]P) induced expression of CYP2N23, bisphenol A (BPA) induced CYP2P18 and CYP2P19, and 4-octylphenol (OP) induced CYP2AD12, but there was no significant response to 4-nonylphenol (NP), implying differential catalytic roles of the enzyme. In this paper, we identify and characterize a CYP2 gene cluster in the mangrove killifish K. marmoratus with differing catalytic roles toward EDCs. Our findings provide insights on the roles of nine co-localized CYP2 genes and their catalytic functions for better understanding of chemical-biological interactions in fish.

Complete mitochondrial genome of Gracilaria changii (Rhodophyta: Gracilariaceae)

Gracilaria changii is an agarophyte that has been used for agar extraction and as salads. Currently, there are only three complete mitochondrial genomes of Gracilaria species available in GenBank. We report here the complete mitogenome of G. changii from Peninsular Malaysia using next-generation sequencing. It had a total length of 25,729 bp, comprising 50 genes (25 protein-coding genes – PCGs, 2 rRNA genes, and 23 tRNA genes). The gene synteny of G. changii was similar to G. salicornia from Peninsular Malaysia with three additional tRNA genes (trnY, trnR, trnS1) in comparison to G. salicornia from Hawaii. G. changii had similar start codons (ATG) for all the PCGs, but three (rps3, rpl20, and secY) had incomplete T-stop codons. Genomic data of G. changii is useful to build a more complete dataset for phylogenetic analyses in Gracilariaceae.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota).

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Comparative mitochondrial and chloroplast genomics of a genetically distinct form of Sargassum contributing to recent "Golden Tides" in the Western Atlantic.

Over the past 5 years, massive accumulations of holopelagic species of the brown macroalga Sargassum in coastal areas of the Caribbean have created “golden tides” that threaten local biodiversity and trigger economic losses associated with beach deterioration and impact on fisheries and tourism. In 2015, the first report identifying the cause of these extreme events implicated a rare form of the holopelagic species Sargassum natans (form VIII). However, since the first mention of S. natans VIII in the 1930s, based solely on morphological characters, no molecular data have confirmed this identification. We generated full‐length mitogenomes and partial chloroplast genomes of all representative holopelagic Sargassum species, S. fluitans III and S. natans I alongside the putatively rare S. natans VIII, to demonstrate small but consistent differences between S. natans I and VIII (7 bp differences out of the 34,727). Our comparative analyses also revealed that both S. natans I and S. natans VIII share a very close phylogenetic relationship with S. fluitans III (94‐ and 96‐bp differences of 34,727). We designed novel primers that amplified regions of the cox2 and cox3 marker genes with consistent polymorphic sites that enabled differentiation between the two S. natans forms (I and VIII) from each other and both from S. fluitans III in over 150 Sargassum samples including those from the 2014 golden tide event. Despite remarkable gene synteny and sequence conservation, the three Sargassum forms differ in morphology, ecology, and distribution patterns, warranting more extensive interrogation of holopelagic Sargassum genomes as a whole.

Characterization and Functional Analysis of Calmodulin and Calmodulin-Like Genes in Fragaria vesca.

Calcium is a universal messenger that is involved in the modulation of diverse developmental and adaptive processes in response to various stimuli. Calmodulin (CaM) and calmodulin-like (CML) proteins are major calcium sensors in all eukaryotes, and they have been extensively investigated for many years in plants and animals. However, little is known about CaMs and CMLs in woodland strawberry (Fragaria vesca). In this study, we performed a genome-wide analysis of the strawberry genome and identified 4 CaM and 36 CML genes. Bioinformatics analyses, including gene structure, phylogenetic tree, synteny and three-dimensional model assessments, revealed the conservation and divergence of FvCaMs and FvCMLs, thus providing insight regarding their functions. In addition, the transcript abundance of four FvCaM genes and the four most related FvCML genes were examined in different tissues and in response to multiple stress and hormone treatments. Moreover, we investigated the subcellular localization of several FvCaMs and FvCMLs, revealing their potential interactions based on the localizations and potential functions. Furthermore, overexpression of five FvCaM and FvCML genes could not induce a hypersensitive response, but four of the five genes could increase resistance to Agrobacterium tumefaciens in Nicotiana benthamiana leaves. This study provides evidence for the biological roles of FvCaM and CML genes, and the results lay the foundation for future functional studies of these genes.

Characterisation of ZBTB46 and DC-SCRIPT/ZNF366 in rainbow trout, transcription factors potentially involved in dendritic cell maturation and activation in fish

ZBTB46 and DC-SCRIPT/ZNF366 are two zinc finger transcription factors that play important roles in regulating differentiation of dendritic cells in mammals. In this study, the ZBTB46 and DC-SCRIPT/ZNF366 homologues were identified in rainbow trout Oncorhynchus mykiss and their expression analysed in vivo and in vitro. As transcription factors, they are well conserved in sequence, genomic organisation and gene synteny. Their expression was differentially modulated by bacterial and viral PAMPs in the monocyte/macrophage-like cell line RTS-11, in primary head kidney (HK) macrophages, and in HK macrophages cultured with IL-4/13A. In the RTS-11 cells and primary HK macrophages, all the ZBTB46 and DC-SCRIPT/ZNF366 homologues were down-regulated by interferon gamma (type II IFN) but unaffected by IFN2 (type I IFN), administered as recombinant proteins to cell cultures. In fish gills, infection with amoebae (Paramoebae perurans) resulted in reduction of ZBTB46 and DC-SCRIPT/ZNF366 expression in Atlantic salmon Salmo salar, whilst infection with Yersinia ruckeri induced gene expression in rainbow trout.

Analysis of the complete plastomes of three species of Membranoptera (Ceramiales, Rhodophyta) from Pacific North America.

Next generation sequence data were generated and used to assemble the complete plastomes of the holotype of Membranoptera weeksiae, the neotype (designated here) of M.tenuis, and a specimen examined by Kylin in making the new combination M.platyphylla. The three plastomes were similar in gene content and length and showed high gene synteny to Calliarthron, Grateloupia, Sporolithon, and Vertebrata. Sequence variation in the plastome coding regions were 0.89% between M.weeksiae and M.tenuis, 5.14% between M.weeksiae and M.platyphylla, and 5.18% between M.tenuis and M.platyphylla. We were unable to decipher the complete mitogenomes of the three species due to low coverage and structural problems; however, we assembled and analyzed, the cytochrome oxidase I, II, and III loci and found that M.weeksiae and M.tenuis differed in sequence by 1.3%, M.weeksiae and M.platyphylla by 8.4%, and M.tenuis and M.platyphylla by 8.1%. Evaluation of standard marker genes indicated that sequences from the rbcL, RuBisCO spacer, and CO1 genes closely approximated the pair-wise genetic distances observed between the plastomes of the three species of Membranoptera. A phylogenetic tree based on rbcL sequences showed that M.tenuis and M.weeksiae were sister taxa. Short rbcL sequences were obtained from type specimens of M.dimorpha, M.multiramosa, and M.edentata and confirmed their conspecificity with M.platyphylla. The data support the recognition of three species of Membranoptera occurring south of Alaska: M.platyphylla, M.tenuis, and M.weeksiae.

Blastocystis Mitochondrial Genomes Appear to Show Multiple Independent Gains and Losses of Start and Stop Codons.

Complete mitochondrion-related organelle (MRO) genomes of several subtypes (STs) of the unicellular stramenopile Blastocystis are presented. Complete conservation of gene content and synteny in gene order is observed across all MRO genomes, comprising 27 protein coding genes, 2 ribosomal RNA genes, and 16 transfer RNA (tRNA) genes. Despite the synteny, differences in the degree of overlap between genes were observed between subtypes and also between isolates within the same subtype. Other notable features include unusual base-pairing mismatches in the predicted secondary structures of some tRNAs. Intriguingly, the rps4 gene in some MRO genomes is missing a start codon and, based on phylogenetic relationships among STs, this loss has happened twice independently. One unidentified open reading frame (orf160) is present in all MRO genomes. However, with the exception of ST4 where the feature has been lost secondarily, orf160 contains variously one or two in-frame stop codons. The overall evidence suggests that both the orf160 and rps4 genes are functional in all STs, but how they are expressed remains unclear.

The landscape of microbial phenotypic traits and associated genes.

Bacteria and Archaea display a variety of phenotypic traits and can adapt to diverse ecological niches. However, systematic annotation of prokaryotic phenotypes is lacking. We have therefore developed ProTraits, a resource containing ∼545 000 novel phenotype inferences, spanning 424 traits assigned to 3046 bacterial and archaeal species. These annotations were assigned by a computational pipeline that associates microbes with phenotypes by text-mining the scientific literature and the broader World Wide Web, while also being able to define novel concepts from unstructured text. Moreover, the ProTraits pipeline assigns phenotypes by drawing extensively on comparative genomics, capturing patterns in gene repertoires, codon usage biases, proteome composition and co-occurrence in metagenomes. Notably, we find that gene synteny is highly predictive of many phenotypes, and highlight examples of gene neighborhoods associated with spore-forming ability. A global analysis of trait interrelatedness outlined clusters in the microbial phenotype network, suggesting common genetic underpinnings. Our extended set of phenotype annotations allows detection of 57 088 high confidence gene-trait links, which recover many known associations involving sporulation, flagella, catalase activity, aerobicity, photosynthesis and other traits. Over 99% of the commonly occurring gene families are involved in genetic interactions conditional on at least one phenotype, suggesting that epistasis has a major role in shaping microbial gene content.

CTDGFinder: A Novel Homology-Based Algorithm for Identifying Closely Spaced Clusters of Tandemly Duplicated Genes.

Closely spaced clusters of tandemly duplicated genes (CTDGs) contribute to the diversity of many phenotypes, including chemosensation, snake venom, and animal body plans. CTDGs have traditionally been identified subjectively as genomic neighborhoods containing several gene duplicates in close proximity; however, CTDGs are often highly variable with respect to gene number, intergenic distance, and synteny. This lack of formal definition hampers the study of CTDG evolutionary dynamics and the discovery of novel CTDGs in the exponentially growing body of genomic data. To address this gap, we developed a novel homology-based algorithm, CTDGFinder, which formalizes and automates the identification of CTDGs by examining the physical distribution of individual members of families of duplicated genes across chromosomes. Application of CTDGFinder accurately identified CTDGs for many well-known gene clusters (e.g., Hox and beta-globin gene clusters) in the human, mouse and 20 other mammalian genomes. Differences between previously annotated gene clusters and our inferred CTDGs were due to the exclusion of nonhomologs that have historically been considered parts of specific gene clusters, the inclusion or absence of genes between the CTDGs and their corresponding gene clusters, and the splitting of certain gene clusters into distinct CTDGs. Examination of human genes showing tissue-specific enhancement of their expression by CTDGFinder identified members of several well-known gene clusters (e.g., cytochrome P450s and olfactory receptors) and revealed that they were unequally distributed across tissues. By formalizing and automating CTDG identification, CTDGFinder will facilitate understanding of CTDG evolutionary dynamics, their functional implications, and how they are associated with phenotypic diversity.

Mitogenomes from type specimens, a genotyping tool for morphologically simple species: ten genomes of agar-producing red algae.

DNA sequences from type specimens provide independent, objective characters that enhance the value of type specimens and permit the correct application of species names to phylogenetic clades and specimens. We provide mitochondrial genomes (mitogenomes) from archival type specimens of ten species in agar-producing red algal genera Gelidium and Pterocladiella. The genomes contain 43–44 genes, ranging in size from 24,910 to 24,970 bp with highly conserved gene synteny. Low Ka/Ks ratios of apocytochrome b and cytochrome oxidase genes support their utility as markers. Phylogenies of mitogenomes and cox1+rbcL sequences clarified classification at the genus and species levels. Three species formerly in Gelidium and Pterocladia are transferred to Pterocladiella: P. media comb. nov., P. musciformis comb. nov., and P. luxurians comb. and stat. nov. Gelidium sinicola is merged with G. coulteri because they share identical cox1 and rbcL sequences. We describe a new species, Gelidium millariana sp. nov., previously identified as G. isabelae from Australia. We demonstrate that mitogenomes from type specimens provide a new tool for typifying species in the Gelidiales and that there is an urgent need for analyzing mitogenomes from type specimens of red algae and other morphologically simple organisms for insight into their nomenclature, taxonomy and evolution.