Mb Of Sequence Research Articles

Cognitive, Behavioral, and Adaptive Profiles in Williams Syndrome With and Without Loss of GTF2IRD2.

The present study aimed to identify the role of GTF2IRD2 in the cognitive, behavioral, and adaptive profile of WS patients. Twelve patients diagnosed with WS participated, four with GTF2IRD2 deletion (atypical WS group), and eight without this deletion (typical WS group). The age range of both groups was 7-18 years old. Each patient's 7q11.23 deletion scope was determined by chromosomal microarray analysis. Cognitive, behavioral, and adaptive abilities were assessed with a battery of neuropsychological tests. Compared with the typical WS group, the atypical WS patients with GTF2IRD2 deletion had more impaired visuospatial abilities and more significant behavioral problems, mainly related to the construct of social cognition. These findings provide new evidence regarding the influence of the GTF2IRD2 gene on the severity of behavioral symptoms of WS related to social cognition and certain visuospatial abilities. (JINS, 2018, 24, 896-904).

Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity.

BackgroundAlthough draft genomes are available for most agronomically important plant species, the majority are incomplete, highly fragmented, and often riddled with assembly and scaffolding errors. These assembly issues hinder advances in tool development for functional genomics and systems biology.FindingsHere we utilized a robust, cost-effective approach to produce high-quality reference genomes. We report a near-complete genome of diploid woodland strawberry (Fragaria vesca) using single-molecule real-time sequencing from Pacific Biosciences (PacBio). This assembly has a contig N50 length of ∼7.9 million base pairs (Mb), representing a ∼300-fold improvement of the previous version. The vast majority (>99.8%) of the assembly was anchored to 7 pseudomolecules using 2 sets of optical maps from Bionano Genomics. We obtained ∼24.96 Mb of sequence not present in the previous version of the F. vesca genome and produced an improved annotation that includes 1496 new genes. Comparative syntenic analyses uncovered numerous, large-scale scaffolding errors present in each chromosome in the previously published version of the F. vesca genome.ConclusionsOur results highlight the need to improve existing short-read based reference genomes. Furthermore, we demonstrate how genome quality impacts commonly used analyses for addressing both fundamental and applied biological questions.

Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers

Genomic and genetic variation among six Italian chicken native breeds (Livornese, Mericanel della Brianza, Milanino, Bionda Piemontese, Bianca di Saluzzo and Siciliana) were studied using single nucleotide polymorphism (SNP) and copy number variants (CNV) as markers. A total of 94 DNA samples genotyped with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix) were used in the analyses. The results showed the genetic and genomic variability occurring among the six Italian chicken breeds. The genetic relationship among animals was established with a principal component analysis. The genetic diversity within breeds was calculated using heterozygosity values (expected and observed) and with Wright’s F-statistics. The individual-based CNV calling, based on log R ratio and B-allele frequency values, was done by the Hidden–Markov Model (HMM) of PennCNV software on autosomes. A hierarchical agglomerative clustering was applied in each population according to the absence or presence of definite CNV regions (CNV were grouped by overlapping of at least 1 bp). The CNV map was built on a total of 1003 CNV found in individual samples, after grouping by overlaps, resulting in 564 unique CNV regions (344 gains, 213 losses and 7 complex), for a total of 9.43 Mb of sequence and 1.03% of the chicken assembly autosome. All the approaches using SNP data showed that the Siciliana breed clearly differentiate from other populations, the Livornese breed separates into two distinct groups according to the feather colour (i.e. white and black) and the Bionda Piemontese and Bianca di Saluzzo breeds are closely related. The genetic variability found using SNP is comparable with that found by other authors in the same breeds using microsatellite markers. The CNV markers analysis clearly confirmed the SNP results.

A New Chicken Genome Assembly Provides Insight into Avian Genome Structure.

The importance of the Gallus gallus (chicken) as a model organism and agricultural animal merits a continuation of sequence assembly improvement efforts. We present a new version of the chicken genome assembly (Gallus_gallus-5.0; GCA_000002315.3), built from combined long single molecule sequencing technology, finished BACs, and improved physical maps. In overall assembled bases, we see a gain of 183 Mb, including 16.4 Mb in placed chromosomes with a corresponding gain in the percentage of intact repeat elements characterized. Of the 1.21 Gb genome, we include three previously missing autosomes, GGA30, 31, and 33, and improve sequence contig length 10-fold over the previous Gallus_gallus-4.0. Despite the significant base representation improvements made, 138 Mb of sequence is not yet located to chromosomes. When annotated for gene content, Gallus_gallus-5.0 shows an increase of 4679 annotated genes (2768 noncoding and 1911 protein-coding) over those in Gallus_gallus-4.0. We also revisited the question of what genes are missing in the avian lineage, as assessed by the highest quality avian genome assembly to date, and found that a large fraction of the original set of missing genes are still absent in sequenced bird species. Finally, our new data support a detailed map of MHC-B, encompassing two segments: one with a highly stable gene copy number and another in which the gene copy number is highly variable. The chicken model has been a critical resource for many other fields of study, and this new reference assembly will substantially further these efforts.

Deep sequencing of 10,000 human genomes

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Anchored pseudo-de novo assembly of human genomes identifies extensive sequence variation from unmapped sequence reads.

The human genome reference (HGR) completion marked the genomics era beginning, yet despite its utility universal application is limited by the small number of individuals used in its development. This is highlighted by the presence of high-quality sequence reads failing to map within the HGR. Sequences failing to map generally represent 2-5% of total reads, which may harbor regions that would enhance our understanding of population variation, evolution, and disease. Alternatively, complete de novo assemblies can be created, but these effectively ignore the groundwork of the HGR. In an effort to find a middle ground, we developed a bioinformatic pipeline that maps paired-end reads to the HGR as separate single reads, exports unmappable reads, de novo assembles these reads per individual and then combines assemblies into a secondary reference assembly used for comparative analysis. Using 45 diverse 1000 Genomes Project individuals, we identified 351,361 contigs covering 195.5Mb of sequence unincorporated in GRCh38. 30,879 contigs are represented in multiple individuals with ~40% showing high sequence complexity. Genomic coordinates were generated for 99.9%, with 52.5% exhibiting high-quality mapping scores. Comparative genomic analyses with archaic humans and primates revealed significant sequence alignments and comparisons with model organism RefSeq gene datasetsidentified novel human genes. If incorporated, these sequences will expand the HGR, but more importantly our data highlight that with this method low coverage (~10-20×) next-generation sequencing can still be used to identify novel unmapped sequences to explore biological functions contributing to human phenotypic variation, disease and functionality for personal genomic medicine.

Chromoanasynthetic Genomic Rearrangement Identified in a N-Ethyl-N-Nitrosourea (ENU) Mutagenesis Screen in Caenorhabditis elegans.

Chromoanasynthesis is a recently discovered phenomenon in humans with congenital diseases that is characterized by complex genomic rearrangements (CGRs) resulting from aberrant repair of catastrophic chromosomal damage. How these CGRs are induced is not known. Here, we describe the structure and function of dpDp667, a causative CGR that emerged from a Caenorhabditis elegans dauer suppressor screen in which animals were treated with the point mutagen N-ethyl-N-nitrosourea (ENU). dpDp667 comprises nearly 3 Mb of sequence on the right arm of the X chromosome, contains three duplications and one triplication, and is devoid of deletions. Sequences from three out of the four breakpoint junctions in dpDp667 reveal microhomologies that are hallmarks of chromoanasynthetic CGRs. Our findings suggest that environmental insults and physiological processes that cause point mutations may give rise to chromoanasynthetic rearrangements associated with congenital disease. The relatively subtle phenotype of animals harboring dpDp667 suggests that the prevalence of CGRs in the genomes of mutant and/or phenotypically unremarkable animals may be grossly underestimated.

SNP Assay Development for Linkage Map Construction, Anchoring Whole-Genome Sequence, and Other Genetic and Genomic Applications in Common Bean.

A total of 992,682 single-nucleotide polymorphisms (SNPs) was identified as ideal for Illumina Infinium II BeadChip design after sequencing a diverse set of 17 common bean (Phaseolus vulgaris L) varieties with the aid of next-generation sequencing technology. From these, two BeadChips each with >5000 SNPs were designed. The BARCBean6K_1 BeadChip was selected for the purpose of optimizing polymorphism among market classes and, when possible, SNPs were targeted to sequence scaffolds in the Phaseolus vulgaris 14× genome assembly with sequence lengths >10 kb. The BARCBean6K_2 BeadChip was designed with the objective of anchoring additional scaffolds and to facilitate orientation of large scaffolds. Analysis of 267 F2 plants from a cross of varieties Stampede × Red Hawk with the two BeadChips resulted in linkage maps with a total of 7040 markers including 7015 SNPs. With the linkage map, a total of 432.3 Mb of sequence from 2766 scaffolds was anchored to create the Phaseolus vulgaris v1.0 assembly, which accounted for approximately 89% of the 487 Mb of available sequence scaffolds of the Phaseolus vulgaris v0.9 assembly. A core set of 6000 SNPs (BARCBean6K_3 BeadChip) with high genotyping quality and polymorphism was selected based on the genotyping of 365 dry bean and 134 snap bean accessions with the BARCBean6K_1 and BARCBean6K_2 BeadChips. The BARCBean6K_3 BeadChip is a useful tool for genetics and genomics research and it is widely used by breeders and geneticists in the United States and abroad.

Structural features of conopeptide genes inferred from partial sequences of the Conus tribblei genome.

The evolvability of venom components (in particular, the gene-encoded peptide toxins) in venomous species serves as an adaptive strategy allowing them to target new prey types or respond to changes in the prey field. The structure, organization, and expression of the venom peptide genes may provide insights into the molecular mechanisms that drive the evolution of such genes. Conus is a particularly interesting group given the high chemical diversity of their venom peptides, and the rapid evolution of the conopeptide-encoding genes. Conus genomes, however, are large and characterized by a high proportion of repetitive sequences. As a result, the structure and organization of conopeptide genes have remained poorly known. In this study, a survey of the genome of Conus tribblei was undertaken to address this gap. A partial assembly of C. tribblei genome was generated; the assembly, though consisting of a large number of fragments, accounted for 2160.5Mb of sequence. A large number of repetitive genomic elements consisting of 642.6Mb of retrotransposable elements, simple repeats, and novel interspersed repeats were observed. We characterized the structural organization and distribution of conotoxin genes in the genome. A significant number of conopeptide genes (estimated to be between 148 and 193) belonging to different superfamilies with complete or nearly complete exon regions were observed, ~60% of which were expressed. The unexpressed conopeptide genes represent hidden but significant conotoxin diversity. The conotoxin genes also differed in the frequency and length of the introns. The interruption of exons by long introns in the conopeptide genes and the presence of repeats in the introns may indicate the importance of introns in facilitating recombination, evolution and diversification of conotoxins. These findings advance our understanding of the structural framework that promotes the gene-level molecular evolution of venom peptides.

Whole-genome sequencing of a Plasmodium vivax isolate from the China-Myanmar border area

Currently, there is a trend of an increasing number of Plasmodium vivaxmalaria cases in China that are imported across its Southeast Asia border, especially in the China-Myanmar border area (CMB). To date, little is known about the genetic diversity of P. vivax in this region. In this paper, we report the first genome sequencing of a P. vivaxisolate (CMB-1) from a vivax malaria patient in CMB. The sequencing data were aligned onto 96.43% of the P. vivax Salvador I reference strain (Sal I) genome with 7.84-fold coverage as well as onto 98.32% of 14 Sal I chromosomes. Using the de novo assembly approach, we generated 8,541 scaffolds and assembled a total of 27.1 Mb of sequence into CMB-1 scaffolds. Furthermore, we identified all 295 known virgenes, which is the largest subtelomeric multigene family in malaria parasites. These results provide an important foundation for further research onP. vivax population genetics.

Analysis of the genomic DNA of the harmful dinoflagellate Prorocentrum minimum: a brief survey focused on the noncoding RNA gene sequences

Dinoflagellate algae are considered some of the most complicated organisms owing to their unusual genomic characteristics and novel gene regulatory mechanisms. Here, we extracted information from the genomic sequences of the marine dinoflagellate Prorocentrum minimum using GS FLX pyrosequencing technology. We obtained 473 Mb of sequences from 1,379,588 reads, which generated 16,599 contigs after assembly. Among the annotated sequences (4902 contigs, 30 %), BLAST analyses showed that 28.95 % (4731 contigs) of the genome fragments were homologous with bacterial sequences and only 0.96 % (156 contigs) had eukaryotic origins. However, analysis of bacterial 16S/23S rDNA sequences in the P. minimum genome revealed that the organism likely acquired this genetic material from symbiotic bacteria, possibly through horizontal gene transfer, thus demonstrating a definitive association between bacteria and P. minimum. Moreover, a specific consensus pattern was defined for dinoflagellate spliced leader (SL) sequences, a well-known genetic marker of dinoflagellates. Further, comparisons of the various genomic proportions with respect to nuclear, mitochondrial, and chloroplast genes revealed unique genomic features that include clusters of noncoding RNA genes and tandem repeats. Taken together, this study details the salient features of noncoding RNA genes in P. minimum and provides further insights into the genome of this and other dinoflagellates.

Integration of genetic and physical maps of the Primula vulgaris S locus and localization by chromosome in situ hybridization.

Summary Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S‐linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location.We have employed a combination of classical genetics and three‐point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization.We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5‐Mb genomic region with 1 Mb of sequence containing 82 S‐linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair.These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

The Y-chromosome point mutation rate in humans

Mutations are the fundamental source of biological variation, and their rate is a crucial parameter for evolutionary and medical studies. Here we used whole-genome sequence data from 753 Icelandic males, grouped into 274 patrilines, to estimate the point mutation rate for 21.3 Mb of male-specific Y chromosome (MSY) sequence, on the basis of 1,365 meioses (47,123 years). The combined mutation rate for 15.2 Mb of X-degenerate (XDG), X-transposed (XTR) and ampliconic excluding palindromes (rAMP) sequence was 8.71 × 10(-10) mutations per position per year (PPPY). We observed a lower rate (P = 0.04) of 7.37 × 10(-10) PPPY for 6.1 Mb of sequence from palindromes (PAL), which was not statistically different from the rate of 7.2 × 10(-10) PPPY for paternally transmitted autosomes. We postulate that the difference between PAL and the other MSY regions may provide an indication of the rate at which nascent autosomal and PAL de novo mutations are repaired as a result of gene conversion.

Genome-wide characteristics of copy number variation in Polish Holstein and Polish Red cattle using SNP genotyping assay.

Copy number variation (CNV), which results from deletions or amplifications of large fragments of genomic DNA, is widespread in mammalian genomes and apart from its potential pathogenic effect it is considered as a source of natural genetic diversity. In cattle populations, this kind of genetic variability remains still insufficiently elucidated and studies focusing on the detection of new structural genomic variants in different cattle populations may contribute to a better understanding of cattle breeds' diversity and genetic basis of production traits. In this study, by using BovineSNP50 assay and cnvPartition algorithm we identified CNVs in two different cattle breeds: Holstein (859 animals) and Polish Red (301). In Holstein cattle we found 648 CNVs which could be reduced to 91 non-redundant variable genomic regions (CNVRs) covering in total 168.6Mb of the genomic sequence. In Polish Red cattle we detected 62 CNVs, localized in 37 variable regions encompassing 22.3Mb of the sequence, corresponding to 0.89% of the autosomal genome. Within the regions we identified 1,192 unique RefSeq genes which are engaged in a variety of biological processes. High concordance of the regions' distribution was found between the studied breeds, however copy number variants seemed to be more common in Holstein cattle. About 26% of the regions described in this study could be classified as newly identified. The results of this study will broaden the knowledge of CNVs in genomes of cattle of different breeds and will provide foundations for further research aiming to identify a relationship between this type of genetic variation and phenotypic traits.

Use of targeted SNP selection for an improved anchoring of the melon (Cucumis melo L.) scaffold genome assembly.

BackgroundThe genome of the melon (Cucumis melo L.) double-haploid line DHL92 was recently sequenced, with 87.5 and 80.8% of the scaffold assembly anchored and oriented to the 12 linkage groups, respectively. However, insufficient marker coverage and a lack of recombination left several large, gene rich scaffolds unanchored, and some anchored scaffolds unoriented. To improve the anchoring and orientation of the melon genome assembly, we used resequencing data between the parental lines of DHL92 to develop a new set of SNP markers from unanchored scaffolds.ResultsA high-resolution genetic map composed of 580 SNPs was used to anchor 354.8 Mb of sequence, contained in 141 scaffolds (average size 2.5 Mb) and corresponding to 98.2% of the scaffold assembly, to the 12 melon chromosomes. Over 325.4 Mb (90%) of the assembly was oriented. The genetic map revealed regions of segregation distortion favoring SC alleles as well as recombination suppression regions coinciding with putative centromere, 45S, and 5S rDNA sites. New chromosome-scale pseudomolecules were created by incorporating to the previous v3.5 version an additional 38.3 Mb of anchored sequence representing 1,837 predicted genes contained in 55 scaffolds. Using fluorescent in situ hybridization (FISH) with BACs that produced chromosome-specific signals, melon chromosomes that correspond to the twelve linkage groups were identified, and a standardized karyotype of melon inbred line T111 was developed.ConclusionsBy utilizing resequencing data and targeted SNP selection combined with a large F2 mapping population, we significantly improved the quantity of anchored and oriented melon scaffold genome assembly. Using genome information combined with FISH mapping provided the first cytogenetic map of an inodorus melon type. With these results it was possible to make inferences on melon chromosome structure by relating zones of recombination suppression to centromeres and 45S and 5S heterochromatic regions. This study represents the first steps towards the integration of the high-resolution genetic and cytogenetic maps with the genomic sequence in melon that will provide more information on genome organization and allow for the improvement of the melon genome draft sequence.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-014-1196-3) contains supplementary material, which is available to authorized users.

Recurrent Loss of Sex Is Associated with Accumulation of Deleterious Mutations in Oenothera

Sexual reproduction is nearly universal among eukaryotes. Theory predicts that the rarity of asexual eukaryotic species is in part caused by accumulation of deleterious mutations and heightened extinction risk associated with suppressed recombination and segregation in asexual species. We tested this prediction with a large data set of 62 transcriptomes from 29 species in the plant genus Oenothera, spanning ten independent transitions between sexual and a functionally asexual genetic system called permanent translocation heterozygosity. Illumina short-read sequencing and de novo transcript assembly yielded an average of 16.4 Mb of sequence per individual. Here, we show that functionally asexual species accumulate more deleterious mutations than sexual species using both population genomic and phylogenetic analysis. At an individual level, asexual species exhibited 1.8 × higher heterozygosity than sexual species. Within species, we detected a higher proportion of nonsynonymous polymorphism relative to synonymous variation within asexual compared with sexual species, indicating reduced efficacy of purifying selection. Asexual species also exhibited a greater proportion of transcripts with premature stop codons. The increased proportion of nonsynonymous mutations was also positively correlated with divergence time between sexual and asexual species, consistent with Muller's ratchet. Between species, we detected repeated increases in the ratio of nonsynonymous to synonymous divergence in asexual species compared with sexually reproducing sister taxa, indicating increased accumulation of deleterious mutations. These results confirm that an important advantage of sex is that it facilitates selection against deleterious alleles, which might help to explain the dearth of extant asexual species.

A SNaPshot of next generation sequencing for forensic SNP analysis.

Forensic phenotyping can provide useful intelligence regarding the biogeographical ancestry (BGA) and externally visible characteristics (EVCs) of the donor of an evidentiary sample. Currently, single nucleotide polymorphism (SNP) based inference of BGA and EVCs is performed most commonly using SNaPshot(®), a single base extension (SBE) assay. However, a single SNaPshot multiplex PCR is limited to 30-40 SNPs. Next generation sequencing (NGS) offers the potential to genotype hundreds to thousands of SNPs from multiple samples in a single experimental run. The PCR multiplexes from five SNaPshot assays (SNPforID 52plex, SNPforID 34plex, Eurasiaplex, IrisPlex and an unpublished BGA assay) were applied to three different DNA template amounts (0.1, 0.2 and 0.3 ng) in three samples (9947A and 007 control DNAs and a male donor). The pooled PCR amplicons containing 136 unique SNPs were sequenced using Life Technologies' Ion Torrent™ PGM system. Approximately 72 Mb of sequence was generated from two 10 Mb Ion 314™ v1 chips. Accurate genotypes were readily obtained from all three template amounts. Of a total of 408 genotypes, 395 (97%) were fully concordant with SNaPshot across all three template amounts. Of those genotypes discordant with SNaPshot, six Ion Torrent sequences (1.5%) were fully concordant with Sanger sequencing across the three template amounts. Seven SNPs (1.7%) were either discordant between template amounts or discordant with Sanger sequencing. Sequence coverage observed in the negative control, and, allele coverage variation for heterozygous genotypes highlights the need to establish a threshold for background levels of sequence output and heterozygous balance. This preliminary study of the Ion Torrent PGM system has demonstrated considerable potential for use in forensic DNA analyses as a low to medium throughput NGS platform using established SNaPshot assays.

Recently published Streptomyces genome sequences.

Recently published Streptomyces genome sequences.

A reference genome for common bean and genome-wide analysis of dual domestications.

Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement.

Palaeosymbiosis revealed by genomic fossils of Wolbachia in a strongyloidean nematode.

Wolbachia are common endosymbionts of terrestrial arthropods, and are also found in nematodes: the animal-parasitic filaria, and the plant-parasite Radopholus similis. Lateral transfer of Wolbachia DNA to the host genome is common. We generated a draft genome sequence for the strongyloidean nematode parasite Dictyocaulus viviparus, the cattle lungworm. In the assembly, we identified nearly 1 Mb of sequence with similarity to Wolbachia. The fragments were unlikely to derive from a live Wolbachia infection: most were short, and the genes were disabled through inactivating mutations. Many fragments were co-assembled with definitively nematode-derived sequence. We found limited evidence of expression of the Wolbachia-derived genes. The D. viviparus Wolbachia genes were most similar to filarial strains and strains from the host-promiscuous clade F. We conclude that D. viviparus was infected by Wolbachia in the past, and that clade F-like symbionts may have been the source of filarial Wolbachia infections.