BAC End Sequences Research Articles

The Melon Zym Locus Conferring Resistance to ZYMV: High Resolution Mapping and Candidate Gene Identification

Zucchini yellow mosaic virus (ZYMV; potyviridae) represents a major pathogen of Cucurbitaceae crops. ZYMV resistance in melon PI 414723 is conditioned by a dominant allele at the Zym locus. This resistant accession restricts viral spread and does not develop mosaic symptoms, but necrosis sometimes develops in response to inoculation. In previous studies, Zym has been mapped to linkage group II of the melon genetic map. In the present study, positional cloning of the locus was undertaken, starting from the CM-AG36 SSR marker at approximately 2 cm distance. We utilized five mapping populations that share the same resistant parent, PI 414723, and analyzed a total of 1630 offspring, to construct a high-resolution genetic map of the Zym locus. Two melon BAC libraries were used for chromosome walking and for developing new markers closer to the resistance gene by BAC-end sequencing. A BAC contig was constructed, and we identified a single BAC clone, from the ZYMV susceptible genotype MR-1, that physically encompasses the resistance gene. A second clone was isolated from another susceptible genotype, WMR 29, and the two clones were fully sequenced and annotated. Additional markers derived from the sequenced region delimited the region to 17.6 kb of a sequence that harbors a NAC-like transcription factor and, depending on the genotype, either two or three R-gene homologs with a CC-NBS-LRR structure. Mapping was confirmed by saturating the map with SNP markers using a single mapping population. The same region was amplified and sequenced also in the ZYMV resistant genotype PI 414723. Because numerous polymorphic sites were noted between genotypes, we could not associate resistance with a specific DNA polymorphism; however, this study enables molecular identification of Zym and paves the way to functional studies of this important locus.

Genomic resources of broomcorn millet: demonstration and application of a high-throughput BAC mapping pipeline

BackgroundWith high-efficient water-use and drought tolerance, broomcorn millet has emerged as a candidate for food security. To promote its research process for molecular breeding and functional research, a comprehensive genome resource is of great importance.ResultsHerein, we constructed a BAC library for broomcorn millet, generated BAC end sequences based on the clone-array pooled shotgun sequencing strategy and Illumina sequencing technology, and integrated BAC clones into genome by a novel pipeline for BAC end profiling. The BAC library consisted of 76,023 clones with an average insert length of 123.48 Kb, covering about 9.9-fold of the 850 Mb genome. Of 9216 clones tested using our pipeline, 8262 clones were mapped on the broomcorn millet cultivar longmi4 genome. These mapped clones covered 308 of the 829 gaps left by the genome. To our knowledge, this is the only BAC resource for broomcorn millet.ConclusionsWe constructed a high-quality BAC libraray for broomcorn millet and designed a novel pipeline for BAC end profiling. BAC clones can be browsed and obtained from our website (http://eightstarsbio.com/gresource/JBrowse-1.16.5/index.html). The high-quality BAC clones mapped on genome in this study will provide a powerful genomic resource for genome gap filling, complex segment sequencing, FISH, functional research and genetic engineering of broomcorn millet.

Long-insert clone experimental evidence for assembly improvement and chimeric chromosomes detection in an allopentaploid beer yeast.

Lager beer is made with the hybrid Saccharomyces pastorianus. Many publicly available S. pastorianus genome assemblies are highly fragmented due to the difficulties of assembling hybrid genomes, such as the presence of homeologous chromosomes from both parental types, and translocations between them. To improve the assembly of a previously sequenced lager yeast hybrid Saccharomyces sp. 790 and elucidate its genome structure, we proposed the use of alternative experimental evidence. We determined the phylogenetic position of Saccharomyces sp. 790 and established it as S. pastorianus 790. Then, we obtained from this yeast a bacterial artificial chromosome (BAC) genomic library with its BAC-end sequences (BESs). To analyze these data, we developed a pipeline (applicable to other assemblies) that classifies BES pairs alignments according to their orientation. For the case of S. pastorianus 790, paired-end BESs alignments validated parts of the assembly and unpaired-end ones suggested contig joins or misassemblies. Importantly, the BACs library was preserved and used for verification experiments. Unpaired-end alignments were used to upgrade the previous assembly and provided an improved detection of translocations. With this, we proposed a genome structure of S. pastorianus 790, which was similar to that of other lager yeasts; however, when we estimated chromosome copy number and experimentally measured its genome size, we discovered that one key difference is the outstanding S. pastorianus 790 ploidy level (allopentaploid). Altogether, our results show the value of combining bioinformatic analyses with experimental data such as long-insert clone information to improve a short-read assembly of a hybrid genome.

Development of a high-density SSR genetic linkage map in sweet potato

Simple sequence repeat (SSR) markers have been proved to be a very powerful tool for quantitative trait locus (QTL) mapping, marker-assisted selection and comparative genomics research in many crop species. However, a high-density SSR genetic linkage map is still lacking because there are only a few SSR markers available in sweet potato. In this study, a total of 2545 simple sequence repeat (SSR) primer pairs, including 1215 genomic SSR (gSSR) primer pairs and 1330 BES-SSR (bSSR) primer pairs designed from the genome sequence and BAC-end sequence of sweet potato, respectively, were screened with sweet potato cultivars Luoxushu 8 and Zhengshu 20 and their randomly sampled two F1 individuals and 571 of them generated polymorphic bands. The selected 571 polymorphic SSR primer pairs and 35 EST-based SSR (eSSR) primer pairs developed at our laboratory were used to genotype 240 F1 individuals derived from a cross between Luoxushu 8 and Zhengshu 20. A double pseudo-test-cross strategy was applied for linkage analysis. The Luoxushu 8 map included 90 linkage groups with 5057 SSR markers and covered 13,299.9 cM with a marker density of 2.6 cM, and the Zhengshu 20 map contained 90 linkage groups with 3009 SSR markers and covered 11,122.9 cM with a marker density of 3.7 cM. Fifteen homologous groups were identified in both parent maps. These are the first SSR linkage maps consisting of the complete 90 linkage groups and 15 homologous groups, which are consistent with the autohexaploid nature of sweetpotato, and are also the linkage maps with the highest SSR marker density reported to date. These results provide a basis for QTL mapping, marker-assisted breeding and comparative genomics research of sweet potato.

Fine mapping of the locus controlling self-incompatibility in European hazelnut

Incompatibility in European hazelnut (Corylus avellana L.) is sporophytic and under the control of a single locus on linkage group 5 between markers G05-510 and AU02-1350. In this study, two rounds of marker development and a population of 192 seedlings with known S-alleles that showed recombination between the flanking markers were used for fine mapping. Using the sequences of random amplified polymorphic DNA and simple sequence repeat (SSR) markers and bacterial artificial chromosome end sequences, 36 contigs from the genome sequence of “Jefferson” hazelnut were identified for pursuit. Di-nucleotide SSR markers in those contigs were developed, characterized, and mapped. This reduced the size to a region of 500 kb that contained the S-locus and 50 predicted genes, in which single-nucleotide polymorphism and additional SSR markers were developed. When the new markers were used in fine mapping, they fully exploited all recombination in the fine mapping population and reduced the region to 193.5 kb containing 18 genes. This 193.5-kb region most likely contains the S1 haplotype. A second region, 2 Mbp away from the first in the “Jefferson” genome (V3), is predicted to represent the S3 haplotype based on SSR marker allele sizes and the ratio of parental reads that align to them. Although they appear side-by-side in the “Jefferson” genome (V3), the mapped markers appear in both sequences in the same order. Gene annotations in the S1 and S3 haplotypes are highly similar and include five probable leucine-rich repeat receptor-like protein kinases with homology to Arabidopsis thaliana genes At1g35710 and three probable receptor-like serine/threonine protein kinases with homology to At5g15080 (PIX7).

Identification of high-copy number long terminal repeat retrotransposons and their expansion in Phalaenopsis orchids

BackgroundTransposable elements (TEs) are fragments of DNA that can insert into new chromosomal locations. They represent a great proportion of eukaryotic genomes. The identification and characterization of TEs facilitates understanding the transpositional activity of TEs with their effects on the orchid genome structure.ResultsWe combined the draft whole-genome sequences of Phalaenopsis equestris with BAC end sequences, Roche 454, and Illumina/Solexa, and identified long terminal repeat (LTR) retrotransposons in these genome sequences by using LTRfinder and classified by using Gepard software. Among the 10 families Gypsy-like retrotransposons, three families Gypsy1, Gypsy2, and Gypsy3, contained the most copies among these predicted elements. In addition, six high-copy retrotransposons were identified according to their reads in the sequenced raw data. The 12-kb Orchid-rt1 contains 18,000 copies representing 220 Mbp of the P. equestris genome. Southern blot and slot blot assays showed that these four retrotransposons Gypsy1, Gypsy2, Gypsy3, and Orchid-rt1 contained high copies in the large-genome-size/large-chromosome species P. violacea and P. bellina. Both Orchid-rt1 and Gypsy1 displayed various ratios of copy number for the LTR sequences versus coding sequences among four Phalaenopsis species, including P. violacea and P. bellina and small-genome-size/small-chromosome P. equestris and P. ahprodite subsp. formosana, which suggests that Orchid-rt1 and Gypsy1 have been through various mutations and homologous recombination events. FISH results showed amplification of Orchid-rt1 in the euchromatin regions among the four Phalaenopsis species. The expression levels of Peq018599 encoding copper transporter 1 is highly upregulated with the insertion of Orchid-rt1, while it is down regulated for Peq009948 and Peq014239 encoding for a 26S proteasome non-ATP regulatory subunit 4 homolog and auxin-responsive factor AUX/IAA-related. In addition, insertion of Orchid-rt1 in these three genes are all in their intron regions.ConclusionOrchid-rt1 and Gypsy1–3 have amplified within Phalaenopsis orchids concomitant with the expanded genome sizes, and Orchid-rt1 and Gypsy1 may have gone through various mutations and homologous recombination events. Insertion of Orchid-rt1 is in the introns and affects gene expression levels.

Improved protocols for BAC insert DNA isolation, BAC end sequencing and FISH for construction of BAC based physical map of genes on the chromosomes.

Bacterial artificial chromosome (BAC) library is an important genomic resource useful in targeted marker development, positional cloning, physical mapping and a substrate for genome sequencing for better understanding the genome organization of a species. The present manuscript elucidates the improvement in protocols for economical and efficient BAC insert DNA isolation, BAC end sequencing and FISH for physical localization on the metaphase chromosome complements. BAC clones of Clarias magur, maintained in 384-well plate format in our laboratory, were used in this study. The protocols gave consistent and efficient results. We use routinely these protocols for BAC insert DNA extraction, generating end sequence data of the clone and constructing DNA probes to hybridize on the metaphase spreads of C. magur using FISH for physical their localization.

Complete Genome Sequence of Sequevar 14M Ralstonia solanacearum Strain HA4-1 Reveals Novel Type III Effectors Acquired Through Horizontal Gene Transfer.

Ralstonia solanacearum, which causes bacterial wilt in a broad range of plants, is considered a “species complex” due to its significant genetic diversity. Recently, we have isolated a new R. solanacearum strain HA4-1 from Hong’an county in Hubei province of China and identified it being phylotype I, sequevar 14M (phylotype I-14M). Interestingly, we found that it can cause various disease symptoms among different potato genotypes and display different pathogenic behavior compared to a phylogenetically related strain, GMI1000. To dissect the pathogenic mechanisms of HA4-1, we sequenced its whole genome by combined sequencing technologies including Illumina HiSeq2000, PacBio RS II, and BAC-end sequencing. Genome assembly results revealed the presence of a conventional chromosome, a megaplasmid as well as a 143 kb plasmid in HA4-1. Comparative genome analysis between HA4-1 and GMI1000 shows high conservation of the general virulence factors such as secretion systems, motility, exopolysaccharides (EPS), and key regulatory factors, but significant variation in the repertoire and structure of type III effectors, which could be the determinants of their differential pathogenesis in certain potato species or genotypes. We have identified two novel type III effectors that were probably acquired through horizontal gene transfer (HGT). These novel R. solanacearum effectors display homology to several YopJ and XopAC family members. We named them as RipBR and RipBS. Notably, the copy of RipBR on the plasmid is a pseudogene, while the other on the megaplasmid is normal. For RipBS, there are three copies located in the megaplasmid and plasmid, respectively. Our results have not only enriched the genome information on R. solanacearum species complex by sequencing the first sequevar 14M strain and the largest plasmid reported in R. solanacearum to date but also revealed the variation in the repertoire of type III effectors. This will greatly contribute to the future studies on the pathogenic evolution, host adaptation, and interaction between R. solanacearum and potato.

A high-density BAC physical map covering the entire MHC region of addax antelope genome

BackgroundThe mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world.ResultsThe high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing.ConclusionsWe report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

Integrating a newly developed BAC-based physical mapping resource for Lolium perenne with a genome-wide association study across a L. perenne European ecotype collection identifies genomic contexts associated with agriculturally important traits.

Background and Aims Lolium perenne (perennial ryegrass) is the most widely cultivated forage and amenity grass species in temperate areas worldwide and there is a need to understand the genetic architectures of key agricultural traits and crop characteristics that deliver wider environmental services. Our aim was to identify genomic regions associated with agriculturally important traits by integrating a bacterial artificial chromosome (BAC)-based physical map with a genome-wide association study (GWAS).MethodsBAC-based physical maps for L. perenne were constructed from ~212 000 high-information-content fingerprints using Fingerprint Contig and Linear Topology Contig software. BAC clones were associated with both BAC-end sequences and a partial minimum tiling path sequence. A panel of 716 L. perenne diploid genotypes from 90 European accessions was assessed in the field over 2 years, and genotyped using a Lolium Infinium SNP array. The GWAS was carried out using a linear mixed model implemented in TASSEL, and extended genomic regions associated with significant markers were identified through integration with the physical map.Key ResultsBetween ~3600 and 7500 physical map contigs were derived, depending on the software and probability thresholds used, and integrated with ~35 k sequenced BAC clones to develop a resource predicted to span the majority of the L. perenne genome. From the GWAS, eight different loci were significantly associated with heading date, plant width, plant biomass and water-soluble carbohydrate accumulation, seven of which could be associated with physical map contigs. This allowed the identification of a number of candidate genes.ConclusionsCombining the physical mapping resource with the GWAS has allowed us to extend the search for candidate genes across larger regions of the L. perenne genome and identified a number of interesting gene model annotations. These physical maps will aid in validating future sequence-based assemblies of the L. perenne genome.

Molecular marker development and primary physical map construction for the tuber shape Ro gene locus in diploid potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important crop worldwide that is consumed fresh and processed in many different ways. The tuber shape is a critical trait in potato breeding, because it directly affects processing quality. Previous studies have demonstrated that the Ro locus on chromosome X is the major locus controlling tuber shape, and the round shape allele, Ro, is dominant over the long shape allele, ro. However, no reports on the cloning or physical mapping of this locus have been published. In this study, whole genome sequencing and bacterial artificial chromosome (BAC) library construction were applied to develop molecular markers and construct a physical map of Ro. Eight round-tuber genotypes were used to construct the BAC library, which contains about 40,000 clones with an average insert size of 85 kb and represents approximately four haploid potato genomes. A total of 66 individuals with extreme round or long tubers were selected to construct a round-tuber shape DNA bulk and a long-tuber shape DNA bulk for WGS. Single-nucleotide polymorphisms (SNPs) identified from WGS data were used to develop nine polymorphic markers linked to Ro. These markers were used to screen the BAC library, and eight positive BAC clones were ultimately identified and then analyzed by restriction enzyme digestion and BAC end sequencing and used to develop markers and identify new positive BAC clones. Finally, a primary physical map of Ro was constructed using a chromosome walking strategy. This map provides an important reference for cloning the potato tuber shape gene and studying tuber shape variation.

Analysis of Genome-Wide Differentiation between Native and Introduced Populations of the Cupped Oysters Crassostrea gigas and Crassostrea angulata.

The Pacific cupped oyster is genetically subdivided into two sister taxa, Crassostrea gigas and Crassostrea angulata, which are in contact in the north-western Pacific. The nature and origin of their genetic and taxonomic differentiation remains controversial due the lack of known reproductive barriers and the high degree of morphologic similarity. In particular, whether the presence of ecological and/or intrinsic isolating mechanisms contributes to species divergence is unknown. The recent co-introduction of both taxa into Europe offers a unique opportunity to test how genetic differentiation is maintained under new environmental and demographic conditions. We generated a pseudochromosome assembly of the Pacific oyster genome using a combination of BAC-end sequencing and scaffold anchoring to a new high-density linkage map. We characterized genome-wide differentiation between C. angulata and C. gigas in both their native and introduced ranges, and showed that gene flow between species has been facilitated by their recent co-introductions in Europe. Nevertheless, patterns of genomic divergence between species remain highly similar in Asia and Europe, suggesting that the environmental transition caused by the co-introduction of the two species did not affect the genomic architecture of their partial reproductive isolation. Increased genetic differentiation was preferentially found in regions of low recombination. Using historical demographic inference, we show that the heterogeneity of differentiation across the genome is well explained by a scenario whereby recent gene flow has eroded past differentiation at different rates across the genome after a period of geographical isolation. Our results thus support the view that low-recombining regions help in maintaining intrinsic genetic differences between the two species.

Genome sequence comparison between Chinese hamster ovary (CHO) DG44 cells and mouse using end sequences of CHO BAC clones based on BAC-FISH results.

Chinese hamster ovary (CHO) cells have frequently been used in biotechnology as a mammalian host cell platform for expressing genes of interest. Previously, we constructed a detailed physical chromosomal map of the CHO DG44 cell line by fluorescence in situ hybridization (FISH) imaging using 303 bacterial artificial chromosome (BAC) clones as hybridization probes (BAC-FISH). BAC-FISH results revealed that the two longest chromosomes were completely paired. However, other chromosomes featured partial deletions or rearrangements. In this study, we determined the end sequences of 303 BAC clones (BAC end sequences), which were used for BAC-FISH probes. Among 606 BAC-end sequences (BESs) (forward and reverse ends), 558 could be determined. We performed a comparison between all determined BESs and mouse genome sequences using NCBI BLAST. Among these 558 BESs, 465 showed high homology to mouse chromosomal sequences. We analyzed the locations of these BACs in chromosomes of the CHO DG44 cell line using a physical chromosomal map. From the obtained results, we investigated the regional similarities among CHO chromosomes (A-T) and mouse chromosomes (1-19 and sex) about 217 BESs (46.7% of 465 high homologous BESs). Twenty-three specific narrow regions in 13 chromosomes of the CHO DG44 cell line showed high homology to mouse chromosomes, but most of other regions did not show significant correlations with the mouse genome. These results contribute to accurate alignments of chromosomes of Chinese hamster and its genome sequence, analysis of chromosomal instability in CHO cells, and the development of target locations for gene and/or genome editing techniques.

Draft Sequencing of the Heterozygous Diploid Genome of Satsuma (Citrus unshiu Marc.) Using a Hybrid Assembly Approach.

Satsuma (Citrus unshiu Marc.) is one of the most abundantly produced mandarin varieties of citrus, known for its seedless fruit production and as a breeding parent of citrus. De novo assembly of the heterozygous diploid genome of Satsuma (“Miyagawa Wase”) was conducted by a hybrid assembly approach using short-read sequences, three mate-pair libraries, and a long-read sequence of PacBio by the PLATANUS assembler. The assembled sequence, with a total size of 359.7 Mb at the N50 length of 386,404 bp, consisted of 20,876 scaffolds. Pseudomolecules of Satsuma constructed by aligning the scaffolds to three genetic maps showed genome-wide synteny to the genomes of Clementine, pummelo, and sweet orange. Gene prediction by modeling with MAKER-P proposed 29,024 genes and 37,970 mRNA; additionally, gene prediction analysis found candidates for novel genes in several biosynthesis pathways for gibberellin and violaxanthin catabolism. BUSCO scores for the assembled scaffold and predicted transcripts, and another analysis by BAC end sequence mapping indicated the assembled genome consistency was close to those of the haploid Clementine, pummel, and sweet orange genomes. The number of repeat elements and long terminal repeat retrotransposon were comparable to those of the seven citrus genomes; this suggested no significant failure in the assembly at the repeat region. A resequencing application using the assembled sequence confirmed that both kunenbo-A and Satsuma are offsprings of Kishu, and Satsuma is a back-crossed offspring of Kishu. These results illustrated the performance of the hybrid assembly approach and its ability to construct an accurate heterozygous diploid genome.

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.

Genetic variation for domestication-related traits revealed in a cultivated rice, Nipponbare (Oryza sativa ssp. japonica) × ancestral rice, O. nivara, mapping population

Oryza nivara is the ancestral species of cultivated rice (Oryza sativa). It has been the source of novel alleles for resistance to biotic and abiotic stresses, as well as yield improvement, lost during the course of domestication. To determine the molecular changes that occurred during domestication, the O. sativa ssp. japonica variety, Nipponbare, from which a reference sequence (RefSeq) was developed, was crossed with the O. nivara accession (IRGC100897), from which BAC-end sequences (BES) were derived. The mapping population composed of 279 F2 progeny lines derived from this cross was phenotyped for 19 traits important to domestication and yield improvement, including basal sheath and culm color, culm angle, days to heading, plant height, seed shattering, flag leaf length and width, panicle type and length, awn length and color, pericarp color, and seed color, length, width, length to width ratio, volume and surface area. The population was genotyped using 95 SSR markers and 114 single nucleotide variation (SNV) markers, selected by comparing the Nipponbare RefSeq and O. nivara BES. At least one major QTL was identified for each trait evaluated, and for 28 of the 46 QTL, the trait increase was attributed to the allele contributed by the O. nivara parent. Candidate genes were identified in 37 of the QTL regions. This study validated SNV markers that can be used for mapping in populations with a wild species parent. In the future, SNVs could be used for marker-assisted selection to incorporate desirable, novel alleles for stress resistance and yield improvement, identified in rice wild species like O. nivara into elite, adapted O. sativa varieties.

Construction and characterization of a bacterial artificial chromosome library for Camellia sinensis

Tea is a popular and natural non-alcoholic beverage, and is produced from fresh leaves of Camellia sinensis. Tea leaves contain many bioactive compounds that have significant health benefits. We constructed a high quality bacterial artificial chromosome (BAC) library by using the fresh petals of C. sinensis “Shuchazao” for genome sequencing and improvement of genomic assembly. BAC library is still a significant tool for studies of functional genomes and preservation of precious genetic resources. The BAC library contains 161,280 clones with an average insert size of 113 kb, which represents approximately 6.2-fold coverage of haploid genome equivalents of C. sinensis. We characterized 20 complete BAC clones and 738 BAC end sequences (BESs) ranging from 105 to 917 bp. In addition, we predicted cis-regulatory elements of LAR (leucoanthocyanidin reductase), TCS (caffeine synthase), and TS (theanine synthetase) involved in tea characteristic metabolite synthesis and identified a larger number of light-responsive cis-acting elements in these three genes. Meanwhile, we analyzed alternative splicing of these three genes. Furthermore, 12 pairs of SSR primers were successfully amplified in tea plant DNA. The tea BAC library was a critical resource to accomplish de novo whole-genome sequencing, accelerate gene discovery and enhance molecular breeding of C. sinensis.

Bacterial artificial chromosome library construction of root-knot nematode resistant pepper genotype HDA149 and identification of clones linked to Me3 resistant locus

Pepper (Capsicum annuum. L.) is a widely cultivated vegetable crop worldwide and has the second largest planting area and the first largest vegetable output and value in China. Pepper root-knot nematode (Meloidogyne spp.) is one of the most serious pests of pepper, which caused huge losses every year. Previous studies showed that the Me3 gene is resistant to a wide range of Meloidogyne species, including M. arenaria, M. javanica, and M. incognita. HDA149, a double haploid pepper genotype, harboring the root-knot nematode resistance gene Me3, was used to construct bacterial artificial chromosome library (BAC) via the vector of CopyControl™ pCC1 in this study. The library consists of 210 200 BAC clones and is equivalent to 5.3 pepper genomes. The average insert size is 95 kb, and most of them are 90–120 kb; but the empty clones are less than 3%. In order to screen the BAC library easily, 550 super pools with 384 BAC clones of each pool were further developed in this study. Specific primers from Me3 gene locus were used for BAC library screening, and more than 20 positive BAC clones were obtained. Then the selected positive BAC clones were analyzed by restriction enzyme digestion, BAC-end sequencing, marker development, and new positive BAC clones exploration, respectively. Finally, the contig with total length of about 300 kb linked to the Me3 locus was constructed based on chromosome walking strategy, which made a solid foundation for the cloning of the important root-knot nematode resistance gene Me3.

In Silico Mining of Conserved miRNAs of Indian Catfish Clarias batrachus (Linnaeus, 1758) from the Contigs, ESTs, and BAC End Sequences.

MicroRNAs are small non-coding RNAs which play significant role in RNA interference. The present work deals with mining of the conserved miRNA and their target genes from the contigs, ESTs, and BAC end sequences of commercially important catfish, Clarias batrachus, from India. A total of 138, 1 and 1 conserved pre-miRNA sequences, were mined from the contigs, ESTs, and BAC end sequences, respectively. The analysis of families of the conserved pre-miRNA revealed conservation of the fish-specific family mir-430 and other important families, such as mir-455, let-7, mir-133, and mir-137. The mir-455 is involved in hypoxia signaling, let-7 family represents potential anti-tumor molecules involved in human cancer therapy, whereas mir-133 and mir-137 have high therapeutic potentials. Using an alternate computational in silico approach, mining of mature miRNAs resulted in identification of 210 mature miRNAs from contigs, 1 from EST, and 2 each from forward as well as reverse BAC end sequences. Target prediction of these putative miRNAs resulted in the identification of 66,758 and 18,747 target genes in C. batrachus and Danio rerio, respectively. Functional annotation of these miRNAs indicated their involvement in diverse biological functions. The findings of the present study can serve as a valuable resource for further functional genomics studies in C. batrachus.

A genome-wide BAC-end sequence survey provides first insights into sweetpotato (Ipomoea batatas (L.) Lam.) genome composition.

BackgroundSweetpotato, Ipomoea batatas (L.) Lam., is an important food crop widely grown in the world. However, little is known about the genome of this species because it is a highly heterozygous hexaploid. Gaining a more in-depth knowledge of sweetpotato genome is therefore necessary and imperative. In this study, the first bacterial artificial chromosome (BAC) library of sweetpotato was constructed. Clones from the BAC library were end-sequenced and analyzed to provide genome-wide information about this species.ResultsThe BAC library contained 240,384 clones with an average insert size of 101 kb and had a 7.93–10.82 × coverage of the genome, and the probability of isolating any single-copy DNA sequence from the library was more than 99%. Both ends of 8310 BAC clones randomly selected from the library were sequenced to generate 11,542 high-quality BAC-end sequences (BESs), with an accumulative length of 7,595,261 bp and an average length of 658 bp. Analysis of the BESs revealed that 12.17% of the sweetpotato genome were known repetitive DNA, including 7.37% long terminal repeat (LTR) retrotransposons, 1.15% Non-LTR retrotransposons and 1.42% Class II DNA transposons etc., 18.31% of the genome were identified as sweetpotato-unique repetitive DNA and 10.00% of the genome were predicted to be coding regions. In total, 3,846 simple sequences repeats (SSRs) were identified, with a density of one SSR per 1.93 kb, from which 288 SSRs primers were designed and tested for length polymorphism using 20 sweetpotato accessions, 173 (60.07%) of them produced polymorphic bands. Sweetpotato BESs had significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum than those of Vitis vinifera, Theobroma cacao and Arabidopsis thaliana.ConclusionsThe first BAC library for sweetpotato has been successfully constructed. The high quality BESs provide first insights into sweetpotato genome composition, and have significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum. These resources as a robust platform will be used in high-resolution mapping, gene cloning, assembly of genome sequences, comparative genomics and evolution for sweetpotato.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3302-1) contains supplementary material, which is available to authorized users.