PipMaker--a web server for aligning two genomic DNA sequences

B52 promotes alternative splicing of Dscam in Chinese mitten crab, Eriocheir sinensis

B52 is a member of the classical serine/arginine (SR)-rich proteins, which are phylogenetically conserved and play significant roles in mRNA maturation, including alternative splicing. In the present study, the docking site, selector sequences and locus control region of the Chinese mitten crab (Eriocheir sinensis) Down syndrome cell adhesion molecule (EsDscam) were identified. Alternative splicing of Dscam is essential to generate different isoforms. We also isolated and characterised the B52 gene from E. sinensis (EsB52). The 876 bp open reading frame of EsB52 encodes a 291 amino acid residue polypeptide, and EsB52 has two RNA recognition motifs (RRMs) at the N-terminus and an arginine/serine-rich domain at the C-terminus. Each RRM contains two degenerate short submotifs, RNP-1 and RNP2. Analysis of tissue distribution revealed that EsB52 mRNA expression was widespread in all tested tissues, and especially high in brain and hemocytes. In hemocytes, EsB52 was upregulated significantly after stimulation with pathogen-associated molecular patterns and bacteria. Furthermore, EsB52 RNAi decreased the number of Ig7 inclusion in mRNA rather than Ig2 or Ig3. Taken together, these findings suggest that EsB52 acts as an alternative splicing activator of EsDscam.

Multiple Conserved Elements Structuring Inverted Repeats in the Mammalian Coat Color-Related Gene Asip

In the agouti signaling gene protein (Asip) of the house mouse (Mus musculus), inverted repeat (IR) arrays are known to exist in a non-coding region adjacent to the ventral-specific promoter region and the accompanying two exons (exons 1A and 1A'), which are around 100 kb upstream from the amino acid coding regions of exons 2, 3, and 4. To determine the gene structure of mammalian Asip and to elucidate trends in its evolution, non-coding sequences of six rodent (mouse, rat, Chinese hamster, squirrel, guinea pig, and naked mole rat) and three non-rodent (rabbit, human, and cow) species were retrieved from databases and compared. Our homology search analyses revealed the presence of three to five highly conserved non-coding elements (CNE). These CNEs were found to form IRs in rodents and lagomorphs. Combinations of IRs were further shown to build symmetric, long IR arrays. Intra- and inter-specific comparisons of the sequences of three universal CNEs showed homogeneity between CNE pairs within species. This implies that certain evolutionary constraints maintained the IR structure in the rodent and rabbit species.

Arid3a regulates nephric tubule regeneration via evolutionarily conserved regeneration signal-response enhancers

Amphibians and fish have the ability to regenerate numerous tissues, whereas mammals have a limited regenerative capacity. Despite numerous developmental genes becoming reactivated during regeneration, an extensive analysis is yet to be performed on whether highly regenerative animals utilize unique cis-regulatory elements for the reactivation of genes during regeneration and how such cis-regulatory elements become activated. Here, we screened regeneration signal-response enhancers at the lhx1 locus using Xenopus and found that the noncoding elements conserved from fish to human function as enhancers in the regenerating nephric tubules. A DNA-binding motif of Arid3a, a component of H3K9me3 demethylases, was commonly found in RSREs. Arid3a binds to RSREs and reduces the H3K9me3 levels. It promotes cell cycle progression and causes the outgrowth of nephric tubules, whereas the conditional knockdown of arid3a using photo-morpholino inhibits regeneration. These results suggest that Arid3a contributes to the regeneration of nephric tubules by decreasing H3K9me3 on RSREs.

Gene fusion of heterophyletic gamma-globin genes in platyrrhine primates

We performed phylogenetic analyses of HBG genes to assess its origin and interspecific variation among primates. Our analyses showed variation in HBG genes copy number ranging from one to three, some of them pseudogenes. For platyrrhines HBG genes, phylogenetic reconstructions of flanking regions recovered orthologous clades with distinct topologies for 5' and 3' flanking regions. The 5' region originated in the common ancestor of platyrrhines but the 3' region had an anthropoid origin. We hypothesize that the platyrrhine HBG genes of 5' and 3' heterophyletic origins arose from subsequent fusions of the (earlier) platyrrhine 5' portion and the (later) anthropoid 3' portion.

Mapping of a responsible region for sex reversal upstream of Sox9 by production of mice with serial deletion in a genomic locus

Sox9 plays critical roles in testis formation. By mapping four familial cases of disorders of sexual development, a 32.5 kb sequence located far upstream of SOX9 was previously identified as being a commonly deleted region and named the XY sex reversal region (XYSR). To narrow down a responsible sequence in XYSR, we generated mutant mice with a series of deletions in XYSR by application of the CRISPR/Cas9 system, using a mixture of sgRNAs targeting several kilobase (kb) intervals in the region. When the whole XYSR corresponding sequence in mice was deleted in XY karyotype individuals, the mutation resulted in female offspring, suggesting that an expression mechanism of SOX9/Sox9 through XYSR is conserved in human and mouse. Male-to-female sex reversal was found in mice with a 4.8 kb deletion. We identified a sequence conserved among humans, mice, and opossum, the deletion of which (783 bp) in mice resulted in male-to-female sex reversal. The sequence includes a recently reported critical gonad enhancer for Sox9. Although it cannot be concluded that the human sequence is responsible for XYSR, it is likely. This method is applicable for fine mapping of responsible sequences for disease-causing deletions especially with regard to rare diseases.

Overview of the Germline and Expressed Repertoires of the TRB Genes in Sus scrofa

The α/β T cell receptor (TR) is a complex heterodimer that recognizes antigenic peptides and binds to major histocompatibility complex (MH) molecules. Both α and β chains are encoded by different genes localized on two distinct chromosomal loci: TRA and TRB. The present study employed the recent release of the swine genome assembly to define the genomic organization of the TRB locus. According to the sequencing data, the pig TRB locus spans approximately 400 kb of genomic DNA and consists of 38 TRBV genes belonging to 24 subgroups located upstream of three in tandem TRBD-J-C clusters, which are followed by a TRBV gene in an inverted transcriptional orientation. Comparative analysis confirms that the general organization of the TRB locus is similar among mammalian species, but the number of germline TRBV genes varies greatly even between species belonging to the same order, determining the diversity and specificity of the immune response. However, sequence analysis of the TRB locus also suggests the presence of blocks of conserved homology in the genomic region across mammals. Furthermore, by analysing a public cDNA collection, we identified the usage pattern of the TRBV, TRBD, and TRBJ genes in the adult pig TRB repertoire, and we noted that the expressed TRBV repertoire seems to be broader and more diverse than the germline repertoire, in line with the presence of a high level of TRBV gene polymorphisms. Because the nucleotide differences seems to be principally concentrated in the CDR2 region, it is reasonable to presume that most T cell β-chain diversity can be related to polymorphisms in pig MH molecules. Domestic pigs represent a valuable animal model as they are even more anatomically, genetically and physiologically similar to humans than are mice. Therefore, present knowledge on the genomic organization of the pig TRB locus allows the collection of increased information on the basic aspects of the porcine immune system and contributes to filling the gaps left by rodent models.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

Evolution and Sequence Diversity of FhuA in Salmonella and Escherichia

The fhuACDB operon, present in a number of Enterobacteriaceae, encodes components essential for the uptake of ferric hydroxamate type siderophores. FhuA acts not only as a transporter for physiologically important chelated ferric iron but also as a receptor for various bacteriophages, toxins, and antibiotics, which are pathogenic to bacterial cells. In this research, fhuA gene distribution and sequence diversity were investigated in Enterobacteriaceae, especially Salmonella and Escherichia Comparative sequence analysis resulted in a fhuA phylogenetic tree that did not match the expected phylogeny of species or trees of the fhuCDB genes. The fhuA sequences showed a unique mosaic clustering pattern. On the other hand, the gene sequences showed high conservation for strains from the same serovar or serotype. In total, six clusters were identified from FhuA proteins in Salmonella and Escherichia, among which typical peptide fragment variations could be defined. Six fragmental insertions/deletions and two substitution fragments were discovered, for which the combination of polymorphism patterns could well classify the different clusters. Structural modeling demonstrated that all the six featured insertions/deletions and one substitution fragment are located at the apexes of the long loops present as part of the FhuA external pocket. These frequently mutated regions are likely under high selection pressure, with bacterial strains balancing escape from phage infection or toxin/antibiotics attack via fhuA gene mutations while maintaining the siderophore uptake activity essential for bacterial survival. The unusual fhuA clustering suggests that high-frequency exchange of fhuA genes has occurred between enterobacterial strains after distinctive species were established.

Diversification of the functions of proglucagon and glucagon receptor genes in fish

The teleost fish-specific genome duplication gave rise to a great number of species inhabiting diverse environments with different access to nutrients and life histories. This event produced duplicated gcg genes, gcga and gcgb, for proglucagon-derived peptides, glucagon and GLP-1 and duplicated gcgr receptor genes, gcgra and gcgrb, which play key roles connecting the consumption of nutrients with glucose metabolism. We conducted a systematic survey of the genomes from 28 species of fish (24 bony (Superclass Osteichthyes), 1 lobe-finned (Class Sarcoperygii), 1 cartilaginous (Superclass Chondrichthyes), and 2 jawless (Superclass Agnatha)) and find that almost all surveyed ray-finned fish contain gcga and gcgb genes with different coding potential and duplicated gcgr genes, gcgra and gcgrb that form two separate clades in the phylogenetic tree consistent with the accepted species phylogeny. All gcgb genes encoded only glucagon and GLP-1 and gcga genes encoded glucagon, GLP-1, and GLP-2, indicating that gcga was subfunctionalized to produce GLP-2. We find a single glp2r, but no glp1r suggesting that duplicated gcgrb was neofunctionalized to bind GLP-1, as demonstrated for the zebrafish gcgrb (Oren et al., 2016). In functional experiments with zebrafish gcgrb and GLP-1 from diverse fish we find that anglerfish GLP-1a, encoded by gcga, is less biologically active than the gcgb anglerfish GLP-1b paralog. But some other fish (zebrafish, salmon, and catfish) gcga GLP-1a display similar biological activities, indicating that the regulation of glucose metabolism by GLP-1 in ray-finned fish is species-specific. Searches of genomes in cartilaginous fish identified a proglucagon gene that encodes a novel GLP-3 peptide in addition to glucagon, GLP-1, and GLP-2, as well as a single gcgr, glp2r, and a new glucagon receptor-like receptor whose identity still needs to be confirmed. The sequence of the shark GLP-1 contained an N-terminal mammalian-like extension that in mammals undergoes a proteolytic cleavage to release biologically active GLP-1. Our results indicate that early in vertebrate evolution diverse regulatory mechanisms emerged for the control of glucose metabolism by proglucagon-derived peptides and their receptors and that in ray-finned fish they included subfunctionalization and neofunctionalization of these genes.

A genome-wide assessment of conserved SNP alleles reveals a panel of regulatory SNPs relevant to the peripheral nerve

Background

Identifying functional non-coding variation is critical for defining the genetic contributions to human disease. While single-nucleotide polymorphisms (SNPs) within cis-acting transcriptional regulatory elements have been implicated in disease pathogenesis, not all cell types have been assessed and functional validations have been limited. In particular, the cells of the peripheral nervous system have been excluded from genome-wide efforts to link non-coding SNPs to altered gene function. Addressing this gap is essential for defining the genetic architecture of diseases that affect the peripheral nerve. We developed a computational pipeline to identify SNPs that affect regulatory function (rSNPs) and evaluated our predictions on a set of 144 regions in Schwann cells, motor neurons, and muscle cells.

Results

We identified 28 regions that display regulatory activity in at least one cell type and 13 SNPs that affect regulatory function. We then tailored our pipeline to one peripheral nerve cell type by incorporating SOX10 ChIP-Seq data; SOX10 is essential for Schwann cells. We prioritized 22 putative SOX10 response elements harboring a SNP and rapidly validated two rSNPs. We then selected one of these elements for further characterization to assess the biological relevance of our approach. Deletion of the element from the genome of cultured Schwann cells—followed by differential gene expression studies—revealed Tubb2b as a candidate target gene. Studying the enhancer in developing mouse embryos revealed activity in SOX10-positive cells including the dorsal root ganglia and melanoblasts.

Conclusions

Our efforts provide insight into the utility of employing strict conservation for rSNP discovery. This strategy, combined with functional analyses, can yield candidate target genes. In support of this, our efforts suggest that investigating the role of Tubb2b in SOX10-positive cells may reveal novel biology within these cell populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4692-z) contains supplementary material, which is available to authorized users.

Egg Case Silk Gene Sequences from Argiope Spiders: Evidence for Multiple Loci and a Loss of Function Between Paralogs

Spiders swath their eggs with silk to protect developing embryos and hatchlings. Egg case silks, like other fibrous spider silks, are primarily composed of proteins called spidroins (spidroin = spider-fibroin). Silks, and thus spidroins, are important throughout the lives of spiders, yet the evolution of spidroin genes has been relatively understudied. Spidroin genes are notoriously difficult to sequence because they are typically very long (≥ 10 kb of coding sequence) and highly repetitive. Here, we investigate the evolution of spider silk genes through long-read sequencing of Bacterial Artificial Chromosome (BAC) clones. We demonstrate that the silver garden spider Argiope argentata has multiple egg case spidroin loci with a loss of function at one locus. We also use degenerate PCR primers to search the genomic DNA of congeneric species and find evidence for multiple egg case spidroin loci in other Argiope spiders. Comparative analyses show that these multiple loci are more similar at the nucleotide level within a species than between species. This pattern is consistent with concerted evolution homogenizing gene copies within a genome. More complicated explanations include convergent evolution or recent independent gene duplications within each species.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

The occurrence of three D-J-C clusters within the dromedary TRB locus highlights a shared evolution in Tylopoda, Ruminantia and Suina

The αβ T cells are important components of the adaptive immune system and can recognize a vast array of peptides presented by MHC molecules. The ability of these T cells to recognize the complex depends on the diversity of the αβ TR, which is generated by a recombination of specific Variable, Diversity and Joining genes for the β chain, and Variable and Joining genes for the α chain. In this study, we analysed the genomic structure and the gene content of the TRB locus in Camelus dromedarius, which is a species belonging to the Tylopoda suborder. The most noteworthy result is the presence of three in tandem TRBD-J-C clusters in the dromedary TRB locus, which is similar to clusters found in sheep, cattle and pigs and suggests a common duplication event occurred prior to the Tylopoda/Ruminantia/Suina divergence. Conversely, a significant contraction of the dromedary TRBV genes, which was previously found in the TRG and TRD loci, was observed with respect to the other artiodactyl species.

Alignment-free sequence comparison: benefits, applications, and tools

Alignment-free sequence analyses have been applied to problems ranging from whole-genome phylogeny to the classification of protein families, identification of horizontally transferred genes, and detection of recombined sequences. The strength of these methods makes them particularly useful for next-generation sequencing data processing and analysis. However, many researchers are unclear about how these methods work, how they compare to alignment-based methods, and what their potential is for use for their research. We address these questions and provide a guide to the currently available alignment-free sequence analysis tools.

Progressive erosion of the Relaxin1 gene in bovids

The relaxin/insulin-like (RLN/INSL) gene family is a group of genes that encode peptide hormones involved in a variety of physiological functions related to reproduction. Previous studies have shown that relaxin plays a key role in widening of the pubic bone during labor and in gamete maturation. Because of these functions, studying the evolution of RLN1, the gene encoding for relaxin, is relevant in livestock species, most of which belong in the group Laurasiatheria, which includes cow, pig, horse, goat, and sheep in addition to bats, cetaceans and carnivores. Experimental evidence suggests that cows do not synthesize relaxin, but respond to it, and sheep apparently have a truncated RLN1 gene. Thus, we made use of genome sequence data to characterize the genomic locus of the RLN1 gene in Laurasiatherian mammals to better understand how cows lost the ability to synthesize this peptide. We found that all ruminants in our study (cow, giraffe, goat, sheep and Tibetan antelope) lack a functional RLN1 gene, and document the progressive loss of RLN1 in the lineage leading to cows. Our analyses indicate that 1 - all ruminants have lost all key regulatory elements upstream of the first exon, 2 - giraffe, goat, sheep and Tibetan antelope have multiple inactivating mutations in the RLN1 pseudogene, and 3 - the cow genome has lost all traces of RLN1. The 5' regulatory sequence plays a key role in activating expression, and the loss of this sequence would impair synthesis of mRNA. Our results suggest that changes in regulatory sequence preceded mutations in coding sequence and highlight the importance of these regions in maintaining proper gene function. In addition, we found that all bovids examined posses copies of the relaxin receptors, which explains why they are able to respond to relaxin despite their inability to produce it.

Comprehensive annotation and evolutionary insights into the canine (Canis lupus familiaris) antigen receptor loci

Dogs are an excellent model for human disease. For example, the treatment of canine lymphoma has been predictive of the human response to that treatment. However, an incomplete picture of canine (Canis lupus familiaris) immunoglobulin (IG) and T cell receptor (TR)-or antigen receptor (AR)-gene loci has restricted their utility. This work advances the annotation of the canine AR loci and looks into breed-specific features of the loci. Bioinformatic analysis of unbiased RNA sequence data was used to complete the annotation of the canine AR genes. This annotation was used to query 107 whole genome sequences from 19 breeds and identified over 5500 alleles across the 550 genes of the seven AR loci: the IG heavy, kappa, and lambda loci; and the TR alpha, beta, gamma, and delta loci. Of note was the discovery that half of the IGK variable (V) genes were located downstream of, and inverted with respect to, the rest of the locus. Analysis of the germline sequences of all the AR V genes identified greater conservation between dog and human than mouse with either. This work brings our understanding of the genetic diversity and expression of AR in dogs to the same completeness as that of mice and men, making it the third species to have all AR loci comprehensively and accurately annotated. The large number of germline sequences serves as a reference for future studies, and has allowed statistically powerful conclusions to be drawn on the pressures that have shaped these loci.

Evolution of the vertebrate insulin receptor substrate (Irs) gene family

BACKGROUND

Insulin receptor substrate (Irs) proteins are essential for insulin signaling as they allow downstream effectors to dock with, and be activated by, the insulin receptor. A family of four Irs proteins have been identified in mice, however the gene for one of these, IRS3, has been pseudogenized in humans. While it is known that the Irs gene family originated in vertebrates, it is not known when it originated and which members are most closely related to each other. A better understanding of the evolution of Irs genes and proteins should provide insight into the regulation of metabolism by insulin.

RESULTS

Multiple genes for Irs proteins were identified in a wide variety of vertebrate species. Phylogenetic and genomic neighborhood analyses indicate that this gene family originated very early in vertebrae evolution. Most Irs genes were duplicated and retained in fish after the fish-specific genome duplication. Irs genes have been lost of various lineages, including Irs3 in primates and birds and Irs1 in most fish. Irs3 and Irs4 experienced an episode of more rapid protein sequence evolution on the ancestral mammalian lineage. Comparisons of the conservation of the proteins sequences among Irs paralogs show that domains involved in binding to the plasma membrane and insulin receptors are most strongly conserved, while divergence has occurred in sequences involved in interacting with downstream effector proteins.

CONCLUSIONS

The Irs gene family originated very early in vertebrate evolution, likely through genome duplications, and in parallel with duplications of other components of the insulin signaling pathway, including insulin and the insulin receptor. While the N-terminal sequences of these proteins are conserved among the paralogs, changes in the C-terminal sequences likely allowed changes in biological function.

Complete Genome Sequence of Staphylococcus lutrae ATCC 700373, a Potential Pathogen Isolated from Deceased Otters

Despite their relevance to human health, not all staphylococcal species have been characterized. As such, the potential zoonotic threats posed by uninvestigated species and their contribution to the staphylococcal pangenome are unclear. Here, we report the complete genome sequence of Staphylococcus lutrae ATCC 700373, a coagulase-positive species isolated from deceased otters.

Genes coding for cyclin-dependent kinase inhibitors are fragile in Xenopus

Cell proliferation is strictly regulated by the dosage balance among cell-cycle regulators such as CDK/cyclin complexes and CDK-Inhibitors. Even in the allotetraploid genome of Xenopus laevis, the dosage balance must be maintained for animals to stay alive, and the duplicated homeologous genes seem to have gradually changed, through evolution, resulting in the best genes for them to thrive. In the Xenopus laevis genome, while homeologous gene pairs of CDKs are fundamentally maintained and a few cyclin genes are amplified, homeologous gene pairs of the important CDK-Inhibitors, CDKn1c and CDKn2a, are deleted from chromosomes L and S. Although losses of CDKn1c and CDKn2a can lead to diseases in humans, their loss in X. laevis does not affect the animals' health. Also, another gene coding CDKn1b is lost besides CDKn1c and CDKn2a in the genome of Xenopus tropicalis. These findings suggest a high resistance of Xenopus to diseases. We also found that CDKn2c.S expression is higher than that of CDKn2c.L, and a conserved noncoding sequence (CNS) of CDKn2c genomic loci on X. laevis chromosome S and X. tropicalis has an enhancement activity in regulating the different expression. These findings together indicate a surprising fragility of CDK inhibitor gene loci in the Xenopus genome in spite of their importance, and may suggest that factors other than CDK-inhibitors decelerate cell-cycling in Xenopus.

Application of a simplified method of chloroplast enrichment to small amounts of tissue for chloroplast genome sequencing

PREMISE OF THE STUDY

High-throughput sequencing of genomic DNA can recover complete chloroplast genome sequences, but the sequence data are usually dominated by sequences from nuclear/mitochondrial genomes. To overcome this deficiency, a simple enrichment method for chloroplast DNA from small amounts of plant tissue was tested for eight plant species including a gymnosperm and various angiosperms.

METHODS

Chloroplasts were enriched using a high-salt isolation buffer without any step gradient procedures, and enriched chloroplast DNA was sequenced by multiplexed high-throughput sequencing.

RESULTS

Using this simple method, significant enrichment of chloroplast DNA-derived reads was attained, allowing deep sequencing of chloroplast genomes. As an example, the chloroplast genome of the conifer Callitris sulcata was assembled, from which polymorphic microsatellite loci were isolated successfully.

DISCUSSION

This chloroplast enrichment method from small amounts of plant tissue will be particularly useful for studies that use sequencers with relatively small throughput and that cannot use large amounts of tissue (e.g., for endangered species).

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology

The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight Bacterial Artificial Chromosome (BAC) clones spanning the horse MHC class II region. The final assembly resulted in a 1,165,328 bp continuous gap free sequence with 35 manually curated genomic loci of which 23 were considered to be functional and 12 to be pseudogenes. In comparison to the MHC class II region in other mammals, the corresponding region in horse shows extraordinary copy number variation and different relative location and directionality of the Eqca-DRB, -DQA, -DQB and -DOB loci. This is the first long-read sequence assembly of the horse MHC class II region with rigorous manual gene annotation, and it will serve as an important resource for association studies of immune-mediated equine diseases and for evolutionary analysis of genetic diversity in this region.

QuickProbs 2: Towards rapid construction of high-quality alignments of large protein families

The ever-increasing size of sequence databases caused by the development of high throughput sequencing, poses to multiple alignment algorithms one of the greatest challenges yet. As we show, well-established techniques employed for increasing alignment quality, i.e., refinement and consistency, are ineffective when large protein families are investigated. We present QuickProbs 2, an algorithm for multiple sequence alignment. Based on probabilistic models, equipped with novel column-oriented refinement and selective consistency, it offers outstanding accuracy. When analysing hundreds of sequences, Quick-Probs 2 is noticeably better than ClustalΩ and MAFFT, the previous leaders for processing numerous protein families. In the case of smaller sets, for which consistency-based methods are the best performing, QuickProbs 2 is also superior to the competitors. Due to low computational requirements of selective consistency and utilization of massively parallel architectures, presented algorithm has similar execution times to ClustalΩ, and is orders of magnitude faster than full consistency approaches, like MSAProbs or PicXAA. All these make QuickProbs 2 an excellent tool for aligning families ranging from few, to hundreds of proteins.

Signatures of Crested Ibis MHC Revealed by Recombination Screening and Short-Reads Assembly Strategy

Whole-genome shotgun (WGS) sequencing has become a routine method in genome research over the past decade. However, the assembly of highly polymorphic regions in WGS projects remains a challenge, especially for large genomes. Employing BAC library constructing, PCR screening and Sanger sequencing, traditional strategy is laborious and expensive, which hampers research on polymorphic genomic regions. As one of the most highly polymorphic regions, the major histocompatibility complex (MHC) plays a central role in the adaptive immunity of all jawed vertebrates. In this study, we introduced an efficient procedure based on recombination screening and short-reads assembly. With this procedure, we constructed a high quality 488-kb region of crested ibis MHC that consists of 3 superscaffolds and contains 50 genes. Our sequence showed comparable quality (97.29% identity) to traditional Sanger assembly, while the workload was reduced almost 7 times. Comparative study revealed distinctive features of crested ibis by exhibiting the COL11A2-BLA-BLB-BRD2 cluster and presenting both ADPRH and odorant receptor (OR) gene in the MHC region. Furthermore, the conservation of the BF-TAP1-TAP2 structure in crested ibis and other vertebrate lineages is interesting in light of the hypothesis that coevolution of functionally related genes in the primordial MHC is responsible for the appearance of the antigen presentation pathways at the birth of the adaptive immune system.

Deletion of conserved sequences in IG-DMR at Dlk1-Gtl2 locus suggests their involvement in expression of paternally expressed genes in mice

Expression regulation of the Dlk1-Dio3 imprinted domain by the intergenic differentially methylated region (IG-DMR) is essential for normal embryonic development in mammals. In this study, we investigated conserved IG-DMR genomic sequences in eutherians to elucidate their role in genomic imprinting of the Dlk1-Dio3 domain. Using a comparative genomics approach, we identified three highly conserved sequences in IG-DMR. To elucidate the functions of these sequences in vivo, we generated mutant mice lacking each of the identified highly conserved sequences using the CRISPR/Cas9 system. Although mutant mice did not exhibit the gross phenotype, deletions of the conserved sequences altered the expression levels of paternally expressed imprinted genes in the mutant embryos without skewing imprinting status. These results suggest that the conserved sequences in IG-DMR are involved in the expression regulation of some of the imprinted genes in the Dlk1-Dio3 domain.

Coevolution of paired receptors in Xenopus carcinoembryonic antigen-related cell adhesion molecule families suggests appropriation as pathogen receptors

Background

In mammals, CEACAM1 and closely related members represent paired receptors with similar extracellular ligand-binding regions and cytoplasmic domains with opposing functions. Human CEACAM1 and CEACAM3 which have inhibitory ITIM/ITSM and activating ITAM-like motifs, respectively, in their cytoplasmic regions are such paired receptors. Various bacterial pathogens bind to CEACAM1 on epithelial and immune cells facilitating both entry into the host and down-regulation of the immune response whereas interaction with granulocyte-specific CEACAM3 leads to their uptake and destruction. It is unclear whether paired CEACAM receptors also exist in other vertebrate clades.

Results

We identified more than 80 ceacam genes in Xenopus tropicalis and X. laevis. They consist of two subgroups containing one or two putative paired receptor pairs each. Analysis of genomic sequences of paired receptors provide evidence that their highly similar ligand binding domains were adjusted by recent gene conversion events. In contrast, selection for diversification is observed among inhibitory receptor orthologs of the two frogs which split some 60 million years ago. The allotetraploid X. laevis arose later by hybridization of two closely related species. Interestingly, despite the conservation of the genomic landscape surrounding the homeologous ceacam loci only one locus resembles the one found in X. tropicalis. From the second X. laevis locus more than 80 % of the ceacam genes were lost including 5 of the 6 paired receptor genes. This suggests that once the gene for one of the paired receptors is lost the remaining gene cluster degrades rapidly probably due to lack of selection pressure exerted by pathogens.

Conclusions

The presence of paired receptors and selection for diversification suggests that also in amphibians CEACAM1-related inhibitory proteins are or were used as pathogen receptors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3279-9) contains supplementary material, which is available to authorized users.

Developmental mechanisms of stripe patterns in rodents

Mammalian colour patterns are among the most recognizable characteristics found in nature and can have a profound impact on fitness. However, little is known about the mechanisms underlying the formation and subsequent evolution of these patterns. Here we show that, in the African striped mouse (Rhabdomys pumilio), periodic dorsal stripes result from underlying differences in melanocyte maturation, which give rise to spatial variation in hair colour. We identify the transcription factor ALX3 as a regulator of this process. In embryonic dorsal skin, patterned expression of Alx3 precedes pigment stripes and acts to directly repress Mitf, a master regulator of melanocyte differentiation, thereby giving rise to light-coloured hair. Moreover, Alx3 is upregulated in the light stripes of chipmunks, which have independently evolved a similar dorsal pattern. Our results show a previously undescribed mechanism for modulating spatial variation in hair colour and provide insights into how phenotypic novelty evolves.

Pig lacks functional NLRC4 and NAIP genes

The NLRC4 inflammasome, which recognizes flagellin and components of the type III secretion system, plays an important role in the clearance of intracellular bacteria. Here, we examined the genomic sequences carrying two genes encoding key components of the NLRC4 inflammasome-NLR family, CARD-containing 4 (NLRC4), and NLR apoptosis inhibitory protein (NAIP)-in pigs. Pigs have a single locus encoding NLRC4 and NAIP. Comparison of the sequences thus obtained with the corresponding regions in humans revealed the deletion of intermediate exons in both pig genes. In addition, the genomic sequences of both pig genes lacked valid open reading frames encoding functional NLRC4 or NAIP protein. Additional pigs representing multiple breeds and wild boars also lacked the exons that we failed to find through genome sequencing. Furthermore, neither the NLRC4 nor the NAIP gene was expressed in pigs. These findings indicate that pigs lack the NLRC4 inflammasome, an important factor involved in monitoring bacterial proteins and contributing to the clearance of intracellular pathogens. These results also suggest that genetic polymorphisms affecting the molecular functions of TLR2, TLR4, TLR5, and other pattern recognition receptors associated with the recognition of bacteria have a more profound influence on disease resistance in pigs than in other species.

Genomic and expression analyses of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) loci reveal a similar basic public γδ repertoire in dolphin and human

BACKGROUND

The bottlenose dolphin (Tursiops truncatus) is a mammal that belongs to the Cetartiodactyla and have lived in marine ecosystems for nearly 60 millions years. Despite its popularity, our knowledge about its adaptive immunity and evolution is very limited. Furthermore, nothing is known about the genomics and evolution of dolphin antigen receptor immunity.

RESULTS

Here we report a evolutionary and expression study of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) genes. We have identified in silico the TRG and TRA/TRD genes and analyzed the relevant mature transcripts in blood and in skin from four subjects. The dolphin TRG locus is the smallest and simplest of all mammalian loci as yet studied. It shows a genomic organization comprising two variable (V1 and V2), three joining (J1, J2 and J3) and a single constant (C), genes. Despite the fragmented nature of the genome assemblies, we deduced the TRA/TRD locus organization, with the recent TRDV1 subgroup genes duplications, as it is expected in artiodactyls. Expression analysis from blood of a subject allowed us to assign unambiguously eight TRAV genes to those annotated in the genomic sequence and to twelve new genes, belonging to five different subgroups. All transcripts were productive and no relevant biases towards TRAV-J rearrangements are observed. Blood and skin from four unrelated subjects expression data provide evidence for an unusual ratio of productive/unproductive transcripts which arise from the TRG V-J gene rearrangement and for a "public" gamma delta TR repertoire. The productive cDNA sequences, shared both in the same and in different individuals, include biases of the TRGV1 and TRGJ2 genes. The high frequency of TRGV1-J2/TRDV1- D1-J4 productive rearrangements in dolphins may represent an interesting oligo-clonal population comparable to that found in human with the TRGV9- JP/TRDV2-D-J T cells and in primates.

CONCLUSIONS

Although the features of the TRG and TRA/TRD loci organization reflect those of the so far examined artiodactyls, genomic results highlight in dolphin an unusually simple TRG locus. The cDNA analysis reveal productive TRA/TRD transcripts and unusual ratios of productive/unproductive TRG transcripts. Comparing multiple different individuals, evidence is found for a "public" gamma delta TCR repertoire thus suggesting that in dolphins as in human the gamma delta TCR repertoire is accompanied by selection for public gamma chain.

Distribution and Evolution of Yersinia Leucine-Rich Repeat Proteins

Leucine-rich repeat (LRR) proteins are widely distributed in bacteria, playing important roles in various protein-protein interaction processes. In Yersinia, the well-characterized type III secreted effector YopM also belongs to the LRR protein family and is encoded by virulence plasmids. However, little has been known about other LRR members encoded by Yersinia genomes or their evolution. In this study, the Yersinia LRR proteins were comprehensively screened, categorized, and compared. The LRR proteins encoded by chromosomes (LRR1 proteins) appeared to be more similar to each other and different from those encoded by plasmids (LRR2 proteins) with regard to repeat-unit length, amino acid composition profile, and gene expression regulation circuits. LRR1 proteins were also different from LRR2 proteins in that the LRR1 proteins contained an E3 ligase domain (NEL domain) in the C-terminal region or an NEL domain-encoding nucleotide relic in flanking genomic sequences. The LRR1 protein-encoding genes (LRR1 genes) varied dramatically and were categorized into 4 subgroups (a to d), with the LRR1a to -c genes evolving from the same ancestor and LRR1d genes evolving from another ancestor. The consensus and ancestor repeat-unit sequences were inferred for different LRR1 protein subgroups by use of a maximum parsimony modeling strategy. Structural modeling disclosed very similar repeat-unit structures between LRR1 and LRR2 proteins despite the different unit lengths and amino acid compositions. Structural constraints may serve as the driving force to explain the observed mutations in the LRR regions. This study suggests that there may be functional variation and lays the foundation for future experiments investigating the functions of the chromosomally encoded LRR proteins of Yersinia.

Evolution of Heat Sensors Drove Shifts in Thermosensation between Xenopus Species Adapted to Different Thermal Niches

Temperature is one of the most critical environmental factors affecting survival, and thus species that inhabit different thermal niches have evolved thermal sensitivities suitable for their respective habitats. During the process of shifting thermal niches, various types of genes expressed in diverse tissues, including those of the peripheral to central nervous systems, are potentially involved in the evolutionary changes in thermosensation. To elucidate the molecular mechanisms behind the evolution of thermosensation, thermal responses were compared between two species of clawed frogs (Xenopus laevis and Xenopus tropicalis) adapted to different thermal environments. X. laevis was much more sensitive to heat stimulation than X. tropicalis at the behavioral and neural levels. The activity and sensitivity of the heat-sensing TRPA1 channel were higher in X. laevis compared with those of X. tropicalis The thermal responses of another heat-sensing channel, TRPV1, also differed between the two Xenopus species. The species differences in Xenopus TRPV1 heat responses were largely determined by three amino acid substitutions located in the first three ankyrin repeat domains, known to be involved in the regulation of rat TRPV1 activity. In addition, Xenopus TRPV1 exhibited drastic species differences in sensitivity to capsaicin, contained in chili peppers, between the two Xenopus species. Another single amino acid substitution within Xenopus TRPV1 is responsible for this species difference, which likely alters the neural and behavioral responses to capsaicin. These combined subtle amino acid substitutions in peripheral thermal sensors potentially serve as a driving force for the evolution of thermal and chemical sensation.

Haplotypes of the D-Amino Acid Oxidase Gene Are Significantly Associated with Schizophrenia and Its Neurocognitive Deficits

D-amino acid oxidase (DAO) has been reported to be associated with schizophrenia. This study aimed to search for genetic variants associated with this gene. The genomic regions of all exons, highly conserved regions of introns, and promoters of this gene were sequenced. Potentially meaningful single-nucleotide polymorphisms (SNPs) obtained from direct sequencing were selected for genotyping in 600 controls and 912 patients with schizophrenia and in a replicated sample consisting of 388 patients with schizophrenia. Genetic associations were examined using single-locus and haplotype association analyses. In single-locus analyses, the frequency of the C allele of a novel SNP rs55944529 located at intron 8 was found to be significantly higher in the original large patient sample (p = 0.016). This allele was associated with a higher level of DAO mRNA expression in the Epstein-Barr virus-transformed lymphocytes. The haplotype distribution of a haplotype block composed of rs11114083-rs2070586-rs2070587-rs55944529 across intron 1 and intron 8 was significantly different between the patients and controls and the haplotype frequencies of AAGC were significantly higher in patients, in both the original (corrected p < 0.0001) and replicated samples (corrected p = 0.0003). The CGTC haplotype was specifically associated with the subgroup with deficits in sustained attention and executive function and the AAGC haplotype was associated with the subgroup without such deficits. The DAO gene was a susceptibility gene for schizophrenia and the genomic region between intron 1 and intron 8 may harbor functional genetic variants, which may influence the mRNA expression of DAO and neurocognitive functions in schizophrenia.

Organization and characteristics of the major histocompatibility complex class II region in the Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis)

Little is known about the major histocompatibility complex (MHC) in the genome of Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) (YFP) or other cetaceans. In this study, a high-quality YFP bacterial artificial chromosome (BAC) library was constructed. We then determined the organization and characterization of YFP MHC class II region by screening the BAC library, followed by sequencing and assembly of positive BAC clones. The YFP MHC class II region consists of two segregated contigs (218,725 bp and 328,435 bp respectively) that include only eight expressed MHC class II genes, three pseudo MHC genes and twelve non-MHC genes. The YFP has fewer MHC class II genes than ruminants, showing locus reduction in DRB, DQA, DQB, and loss of DY. In addition, phylogenic and evolutionary analyses indicated that the DRB, DQA and DQB genes might have undergone birth-and-death evolution, whereas the DQB gene might have evolved under positive selection in cetaceans. These findings provide an essential foundation for future work, such as estimating MHC genetic variation in the YFP or other cetaceans. This work is the first report on the MHC class II region in cetaceans and offers valuable information for understanding the evolution of MHC genome in cetaceans.

Contrasting Patterns in the Evolution of Vertebrate MLX Interacting Protein (MLXIP) and MLX Interacting Protein-Like (MLXIPL) Genes

ChREBP and MondoA are glucose-sensitive transcription factors that regulate aspects of energy metabolism. Here we performed a phylogenomic analysis of Mlxip (encoding MondoA) and Mlxipl (encoding ChREBP) genes across vertebrates. Analysis of extant Mlxip and Mlxipl genes suggests that the most recent common ancestor of these genes was composed of 17 coding exons. Single copy genes encoding both ChREBP and MondoA, along with their interacting partner Mlx, were found in diverse vertebrate genomes, including fish that have experienced a genome duplication. This observation suggests that a single Mlx gene has been retained to maintain coordinate regulation of ChREBP and MondoA. The ChREBP-β isoform, the more potent and constitutively active isoform, appeared with the evolution of tetrapods and is absent from the Mlxipl genes of fish. Evaluation of the conservation of ChREBP and MondoA sequences demonstrate that MondoA is better conserved and potentially mediates more ancient function in glucose metabolism.

Draft Genome Sequence of Lactococcus lactis subsp. lactis bv. diacetylactis CRL264, a Citrate-Fermenting Strain

We report the draft genome sequence of Lactococcus lactis subsp. lactis bv. diacetylactis CRL264, a natural strain isolated from artisanal cheese from northwest Argentina. L. lactis subsp. lactis bv. diacetylactis is one of the most important microorganisms used as starter culture around the world. The CRL264 strain constitutes a model microorganism in the studies on the generation of aroma compounds (diacetyl, acetoin, and 2,3-butanediol) by lactic acid bacteria. Our genome analysis shows similar genetic organization to other available genomes of L. lactis bv. diacetylactis strains.

Draft Genome Sequences of Four Enterococcus faecium Strains Isolated from Argentine Cheese

We report the draft genome sequences of four Enterococcus faecium strains isolated from Argentine regional cheeses. These strains were selected based on their technological properties, i.e., their ability to produce aroma compounds (diacetyl, acetoin, and 2,3-butanediol) from citrate. The goal of our study is to provide further genetic evidence for the rational selection of enterococci strains based on their pheno- and genotype in order to be used in cheese production.

Gene structure variation in segmental duplication block C of human chromosome 7q 11.23 during primate evolution

Segmental duplication, or low-copy repeat (LCR) event, occurs during primate evolution and is an important source of genomic diversity, including gain or loss of gene function. The human chromosome 7q 11.23 is related to the William-Beuren syndrome and contains large region-specific LCRs composed of blocks A, B, and C that have different copy numbers in humans and different primates. We analyzed the structure of POM121, NSUN5, FKBP6, and TRIM50 genes in the LCRs of block C. Based on computational analysis, POM121B created by a segmental duplication acquired a new exonic region, whereas NSUN5B (NSUN5C) showed structural variation by integration of HERV-K LTR after duplication from the original NSUN5 gene. The TRIM50 gene originally consists of seven exons, whereas the duplicated TRIM73 and TRIM74 genes present five exons because of homologous recombination-mediated deletion. In addition, independent duplication events of the FKBP6 gene generated two pseudogenes at different genomic locations. In summary, these clustered genes are created by segmental duplication, indicating that they show dynamic evolutionary events, leading to structure variation in the primate genome.

Single fluorophore melting curve analysis for detection of hypervirulent Clostridium difficile

This study demonstrates a novel detection assay able to identify and subtype strains of Clostridium difficile. Primers carefully designed for melting curve analysis amplify DNA from three C. difficile genes, tcdB, tcdC and cdtB, during quantitative (q)PCR. The tcdB gene allows for confirmation of organism presence, whilst the tcdC and cdtB genes allow for differentiation of virulence status, as deletions in the tcdC gene and the concurrent presence of the cdtB gene, which produces binary toxin, are associated with hypervirulence. Following qPCR, subtyping is then achieved by automated, inline melting curve analysis using only a single intercalating dye and verified by microchip electrophoresis. This assay represents a novel means of distinguishing between toxigenic and hypervirulent C. difficile strains NAP1/027/BI and 078 ribotype, which are highly prevalent hypervirulent strains in humans. This methodology can help rapidly detect and identify C. difficile strains that impose a significant health and economic burden in hospitals and other healthcare settings.

Rapid evolutionary divergence of Gossypium barbadense and G. hirsutum mitochondrial genomes

Background

The mitochondrial genome from upland cotton, G. hirsutum, was previously sequenced. To elucidate the evolution of mitochondrial genomic diversity within a single genus, we sequenced the mitochondrial genome from Sea Island cotton (Gossypium barbadense L.).

Methods

Mitochondrial DNA from week-old etiolated seedlings was extracted from isolated organelles using discontinuous sucrose density gradient method. Mitochondrial genome was sequenced with Solexa using paired-end, 90 bp read. The clean reads were assembled into contigs using ABySS and finished via additional fosmid and BAC sequencing. Finally, the genome was annotated and analyzed using different softwares.

Results

The G. barbadense (Sea Island cotton) mitochondrial genome was fully sequenced (677,434-bp) and compared to the mitogenome of upland cotton. The G. barbadense mitochondrial DNA contains seven more genes than that of upland cotton, with a total of 40 protein coding genes (excluding possible pseudogenes), 6 rRNA genes, and 29 tRNA genes. Of these 75 genes, atp1, mttB, nad4, nad9, rrn5, rrn18, and trnD(GTC)-cp were each represented by two identical copies. A single 64 kb repeat was largely responsible for the 9 % difference in genome size between the two mtDNAs. Comparison of genome structures between the two mitochondrial genomes revealed 8 rearranged syntenic regions and several large repeats. The largest repeat was missing from the master chromosome in G. hirsutum. Both mitochondrial genomes contain a duplicated copy of rps3 (rps3-2) in conjunction with a duplication of repeated sequences. Phylogenetic and divergence considerations suggest that a 544-bp fragment of rps3 was transferred to the nuclear genome shortly after divergence of the A- and D- genome diploid cottons.

Conclusion

These results highlight the insights to the evolution of structural variation between Sea Island and upland cotton mitochondrial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1988-0) contains supplementary material, which is available to authorized users.

Sheep (Ovis aries) T cell receptor alpha (TRA) and delta (TRD) genes and genomic organization of the TRA/TRD locus

Background

In mammals, T cells develop along two discrete pathways characterized by expression of either the αβ or the γδ T cell receptors. Human and mouse display a low peripheral blood γδ T cell percentage ("γδ low species") while sheep, bovine and pig accounts for a high proportion of γδ T lymphocytes ("γδ high species"). While the T cell receptor alpha (TRA) and delta (TRD) genes and the genomic organization of the TRA/TRD locus has been determined in human and mouse, this information is still poorly known in artiodactyl species, such as sheep.

Results

The analysis of the current Ovis aries whole genome assembly, Oar_v3.1, revealed that, as in the other mammalian species, the sheep TRD locus is nested within the TRA locus. In the most 5’ part the TRA/TRD locus contains TRAV genes which are intermingled with TRDV genes, then TRD genes which include seven TRDD, four TRDJ genes, one TRDC and a single TRDV gene with an inverted transcriptional orientation, and finally in the most 3’ part, the TRA locus is completed by 61 TRAJ genes and one TRAC gene.

Comparative sequence and analysis and annotation led to the identification of 66 TRAV genes assigned to 34 TRAV subgroups and 25 TRDV genes belonging to the TRDV1 subgroup, while one gene was found for each TRDV2, TRDV3 and TRDV4 subgroups. Multiple duplication events within several TRAV subgroups have generated the sheep TRAV germline repertoire, which is substantially larger than the human one. A significant proportion of these TRAV gene duplications seems to have occurred simultaneously with the amplification of the TRDV1 subgroup genes. This dynamic of expansion has also generated novel multigene subgroups, which are species-specific. Ovis aries TRA and TRD genes identified in this study were assigned IMGT definitive or temporary names and were approved by the IMGT/WHO-IUIS nomenclature committee.

The completeness of the genome assembly in the 3' part of the locus has allowed us to interpret rearranged CDR3 of cDNA from both TRA and TRD chain repertoires. The involvement of one up to four TRDD genes into a single transcript makes the potential sheep TRD chain much larger than any known TR chain repertoire.

Conclusions

The sheep genome, as the bovine genome, contains a large and diverse repertoire of TRA and TRD genes when compared to the “γδ T cell low” species genomes. The composition and length of the rearranged CDR3 in TRD V-delta domains influence the three-dimensional configuration of the antigen-combining site thus suggesting that in ruminants, γδ T cells play a more important and specific role in immune recognition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1790-z) contains supplementary material, which is available to authorized users.

Genomic approaches to studying human-specific developmental traits

Changes in developmental regulatory programs drive both disease and phenotypic differences among species. Linking human-specific traits to alterations in development is challenging, because we have lacked the tools to assay and manipulate regulatory networks in human and primate embryonic cells. This field was transformed by the sequencing of hundreds of genomes – human and non-human – that can be compared to discover the regulatory machinery of genes involved in human development. This approach has identified thousands of human-specific genome alterations in developmental genes and their regulatory regions. With recent advances in stem cell techniques, genome engineering, and genomics, we can now test these sequences for effects on developmental gene regulation and downstream phenotypes in human cells and tissues.

A nonclassical MHC class I U lineage locus in zebrafish with a null haplotypic variant

Three sequence lineages of MHC class I genes have been described in zebrafish (Danio rerio): U, Z, and L. The U lineage genes encoded on zebrafish chromosome 19 are predicted to provide the classical function of antigen presentation. This MHC class I locus displays significant haplotypic variation and is the only MHC class I locus in zebrafish that shares conserved synteny with the core mammalian MHC. Here, we describe two MHC class I U lineage genes, mhc1ula and mhc1uma, that map to chromosome 22. Unlike the U lineage proteins encoded on chromosome 19, Ula and Uma likely play a nonclassical role as they lack conservation of key peptide binding residues, display limited polymorphic variation, and exhibit tissue-specific expression. We also describe a null haplotype at this chromosome 22 locus in which the mhc1ula and mhc1uma genes are absent due to a ~30 kb deletion with no other MHC class I sequences present. Functional and non-functional transcripts of mhc1ula and mhc1uma were identified; however, mhc1uma transcripts were often not amplified or amplified at low levels from individuals possessing an apparently bona fide gene. These distinct U lineage genes may be restricted to the superorder Ostariophysi as similar sequences only could be identified from the blind cavefish (Astyanax mexicanus), fathead minnow (Pimephales promelas), goldfish (Carassius auratus), and grass carp (Ctenopharyngodon idella).

Evolutionary Analyses and Natural Selection of Betaine-Homocysteine S-Methyltransferase (BHMT) and BHMT2 Genes

Betaine-homocysteine S-methyltransferase (BHMT) and BHMT2 convert homocysteine to methionine using betaine and S-methylmethionine, respectively, as methyl donor substrates. Increased levels of homocysteine in blood are associated with cardiovascular disease. Given their role in human health and nutrition, we identified BHMT and BHMT2 genes and proteins from 38 species of deuterostomes including human and non-human primates. We aligned the genes to look for signatures of selection, to infer evolutionary rates and events across lineages, and to identify the evolutionary timing of a gene duplication event that gave rise to two genes, BHMT and BHMT2. We found that BHMT was present in the genomes of the sea urchin, amphibians, reptiles, birds and mammals; BHMT2 was present only across mammals. BHMT and BHMT2 were present in tandem in the genomes of all monotreme, marsupial and placental species examined. Evolutionary rates were accelerated for BHMT2 relative to BHMT. Selective pressure varied across lineages, with the highest dN/dS ratios for BHMT and BHMT2 occurring immediately following the gene duplication event, as determined using GA Branch analysis. Nine codons were found to display signatures suggestive of positive selection; these contribute to the enzymatic or oligomerization domains, suggesting involvement in enzyme function. Gene duplication likely occurred after the divergence of mammals from other vertebrates but prior to the divergence of extant mammalian subclasses, followed by two deletions in BHMT2 that affect oligomerization and methyl donor specificity. The faster evolutionary rate of BHMT2 overall suggests that selective constraints were reduced relative to BHMT. The dN/dS ratios in both BHMT and BHMT2 was highest following the gene duplication, suggesting that purifying selection played a lesser role as the two paralogs diverged in function.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.

A Vaccine Approach for the Prevention of Infections by Multidrug-resistant Enterococcus faecium

The incidence of multidrug-resistant Enterococcus faecium hospital infections has been steadily increasing. With the goal of discovering new vaccine antigens, we systematically fractionated and purified four distinct surface carbohydrates from E. faecium endocarditis isolate Tx16, shown previously to be resistant to phagocytosis in the presence of human serum. The two most abundant polysaccharides consist of novel branched heteroglycan repeating units that include signature sugars altruronic acid and legionaminic acid, respectively. A minor high molecular weight polysaccharide component was recognized as the fructose homopolymer levan, and a glucosylated lipoteichoic acid (LTA) was identified in a micellar fraction. The polysaccharides were conjugated to the CRM197 carrier protein, and the resulting glycoconjugates were used to immunize rabbits. Rabbit immune sera were evaluated for their ability to kill Tx16 in opsonophagocytic assays and in a mouse passive protection infection model. Although antibodies raised against levan failed to mediate opsonophagocytic killing, the other glycoconjugates induced effective opsonic antibodies, with the altruronic acid-containing polysaccharide antisera showing the greatest opsonophagocytic assay activity. Antibodies directed against either novel heteroglycan or the LTA reduced bacterial load in mouse liver or kidney tissue. To assess antigen prevalence, we screened a diverse collection of blood isolates (n = 101) with antibodies to the polysaccharides. LTA was detected on the surface of 80% of the strains, and antigens recognized by antibodies to the two major heteroglycans were co-expressed on 63% of these clinical isolates. Collectively, these results represent the first steps toward identifying components of a glycoconjugate vaccine to prevent E. faecium infection.

Evolution of the Vertebrate Resistin Gene Family

Resistin (encoded by Retn) was previously identified in rodents as a hormone associated with diabetes; however human resistin is instead linked to inflammation. Resistin is a member of a small gene family that includes the resistin-like peptides (encoded by Retnl genes) in mammals. Genomic searches of available genome sequences of diverse vertebrates and phylogenetic analyses were conducted to determine the size and origin of the resistin-like gene family. Genes encoding peptides similar to resistin were found in Mammalia, Sauria, Amphibia, and Actinistia (coelacanth, a lobe-finned fish), but not in Aves or fish from Actinopterygii, Chondrichthyes, or Agnatha. Retnl originated by duplication and transposition from Retn on the early mammalian lineage after divergence of the platypus, but before the placental and marsupial mammal divergence. The resistin-like gene family illustrates an instance where the locus of origin of duplicated genes can be identified, with Retn continuing to reside at this location. Mammalian species typically have a single copy Retn gene, but are much more variable in their numbers of Retnl genes, ranging from 0 to 9. Since Retn is located at the locus of origin, thus likely retained the ancestral expression pattern, largely maintained its copy number, and did not display accelerated evolution, we suggest that it is more likely to have maintained an ancestral function, while Retnl, which transposed to a new location, displays accelerated evolution, and shows greater variability in gene number, including gene loss, likely evolved new, but potentially lineage-specific, functions.

YOC, A new strategy for pairwise alignment of collinear genomes

Background

Comparing and aligning genomes is a key step in analyzing closely related genomes. Despite the development of many genome aligners in the last 15 years, the problem is not yet fully resolved, even when aligning closely related bacterial genomes of the same species. In addition, no procedures are available to assess the quality of genome alignments or to compare genome aligners.

Results

We designed an original method for pairwise genome alignment, named YOC, which employs a highly sensitive similarity detection method together with a recent collinear chaining strategy that allows overlaps. YOC improves the reliability of collinear genome alignments, while preserving or even improving sensitivity. We also propose an original qualitative evaluation criterion for measuring the relevance of genome alignments. We used this criterion to compare and benchmark YOC with five recent genome aligners on large bacterial genome datasets, and showed it is suitable for identifying the specificities and the potential flaws of their underlying strategies.

Conclusions

The YOC prototype is available at https://github.com/ruricaru/YOC. It has several advantages over existing genome aligners: (1) it is based on a simplified two phase alignment strategy, (2) it is easy to parameterize, (3) it produces reliable genome alignments, which are easier to analyze and to use.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0530-3) contains supplementary material, which is available to authorized users.