Long‐read sequence capture of the haemoglobin gene clusters across codfish species

Combining high‐throughput sequencing with targeted sequence capture has become an attractive tool to study specific genomic regions of interest. Most studies have so far focused on the exome using short‐read technology. These approaches are not designed to capture intergenic regions needed to reconstruct genomic organization, including regulatory regions and gene synteny. Here, we demonstrate the power of combining targeted sequence capture with long‐read sequencing technology for comparative genomic analyses of the haemoglobin (Hb) gene clusters across eight species separated by up to 70 million years. Guided by the reference genome assembly of the Atlantic cod (Gadus morhua) together with genome information from draft assemblies of selected codfishes, we designed probes covering the two Hb gene clusters. Use of custom‐made barcodes combined with PacBio RSII sequencing led to highly continuous assemblies of the LA (~100 kb) and MN (~200 kb) clusters, which include syntenic regions of coding and intergenic sequences. Our results revealed an overall conserved genomic organization of the Hb genes within this lineage, yet with several, lineage‐specific gene duplications. Moreover, for some of the species examined, we identified amino acid substitutions at two sites in the Hbb1 gene as well as length polymorphisms in its regulatory region, which has previously been linked to temperature adaptation in Atlantic cod populations. This study highlights the use of targeted long‐read capture as a versatile approach for comparative genomic studies by generation of a cross‐species genomic resource elucidating the evolutionary history of the Hb gene family across the highly divergent group of codfishes.     

The ergot alkaloid gene cluster: Functional analyses and evolutionary aspects

Ergot alkaloids and their derivatives have been traditionally used as therapeutic agents in migraine, blood pressure regulation and help in childbirth and abortion. Their production in submerse culture is a long established biotechnological process. Ergot alkaloids are produced mainly by members of the genus Claviceps, with Claviceps purpurea as best investigated species concerning the biochemistry of ergot alkaloid synthesis (EAS). Genes encoding enzymes involved in EAS have been shown to be clustered; functional analyses of EAS cluster genes have allowed to assign specific functions to several gene products. Various Claviceps species differ with respect to their host specificity and their alkaloid content; comparison of the ergot alkaloid clusters in these species (and of clavine alkaloid clusters in other genera) yields interesting insights into the evolution of cluster structure. This review focuses on recently published and also yet unpublished data on the structure and evolution of the EAS gene cluster and on the function and regulation of cluster genes. These analyses have also significant biotechnological implications: the characterization of non-ribosomal peptide synthetases (NRPS) involved in the synthesis of the peptide moiety of ergopeptines opened interesting perspectives for the synthesis of ergot alkaloids; on the other hand, defined mutants could be generated producing interesting intermediates or only single peptide alkaloids (instead of the alkaloid mixtures usually produced by industrial strains).     

Taxol from fungal endophytes and the issue of biodiversity

Fungi represent one of the most understudied and diverse group of organisms. Commonly, these organisms make associations with higher life forms and may proceed to biochemically mimic the host organism. An excellent example of this is the anticancer drug, taxol, which had been previously supposed to occur only in the plant genusTaxus (yew). However, taxol has been reported in a novel endophytic fungus—Taxomyces andreanae, but also has been demonstrated to occur in a number of unrelated fungal endophytes includingPestalotia, Pestalotiopsis, Fusarium, Alternaria, Pithomyces, Monochaetia and others. Thus, this report presents information on the presence of taxol among disparate fungal genera, and uses these observations as an additional argument to support efforts to study fungal endophytes and preserve their associated host plants.     

Recent evolutionary history of the bluethroat (Luscinia svecica) across Eurasia

We analysed mitochondrial DNA (mtDNA) sequences from 154 bluethroats (Luscinia svecica) sampled at 21 sites throughout much of their Eurasian range. A previously reported, single base-pair mtDNA difference between L. s. svecica and L. s. namnetum was inconsistent upon expanded geographical sampling. A significant FST value (0.29) and an isolation-by-distance effect show the existence of geographical differentiation. Phylogenetic analysis of haplotypes revealed northern and southern groups, although lineage sorting is incomplete. There was no geographical structure to the haplotype tree within groups, and currently recognized subspecies were not supported. A minimum evolution tree based on pairwise mtDNA genetic distances among average samples showed the same two broadly distributed northern and southern groups. These groups abut in the centre of the latitudinal range, and were possibly isolated by forest that developed and spread westward over the last 15 000 years. Pairwise FST values averaged 0.16 in the southern group, 0.04 in the northern group, and 0.42 between groups. Mismatch distributions suggested population growth in each group, with that in the south being more recent. In the northern group, the geographical pattern in tau suggested northward and eastward expansion. Analysis of nucleotide diversity suggested westward expansion in the southern group. The northern group had higher nucleotide diversity than the southern group, consistent with a larger current population size in the north. Given the significant FST, incompletely sorted haplotype tree, and broadly patterned minimum evolution tree, L. svevica appears to represent a species at an intermediate stage of differentiation between panmixia and reciprocal monophyly.     

Support Vector Machine-based method for predicting subcellular localization of mycobacterial proteins using evolutionary information and motifs

Background

Annotations that describe the function of sequences are enormously important to researchers during laboratory investigations and when making computational inferences. However, there has been little investigation into the data quality of sequence function annotations. Here we have developed a new method of estimating the error rate of curated sequence annotations, and applied this to the Gene Ontology (GO) sequence database (GOSeqLite). This method involved artificially adding errors to sequence annotations at known rates, and used regression to model the impact on the precision of annotations based on BLAST matched sequences.

Results

We estimated the error rate of curated GO sequence annotations in the GOSeqLite database (March 2006) at between 28% and 30%. Annotations made without use of sequence similarity based methods (non-ISS) had an estimated error rate of between 13% and 18%. Annotations made with the use of sequence similarity methodology (ISS) had an estimated error rate of 49%.

Conclusion

While the overall error rate is reasonably low, it would be prudent to treat all ISS annotations with caution. Electronic annotators that use ISS annotations as the basis of predictions are likely to have higher false prediction rates, and for this reason designers of these systems should consider avoiding ISS annotations where possible. Electronic annotators that use ISS annotations to make predictions should be viewed sceptically. We recommend that curators thoroughly review ISS annotations before accepting them as valid. Overall, users of curated sequence annotations from the GO database should feel assured that they are using a comparatively high quality source of information.     

Diversity and evolutionary history of the MHC DQA gene in leporids

The European rabbit (Oryctolagus cuniculus) is used as a model for many human diseases, yet comparatively little is known of its genetics, particularly at important loci such as the major histocompatibility complex (MHC). This study investigated genetic diversity and evolutionary history of the DQA gene in a range of leporid species by analysing coding sequence diversity of exon 2 and intron 2 in 53 individuals of 16 different species. Fifty leporid DQA alleles were detected, including 13 novel European rabbit alleles. In the rabbit, the highest levels of diversity were observed in wild rabbits from Portugal, with wild rabbits from England and domestic rabbits showing less diversity. Within the sample, several recombination events were detected and trans-specific evolution of alleles was evidenced, both being general characteristics of mammalian MHC genes. Positive selection is implicated as operating on six codons within exon 2, which are also subject to positive selection in other mammals. Some of these positions are putative antigen recognition sites and underline the importance of pathogen-driven selection on these MHC genes.     

Inferring the evolutionary history of vibrios by means of multilocus sequence analysis

We performed the first broad study aiming at the reconstruction of the evolutionary history of vibrios by means of multilocus sequence analysis of nine genes. Overall, 14 distinct clades were recognized using the SplitsTree decomposition method. Some of these clades may correspond to families, e.g., the clades Salinivibrio and Photobacteria, while other clades, e.g., Splendidus and Harveyi, correspond to genera. The common ancestor of all vibrios was estimated to have been present 600 million years ago. We can define species of vibrios as groups of strains that share >95% gene sequence similarity and >99.4% amino acid identity based on the eight protein-coding housekeeping genes. The gene sequence data were used to refine the standard online electronic taxonomic scheme for vibrios (http://www.taxvibrio.lncc.br).     

Phosphoglucose isomerase gene duplication in the bony fishes: An evolutionary history

Electrophoretic patterns of phosphoglucose isomerase (PGI) in bony fishes provide strong evidence for a model of genetic control by two independent structural gene loci, most likely resulting from a gene duplication. This model is confirmed by a comparison of certain kinetic and molecular properties of the PGI homodimers (PGI-1 and PGI-2) isolated from extracts of the teleost Astyanax mexicanus. In addition, in most higher teleosts examined, the PGI enzymes show a regular pattern of tissue distribution, with PGI-2 predominant in muscle, the heterodimer often strongest in the heart, and PGI-1 predominant in liver and other organs. An examination of 53 species of bony fishes belonging to 38 families indicates a widespread occurrence of duplicate PGI loci and an early origin of the gene duplication, perhaps in the Leptolepiformes. The apparent presence of three PGI loci in trout and goldfish exemplifies how new loci can be incorporated into the genome through polyploidization.This research was supported in part by a NSF graduate traineeship to J.C.A., by the Clayton Foundation for Research in Biochemistry (G.B.K.), by NSF Grant GB-15644 and NIH Grant GM-15769 to Robert K. Selander, and by contract AT(38-1)-310 between the University of Georgia and the U.S. Atomic Energy Commission.     

An evolutionary history of the selectin gene cluster in humans

Molecules involved in leukocyte trafficking have a central role in the development of inflammatory and immune responses. We performed F(ST) analysis of the selectin cluster, as well as of SELPLG, ICAM1 and VCAM1. Peaks of significantly high population genetic differentiation were restricted to two regions in SELP and one in SELPLG. Resequencing data indicated that the region covering SELP exons 11-13 displays high nucleotide diversity in Africans and Europeans (CEU), and a high level of within-species diversity compared with inter-specific divergence. Analysis of inferred haplotypes revealed a complex phylogeny with two deeply separated clades that coalesce at ~3.5 million years (MY) plus a minor clade with a TMRCA (time to the most recent common ancestor) of ~2.2 MY. A splicing assay indicated no haplotype-specific effect on SELP exon 14 inclusion. These data are consistent with a model of multiallelic balancing selection; single-nucleotide polymorphism analysis indicated that the Val640Leu variant represents a likely selection target. In populations of Asian ancestry a distinct haplotype, possibly carrying regulatory variants, has been driven to high frequency by positive selection. No deviation from neutrality was observed for the SELPLG region. Resequencing of SELP in chimpanzees revealed a haplotype phylogeny with extremely deep basal branches, suggesting either long-standing balancing selection or ancestral population structure. Thus, SELP has experienced a complex selective history, possibly as a result of local adaptation. Variants in the gene have been associated with autoimmune and cardiovascular diseases. Association studies would benefit from both taking the complex SELP haplotype structure into account and from analysis of possible regulatory variants in the gene.     

Definition of the tempo of sequence diversity across an alignment and automatic identification of sequence motifs: Application to protein homologous families and superfamilies

It is often possible to identify sequence motifs that characterize a protein family in terms of its fold and/or function from aligned protein sequences. Such motifs can be used to search for new family members. Partitioning of sequence alignments into regions of similar amino acid variability is usually done by hand. Here, I present a completely automatic method for this purpose: one that is guaranteed to produce globally optimal solutions at all levels of partition granularity. The method is used to compare the tempo of sequence diversity across reliable three-dimensional (3D) structure-based alignments of 209 protein families (HOMSTRAD) and that for 69 superfamilies (CAMPASS). (The mean alignment length for HOMSTRAD and CAMPASS are very similar.) Surprisingly, the optimal segmentation distributions for the closely related proteins and distantly related ones are found to be very similar. Also, optimal segmentation identifies an unusual protein superfamily. Finally, protein 3D structure clues from the tempo of sequence diversity across alignments are examined. The method is general, and could be applied to any area of comparative biological sequence and 3D structure analysis where the constraint of the inherent linear organization of the data imposes an ordering on the set of objects to be clustered.     

The alpha-amylase gene in Drosophila melanogaster: nucleotide sequence, gene structure and expression motifs.

We present the complete nucleotide sequence of a Drosophila alpha-amylase gene and its flanking regions, as determined by cDNA and genomic sequence analysis. This gene, unlike its mammalian counterparts, contains no introns. Nevertheless the insect and mammalian genes share extensive nucleotide similarity and the insect protein contains the four amino acid sequence blocks common to all alpha-amylases. In Drosophila melanogaster, there are two closely-linked copies of the alpha-amylase gene and they are divergently transcribed. In the 5'-regions of the two gene-copies we find high sequence divergence, yet the typical eukaryotic gene expression motifs have been maintained. The 5'-terminus of the alpha-amylase mRNA, as determined by primer extension analysis, maps to a characteristic Drosophila sequence motif. Additional conserved elements upstream of both genes may also be involved in amylase gene expression which is known to be under complex controls that include glucose repression.     

The evolutionary history of human and chimpanzee Y-chromosome gene loss

Recent studies have suggested that gene gain and loss may contribute significantly to the divergence between humans and chimpanzees. Initial comparisons of the human and chimpanzee Y-chromosomes indicate that chimpanzees have a disproportionate loss of Y-chromosome genes, which may have implications for the adaptive evolution of sex-specific as well as reproductive traits, especially because one of the genes lost in chimpanzees is critically involved in spermatogenesis in humans. Here we have characterized Y-chromosome sequences in gorilla, bonobo, and several chimpanzee subspecies for 7 chimpanzee gene-disruptive mutations. Our analyses show that 6 of these gene-disruptive mutations predate chimpanzee-bonobo divergence at approximately 1.8 MYA, which indicates significant Y-chromosome change in the chimpanzee lineage relatively early in the evolutionary divergence of humans and chimpanzees.     

Comparison of lantibiotic gene clusters and encoded proteins

Lantibiotics form a group of modified peptides with unique structures, containing post-translationally modified amino acids such as dehydrated and lanthionine residues. In the gram-positive bacteria that secrete these lantibiotics, the gene clusters flanking the structural genes for various linear (type A) lantibiotics have recently been characterized. The best studied representatives are those of nisin (nis), subtilin (spa), epidermin (epi), Pep5 (pep), cytolysin (cyl), lactocin S (las) and lacticin 481 (lct). Comparison of the lantibiotic gene clusters shows that they contain conserved genes that probably encode similar functions.The nis, spa, epi and pep clusters contain lanB and lanC genes that are presumed to code for two types of enzymes that have been implicated in the modification reactions characteristic of all lantibiotics, i.e. dehydration and thio-ether ring formation. The cyl, las and lct gene clusters have no homologue of the lanB gene, but they do contain a much larger lanM gene that is the lanC gene homologue. Most lantibiotic gene clusters contain a lanP gene encoding a serine protease that is presumably involved in the proteolytic processing of the prelantibiotics. All clusters contain a lanT gene encoding and ABC transporter likely to be involved in the export of (precursors of) the lantibiotics. The lanE, lanF and lanG genes in the nis, spa and epi clusters encode another transport system that is possibly involved in self-protection. In the nisin and subtilin gene clusters two tandem genes, lanR and lanK, have been located that code for a two-component regulatory system.Finally, non-homologous genes are found in some lantibiotic gene clusters. The nisl and spal genes encode lipoproteins that are involved in immunity, the pepI gene encodes a membrane-located immunity protein, and epiD encodes an enzyme involved in a post-translational modification found only in the C-terminus of epidermin. Several genes of unknown function are also found in the las gene cluster.A database has been assembled for all putative gene products of type A lantibiotic gene clusters. Database searches, multiple sequence alignment and secondary structure prediction have been used to identify conserved sequence segments in the LanB, LanC, LanE, LanF, LanG, LanK, LanM, LanP, LanR and LanT gene products that may be essential for structure and function. This database allows for a rapid screening of newly determined sequences in lantibiotic gene clusters.     

The floral genome: an evolutionary history of gene duplication and shifting patterns of gene expression

Through multifaceted genome-scale research involving phylogenomics, targeted gene surveys, and gene expression analyses in diverse basal lineages of angiosperms, our studies provide insights into the most recent common ancestor of all extant flowering plants. MADS-box gene duplications have played an important role in the origin and diversification of angiosperms. Furthermore, early angiosperms possessed a diverse tool kit of floral genes and exhibited developmental 'flexibility', with broader patterns of expression of key floral organ identity genes than are found in eudicots. In particular, homologs of B-function MADS-box genes are more broadly expressed across the floral meristem in basal lineages. These results prompted formulation of the 'fading borders' model, which states that the gradual transitions in floral organ morphology observed in some basal angiosperms (e.g. Amborella) result from a gradient in the level of expression of floral organ identity genes across the developing floral meristem.