Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.     

用于目标序列染色体定位的镜鲤BAC-FISH体系构建

为了构建用于镜鲤(Cyprinus carpio var. specularis)特定基因组序列染色体定位的实验体系, 在细菌人工染色体(Bacterial Artificial Chromosome, BAC)文库筛选池中对已知短序列基因组片段进行PCR扩增, 筛选出包含目标序列的BAC克隆, 提取BAC质粒进行缺刻平移标记制备探针, 开展荧光原位杂交(Fluorescence in situ hybridization, FISH)实验。通过对染色体片前处理、BAC质粒探针制备、C0t-1 DNA封闭基因组重复序列、预杂交、荧光染料选择、信号放大等一系列实验条件和方法的探索优化, 成功实现了目标序列在镜鲤有丝分裂中期染色体上的定位。定位对象既包括在染色体上有单一位点的序列, 如斑马鱼微卫星标记Z6884和Z4268, 也包括在染色体上有多个位点的重复序列, 如黄河鲤性别相关标记CCmf1。来自斑马鱼同一条染色体上的两个微卫星标记被分别定位于镜鲤不同染色体上, 为鲤鱼染色体数目加倍的进化假设提供了一项直接实验证据, 同时将现有遗传连锁图谱与染色体对应起来, 可作为染色体识别和细胞遗传学图谱构建的依据。黄河鲤性别相关重复序列被定位于不少于四条染色体上, 为性别决定相关基因的筛查提供了研究线索。这一BAC-FISH实验体系将成为鲤细胞遗传学图谱构建、基因组进化和比较基因组学研究中的重要研究工具。       

Development of Y chromosome intraspecific polymorphic markers in the Felidae

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.     

Construction of a BAC library and generation of BAC end sequence-tagged connectors for genome sequencing of the African malaria mosquito Anopheles gambiae

A Bacterial Artificial Chromosome (BAC) genomic DNA library of Anopheles gambiae, the major human malaria vector in sub-Saharan Africa, was constructed and characterized. This library (ND-TAM) is composed of 30,720 BAC clones in eighty 384-well plates. The estimated average insert size of the library is 133 kb, with an overall genome coverage of approximately 14-fold. The ends of approximately two-thirds of the clones in the library were sequenced, yielding 32,340 pair-mate ends. A statistical analysis (G-test) of the results of PCR screening of the library indicated a random distribution of BACs in the genome, although one gap encompassing the white locus on the X-chromosome was identified. Furthermore, combined with another previously constructed BAC library (ND-1), ~2,000 BACs have been physically mapped by polytene chromosomal in situ hybridization. These BAC end pair mates and physically mapped BACs have been useful for both the assembly of a fully sequenced A. gambiae genome and for linking the assembled sequence to the three polytene chromosomes. This ND-TAM library is now publicly available at both http://www.malaria.mr4.org/mr4pages/index.html/ and http://hbz.tamu.edu/, providing a valuable resource to the mosquito research community.     

Satellite DNA from the Y chromosome of the malaria vector Anopheles gambiae

Satellite DNA is an enigmatic component of genomic DNA with unclear function that has been regarded as "junk." Yet, persistence of these tandem highly repetitive sequences in heterochromatic regions of most eukaryotic chromosomes attests to their importance in the genome. We explored the Anopheles gambiae genome for the presence of satellite repeats and identified 12 novel satellite DNA families. Certain families were found in close juxtaposition within the genome. Six satellites, falling into two evolutionarily linked groups, were investigated in detail. Four of them were experimentally confirmed to be linked to the Y chromosome, whereas their relatives occupy centromeric regions of either the X chromosome or the autosomes. A complex evolutionary pattern was revealed among the AgY477-like satellites, suggesting their rapid turnover in the A. gambiae complex and, potentially, recombination between sex chromosomes. The substitution pattern suggested rolling circle replication as an array expansion mechanism in the Y-linked 53-bp satellite families. Despite residing in different portions of the genome, the 53-bp satellites share the same monomer lengths, apparently maintained by molecular drive or structural constraints. Potential functional centromeric DNA structures, consisting of twofold dyad symmetries flanked by a common sequence motif, have been identified in both satellite groups.     

Integrated Genetic Map of Anopheles Gambiae: Use of Rapd Polymorphisms for Genetic, Cytogenetic and Sts Landmarks

Randomly amplified polymorphic DNA (RAPD) markers have been integrated in the genetic and cytogenetic maps of the malaria vector mosquito, Anopheles gambiae. Fifteen of these markers were mapped by recombination, relative to microsatellite markers that had been mapped previously. Thirty-four gel-purified RAPD bands were cloned and sequenced, generating sequence tagged sites (STSs) that can be used as entry points to the A. gambiae genome. Thirty one of these STSs were localized on nurse cell polytene chromosomes through their unique hybridization signal in in situ hybridization experiments. Five STSs map close to the breakpoints of polymorphic inversions, which are notable features of the Anopheles genome. The usefulness and limitations of this integrated mosquito map are discussed.     

Construction of a sequence-tagged high-density genetic map of papaya for comparative structural and evolutionary genomics in brassicales

A high-density genetic map of papaya (Carica papaya L.) was constructed using microsatellite markers derived from BAC end sequences and whole-genome shot gun sequences. Fifty-four F(2) plants derived from varieties AU9 and SunUp were used for linkage mapping. A total of 707 markers, including 706 microsatellite loci and the morphological marker fruit flesh color, were mapped into nine major and three minor linkage groups. The resulting map spanned 1069.9 cM with an average distance of 1.5 cM between adjacent markers. This sequence-based microsatellite map resolved the very large linkage group 2 (LG 2) of the previous high-density map using amplified fragment length polymorphism markers. The nine major LGs of our map represent papaya's haploid nine chromosomes with LG 1 of the sex chromosome being the largest. This map validates the suppression of recombination at the male-specific region of the Y chromosome (MSY) mapped on LG 1 and at potential centromeric regions of other LGs. Segregation distortion was detected in a large region on LG 1 surrounding the MSY region due to the abortion of the YY genotype and in a region of LG6 due to an unknown cause. This high-density sequence-tagged genetic map is being used to integrate genetic and physical maps and to assign genome sequence scaffolds to papaya chromosomes. It provides a framework for comparative structural and evolutional genomic research in the order Brassicales.     

Molecular characterization of a bovine Y-specific DNA sequence conserved in taurine and zebu breeds.

The identification of new bovine male-specific DNA sequences is of great interest because the bovine Y chromosome remains poorly characterized in terms of physical and genetic maps. Since taurine and zebu Y chromosomes are structurally different, the identification of Y-specific sequences present in both sub-species is particularly important: these sequences are of evolutionary significance and can be broadly used for embryo sexing. In this work, we initially used the random amplified polymorphic DNA (RAPD) technique to search for male-specific sequences present as monomorphic markers in genomic DNA from zebu and taurine bulls. A male-specific RAPD band was found to be present and highly conserved in both sub-species, as demonstrated by Southern blotting, fluorescent in situ hybridization (FISH) and DNA sequencing. In a previous work, a pair of primers derived from this marker was successfully used in taurine and zebu embryo sexing.     

Functional features of a single chromosome arm in wheat (1AL) determined from its structure

Bread wheat (Triticum aestivum L.) is one of the most important crops globally and a high priority for genetic improvement, but its large and complex genome has been seen as intractable to whole genome sequencing. Isolation of individual wheat chromosome arms has facilitated large-scale sequence analyses. However, so far there is no such survey of sequences from the A genome of wheat. Greater understanding of an A chromosome could facilitate wheat improvement and future sequencing of the entire genome. We have constructed BAC library from the long arm of T. aestivum chromosome 1A (1AL) and obtained BAC end sequences from 7,470 clones encompassing the arm. We obtained 13,445 (89.99%) useful sequences with a cumulative length of 7.57 Mb, representing 1.43% of 1AL and about 0.14% of the entire A genome. The GC content of the sequences was 44.7%, and 90% of the chromosome was estimated to comprise repeat sequences, while just over 1% encoded expressed genes. From the sequence data, we identified a large number of sites suitable for development of molecular markers (362 SSR and 6,948 ISBP) which will have utility for mapping this chromosome and for marker assisted breeding. From 44 putative ISBP markers tested 23 (52.3%) were found to be useful. The BAC end sequence data also enabled the identification of genes and syntenic blocks specific to chromosome 1AL, suggesting regions of particular functional interest and targets for future research.     

Accumulation of chloroplast DNA sequences on the Y chromosome of Silene latifolia

Silene latifolia is a model dioecious plant with heteromorphic sex chromosomes. The Y chromosome is the largest in this species. Theoretical models propose an accumulation of repetitive DNA sequences in non-recombining parts of the Y chromosome. In this study, we isolated a BAC7H5 clone preferentially hybridizing to the Y chromosome of S. latifolia. Sequence analysis revealed that this BAC7H5 contains part of the chloroplast genome, indicating that these chloroplast sequences have accumulated on the Y chromosome and also may contribute to its large size. We constructed Y chromosome- and X chromosome-specific libraries and screened them to find Y- and/or X-linked copies of chloroplast sequences. Sequence analysis revealed higher divergence of a non-genic region of the chloroplast sequences located on the Y chromosome while genic regions tested showed only very low (max 0.9%) divergence from their chloroplast homologues.     

Comparative genome analysis of the yellow fever mosquito Aedes aegypti with Drosophila melanogaster and the malaria vector mosquito Anopheles gambiae

An in silico comparative genomics approach was used to identify putative orthologs to genetically mapped genes from the mosquito, Aedes aegypti, in the Drosophila melanogaster and Anopheles gambiae genome databases. Comparative chromosome positions of 73 D. melanogaster orthologs indicated significant deviations from a random distribution across each of the five A. aegypti chromosomal regions, suggesting that some ancestral chromosome elements have been conserved. However, the two genomes also reflect extensive reshuffling within and between chromosomal regions. Comparative chromosome positions of A. gambiae orthologs indicate unequivocally that A. aegypti chromosome regions share extensive homology to the five A. gambiae chromosome arms. Whole-arm or near-whole-arm homology was contradicted with only two genes among the 75 A. aegypti genes for which orthologs to A. gambiae were identified. The two genomes contain large conserved chromosome segments that generally correspond to break/fusion events and a reciprocal translocation with extensive paracentric inversions evident within. Only very tightly linked genes are likely to retain conserved linear orders within chromosome segments. The D. melanogaster and A. gambiae genome databases therefore offer limited potential for comparative positional gene determinations among even closely related dipterans, indicating the necessity for additional genome sequencing projects with other dipteran species.     

Generation of a 5.5-Mb BAC/PAC contig of pig chromosome 6q1.2 and its integration with existing RH, genetic and comparative maps

We generated a sequence-ready BAC/PAC contig spanning approximately 5.5 Mb on porcine chromosome 6q1.2, which represents a very gene-rich genome region. STS content mapping was used as the main strategy for the assembly of the contig and a total of 6 microsatellite markers, 53 gene-related STS and 116 STS corresponding to BAC and PAC end sequences were analyzed. The contig comprises 316 BAC and PAC clones covering the region between the genes GPI and LIPE. The correct contig assembly was verified by RH-mapping of STS markers and comparative mapping of BAC/PAC end sequences using BLAST searches. The use of microsatellite primer pairs allowed the integration of the physical maps with the genetic map of this region. Comparative mapping of the porcine BAC/PAC contig with respect to the gene-rich region on the human chromosome 19q13.1 map revealed a completely conserved gene order of this segment, however, physical distances differ somewhat between HSA19q13.1 and SSC6q1.2. Three major differences in DNA content between human and pig are found in two large intergenic regions and in one region of a clustered gene family, respectively. While there is a complete conservation of gene order between pig and human, the comparative analysis with respect to the rodent species mouse and rat shows one breakpoint where a genome segment is inverted.     

BAC-FISH in wheat identifies chromosome landmarks consisting of different types of transposable elements

Fluorescence in situ hybridization (FISH) has been widely used in the physical mapping of genes and chromosome landmarks in plants and animals. Bacterial artificial chromosomes (BACs) contain large inserts making them amenable for FISH mapping. We used BAC-FISH to study genome organization and evolution in hexaploid wheat and its relatives. We selected 56 restriction fragment length polymorphism (RFLP) locus-specific BAC clones from libraries of Aegilops tauschii (the D-genome donor of hexaploid wheat) and A-genome diploid Triticum monococcum. Different types of repetitive sequences were identified using BAC-FISH. Two BAC clones gave FISH patterns similar to the repetitive DNA family pSc119; one BAC clone gave a FISH pattern similar to the repetitive DNA family pAs1. In addition, we identified several novel classes of repetitive sequences: one BAC clone hybridized to the centromeric regions of wheat and other cereal species, except rice; one BAC clone hybridized to all subtelomeric chromosome regions in wheat, rye, barley and oat; one BAC clone contained a localized tandem repeat and hybridized to five D-genome chromosome pairs in wheat; and four BAC clones hybridized only to a proximal region in the long arm of chromosome 4A of hexaploid wheat. These repeats are valuable markers for defined chromosome regions and can also be used for chromosome identification. Sequencing results revealed that all these repeats are transposable elements (TEs), indicating the important role of TEs, especially retrotransposons, in genome evolution of wheat.Communicated by P.B. Moens     

Construction of random sheared fosmid library from Chinese cabbage and its use for Brassica rapa genome sequencing project

As a part of the Multinational Genome Sequencing Project of Brassica rapa, linkage group R9 and R3 were sequenced using a bacterial artificial chromosome (BAC) by BAC strategy. The current physical contigs are expected to cover approximately 90% euchromatins of both chromosomes. As the project progresses, BAC selection for sequence extension becomes more limited because BAC libraries are restriction enzyme-specific. To support the project, a random sheared fosmid library was constructed. The library consists of 97536 clones with average insert size of approximately 40 kb corresponding to seven genome equivalents, assuming a Chinese cabbage genome size of 550 Mb. The library was screened with primers designed at the end of sequences of nine points of scaffold gaps where BAC clones cannot be selected to extend the physical contigs. The selected positive clones were end-sequenced to check the overlap between the fosmid clones and the adjacent BAC clones.Nine fosmid clones were selected and fully sequenced. The sequences revealed two completed gap filling and seven sequence extensions, which can be used for further selection of BAC clones confirming that the fosmid library will facilitate the sequence completion of B. rapa.     

The molecular basis of insecticide resistance in mosquitoes

Insecticide resistance is an inherited characteristic involving changes in one or more insect gene. The molecular basis of these changes are only now being fully determined, aided by the availability of the Drosophila melanogaster and Anopheles gambiae genome sequences. This paper reviews what is currently known about insecticide resistance conferred by metabolic or target site changes in mosquitoes.     

甘蓝型油菜A基因组候选BAC克隆的筛选及其插入片段大小分析

选择甘蓝型油菜A基因组10个连锁群上特有的85个SSR分子标记,合成其引物序列,采用四维PCR法筛选甘蓝型油菜-新疆野生油菜二体异附加系BAC文库,成功筛选到甘蓝型油菜A基因组39个BAC克隆,其插入片段介于50～300 kb之间,平均为120 kb。甘蓝型油菜A基因组10个连锁群BAC克隆的获得,对后续开展甘蓝型油菜A基因组染色体识别、基因染色体定位、遗传距离与物理距离间关系分析等均具有重要价值。     

Sex determination in cattle based on simultaneous amplification of a new male-specific DNA sequence and an autosomal locus using the same primers

A PCR-based method for sex determination of bovine DNA samples and embryo biopsies is presented. Using only one primer pair both the male-specific sequence FBNY (127 bp) and a sex-independent control PCR-fragment, the microsatellite marker FBN17 (136-140 bp) are generated in the same PCR reaction. Synteny mapping assigned the male-specific sequence to bovine chromosome Y (BTA Y), whereas FBN17 was mapped to bovine chromosome 2. Localisation of FBNY on BTA Y was confirmed by fluorescence in hybridisation of two BAC clones containing the male-specific sequence. There was no amplification of the male-specific target sequence FBNY in sheep, pig, goat, mice, man, and several wild species of the tribe Bovini. The bovine male-specific fragment was detected in dilutions containing as little as 10 pg genomic DNA and in blastomeres from embryo biopsies. The PCR assay presented here does require neither restriction endonuclease digestion of the PCR product nor additional nested PCR steps. Owing to the advantage of parallel amplification of the autosomal locus FBN17 no additional control fragment is necessary to detect PCR failure. The results of sex determination in embryo biopsies using FBNY were in agreement with the outcome from a reference assay used in commercial breeding programs.     

Gene Finding on the Y: Fruitful Strategy in Drosophila does not Deliver in Anopheles

The Anopheles gambiae genome project yielded almost complete sequences for the autosomes and for a large part of the X chromosome, however, no information for the Y chromosome was obtained. Yet, by design, fragmented Y chromosome sequences should be present in the resulting assembly. Here we report the search for Anopheles Y chromosome genes using a strategy successfully applied for identification of Y genes in Drosophila. A complete set of the unmapped scaffolds was targeted in a broad TBLASTN search using both A. gambiae predicted genes and all proteins from nr database as query sequences. After filtering of the BLAST report, we selected 181 scaffolds possibly containing fragments of Y chromosome genes to experimentally test their Y-linkage. Surprisingly, none of the tested sequences appeared to originate from the Y chromosome. Several factors could account for the failure to detect Y genes, including their different organization in A. gambiae compared to Drosophila and the suboptimal quality of the assembly and annotation of the Anopheles genome. Regardless of the cause, our results illuminate problems associated with the genome analysis of outbred organisms.