Physical mapping in large genomes: accelerating anchoring of BAC contigs to genetic maps through in silico analysis

Anchored physical maps represent essential frameworks for map-based cloning, comparative genomics studies, and genome sequencing projects. High throughput anchoring can be achieved by polymerase chain reaction (PCR) screening of bacterial artificial chromosome (BAC) library pools with molecular markers. However, for large genomes such as wheat, the development of high dimension pools and the number of reactions that need to be performed can be extremely large making the screening laborious and costly. To improve the cost efficiency of anchoring in such large genomes, we have developed a new software named Elephant (electronic physical map anchoring tool) that combines BAC contig information generated by FingerPrinted Contig with results of BAC library pools screening to identify BAC addresses with a minimal amount of PCR reactions. Elephant was evaluated during the construction of a physical map of chromosome 3B of hexaploid wheat. Results show that a one dimensional pool screening can be sufficient to anchor a BAC contig while reducing the number of PCR by 384-fold thereby demonstrating that Elephant is an efficient and cost-effective tool to support physical mapping in large genomes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. E. Paux and F. Legeai contributed equally to this work.     

SALMFamide2 and serotonin immunoreactivity in the nervous system of some acoels (Xenacoelomorpha)

Acoel worms are simple, often microscopic animals with direct development, a multiciliated epidermis, a statocyst, and a digestive parenchyma instead of a gut epithelium. Morphological characters of acoels have been notoriously difficult to interpret due to their relative scarcity. The nervous system is one of the most accessible and widely used comparative features in acoels, which have a so‐called commissural brain without capsule and several major longitudinal neurite bundles. Here, we use the selective binding properties of a neuropeptide antibody raised in echinoderms (SALMFamide2, or S2), and a commercial antibody against serotonin (5‐HT) to provide additional characters of the acoel nervous system. We have prepared whole‐mount immunofluorescent stainings of three acoel species: Symsagittifera psammophila (Convolutidae), Aphanostoma pisae, and the model acoel Isodiametra pulchra (both Isodiametridae). The commissural brain of all three acoels is delimited anteriorly by the ventral anterior commissure, and posteriorly by the dorsal posterior commissure. The dorsal anterior commissure is situated between the ventral anterior commissure and the dorsal posterior commissure, while the statocyst lies between dorsal anterior and dorsal posterior commissure. S2 and serotonin do not co‐localise, and they follow similar patterns to each other within an animal. In particular, S2, but not 5‐HT, stains a prominent commissure posterior to the main (dorsal) posterior commissure. We have for the first time observed a closed posterior loop of the main neurite bundles in S. psammophila for both the amidergic and the serotonergic nervous system. In I. pulchra, the lateral neurite bundles also form a posterior loop in our serotonergic nervous system stainings.     

Molecular characterization and genomic organization of low molecular weight glutenin subunit genes at the Glu-3 loci in hexaploid wheat (Triticum aestivum L.)

In this study, we report on the molecular characterization and genomic organization of the low molecular weight glutenin subunit (LMW-GS) gene family in hexaploid wheat (Triticum aestivum L.). Eighty-two positive BAC clones were identified to contain LMW-GS genes from the hexaploid wheat ‘Glenlea’ BAC library via filter hybridization and PCR validation. Twelve unique LMW glutenin genes and seven pseudogenes were isolated from these positive BAC clones by primer-template mismatch PCR and subsequent primer walking using hemi-nested touchdown PCR. These genes were sequenced and each consisted of a single-open reading frame (ORF) and untranslated 5′ and 3′ flanking regions. All 12 LMW glutenin subunits contained eight cysteine residues. The LMW-m-type subunits are the most abundant in hexaploid wheat. Of the 12 LMW-GS, 1, 2 and 9 are i-type, s-type and m-type, respectively. The phylogenetic analysis suggested that the LMW-i type gene showed greater differences to LMW-s and LMW-m-type genes, which, in turn, were more closely related to one another. On the basis of their N-terminal sequences, they were classified into nine groups. Fingerprinting of the 82 BAC clones indicated 30 BAC clones assembled into eight contigs, while the remaining clones were singletons. BAC end sequencing of the 82 clones revealed that long terminal repeat (LTR) retrotransposons were abundant in the Glu-3 regions. The average physical distance between two adjacent LMW-GS genes was estimated to be 81 kb. Most of LMW-GS genes are located in the d-genome, suggesting that the Glu-D3 locus is much larger than the Glu-B3 locus and Glu-A3 locus. Alignments of sequences indicated that the same type (starting with the same N-terminal sequence) LMW-GS genes were highly conserved in the homologous genomes between hexaploid wheat and its donors such as durum wheat and T. tauschii. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular cloning of zebra finch W chromosome repetitive sequences: evolution of the avian W chromosome

The zebra finch (Taeniopygia guttata) has a large Z chromosome and highly condensed W chromosome. We used the random amplified polymorphic DNA (RAPD) polymerase chain reaction (PCR) technique to isolate female-specific sequences ZBM1 and ZBM2. Southern blot hybridization to male and female zebra finch genomic DNA suggested that these sequences were located on the W chromosome, although homologous sequences appeared to be autosomal or Z-linked. Fluorescent in situ hybridization (FISH) using bacterial artificial chromosome (BAC) clones corresponding to ZBM sequences showed hybridization to the whole W chromosome, suggesting that the BACs encode sequences that are repeated across the entire W chromosome. Based on the sequencing of a ZBM repetitive sequence and Z chromosome derived BAC clones, we demonstrate a random distribution of repeat sequences that are specific to the W chromosome or encoded by both Z and W. The positions of ZW-common repeat sequences mapped to a noncoding region of a Z chromosome BAC clone containing the CHD1Z gene. The apparent lineage-specificity of W chromosome repeat sequences in passerines and galliform birds suggest that the W chromosome had not differentiated well from the Z at the time of divergence of these lineages. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development of a BAC library for yellow-poplar (Liriodendron tulipifera) and the identification of genes associated with flower development and lignin biosynthesis

Liriodendron tulipifera L., a member of the Magnoliaceae, occupies an important phylogenetic position as a basal angiosperm that has retained numerous putatively ancestral morphological characters, and thus has often been used in studies of the evolution of flowering plants and of specific gene families. However, genomic resources for these early branching angiosperm lineages are very limited. In this study, we describe the construction of a large-insert bacterial artificial chromosome (BAC) library from L. tulipifera. Flow cytometry estimates that this nuclear genome is approximately 1,802 Mbp per haploid genome (±16 SD). The BAC library contains 73,728 clones, a 4.8-fold genome coverage, with an average insert size of 117 kb, a chloroplast DNA content of 0.2%, and little to no bacterial sequences nor empty vector content clones. As a test of the utility of this BAC library, we screened the library with six single/low-copy genic probes. We obtained at least two positive clones for each gene and confirmed the clones by DNA sequencing. A total of 182 paired end sequences were obtained from 96 of the BAC clones. Using BLAST searches, we found that 25% of the BAC end sequences were similar to DNA sequences in GenBank. Of these, 68% shared sequence with transposable elements and 25% with genes from other taxa. This result closely reflected the content of random sequences obtained from a small insert genomic library for L. tulipifera, indicating that the BAC library construction process was not biased. The first genomic DNA sequences for Liriodendron genes are also reported. All the Liriodendron genomic sequences described in this paper have been deposited in the GenBank data library. The end sequences from shotgun genomic clones and BAC clones are under accession DU169330–DU169684. Partial sequences of Gigantea, Frigida, LEAFY, cinnamyl alcohol dehydrogenase, 4-coumarate:CoA ligase, and phenylalanine ammonia-lyase genes are under accession DQ223429–DQ223434. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Physical analysis of the complex rye (<Emphasis Type="Italic">Secale cereale</Emphasis> L.) <Emphasis Type="Italic">Alt4</Emphasis> aluminium (aluminum) tolerance locus using a whole-genome BAC library of rye cv. Blanco

Rye is a diploid crop species with many outstanding qualities, and is important as a source of new traits for wheat and triticale improvement. Rye is highly tolerant of aluminum (Al) toxicity, and possesses a complex structure at the Alt4 Al tolerance locus not found at the corresponding locus in wheat. Here we describe a BAC library of rye cv. Blanco, representing a valuable resource for rye molecular genetic studies, and assess the library’s suitability for investigating Al tolerance genes. The library provides 6 × genome coverage of the 8.1 Gb rye genome, has an average insert size of 131 kb, and contains only ~2% of empty or organelle-derived clones. Genetic analysis attributed the Al tolerance of Blanco to the Alt4 locus on the short arm of chromosome 7R, and revealed the presence of multiple allelic variants (haplotypes) of the Alt4 locus in the BAC library. BAC clones containing ALMT1 gene clusters from several Alt4 haplotypes were identified, and will provide useful starting points for exploring the basis for the structural variability and functional specialization of ALMT1 genes at this locus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome organization of Bowman–Birk inhibitor in common bean (Phaseolus vulgaris L.)

A BAC library from common bean has been used in order to isolate the entire multigene Bowman–Birk serine protease inhibitor family and to study its genome organization. Using a previously isolated trypsin/chymotrypsin inhibitor nucleotide sequence as probe, two positive BAC clones were identified. The P2B8 BAC clone, of about 135 kbp and containing the complete BBI family, was chosen and partially sequenced. Our results confirm that a small multigene family codes for three double-headed inhibitors named: tc-BBI-1, tc-BBI-2 and et-BBI. They contain the binding loop trypsin/chymotrypsin (tc-BBI-1 and tc-BBI-2) and the elastase/trypsin one (et-BBI), respectively. Genes coding for tc-BBI-1 and et-BBI, were found to be very close to each other and arranged in a head to head fashion. Southern blot hybridisation on genomic DNA digested with PstI enzyme suggests that all three genes are present in a fragment of 19 kbp. Northern blot analyses on RNA isolated from various common bean organs showed that the expression of tc-BBI-1 and et-BBI was restricted to the developing cotyledons. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic and physical mapping of a high recombination region on chromosome 7H(1) in barley

Approaches utilizing microlinearity between related species allow for the identification of syntenous regions and orthologous genes. Within the barley Chromosome 7H(1) is a region of high recombination flanked by molecular markers cMWG703 and MWG836. We present the constructed physical contigs linked to molecular markers across this region using bacterial artificial chromosomes (BAC) from the cultivar Morex. Barley expressed sequence tags (EST), identified by homology to rice chromosome 6 between the rice molecular markers C425A and S1434, corresponded to the barley syntenous region of Chromosome 7H(1) Bins 2–5 between molecular markers cMWG703-MWG836. Two hundred and thirteen ESTs were genetically mapped yielding 267 loci of which 101 were within the target high recombination region while 166 loci mapped elsewhere. The 101 loci were joined by 43 other genetic markers resulting in a highly saturated genetic map. In order to develop a physical map of the region, ESTs and all other molecular markers were used to identify Morex BAC clones. Seventy-four BAC contigs were formed containing 2–102 clones each with an average of 19 and a median of 13 BAC clones per contig. Comparison of the BAC contigs, generated here, with the Barley Physical Mapping Database contigs, resulted in additional overlaps and a reduction of the contig number to 56. Within cMWG703-MWG836 are 24 agriculturally important traits including the seedling spot blotch resistance locus, Rcs5. Genetic and physical analysis of this region and comparison to rice indicated an inversion distal of the Rcs5 locus. Three BAC clone contigs spanning the Rcs5 locus were identified. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The involvement of <Emphasis Type="Italic">engrailed</Emphasis> and <Emphasis Type="Italic">wingless</Emphasis> during segmentation in the onychophoran <Emphasis Type="Italic">Euperipatoides kanangrensis</Emphasis> (Peripatopsidae: Onychophora) (Reid 1996)

As the putative sister group to the arthropods, onychophorans can provide insight into ancestral developmental mechanisms in the panarthropod clade. Here, we examine the expression during segmentation of orthologues of wingless (Wnt1) and engrailed, two genes that play a key role in defining segment boundaries in Drosophila and that appear to play a role in segmentation in many other arthropods. Both are expressed in segmentally reiterated stripes in all forming segments except the first (brain) segment, which only shows an engrailed stripe. Engrailed is expressed before segments are morphologically visible and is expressed in both mesoderm and ectoderm. Segmental wingless expression is not detectable until after mesodermal somites are clearly distinct. Early engrailed expression lies in and extends to both sides of the furrow that first demarcates segments in the ectoderm, but is largely restricted to the posterior part of somites. Wingless expression lies immediately anterior to engrailed expression, as it does in many arthropods, but there is no precise cellular boundary between the two expression domains analogous to the overt parasegment boundary seen in Drosophila. Engrailed stripes extend along the posterior part of each limb bud, including the antenna, while wingless is restricted to the distal tip of the limbs and the neurectoderm basal to the limbs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A fine physical map of the rice chromosome 5

A fine physical map of the rice (Oryza sativa spp. Japonica var. Nipponbare) chromosome 5 with bacterial artificial chromosome (BAC) and PI-derived artificial chromosome (PAC) clones was constructed through integration of 280 sequenced BAC/PAC clones and 232 sequence tagged site/expressed sequence tag markers with the use of fingerprinted contig data of the Nipponbare genome. This map consists of five contigs covering 99% of the estimated chromosome size (30.08 Mb). The four physical gaps were estimated at 30 and 20 kb for gaps 1–3 and gap 4, respectively. We have submitted 42.2-Mb sequences with 29.8 Mb of nonoverlapping sequences to public databases. BAC clones corresponding to telomere and centromere regions were confirmed by BAC-fluorescence in situ hybridization (FISH) on a pachytene chromosome. The genetically centromeric region at 54.6 cM was covered by a minimum tiling path spanning 2.1 Mb with no physical gaps. The precise position of the centromere was revealed by using three overlapping BAC/PACs for ~150 kb. In addition, FISH results revealed uneven chromatin condensation around the centromeric region at the pachytene stage. This map is of use for positional cloning and further characterization of the rice functional genomics. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. Chia-Hsiung Cheng and Mei-Chu Chung have equal contributions.     

BAC transgenic zebrafish reveal hypothalamic enhancer activity around obesity associated SNP rs9939609 within the human FTO gene

Single Nucleotide Polymorphisms in FTO intron 1 have been associated with obesity risk, leading to the hypothesis that FTO is the obesity‐related gene. However, other studies have shown that the FTO gene is part of the regulatory domain of the neighboring IRX3 gene and that enhancers in FTO intron 1 regulate IRX3. While Irx3 activity was shown to be necessary in the hypothalamus for the metabolic function of Irx3 in mouse, no enhancers with hypothalamic activity have been demonstrated in the risk‐associated region within FTO. In order to identify potential enhancers at the human FTO locus in vivo, we tested regulatory activity in FTO intron 1 using BAC transgenesis in zebrafish. A minimal gata2 promoter‐GFP cassette was inserted 1.3 kb upstream of the obesity associated SNP rs9939609 in a human FTO BAC plasmid. In addition to the previously identified expression domains in notochord and kidney, human FTO BAC:GFP transgenic zebrafish larvae expressed GFP in the ventral posterior tuberculum, the posterior hypothalamus and the anterior brainstem, which are also expression domains of zebrafish irx3a. In contrast, an in‐frame insertion of a GFP cassette at the FTO start codon resulted in weak ubiquitous GFP expression indicating that the promoter of FTO does likely not react to enhancers located in the obesity risk‐associated region. genesis 53:640–651, 2015. © 2015 The Authors. genesis Published by Wiley Periodicals, Inc.     

A 4D-microscopic analysis of the germ band in the isopod crustacean Porcellio scaber (Malacostraca, Peracarida)—developmental and phylogenetic implications

Malacostracan crustaceans have evolved a conserved stereotyped cell division pattern in the post-naupliar germ band. This cleavage pattern is unique in arthropods investigated so far, and allows a combined analysis of gene expression and cell lineage during segmentation and organ development at the level of individual cells. To investigate the cell lineage in the germ band of the isopod Porcellio scaber, we used a 4D-microscopy system, which enables us to analyse every cell event in the living embryo. The study was combined with the analysis of the expression of the gene engrailed (en) at different stages of germ band formation. Our findings confirm the results of earlier investigations of the cell division pattern in the posterior part of the isopod germ band. Furthermore, we can show that in the anterior region, in contrast to the posterior part, cleavage directions are variable and cell sorting takes place—similar to other arthropod germ bands. Additionally, the gene expression pattern of en in this region is not as regular as in the post-naupliar germ band, and only later becomes regulated into its characteristic stripe pattern. The comparison of the cell lineage of P. scaber with that of other malacostracan crustaceans shows an enhancement in the velocity of cell divisions relative to the arrangement of these cells in rows in the isopod germ band. The striking similarity of the formation of the genealogical units in the anterior part suggests a sister group relationship between the peracarid taxa Tanaidacea and Isopoda.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Cleavage pattern,gastrulation, and neurulation in the appendicularian, <Emphasis Type="Italic">Oikopleura dioica</Emphasis>

The appendicularian, Oikopleura dioica is a chordate. Its life cycle is extremely short—approximately 5 days—and its tadpole shape with a beating tail is retained throughout entire life. The tadpole hatches after 3 h of development at 20°C. Here, we describe the cleavage pattern and morphogenetic cell movements during gastrulation and neurulation. Cleavage showed an invariant pattern. It is basically bilateral but also shows various minor left–right asymmetries starting from the four-cell stage. We observed two rounds of unequal cleavage of the posterior-vegetal B-line cells at the posterior pole. The nature of the unequal cleavages is reminiscent of those in ascidian embryos and suggests the presence of a centrosome-attracting body, a special subcellular structure at the posterior pole. The representation of the cell division pattern in this report will aid the identification of each cell, a prerequisite for clarifying the gene expression patterns in early embryos. Gastrulation started as early as the 32-cell stage and progressed in three phases. By the end of the second phase at the 64-cell stage, every vegetal cell had ingressed into the embryo, and animal cells had covered the entire embryo by epiboly. There was no archenteron formation. In the anterior region, eight A-line cells were aligned as a 2 × 4 array along the anterior–posterior axis and become internalized during the 64-cell stage. This process was considered to correspond to neurulation. The simple and accelerated development of Oikopleura, nevertheless giving rise to a conserved chordate body plan, is advantageous for studying developmental mechanisms using molecular and genetic approaches and makes this animal the simplest model organism in the phylum Chordata. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Fine mapping and DNA fiber FISH analysis locates the tobamovirus resistance gene <Emphasis Type="Italic">L</Emphasis><Superscript>3</Superscript> of <Emphasis Type="Italic">Capsicum chinense</Emphasis> in a 400-kb region of R-like genes cluster embedded in highly repetitive sequences

The tobamovirus resistance gene L ³ of Capsicum chinense was mapped using an intra-specific F2 population (2,016 individuals) of Capsicum annuum cultivars, into one of which had been introduced the C. chinense L ³ gene, and an inter-specific F2 population (3,391 individuals) between C. chinense and Capsicum frutescence. Analysis of a BAC library with an AFLP marker closely linked to L ³-resistance revealed the presence of homologs of the tomato disease resistance gene I2. Partial or full-length coding sequences were cloned by degenerate PCR from 35 different pepper I2 homologs and 17 genetic markers were generated in the inter-specific combination. The L ³ gene was mapped between I2 homolog marker IH1-04 and BAC-end marker 189D23M, and located within a region encompassing two different BAC contigs consisting of four and one clones, respectively. DNA fiber FISH analysis revealed that these two contigs are separated from each other by about 30 kb. DNA fiber FISH results and Southern blotting of the BAC clones suggested that the L ³ locus-containing region is rich in highly repetitive sequences. Southern blot analysis indicated that the two BAC contigs contain more than ten copies of the I2 homologs. In contrast to the inter-specific F2 population, no recombinant progeny were identified to have a crossover point within two BAC contigs consisting of seven and two clones in the intra-specific F2 population. Moreover, distribution of the crossover points differed between the two populations, suggesting linkage disequilibrium in the region containing the L locus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. R. Tomita and J. Murai contributed equally to this work.     

Mapping of genome-wide resistance gene analogs (RGAs) in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Isolation and mapping of genome-wide resistance (R) gene analogs (RGAs) is of importance in identifying candidate(s) for a particular resistance gene/QTL. Here we reported our result in mapping totally 228 genome-wide RGAs in maize. By developing RGA-tagged markers and subsequent genotyping a population consisting of 294 recombinant inbred lines (RILs), 67 RGAs were genetically mapped on maize genome. Meanwhile, in silico mapping was conducted to anchor 113 RGAs by comparing all 228 RGAs to those anchored EST and BAC/BAC-end sequences via tblastx search (E-value < 10⁻²⁰). All RGAs from different mapping efforts were integrated into the existing SSR linkage map. After accounting for redundancy, the resultant RGA linkage map was composed of 153 RGAs that were mapped onto 172 loci on maize genome, and the mapped RGAs accounted for approximate three quarters of the genome-wide RGAs in maize. The extensive co-localizations were observed between mapped RGAs and resistance gene/QTL loci, implying the usefulness of this RGA linkage map in R gene cloning via candidate gene approach. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Wenkai Xiao, Jing Zhao and Shengci Fan have contributed equally to this research.     

Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana

We sequenced five BAC clones of Brassica oleracea doubled haploid ‘Early Big' broccoli containing major genes in the aliphatic glucosinolate pathway, and comparatively analyzed them with similar sequences in A. thaliana and B. rapa. Additionally, we included in the analysis published sequences from three other B. oleracea BAC clones and a contig of this species corresponding to segments in A. thaliana chromosomes IV and V. A total of 2,946 kb of B. oleracea, 1,069 kb of B. rapa sequence and 2,607 kb of A. thaliana sequence were compared and analyzed. We found conserved collinearity for gene order and content restricted to specific chromosomal segments, but breaks in collinearity were frequent resulting in gene absence likely not due to gene loss but rearrangements. B. oleracea has the lowest gene density of the three species, followed by B. rapa. The genome expansion of the Brassica species, B. oleracea in particular, is due to larger introns and gene spacers resulting from frequent insertion of DNA transposons and retrotransposons. These findings are discussed in relation to the possible origin and evolution of the Brassica genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The endothelial-specific regulatory mutation, <Emphasis Type="Italic">Mvwf1</Emphasis>, is a common mouse founder allele

Mvwf1 is a cis-regulatory mutation previously identified in the RIIIS/J mouse strain that causes a unique tissue-specific switch in the expression of an N-acetylgalactosaminyltransferase, B4GALNT2, from intestinal epithelium to vascular endothelium. Vascular B4galnt2 expression results in aberrant glycosylation of von Willebrand Factor (VWF) and accelerated VWF clearance from plasma. We now report that 13 inbred mouse strains share the Mvwf1 tissue-specific switch and low VWF phenotype, including five wild-derived strains. Genomic sequencing identified a highly conserved 97-kb Mvwf1 haplotype block shared by these strains that encompasses a 30-kb region of high nucleotide sequence divergence from C57BL6/J flanking B4galnt2 exon 1. The analysis of a series of bacterial artificial chromosome (BAC) transgenes containing B4galnt2 derived from the RIIIS/J or C57BL6/J inbred mouse strains demonstrates that the corresponding sequences are sufficient to confer the vessel (RIIIS/J) or intestine (C57BL6/J)-specific expression patterns. Taken together, our data suggest that the region responsible for the Mvwf1 regulatory switch lies within an approximately 30-kb genomic interval upstream of the B4galnt2 gene. The observation that Mvwf1 is present in multiple wild-derived strains suggests that this locus may be retained in wild mouse populations due to positive selection. Similar selective pressures could contribute to the high prevalence of von Willebrand disease in humans. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

<Emphasis Type="Italic">kitb</Emphasis>, a second zebrafish ortholog of mouse <Emphasis Type="Italic">Kit</Emphasis>

The large numbers of duplicated pairs of genes in zebrafish compared to their mammalian counterparts has lead to the notion that expression of zebrafish co-orthologous pairs in some cases can together describe the expression of their mammalian counterpart. Here, we explore this notion by identification and analysis of a second zebrafish ortholog of the mammalian Kit receptor tyrosine kinase (kitb). We show that in embryos, kitb is expressed in a non-overlapping pattern to that of kita, in the anterior ventral mesoderm, Rohon-beardRohon–Beard neurons, the otic vesicle, and trigeminal ganglia. The expression pattern of kita and kitb in zebrafish together approximates that of Kit in mouse, with the exception that neither zebrafish kit gene is expressed in primordial germ cells, a site of kit expression in the mouse embryo. In addition, zebrafish kita is expressed in a site of zebrafish primitive hematopoiesis but not required for blood development, and we fail to detect kitb expression in sites of zebrafish hematopoiesis. Thus, the expression and function of zebrafish kit genes cannot be described as a simple partition of the expression and function of mouse Kit. We discuss the possibility that these unaccounted for expression domains and functions are derived from more ancestral gene duplications and partitioning instead of the relatively recent teleost teleost-specific duplication. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

Evolution, expression and effectiveness in a cluster of novel bovine β-defensins

The anti-microbial peptides β-defensins constitute a large family of innate immune effector molecules, conserved across a wide species range. In this paper, we describe a systematic search of the sequenced bovine genome to characterise this extensive gene family in Bos taurus, providing an insight into the pattern of conservation of β-defensin genes between species. We have sequenced a sub-set of these newly discovered bovine β-defensin genes and also report expression data for these genes across a range of tissues. We have synthesised the peptide product of one of these genes, bovine β-defensin 123, and found it to be a potent inhibitor of several pathogenic microbes, particularly Escherichia coli and Listeria monocytogenes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.