An improved approach for construction of bacterial artificial chromosome libraries

Independent expansion, selection and hypervariability of the TBC1D3 gene family in humans

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on chromosome 17. We find that most humans vary along two TBC1D3 clusters where human haplotypes are highly variable in copy number, differing by as many as 20 copies, and structure (structural heterozygosity 90%). We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Lastly, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL. These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

Autophosphorylation of the KaiC-like protein ArlH inhibits oligomerization and interaction with ArlI, the motor ATPase of the archaellum

Motile archaea are propelled by the archaellum, whose motor complex consists of the membrane protein ArlJ, the ATPase ArlI, and the ATP-binding protein ArlH. Despite its essential function and the existence of structural and biochemical data on ArlH, the role of ArlH in archaellum assembly and function remains elusive. ArlH is a structural homolog of KaiC, the central component of the cyanobacterial circadian clock. Since autophosphorylation and dephosphorylation of KaiC are central properties for the function of KaiC, we asked whether autophosphorylation is also a property of ArlH proteins. We observed that both ArlH from the euryarchaeon Pyrococcus furiosus (PfArlH) and from the crenarchaeon Sulfolobus acidocaldarius (SaArlH) have autophosphorylation activity. Using a combination of single-molecule fluorescence measurements and biochemical assays, we show that autophosphorylation of ArlH is closely linked to its oligomeric state when bound to hexameric ArlI. These experiments also strongly suggest that ArlH is a hexamer in its ArlI-bound state. Mutagenesis of the putative catalytic residue (Glu-57 in SaArlH) in ArlH results in a reduced autophosphorylation activity and abolished archaellation and motility in S. acidocaldarius, indicating that optimum phosphorylation activity of ArlH is essential for archaellation and motility.

High-Throughput Screening for CEBPD-Modulating Compounds in THP-1-Derived Reporter Macrophages Identifies Anti-Inflammatory HDAC and BET Inhibitors

This study aimed to identify alternative anti-inflammatory compounds that modulate the activity of a relevant transcription factor, CCAAT/enhancer binding protein delta (C/EBPδ). C/EBPδ is a master regulator of inflammatory responses in macrophages (Mϕ) and is mainly regulated at the level of CEBPD gene transcription initiation. To screen for CEBPD-modulating compounds, we generated a THP-1-derived reporter cell line stably expressing secreted alkaline phosphatase (SEAP) under control of the defined CEBPD promoter (CEBPD::SEAP). A high-throughput screening of LOPAC^®1280 and ENZO^®774 libraries on LPS- and IFN-γ-activated THP-1 reporter Mϕ identified four epigenetically active hits: two bromodomain and extraterminal domain (BET) inhibitors, I-BET151 and Ro 11-1464, as well as two histone deacetylase (HDAC) inhibitors, SAHA and TSA. All four hits markedly and reproducibly upregulated SEAP secretion and CEBPD::SEAP mRNA expression, confirming screening assay reliability. Whereas BET inhibitors also upregulated the mRNA expression of the endogenous CEBPD, HDAC inhibitors completely abolished it. All hits displayed anti-inflammatory activity through the suppression of IL-6 and CCL2 gene expression. However, I-BET151 and HDAC inhibitors simultaneously upregulated the mRNA expression of pro-inflammatory IL-1ß. The modulation of CEBPD gene expression shown in this study contributes to our understanding of inflammatory responses in Mϕ and may offer an approach to therapy for inflammation-driven disorders.

Cytogenetic Characterization of a Small Evolutionary Rearrangement Involving Chromosomes BTA21 and OAR18

Both cattle (Bos taurus) and sheep (Ovis aries) belong to the Bovidae family but to different subfamilies, Bovinae and Caprinae, respectively. From a chromosomal point of view, apart from the already known centric fusions (that occurred during the evolutionary process in the Bovidae family) and the small differences in the chromosome classification, the 2 karyotypes are very similar in banding patterns. In this study, the combination of bioinformatics techniques and physical mapping of DNA markers enabled the identification of a micro-rearrangement, a small inversion involving bovine chromosome 21 (BTA21) and the corresponding sheep chromosome 18 (OAR18). The aim of this study was the cytogenetic characterization of this difference in genomic assemblies between cattle and sheep in this single chromosome region. To verify the inversion in FISH experiments, we used the BACs 442H08 and 222H03 from the INRA library and BACs 134H22 and 436P08 from the sheep-specific CHORI library. The results confirmed the presence of the inverted fragment in sheep compared to the cattle genome. Genomic rearrangements may have consequences depending on their influence on gene activity, but in this case no gene or transcribed DNA portion seemed to be involved. In conclusion, we showed for the first time, concerning autosomes, that besides the already known centric fusions also other differences exist between the bovine and sheep karyotypes. Furthermore, we demonstrated that the combination of a bioinformatics approach and physical mapping is a valid tool for the identification of currently unknown rearrangements between related species.

Prospect and challenge of detecting dynamic gene copy number increases in stem cells by whole genome sequencing

Abstract

Gene amplification is an evolutionarily well-conserved and highly efficient mechanism to increase the amount of specific proteins. In humans, gene amplification is a hallmark of cancer and has recently been found during stem cell differentiation. Amplifications in stem cells are restricted to specific tissue areas and time windows, rendering their detection difficult. Here, we report on the performance of deep WGS sequencing (average 82-fold depth of coverage) on the BGISEQ with nanoball technology to detect amplifications in human mesenchymal and neural stem cells. As reference technology, we applied array-based comparative genomic hybridization (aCGH), fluorescence in situ hybridization (FISH), and qPCR. Using different in silico strategies for amplification detection, we analyzed the potential of WGS for amplification detection. Our results provide evidence that WGS accurately identifies changes of the copy number profiles in human stem cell differentiation. However, the identified changes are not in all cases consistent between WGS and aCGH. The results between WGS and the validation by qPCR were concordant in 83.3% of all tested 36 cases. In sum, both genome-wide techniques, aCGH and WGS, have unique advantages and specific challenges, calling for locus-specific confirmation by the low-throughput approaches qPCR or FISH.

Key messages

WGS allows for the identification of dynamic copy number changes in human stem cells.

Less stringent threshold setting is crucial for detection of copy number increase.

Broad knowledge of dynamic copy number is pivotal to estimate stem cell capabilities.

Electronic supplementary material

The online version of this article (10.1007/s00109-019-01792-y) contains supplementary material, which is available to authorized users.

Complete sequencing of the SMN2 gene in SMA patients detects SMN gene deletion junctions and variants in SMN2 that modify the SMA phenotype

Spinal muscular atrophy (SMA) is a progressive motor neuron disease caused by loss or mutation of the survival motor neuron 1 (SMN1) gene and retention of SMN2. We performed targeted capture and sequencing of the SMN2, CFTR, and PLS3 genes in 217 SMA patients. We identified a 6.3 kilobase deletion that occurred in both SMN1 and SMN2 (SMN1/2) and removed exons 7 and 8. The deletion junction was flanked by a 21 bp repeat that occurred 15 times in the SMN1/2 gene. We screened for its presence in 466 individuals with the known SMN1 and SMN2 copy numbers. In individuals with 1 SMN1 and 0 SMN2 copies, the deletion occurred in 63% of cases. We modeled the deletion junction frequency and determined that the deletion occurred in both SMN1 and SMN2. We have identified the first deletion junction where the deletion removes exons 7 and 8 of SMN1/2. As it occurred in SMN1, it is a pathogenic mutation. We called variants in the PLS3 and SMN2 genes, and tested for association with mild or severe exception patients. The variants A-44G, A-549G, and C-1897T in intron 6 of SMN2 were significantly associated with mild exception patients, but no PLS3 variants correlated with severity. The variants occurred in 14 out of 58 of our mild exception patients, indicating that mild exception patients with an intact SMN2 gene and without modifying variants occur. This sample set can be used in the association analysis of candidate genes outside of SMN2 that modify the SMA phenotype.

Sequencing through thick and thin: Historiographical and philosophical implications

DNA sequencing has been characterised by scholars and life scientists as an example of 'big', 'fast' and 'automated' science in biology. This paper argues, however, that these characterisations are a product of a particular interpretation of what sequencing is, what I call 'thin sequencing'. The 'thin sequencing' perspective focuses on the determination of the order of bases in a particular stretch of DNA. Based upon my research on the pig genome mapping and sequencing projects, I provide an alternative 'thick sequencing' perspective, which also includes a number of practices that enable the sequence to travel across and be used in wider communities. If we take sequencing in the thin manner to be an event demarcated by the determination of sequences in automated sequencing machines and computers, this has consequences for the historical analysis of sequencing projects, as it focuses attention on those parts of the work of sequencing that are more centralised, fast (and accelerating) and automated. I argue instead that sequencing can be interpreted as a more open-ended process including activities such as the generation of a minimum tile path or annotation, and detail the historiographical and philosophical consequences of this move.

Heterologous expression-facilitated natural products' discovery in actinomycetes

Actinomycetes produce many of the drugs essential for human and animal health as well as crop protection. Genome sequencing projects launched over the past two decades reveal dozens of cryptic natural product biosynthetic gene clusters in each actinomycete genome that are not expressed under regular laboratory conditions. This so-called 'chemical dark matter' represents a potentially rich untapped resource for drug discovery in the genomic era. Through improved understanding of natural product biosynthetic logic coupled with the development of bioinformatic and genetic tools, we are increasingly able to access this 'dark matter' using a wide variety of strategies with downstream potential application in drug development. In this review, we discuss recent research progress in the field of cloning of natural product biosynthetic gene clusters and their heterologous expression in validating the potential of this methodology to drive next-generation drug discovery.

Related Endogenous Retrovirus-K Elements Harbor Distinct Protease Active Site Motifs

Background: Endogenous retrovirus-K is a group of related genomic elements descending from retroviral infections in human ancestors. HML2 is the clade of these viruses which contains the most intact provirus copies. These elements can be transcribed and translated in healthy and diseased tissues, and some of them produce active retroviral enzymes, such as protease. Retroviral gene products, including protease, contribute to illness in exogenous retroviral infections. There are ongoing efforts to test anti-retroviral regimens against endogenous retroviruses. Herein, we examine the potential activity and diversity of human endogenous retrovirus-K proteases, and their potential for impact on immunity and human disease.

Results: Sequences similar to the endogenous retrovirus-K HML2 protease and reverse transcriptase were identified in the human genome, classified by phylogenetic inference and compared to Repbase reference sequences. The topologies of trees inferred from protease and reverse transcriptase sequences were similar and agreed with the classification using reference sequences. Surprisingly, only 62/480 protease sequences identified by BLAST were classified as HML2; the remainder were classified as other HML groups, with the majority (216) classified as HML3. Variation in functionally significant protease motifs was explored, and two major active site variants were identified – the DTGAD variant is common in all groups, but the DTGVD motif appears limited to HML3, HML5, and HML6. Furthermore, distinct RNA expression patterns of protease variants are seen in disease states, such as amyotrophic lateral sclerosis, breast cancer, and prostate cancer.

Conclusion: Transcribed ERVK proteases exhibit a diversity which could impact immunity and inhibitor-based treatments, and these facets should be considered when designing therapeutic regimens.

A new targeted capture method using bacterial artificial chromosome (BAC) libraries as baits for sequencing relatively large genes

To analyze a specific genome region using next-generation sequencing technologies, the enrichment of DNA libraries with targeted capture methods has been standardized. For enrichment of mitochondrial genome, a previous study developed an original targeted capture method that use baits constructed from long-range polymerase chain reaction (PCR) amplicons, common laboratory reagents, and equipment. In this study, a new targeted capture method is presented, that of bacterial artificial chromosome (BAC) double capture (BDC), modifying the previous method, but using BAC libraries as baits for sequencing a relatively large gene. We applied the BDC approach for the 214 kb autosomal region, ring finger protein 213, which is the susceptibility gene of moyamoya disease (MMD). To evaluate the reliability of BDC, cost and data quality were compared with those of a commercial kit. While the ratio of duplicate reads was higher, the cost was less than that of the commercial kit. The data quality was sufficiently the same as that of the kit. Thus, BDC can be an easy, low-cost, and useful method for analyzing individual genome regions with substantial length.

Specific amplifications and copy number decreases during human neural stem cells differentiation towards astrocytes, neurons and oligodendrocytes

There is growing evidence that gene amplifications are an attribute of normal cells during development and differentiation. During neural progenitor cell differentiation half of the genome is involved in amplification process. To answer the question how specific amplifications occur at different stages and in different lineages of differentiation we analyzed the genes CDK4, MDM2, EGFR, GINS2, GFAP, TP53, DDB1 and MDM4 in human neural stem cells that were induced to differentiate towards astrocytes, neurons and oligodendrocytes. We found specific amplification pattern for each of the eight analyzed genes both in undifferentiated neural stem and progenitor cells and in cells that were induced for differentiation. Different amplification patterns were also found between adherently grown neural stem cells and cells that were grown as spheres. The most frequently amplified genes were MDM2 and CDK4 with the latter amplified in all three lineages at all analyzed stages. Amplification of the analyzed genes was also found in four glioma stem-like cells. The combined amplification data of stem cells and of tumor stem cells can help to define cell populations at the origin of the tumor. Furthermore, we detected a decrease of gene copies at specific differentiation stages most frequently for MDM4. This study shows specific amplification pattern in defined stem cell populations within specific time windows during differentiation processes indicating that amplifications occur in an orderly sequence during the differentiation of human neural stem and progenitor cells.

Gene amplification in mesenchymal stem cells and during differentiation towards adipocytes or osteoblasts

Gene amplifications are an attribute of tumor cells and have for long time been overlooked in normal cells. A growing number of investigations describe gene amplifications in normal mammalian cells during development and differentiation. Possibly, tumor cells have rescued the gene amplification mechanism as a physiological attribute of stem cells. Here, we investigated human mesenchymal stem cells (hMSCs) for gene amplification using array-CGH, single cell fluorescence in situ hybridization and qPCR. Gene amplifications were detected in mesenchymal stem cells and in mesenchymal stem cells during differentiation towards adipocytes and osteoblasts. Undifferentiated hMSCs harbor 12 amplified chromosomal regions, hMSCs that differentiated towards adipocytes 18 amplified chromosome regions, and hMSCs that differentiate towards osteoblasts 19 amplified regions. Specifically, hMSCs that differentiated towards adipocytes or osteoblasts harbor CDK4 and MDM2 amplifications both of which frequently occur in osteosarcoma and liposarcoma that are both of same cell origin. Beside the amplifications, we identified 36 under-replicated regions in undifferentiated and in differentiating hMSC cells.

Sub genome anchored physical frameworks of the allotetraploid Upland cotton (Gossypium hirsutum L.) genome, and an approach toward reference-grade assemblies of polyploids

Like those of many agricultural crops, the cultivated cotton is an allotetraploid and has a large genome (~2.5 gigabase pairs). The two sub genomes, A and D, are highly similar but unequally sized and repeat-rich, which pose significant challenges for accurate genome reconstruction using standard approaches. Here we report the development of BAC libraries, sub genome specific physical maps, and a new-generation sequencing approach that will lead to a reference-grade genome assembly for Upland cotton. Three BAC libraries were constructed, fingerprinted, and integrated with BAC-end sequences (BES) to produce a de novo whole-genome physical map. The BAC map was partitioned by sub genomes through alignment to the diploid progenitor D-genome reference sequence with densely spaced BES anchor points and computational filtering. The physical maps were validated with FISH and genetic mapping of SNP markers derived from BES. Two pairs of homeologous chromosomes, A11/D11 and A12/D12, were used to assess multiplex sequencing approaches for completeness and scalability. The results represent the first sub genome anchored physical maps of Upland cotton, and a new-generation approach to the whole-genome sequencing, which will lead to the reference-grade assembly of allopolyploid cotton and serve as a general strategy for sequencing other polyploid 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.

β-mercaptoethanol assists efficient construction of sperm bacterial artificial chromosome library

Bacterial artificial chromosome (BAC) library plays a critical role in the strategic research in genomics. Sperm is known as a good source for BAC library construction. However, preparation of intact DNA from the highly condensed sperm nuclei is not easy. Here we developed and validated an efficient DNA extraction strategy for BAC library construction from sperm embedded in agarose plugs. The protocol used a combination of lauroylsarcosine, proteinase K and β-mercaptoethanol (a reducing agent of nucleus) In comparison with the normal protocol without reducing agents, β-mercaptoethanol released high-molecular-weight DNA from the protamines which permit DNA to be packed very densely within the spermatozoan nucleus, without damaging DNA. Extracted DNA by this method was readily digested by restriction enzymes and ideal for BAC library construction.

Transient RUNX1 Expression during Early Mesendodermal Differentiation of hESCs Promotes Epithelial to Mesenchymal Transition through TGFB2 Signaling

The transition of human embryonic stem cells (hESCs) from pluripotency to lineage commitment is not fully understood, and a role for phenotypic transcription factors in the initial stages of hESC differentiation remains to be explored. From a screen of candidate factors, we found that RUNX1 is selectively and transiently upregulated early in hESC differentiation to mesendodermal lineages. Transcriptome profiling and functional analyses upon RUNX1 depletion established a role for RUNX1 in promoting cell motility. In parallel, we discovered a loss of repression for several epithelial genes, indicating that loss of RUNX1 impaired an epithelial to mesenchymal transition during differentiation. Cell biological and biochemical approaches revealed that RUNX1 depletion specifically compromised TGFB2 signaling. Both the decrease in motility and deregulated epithelial marker expression upon RUNX1 depletion were rescued by reintroduction of TGFB2, but not TGFB1. These findings identify roles for RUNX1-TGFB2 signaling in early events of mesendodermal lineage commitment.

•

RUNX1 is transiently upregulated during early mesendoderm differentiation of hESCs

•

RUNX1 promotes motility and the EMT process during mesendodermal differentiation

•

RUNX1 knockdown specifically inhibits TGFB2 signaling

•

Reintroduction of TGFB2, but not TGFB1, rescues the phenotype of RUNX1 depletion

Focus on Kir7.1: physiology and channelopathy

Genetic studies have linked alterations in Kir7.1 channel to diverse pathologies. We summarize functional relevance of Kir7.1 channel in retinal pigment epithelium (RPE), regulation of channel function by various cytoplasmic metabolites, and mutations that cause channelopathies. At the apical membrane of RPE, K(+) channels contribute to subretinal K(+) homeostasis and support Na(+)/K(+) pump and Na(+)-K(+)-2Cl(-) cotransporter function by providing a pathway for K(+) secretion. Electrophysiological studies have established that barium- and cesium-sensitive inwardly rectifying K(+) (Kir) channels make up a major component of the RPE apical membrane K(+) conductance. Native human RPE expresses transcripts for Kir1.1, Kir2.1, Kir2.2, Kir3.1, Kir3.4, Kir4.2, and Kir6.1, albeit at levels at least 50-fold lower than Kir7.1. Kir7.1 is structurally similar to other Kir channels, consisting of 2 trans-membrane domains, a pore-forming loop that contains the selectivity filter, and 2 cytoplasmic polar tails. Within the cytoplasmic structure, clusters of amino acid sequences form regulatory domains that interact with cellular metabolites and control the opening and closing of the channel. Recent evidence indicated that intrinsic sequence motifs present in Kir7.1 control surface expression. Mutant Kir7.1 channels are associated with inherited eye pathologies such as Snowflake Vitreoretinal Degeneration (SVD) and Lebers Congenital Amaurosis (LCA16). Based on the current evidence, mutations implicated in channelopathies have the potential to be used for genetic testing to diagnose blindness due to Kir7.1.

Construction and characterization of three wheat bacterial artificial chromosome libraries

We have constructed three bacterial artificial chromosome (BAC) libraries of wheat cultivar Triticum aestivum Wangshuibai, germplasms T. monococcum TA2026 and TA2033. A total of 1,233,792,170,880 and 263,040 clones were picked and arrayed in 384-well plates. On the basis of genome sizes of 16.8 Gb for hexaploid wheat and 5.6 Gb for diploid wheat, the three libraries represented 9.05-, 2.60-, and 3.71-fold coverage of the haploid genomes, respectively. An improved descending pooling system for BAC libraries screening was established. This improved strategy can save 80% of the time and 68% of polymerase chain reaction (PCR) with the same successful rate as the universal 6D pooling strategy.

The first report of a Pelecaniformes defensin cluster: characterization of β-defensin genes in the crested ibis based on BAC libraries

Defensins play a key role in the innate immunity of various organisms. Detailed genomic studies of the defensin cluster have only been reported in a limited number of birds. Herein, we present the first characterization of defensins in a Pelecaniformes species, the crested ibis (Nipponia nippon), which is one of the most endangered birds in the world. We constructed bacterial artificial chromosome libraries, including a 4D-PCR library and a reverse-4D library, which provide at least 40 equivalents of this rare bird's genome. A cluster including 14 β-defensin loci within 129 kb was assigned to chromosome 3 by FISH, and one gene duplication of AvBD1 was found. The ibis defensin genes are characterized by multiform gene organization ranging from two to four exons through extensive exon fusion. Splicing signal variations and alternative splice variants were also found. Comparative analysis of four bird species identified one common and multiple species-specific duplications, which might be associated with high GC content. Evolutionary analysis revealed birth-and-death mode and purifying selection for avian defensin evolution, resulting in different defensin gene numbers among bird species and functional conservation within orthologous genes, respectively. Additionally, we propose various directions for further research on genetic conservation in the crested ibis.

Sphingosine-1-phosphate lyase downregulation promotes colon carcinogenesis through STAT3-activated microRNAs

Growing evidence supports a link between inflammation and cancer; however, mediators of the transition between inflammation and carcinogenesis remain incompletely understood. Sphingosine-1-phosphate (S1P) lyase (SPL) irreversibly degrades the bioactive sphingolipid S1P and is highly expressed in enterocytes but downregulated in colon cancer. Here, we investigated the role of SPL in colitis-associated cancer (CAC). We generated mice with intestinal epithelium-specific Sgpl1 deletion and chemically induced colitis and tumor formation in these animals. Compared with control animals, mice lacking intestinal SPL exhibited greater disease activity, colon shortening, cytokine levels, S1P accumulation, tumors, STAT3 activation, STAT3-activated microRNAs (miRNAs), and suppression of miR-targeted anti-oncogene products. This phenotype was attenuated by STAT3 inhibition. In fibroblasts, silencing SPL promoted tumorigenic transformation through a pathway involving extracellular transport of S1P through S1P transporter spinster homolog 2 (SPNS2), S1P receptor activation, JAK2/STAT3-dependent miR-181b-1 induction, and silencing of miR-181b-1 target cylindromatosis (CYLD). Colon biopsies from patients with inflammatory bowel disease revealed enhanced S1P and STAT3 signaling. In mice with chemical-induced CAC, oral administration of plant-type sphingolipids called sphingadienes increased colonic SPL levels and reduced S1P levels, STAT3 signaling, cytokine levels, and tumorigenesis, indicating that SPL prevents transformation and carcinogenesis. Together, our results suggest that dietary sphingolipids can augment or prevent colon cancer, depending upon whether they are metabolized to S1P or promote S1P metabolism through the actions of SPL.

Recombineering linear BACs

Recombineering is a powerful genetic engineering technique based on homologous recombination that can be used to accurately modify DNA independent of its sequence or size. One novel application of recombineering is the assembly of linear BACs in E. coli that can replicate autonomously as linear plasmids. A circular BAC is inserted with a short telomeric sequence from phage N15, which is subsequently cut and rejoined by the phage protelomerase enzyme to generate a linear BAC with terminal hairpin telomeres. Telomere-capped linear BACs are protected against exonuclease attack both in vitro and in vivo in E. coli cells and can replicate stably. Here we describe step-by-step protocols to linearize any BAC clone by recombineering, including inserting and screening for presence of the N15 telomeric sequence, linearizing BACs in vivo in E. coli, extracting linear BACs, and verifying the presence of hairpin telomere structures. Linear BACs may be useful for functional expression of genomic loci in cells, maintenance of linear viral genomes in their natural conformation, and for constructing innovative artificial chromosome structures for applications in mammalian and plant cells.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

Harnessing mobile genetic elements to explore gene regulation

Sequences that regulate expression of a gene in cis but are located at large distances along the DNA from the gene, as found with most developmentally regulated genes in higher vertebrates, are difficult to identify if those sequences are not conserved across species. Mutating suspected gene-regulatory sequences to alter expression then becomes a hit-or-miss affair. The relaxed specificity of transposon insertions offers an opportunity to develop alternate strategies, to scan in an unbiased manner, pieces of chromosomal DNA cloned in BACs for transcription enhancing elements. This article illustrates how insertions of Tn10 with enhancer-traps into BAC DNA containing the gene, and its germ-line expression in zebrafish, have identified distal regulatory elements functionally. Transposition of Tn10 first introduces the enhancer-trap with a loxP site randomly into BAC DNA. Cre-recombination between the inserted loxP and the loxP endogenous to a BAC-end positions the enhancer-trap to the newly created truncated end of BAC DNA. The procedure generates a library of integration-ready enhancer-trap BACs with progressive truncations from an end in a single experiment. Individual enhancer-trap BACs from the library can be evaluated functionally in zebrafish or mice. Furthermore, the ability to readily alter sequences in a small transposon plasmid containing a regulatory domain of the gene allows re-introduction of altered parts of a BAC back into itself. It serves as a useful strategy to functionally dissect multiple discontinuous regulatory domains of a gene quickly. These methodologies have been successfully used in identifying novel regulatory domains of the Amyloid Precursor Protein (appb) gene in zebrafish, and provided important clues for regulation of the gene in humans.

Finished sequence and assembly of the DUF1220-rich 1q21 region using a haploid human genome

Background

Although the reference human genome sequence was declared finished in 2003, some regions of the genome remain incomplete due to their complex architecture. One such region, 1q21.1-q21.2, is of increasing interest due to its relevance to human disease and evolution. Elucidation of the exact variants behind these associations has been hampered by the repetitive nature of the region and its incomplete assembly. This region also contains 238 of the 270 human DUF1220 protein domains, which are implicated in human brain evolution and neurodevelopment. Additionally, examinations of this protein domain have been challenging due to the incomplete 1q21 build. To address these problems, a single-haplotype hydatidiform mole BAC library (CHORI-17) was used to produce the first complete sequence of the 1q21.1-q21.2 region.

Results

We found and addressed several inaccuracies in the GRCh37sequence of the 1q21 region on large and small scales, including genomic rearrangements and inversions, and incorrect gene copy number estimates and assemblies. The DUF1220-encoding NBPF genes required the most corrections, with 3 genes removed, 2 genes reassigned to the 1p11.2 region, 8 genes requiring assembly corrections for DUF1220 domains (~91 DUF1220 domains were misassigned), and multiple instances of nucleotide changes that reassigned the domain to a different DUF1220 subtype. These corrections resulted in an overall increase in DUF1220 copy number, yielding a haploid total of 289 copies. Approximately 20 of these new DUF1220 copies were the result of a segmental duplication from 1q21.2 to 1p11.2 that included two NBPF genes. Interestingly, this duplication may have been the catalyst for the evolutionarily important human lineage-specific chromosome 1 pericentric inversion.

Conclusions

Through the hydatidiform mole genome sequencing effort, the 1q21.1-q21.2 region is complete and misassemblies involving inter- and intra-region duplications have been resolved. The availability of this single haploid sequence path will aid in the investigation of many genetic diseases linked to 1q21, including several associated with DUF1220 copy number variations. Finally, the corrected sequence identified a recent segmental duplication that added 20 additional DUF1220 copies to the human genome, and may have facilitated the chromosome 1 pericentric inversion that is among the most notable human-specific genomic landmarks.

Construction and characterization of a bacterial artificial chromosome library for the hexaploid wheat line 92R137

For map-based cloning of genes conferring important traits in the hexaploid wheat line 92R137, a bacterial artificial chromosome (BAC) library, including two sublibraries, was constructed using the genomic DNA of 92R137 digested with restriction enzymes HindIII and BamHI. The BAC library was composed of total 765,696 clones, of which 390,144 were from the HindIII digestion and 375,552 from the BamHI digestion. Through pulsed-field gel electrophoresis (PFGE) analysis of 453 clones randomly selected from the HindIII sublibrary and 573 clones from the BamHI sublibrary, the average insert sizes were estimated as 129 and 113 kb, respectively. Thus, the HindIII sublibrary was estimated to have a 3.01-fold coverage and the BamHI sublibrary a 2.53-fold coverage based on the estimated hexaploid wheat genome size of 16,700 Mb. The 765,696 clones were arrayed in 1,994 384-well plates. All clones were also arranged into plate pools and further arranged into 5-dimensional (5D) pools. The probability of identifying a clone corresponding to any wheat DNA sequence (such as gene Yr26 for stripe rust resistance) from the library was estimated to be more than 99.6%. Through polymerase chain reaction screening the 5D pools with Xwe173, a marker tightly linked to Yr26, six BAC clones were successfully obtained. These results demonstrate that the BAC library is a valuable genomic resource for positional cloning of Yr26 and other genes of interest.

Genome sequence of the tsetse fly (Glossina morsitans): vector of African trypanosomiasis

Tsetse flies are the sole vectors of human African trypanosomiasis throughout sub-Saharan Africa. Both sexes of adult tsetse feed exclusively on blood and contribute to disease transmission. Notable differences between tsetse and other disease vectors include obligate microbial symbioses, viviparous reproduction, and lactation. Here, we describe the sequence and annotation of the 366-megabase Glossina morsitans morsitans genome. Analysis of the genome and the 12,308 predicted protein-encoding genes led to multiple discoveries, including chromosomal integrations of bacterial (Wolbachia) genome sequences, a family of lactation-specific proteins, reduced complement of host pathogen recognition proteins, and reduced olfaction/chemosensory associated genes. These genome data provide a foundation for research into trypanosomiasis prevention and yield important insights with broad implications for multiple aspects of tsetse biology.

Construction and analysis of Siberian tiger bacterial artificial chromosome library with approximately 6.5-fold genome equivalent coverage

Bacterial artificial chromosome (BAC) libraries are extremely valuable for the genome-wide genetic dissection of complex organisms. The Siberian tiger, one of the most well-known wild primitive carnivores in China, is an endangered animal. In order to promote research on its genome, a high-redundancy BAC library of the Siberian tiger was constructed and characterized. The library is divided into two sub-libraries prepared from blood cells and two sub-libraries prepared from fibroblasts. This BAC library contains 153,600 individually archived clones; for PCR-based screening of the library, BACs were placed into 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 116.5 kb, representing approximately 6.46 genome equivalents of the haploid genome and affording a 98.86% statistical probability of obtaining at least one clone containing a unique DNA sequence. Screening the library with 19 microsatellite markers and a SRY sequence revealed that each of these markers were present in the library; the average number of positive clones per marker was 6.74 (range 2 to 12), consistent with 6.46 coverage of the tiger genome. Additionally, we identified 72 microsatellite markers that could potentially be used as genetic markers. This BAC library will serve as a valuable resource for physical mapping, comparative genomic study and large-scale genome sequencing in the tiger.

Assembly of eukaryotic algal chromosomes in yeast

BACKGROUND

Synthetic genomic approaches offer unique opportunities to use powerful yeast and Escherichia coli genetic systems to assemble and modify chromosome-sized molecules before returning the modified DNA to the target host. For example, the entire 1 Mb Mycoplasma mycoides chromosome can be stably maintained and manipulated in yeast before being transplanted back into recipient cells. We have previously demonstrated that cloning in yeast of large (> ~ 150 kb), high G + C (55%) prokaryotic DNA fragments was improved by addition of yeast replication origins every ~100 kb. Conversely, low G + C DNA is stable (up to at least 1.8 Mb) without adding supplemental yeast origins. It has not been previously tested whether addition of yeast replication origins similarly improves the yeast-based cloning of large (>150 kb) eukaryotic DNA with moderate G + C content. The model diatom Phaeodactylum tricornutum has an average G + C content of 48% and a 27.4 Mb genome sequence that has been assembled into chromosome-sized scaffolds making it an ideal test case for assembly and maintenance of eukaryotic chromosomes in yeast.

RESULTS

We present a modified chromosome assembly technique in which eukaryotic chromosomes as large as ~500 kb can be assembled from cloned ~100 kb fragments. We used this technique to clone fragments spanning P. tricornutum chromosomes 25 and 26 and to assemble these fragments into single, chromosome-sized molecules. We found that addition of yeast replication origins improved the cloning, assembly, and maintenance of the large chromosomes in yeast. Furthermore, purification of the fragments to be assembled by electroelution greatly increased assembly efficiency.

CONCLUSIONS

Entire eukaryotic chromosomes can be successfully cloned, maintained, and manipulated in yeast. These results highlight the improvement in assembly and maintenance afforded by including yeast replication origins in eukaryotic DNA with moderate G + C content (48%). They also highlight the increased efficiency of assembly that can be achieved by purifying fragments before assembly.

Enhancing metagenomics investigations of microbial interactions with biofilm technology

Investigations of microbial ecology and diversity have been greatly enhanced by the application of culture-independent techniques. One such approach, metagenomics, involves sample collections from soil, water, and other environments. Extracted nucleic acids from bulk environmental samples are sequenced and analyzed, which allows microbial interactions to be inferred on the basis of bioinformatics calculations. In most environments, microbial interactions occur predominately in surface-adherent, biofilm communities. In this review, we address metagenomics sampling and biofilm biology, and propose an experimental strategy whereby the resolving power of metagenomics can be enhanced by incorporating a biofilm-enrichment step during sample acquisition.

Identifying Distal cis-acting Gene-Regulatory Sequences by Expressing BACs Functionalized with loxP-Tn10 Transposons in Zebrafish

Bacterial Artificial Chromosomes (BACs) are large pieces of DNA from the chromosomes of organisms propagated faithfully in bacteria as large extra-chromosomal plasmids. Expression of genes contained in BACs can be monitored after functionalizing the BAC DNA with reporter genes and other sequences that allow stable maintenance and propagation of the DNA in the new host organism. The DNA in BACs can be altered within its bacterial host in several ways. Here we discuss one such approach, using Tn10 mini-transposons, to introduce exogenous sequences into BACs for a variety of purposes. The largely random insertions of Tn10 transposons carrying lox sites have been used to position mammalian cell-selectable antibiotic resistance genes, enhancer-traps and inverted repeat ends of the vertebrate transposon Tol2 precisely at the ends of the genomic DNA insert in BACs. These modified BACs are suitable for expression in zebrafish or mouse, and have been used to functionally identify important long-range gene regulatory sequences in both species. Enhancer-trapping using BACs should prove uniquely useful in analyzing multiple discontinuous DNA domains that act in concert to regulate expression of a gene, and is not limited by genome accessibility issues of traditional enhancer-trapping methods.

The RNase II/RNB family of exoribonucleases: putting the 'Dis' in disease

Important findings over the last years have shed new light onto the mechanistic details of RNA degradation by members of the RNase II/RNB family of exoribonucleases. Members of this family have been shown to be involved in growth, normal chloroplast biogenesis, mitotic control and cancer. Recently, different publications have linked human orthologs (Dis3 and Dis3L2) to important human diseases. This article describes the structural and biochemical characteristics of members of this family of enzymes, and the physiological implications that relate them with disease.

Construction and preliminary characterization analysis of Wuzhishan miniature pig bacterial artificial chromosome library with approximately 8-fold genome equivalent coverage

Bacterial artificial chromosome (BAC) libraries have been invaluable tools for the genome-wide genetic dissection of complex organisms. Here, we report the construction and characterization of a high-redundancy BAC library from a very valuable pig breed in China, Wuzhishan miniature pig (Sus scrofa), using its blood cells and fibroblasts, respectively. The library contains approximately 153,600 clones ordered in 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 152.3 kb, representing approximately 7.68 genome equivalents of the porcine haploid genome and a 99.93% statistical probability of obtaining at least one clone containing a unique DNA sequence in the library. 19 pairs of microsatellite marker primers covering porcine chromosomes were used for screening the BAC library, which showed that each of these markers was positive in the library; the positive clone number was 2 to 9, and the average number was 7.89, which was consistent with 7.68-fold coverage of the porcine genome. And there were no significant differences of genomic BAC library from blood cells and fibroblast cells. Therefore, we identified 19 microsatellite markers that could potentially be used as genetic markers. As a result, this BAC library will serve as a valuable resource for gene identification, physical mapping, and comparative genomics and large-scale genome sequencing in the porcine.

Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution

Lampreys are representatives of an ancient vertebrate lineage that diverged from our own ∼500 million years ago. By virtue of this deeply shared ancestry, the sea lamprey (P. marinus) genome is uniquely poised to provide insight into the ancestry of vertebrate genomes and the underlying principles of vertebrate biology. Here, we present the first lamprey whole-genome sequence and assembly. We note challenges faced owing to its high content of repetitive elements and GC bases, as well as the absence of broad-scale sequence information from closely related species. Analyses of the assembly indicate that two whole-genome duplications likely occurred before the divergence of ancestral lamprey and gnathostome lineages. Moreover, the results help define key evolutionary events within vertebrate lineages, including the origin of myelin-associated proteins and the development of appendages. The lamprey genome provides an important resource for reconstructing vertebrate origins and the evolutionary events that have shaped the genomes of extant organisms.

A PKS/NRPS/FAS hybrid gene cluster from Serratia plymuthica RVH1 encoding the biosynthesis of three broad spectrum, zeamine-related antibiotics

Serratia plymuthica strain RVH1, initially isolated from an industrial food processing environment, displays potent antimicrobial activity towards a broad spectrum of Gram-positive and Gram-negative bacterial pathogens. Isolation and subsequent structure determination of bioactive molecules led to the identification of two polyamino antibiotics with the same molecular structure as zeamine and zeamine II as well as a third, closely related analogue, designated zeamine I. The gene cluster encoding the biosynthesis of the zeamine antibiotics was cloned and sequenced and shown to encode FAS, PKS as well as NRPS related enzymes in addition to putative tailoring and export enzymes. Interestingly, several genes show strong homology to the pfa cluster of genes involved in the biosynthesis of long chain polyunsaturated fatty acids in marine bacteria. We postulate that a mixed FAS/PKS and a hybrid NRPS/PKS assembly line each synthesize parts of the backbone that are linked together post-assembly in the case of zeamine and zeamine I. This interaction reflects a unique interplay between secondary lipid and secondary metabolite biosynthesis. Most likely, the zeamine antibiotics are produced as prodrugs that undergo activation in which a nonribosomal peptide sequence is cleaved off.

A physical map of a BAC clone contig covering the entire autosome insertion between ovine MHC Class IIa and IIb

Background

The ovine Major Histocompatibility Complex (MHC) harbors genes involved in overall resistance/susceptibility of the host to infectious diseases. Compared to human and mouse, the ovine MHC is interrupted by a large piece of autosome insertion via a hypothetical chromosome inversion that constitutes ~25% of ovine chromosome 20. The evolutionary consequence of such an inversion and an insertion (inversion/insertion) in relation to MHC function remains unknown. We previously constructed a BAC clone physical map for the ovine MHC exclusive of the insertion region. Here we report the construction of a high-density physical map covering the autosome insertion in order to address the question of what the inversion/insertion had to do with ruminants during the MHC evolution.

Results

A total of 119 pairs of comparative bovine oligo primers were utilized to screen an ovine BAC library for positive clones and the orders and overlapping relationships of the identified clones were determined by DNA fingerprinting, BAC-end sequencing, and sequence-specific PCR. A total of 368 positive BAC clones were identified and 108 of the effective clones were ordered into an overlapping BAC contig to cover the consensus region between ovine MHC class IIa and IIb. Therefore, a continuous physical map covering the entire ovine autosome inversion/insertion region was successfully constructed. The map confirmed the bovine sequence assembly for the same homologous region. The DNA sequences of 185 BAC-ends have been deposited into NCBI database with the access numbers HR309252 through HR309068, corresponding to dbGSS ID 30164010 through 30163826.

Conclusions

We have constructed a high-density BAC clone physical map for the ovine autosome inversion/insertion between the MHC class IIa and IIb. The entire ovine MHC region is now fully covered by a continuous BAC clone contig. The physical map we generated will facilitate MHC functional studies in the ovine, as well as the comparative MHC evolution in ruminants.

A multiple gene complex on rice chromosome 4 is involved in durable resistance to rice blast

Quantitative trait loci (QTLs) for resistance to rice blast offer a potential source of durable disease resistance in rice. However, few QTLs have been validated in progeny testing, on account of their small phenotypic effects. To understand the genetic basis for QTL-mediated resistance to blast, we dissected a resistance QTL, qBR4-2, using advanced backcross progeny derived from a chromosome segment substitution line in which a 30- to 34-Mb region of chromosome 4 from the resistant cultivar Owarihatamochi was substituted into the genetic background of the highly susceptible Aichiasahi. The analysis resolved qBR4-2 into three loci, designated qBR4-2a, qBR4-2b, and qBR4-2c. The sequences of qBR4-2a and qBR4-2b, which lie 181 kb apart from each other and measure, 113 and 32 kb, respectively, appear to encode proteins with a putative nucleotide-binding site (NBS) and leucine-rich repeats (LRRs). Sequence analysis of the donor allele of qBR4-2a, the region with the largest effect among the three, revealed sequence variations in the NBS-LRR region. The effect of qBR4-2c was smallest among the three, but its combination with the donor alleles of qBR4-2a and qBR4-2b significantly enhanced blast resistance. qBR4-2 comprises three tightly linked QTLs that control blast resistance in a complex manner, and thus gene pyramiding or haplotype selection is the recommended strategy for improving QTL-mediated resistance to blast disease through the use of this chromosomal region.

Evolution of human-specific neural SRGAP2 genes by incomplete segmental duplication

Gene duplication is an important source of phenotypic change and adaptive evolution. We leverage a haploid hydatidiform mole to identify highly identical sequences missing from the reference genome, confirming that the cortical development gene Slit-Robo Rho GTPase-activating protein 2 (SRGAP2) duplicated three times exclusively in humans. We show that the promoter and first nine exons of SRGAP2 duplicated from 1q32.1 (SRGAP2A) to 1q21.1 (SRGAP2B) ∼3.4 million years ago (mya). Two larger duplications later copied SRGAP2B to chromosome 1p12 (SRGAP2C) and to proximal 1q21.1 (SRGAP2D) ∼2.4 and ∼1 mya, respectively. Sequence and expression analyses show that SRGAP2C is the most likely duplicate to encode a functional protein and is among the most fixed human-specific duplicate genes. Our data suggest a mechanism where incomplete duplication created a novel gene function-antagonizing parental SRGAP2 function-immediately "at birth" 2-3 mya, which is a time corresponding to the transition from Australopithecus to Homo and the beginning of neocortex expansion.

Antigen-presenting genes and genomic copy number variations in the Tasmanian devil MHC

Background

The Tasmanian devil (Sarcophilus harrisii) is currently under threat of extinction due to an unusual fatal contagious cancer called Devil Facial Tumour Disease (DFTD). DFTD is caused by a clonal tumour cell line that is transmitted between unrelated individuals as an allograft without triggering immune rejection due to low levels of Major Histocompatibility Complex (MHC) diversity in Tasmanian devils.

Results

Here we report the characterization of the genomic regions encompassing MHC Class I and Class II genes in the Tasmanian devil. Four genomic regions approximately 960 kb in length were assembled and annotated using BAC contigs and physically mapped to devil Chromosome 4q. 34 genes and pseudogenes were identified, including five Class I and four Class II loci. Interestingly, when two haplotypes from two individuals were compared, three genomic copy number variants with sizes ranging from 1.6 to 17 kb were observed within the classical Class I gene region. One deletion is particularly important as it turns a Class Ia gene into a pseudogene in one of the haplotypes. This deletion explains the previously observed variation in the Class I allelic number between individuals. The frequency of this deletion is highest in the northwestern devil population and lowest in southeastern areas.

Conclusions

The third sequenced marsupial MHC provides insights into the evolution of this dynamic genomic region among the diverse marsupial species. The two sequenced devil MHC haplotypes revealed three copy number variations that are likely to significantly affect immune response and suggest that future work should focus on the role of copy number variations in disease susceptibility in this species.

Genomic restructuring in the Tasmanian devil facial tumour: chromosome painting and gene mapping provide clues to evolution of a transmissible tumour

Devil facial tumour disease (DFTD) is a fatal, transmissible malignancy that threatens the world's largest marsupial carnivore, the Tasmanian devil, with extinction. First recognised in 1996, DFTD has had a catastrophic effect on wild devil numbers, and intense research efforts to understand and contain the disease have since demonstrated that the tumour is a clonal cell line transmitted by allograft. We used chromosome painting and gene mapping to deconstruct the DFTD karyotype and determine the chromosome and gene rearrangements involved in carcinogenesis. Chromosome painting on three different DFTD tumour strains determined the origins of marker chromosomes and provided a general overview of the rearrangement in DFTD karyotypes. Mapping of 105 BAC clones by fluorescence in situ hybridisation provided a finer level of resolution of genome rearrangements in DFTD strains. Our findings demonstrate that only limited regions of the genome, mainly chromosomes 1 and X, are rearranged in DFTD. Regions rearranged in DFTD are also highly rearranged between different marsupials. Differences between strains are limited, reflecting the unusually stable nature of DFTD. Finally, our detailed maps of both the devil and tumour karyotypes provide a physical framework for future genomic investigations into DFTD.

The world's largest carnivorous marsupial, the Tasmanian devil, is threatened with extinction due to the emergence of devil facial tumour disease (DFTD), a fatal transmissible tumour. Critical loss of genetic diversity has rendered the devil vulnerable to transmission of tumour cells by grafting or transplanting the cells while biting and jaw wrestling. Initial studies of DFTD tumours revealed rearrangements among tumour chromosomes, with several missing chromosomes and four additional marker chromosomes of unknown origin. Since then, new strains of the disease have emerged and appear to be derived from the original strain. With no prior information available regarding the location of genes on normal devil chromosomes, a necessary first step towards characterisation of chromosome rearrangements in DFTD was to construct a map of the normal devil genome. This enabled us to elucidate the chromosome rearrangements in three DFTD strains. In doing so we determined the origin of the marker chromosomes and compared the three strains to determine which areas of the genome are involved in ongoing tumour evolution. Interestingly, rearrangements between strains are limited to particular genomic regions, demonstrating the unusual stability of this unique cancer. This study is therefore an important first step towards understanding the genetics of DFTD.

Construction of BIBAC and BAC libraries from a variety of organisms for advanced genomics research

Large-insert BAC (bacterial artificial chromosome) and BIBAC (binary BAC) libraries are essential for modern genomics research for all organisms. We helped pioneer the BAC and BIBAC technologies, and by using them we have constructed hundreds of BAC and BIBAC libraries for different species of plants, animals, marine animals, insects, algae and microbes. These libraries have been used globally for different aspects of genomics research. Here we describe the procedure with the latest improvements that we have made and used for construction of BIBAC libraries. The procedure includes the preparation of BIBAC vectors, the preparation of clonable fragments of the desired size from the source DNA, the construction and transformation of BIBACs and, finally, the characterization and assembly of BIBAC libraries. We also specify the modifications necessary for construction of BAC libraries using the protocol. The entire protocol takes ∼7 d.

Evolution and function of the globin intergenic regulatory regions of the antarctic dragonfishes (Notothenioidei: Bathydraconidae)

As the Southern Ocean cooled to -1.8 °C over the past 40 My, the teleostean clade Notothenioidei diversified and, under reduced selection pressure for an oxygen-transporting apparatus, became less reliant on hemoglobin and red blood cells. At the extreme of this trend, the crown group of Antarctic icefishes (Channichthyidae) lost both components of oxygen transport. Under the decreased selection scenario, we hypothesized that the Antarctic dragonfishes (Bathydraconidae, the red-blooded sister clade to the icefishes) evolved lower blood hemoglobin concentrations because their globin gene complexes (α- and β-globin gene pairs linked by a regulatory intergene) transcribe globin mRNAs less effectively than those of basal notothenioids (e.g., the Nototheniidae [notothens]). To test our hypothesis, we 1) sequenced the α/β-intergenes of the adult globin complexes of three notothen and eight dragonfish species and 2) measured globin transcript levels in representative species from each group. The typical nototheniid intergene was ∼3-4 kb in length. The bathydraconid intergenes resolved into three subclasses (long [3.8 kb], intermediate [3.0 kb], and short [1.5-2.3 kb]) that corresponded to the three subclades proposed for the taxon. Although they varied in length due to indels, the three notothen and eight dragonfish intergenes contained a conserved ∼90-nt element that we have previously shown to be required for globin gene transcription. Using the quantitative polymerase chain reaction, we found that globin mRNA levels in red cells from one notothen species and from one species of each dragonfish subclade were equivalent statistically. Thus, our results indicate that the bathydraconids have evolved adult globin loci whose regulatory intergenes tend to be shorter than those of the more basal nototheniids yet are equivalent in transcriptional efficacy. Their low blood hemoglobin concentrations are most likely due to reduction in hematocrit.

A pair of new BAC and BIBAC vectors that facilitate BAC/BIBAC library construction and intact large genomic DNA insert exchange

BACKGROUND

Large-insert BAC and BIBAC libraries are important tools for structural and functional genomics studies of eukaryotic genomes. To facilitate the construction of BAC and BIBAC libraries and the transfer of complete large BAC inserts into BIBAC vectors, which is desired in positional cloning, we developed a pair of new BAC and BIBAC vectors.

RESULTS

The new BAC vector pIndigoBAC536-S and the new BIBAC vector BIBAC-S have the following features: 1) both contain two 18-bp non-palindromic I-SceI sites in an inverted orientation at positions that flank an identical DNA fragment containing the lacZ selection marker and the cloning site. Large DNA inserts can be excised from the vectors as single fragments by cutting with I-SceI, allowing the inserts to be easily sized. More importantly, because the two vectors contain different antibiotic resistance genes for transformant selection and produce the same non-complementary 3' protruding ATAA ends by I-SceI that suppress self- and inter-ligations, the exchange of intact large genomic DNA inserts between the BAC and BIBAC vectors is straightforward; 2) both were constructed as high-copy composite vectors. Reliable linearized and dephosphorylated original low-copy pIndigoBAC536-S and BIBAC-S vectors that are ready for library construction can be prepared from the high-copy composite vectors pHZAUBAC1 and pHZAUBIBAC1, respectively, without the need for additional preparation steps or special reagents, thus simplifying the construction of BAC and BIBAC libraries. BIBAC clones constructed with the new BIBAC-S vector are stable in both E. coli and Agrobacterium. The vectors can be accessed through our website http://GResource.hzau.edu.cn.

CONCLUSIONS

The two new vectors and their respective high-copy composite vectors can largely facilitate the construction and characterization of BAC and BIBAC libraries. The transfer of complete large genomic DNA inserts from one vector to the other is made straightforward.

Construction of an American mink bacterial artificial chromosome (BAC) library and sequencing candidate genes important for the fur industry

Background

Bacterial artificial chromosome (BAC) libraries continue to be invaluable tools for the genomic analysis of complex organisms. Complemented by the newly and fast growing deep sequencing technologies, they provide an excellent source of information in genomics projects.

Results

Here, we report the construction and characterization of the CHORI-231 BAC library constructed from a Danish-farmed, male American mink (Neovison vison). The library contains approximately 165,888 clones with an average insert size of 170 kb, representing approximately 10-fold coverage. High-density filters, each consisting of 18,432 clones spotted in duplicate, have been produced for hybridization screening and are publicly available. Overgo probes derived from expressed sequence tags (ESTs), representing 21 candidate genes for traits important for the mink industry, were used to screen the BAC library. These included candidate genes for coat coloring, hair growth and length, coarseness, and some receptors potentially involved in viral diseases in mink. The extensive screening yielded positive results for 19 of these genes. Thirty-five clones corresponding to 19 genes were sequenced using 454 Roche, and large contigs (184 kb in average) were assembled. Knowing the complete sequences of these candidate genes will enable confirmation of the association with a phenotype and the finding of causative mutations for the targeted phenotypes.

Additionally, 1577 BAC clones were end sequenced; 2505 BAC end sequences (80% of BACs) were obtained. An excess of 2 Mb has been analyzed, thus giving a snapshot of the mink genome.

Conclusions

The availability of the CHORI-321 American mink BAC library will aid in identification of genes and genomic regions of interest. We have demonstrated how the library can be used to identify specific genes of interest, develop genetic markers, and for BAC end sequencing and deep sequencing of selected clones. To our knowledge, this is the first report of 454 sequencing of selected BAC clones in mammals and re-assures the suitability of this technique for obtaining the sequence information of genes of interest in small genomics projects. The BAC end sequences described in this paper have been deposited in the GenBank data library [HN339419-HN341884, HN604664-HN604702]. The 454 produced contigs derived from selected clones are deposited with reference numbers [GenBank: JF288166-JF288183 &JF310744].

Generating libraries of iTol2-end insertions at BAC ends using loxP and lox511 Tn10 transposons

BACKGROUND

Bacterial Artificial Chromosomes (BACs) have been widely used as transgenes in vertebrate model systems such as mice and zebrafish, for a variety of studies. BAC transgenesis has been a powerful tool to study the function of the genome, and gene regulation by distal cis-regulatory elements. Recently, BAC transgenesis in both mice and zebrafish was further facilitated by development of the transposon-mediated method using the Tol2 element. Tol2 ends, in the inverted orientation and flanking a 1 kb spacer DNA (iTol2), were introduced into the BAC DNA within the bacterial host using recombination of homologous sequences. Here we describe experiments designed to determine if a simpler and more flexible system could modify BACs so that they would be suitable for transgenesis into zebrafish or mouse embryos using the Tol2 transposase.

RESULTS

A new technique was developed to introduce recognition sequences for the Tol2 transposase into BACs in E. coli using the Tn10 transposon vector system. We constructed pTnloxP-iTol2kan and pTnlox511-iTol2kan to introduce the loxP or lox511 site and iTol2 cassette, containing the Tol2 cis-sequences in the inverted orientation, into BACs that have loxP and lox511 sites flanking genomic DNA inserts by Tn10-mediated transposition. The procedure enables rapid generation of a large collection of BACs ready for transgenesis with the iTol2 cassette at the new end of a progressively truncated genomic insert via lox-Cre recombination. The iTol2 ends are efficiently recognized by the Tol2 transposase, and the BACs readily integrate into zebrafish chromosomes.

CONCLUSION

The new technology described here can rapidly introduce iTol2 ends at a BAC end of choice, and simultaneously generate a large collection of BACs with progressive deletions of the genomic DNA from that end in a single experiment. This procedure should be applicable to a wider variety of BACs containing lox sites flanking the genomic DNA insert, including those with sequence repeats. The libraries of iTol2 inserted BACs with truncations from an end should facilitate studies on the impact of distal cis-regulatory sequences on gene function, as well as standard BAC transgenesis with precisely trimmed genes in zebrafish or mouse embryos using Tol2 transposition.

Construction of Japanese BAC library Yamato-2 (JY2): a set of 330K clone resources of damage-minimized DNA taken from a genetically established Japanese individual

A bacterial artificial chromosome (BAC) library referred to as Yamato-2 (JY2), was constructed from a Japanese individual and contained 330,000 clones. Library construction was based on 2 concepts: Japanese pedigree and non-immortalization. Genomic DNA was extracted from white blood cells from umbilical cord blood of a Japanese male individual. Four traits of the sample, (1) amelogenin DNA, (2) short tandem repeat (STR), (3) mitochondrial DNA (mtDNA), and (4) HLA-allele typing, were investigated to verify attribution of the donor. One of the samples with quite good Japanese characteristics was named JY2 and used as a resource for construction of a BAC library. Amelogenin DNA indicated male. STR indicated Mongoloid. MtDNA suggested haplogroup B, which is different from any other diploid whose sequence has been reported. The HLA gene was classified into east-Asian specific haplotype. These results revealed that JY2 was obtained from a Japanese male. We sequenced both ends of 185,012 BAC clones. By using the BLAST search, BAC end sequences (BESs) were mapped on the human reference sequence provided by NCBI. Inserts of individual BAC clones were mapped with both ends properly placed. As a result, 103,647 BAC clones were successfully mapped. The average insert size of BAC calculated from the mapping information was 130 kb. Coverage and redundancy of the reference sequence by successfully mapped BAC clones were 96.4% and 3.9-fold, respectively. This library will be especially suitable as a Japanese standard genome resource. The availability of an accurate library is indispensable for diagnostics or drug-design based on genome information, and JY2 will provide an accurate sequence of the Japanese genome as an important addition to the human genome.

BAC library resources for map-based cloning and physical map construction in barley (Hordeum vulgare L.)

BACKGROUND

Although second generation sequencing (2GS) technologies allow re-sequencing of previously gold-standard-sequenced genomes, whole genome shotgun sequencing and de novo assembly of large and complex eukaryotic genomes is still difficult. Availability of a genome-wide physical map is therefore still a prerequisite for whole genome sequencing for genomes like barley. To start such an endeavor, large insert genomic libraries, i.e. Bacterial Artificial Chromosome (BAC) libraries, which are unbiased and representing deep haploid genome coverage, need to be ready in place.

RESULT

Five new BAC libraries were constructed for barley (Hordeum vulgare L.) cultivar Morex. These libraries were constructed in different cloning sites (HindIII, EcoRI, MboI and BstXI) of the respective vectors. In order to enhance unbiased genome representation and to minimize the number of gaps between BAC contigs, which are often due to uneven distribution of restriction sites, a mechanically sheared library was also generated. The new BAC libraries were fully characterized in depth by scrutinizing the major quality parameters such as average insert size, degree of contamination (plate wide, neighboring, and chloroplast), empty wells and off-scale clones (clones with <30 or >250 fragments). Additionally a set of gene-based probes were hybridized to high density BAC filters and showed that genome coverage of each library is between 2.4 and 6.6 X.

CONCLUSION

BAC libraries representing >20 haploid genomes are available as a new resource to the barley research community. Systematic utilization of these libraries in high-throughput BAC fingerprinting should allow developing a genome-wide physical map for the barley genome, which will be instrumental for map-based gene isolation and genome sequencing.