Mating barriers between genetically divergent strains of the parasitic nematode Haemonchus contortus suggest incipient speciation

Haemonchus contortus, in common with many nematode species, has extremely high levels of genetic variation within and between field populations derived from distant geographical locations. MHco10(CAVR), MHco3(ISE) and MHco4(WRS) are genetically divergent H. contortus strains, originally derived from Australia, Kenya and South Africa, respectively, that have been maintained by numerous rounds of in vivo experimental infection of sheep. In order to explore potential pre-zygotic competition or post-zygotic incompatibility between the strains, we have investigated the ability of MHco10(CAVR) to interbreed with either MHco3(ISE) or MHco4(WRS) during dual strain co-infections. Sheep were experimentally co-infected with 4000 infective larvae (L₃) per os of the MHco10(CAVR) strain and an equal number of either the MHco3(ISE) or the MHco4(WRS) strain L₃. The adult worm establishement rates and the proportions of F₁ progeny resulting from intra- and inter-strain mating events were determined by admixture analysis of microsatellite multi-locus genotypes. Although there was no difference in adult worm establishment rates, the proportions of F₁ progeny of both the MHco10(CAVR) × MHco3(ISE) and MHco10(CAVR) × MHco4(WRS) dual strain co-infections departed from Mendelian expectations. The proportions of inter-strain hybrid F₁ progeny were lower than the expected 50%, suggesting either pre-zygotic competition or post-zygotic incompatibility between the co-infecting strains. To investigate this further, both eggs and hatched L₁ of broods from single adult female worms recovered from each dual co-infection were genotyped. Unhatched eggs from the broods revealed no inter-strain hybrid genotype deficit, suggesting there is no pre-zygotic competition between the strains. In contrast, there was a deficit in L₁ inter-strain hybrid genotypes in the broods derived from MHco3(ISE) or MHco4(WRS) maternal parents, but not from MHco10(CAVR) maternal parents. This suggests that hybrid progeny of MHco10(CAVR) paternal parents have reduced post-zygotic development and/or viability consistent with incipient speciation of the MHco10(CAVR) strain. The presence of mating barriers between allopatric H. contortus strains has important implications for parasite ecology, including the ability of newly introduced anthelmintic-resistant parasite populations to compete and interbreed with populations already established in a region.     

A method for single pair mating in an obligate parasitic nematode

Parasitic nematode species have extremely high levels of genetic diversity, presenting a number of experimental challenges for genomic and genetic work. Consequently, there is a need to develop inbred laboratory strains with reduced levels of polymorphism. The most efficient approach to inbred line development is single pair mating, but this is challenging for obligate parasites where the adult sexual reproductive stages are inside the host, and thus difficult to experimentally manipulate. This paper describes a successful approach to single pair mating of a parasitic nematode, Haemonchus contortus. The method allows for polyandrous mating behaviour and involves the surgical transplantation of a single adult male worm with multiple immature adult females directly into the sheep abomasum. We used a panel of microsatellite markers to monitor and validate the single pair mating crosses and to ensure that the genotypes of progeny and subsequent filial generations were consistent with those expected from a mating between a single female parent of known genotype and a single male parent of unknown genotype. We have established two inbred lines that both show a significant overall reduction in genetic diversity based on microsatellite genotyping and genome-wide single nucleotide polymorphism. There was an approximately 50% reduction in heterozygous SNP sites across the genome in the MHco3.N1 line compared with the MoHco3(ISE) parental strain. The MHco3.N1 inbred line has subsequently been used to provide DNA template for whole genome sequencing of H. contortus. This work provides proof of concept and methodologies for forward genetic analysis of obligate parasitic nematodes.     

Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans

Background

Understanding how DNA sequence polymorphism relates to variation in gene expression is essential to connecting genotypic differences with phenotypic differences among individuals. Addressing this question requires linking population genomic data with gene expression variation.

Results

Using whole genome expression data and recent light shotgun genome sequencing of six Drosophila simulans genotypes, we assessed the relationship between expression variation in males and females and nucleotide polymorphism across thousands of loci. By examining sequence polymorphism in gene features, such as untranslated regions and introns, we find that genes showing greater variation in gene expression between genotypes also have higher levels of sequence polymorphism in many gene features. Accordingly, X-linked genes, which have lower sequence polymorphism levels than autosomal genes, also show less expression variation than autosomal genes. We also find that sex-specifically expressed genes show higher local levels of polymorphism and divergence than both sex-biased and unbiased genes, and that they appear to have simpler regulatory regions.

Conclusion

The gene-feature-based analyses and the X-to-autosome comparisons suggest that sequence polymorphism in cis-acting elements is an important determinant of expression variation. However, this relationship varies among the different categories of sex-biased expression, and trans factors might contribute more to male-specific gene expression than cis effects. Our analysis of sex-specific gene expression also shows that female-specific genes have been overlooked in analyses that only point to male-biased genes as having unusual patterns of evolution and that studies of sexually dimorphic traits need to recognize that the relationship between genetic and expression variation at these traits is different from the genome as a whole.     

Pathogenesis gene families in the common minimal genome of Staphylococcus aureus are hypervariable

Staphylococcus aureus is a versatile pathogen that shows high levels of inter-strain genetic variability and positive evolution in certain pathogenesis-related genes. Apart from gene content differences, variability in shared genes may affect pathogenicity. Studying such variability requires that the common minimal genome (CMG) be identified. In this study, we have surveyed the CMG of S. aureus with respect to variability amongst orthologous family members, and determined that genes involved in pathogenesis preferentially accumulate variations. A negative correlation between variability of genes and their evolution was found, suggesting a preservation of host-specific function while exhibiting sequence diversity. Variation in key pathogenesis genes in S. aureus might predispose them to functional modulation, thereby playing an important role in evasion of host immunity.     

Gene Expression Variability within and between Human Populations and Implications toward Disease Susceptibility

Variations in gene expression level might lead to phenotypic diversity across individuals or populations. Although many human genes are found to have differential mRNA levels between populations, the extent of gene expression that could vary within and between populations largely remains elusive. To investigate the dynamic range of gene expression, we analyzed the expression variability of ∼18, 000 human genes across individuals within HapMap populations. Although ∼20% of human genes show differentiated mRNA levels between populations, our results show that expression variability of most human genes in one population is not significantly deviant from another population, except for a small fraction that do show substantially higher expression variability in a particular population. By associating expression variability with sequence polymorphism, intriguingly, we found SNPs in the untranslated regions (5′ and 3′UTRs) of these variable genes show consistently elevated population heterozygosity. We performed differential expression analysis on a genome-wide scale, and found substantially reduced expression variability for a large number of genes, prohibiting them from being differentially expressed between populations. Functional analysis revealed that genes with the greatest within-population expression variability are significantly enriched for chemokine signaling in HIV-1 infection, and for HIV-interacting proteins that control viral entry, replication, and propagation. This observation combined with the finding that known human HIV host factors show substantially elevated expression variability, collectively suggest that gene expression variability might explain differential HIV susceptibility across individuals.     

Analysis of Multiple Brachyspira hyodysenteriae Genomes Confirms That the Species Is Relatively Conserved but Has Potentially Important Strain Variation

The intestinal spirochete Brachyspira hyodysenteriae is an important pathogen in swine, causing mucohemorrhagic colitis in a disease known as swine dysentery. Based on the detection of significant linkage disequilibrium in multilocus sequence data, the species is considered to be clonal. An analysis of the genome sequence of Western Australian B. hyodysenteriae strain WA1 has been published, and in the current study 19 further strains from countries around the world were sequenced with Illumina technology. The genomes were assembled and aligned to over 97.5% of the reference WA1 genome at a percentage sequence identity better than 80%. Strain regions not aligned to the reference ranged between 0.2 and 2.5%. Clustering of the strain genes found on average 2,354 (88%) core genes, 255 (8.6%) ancillary genes and 77 (2.9%) unique genes per strain. Depending on the strain the proportion of genes with 100% sequence identity to WA1 ranged from 85% to 20%. The result is a global comparative genomic analysis of B. hyodysenteriae genomes revealing potential differential phenotypic markers for numerous strains. Despite the differences found, the genomes were less varied than those of the related pathogenic species Brachyspira pilosicoli, and the analysis supports the clonal nature of the species. From this study, a public genome resource has been created that will serve as a repository for further genetic and phenotypic studies of these important porcine bacteria. This is the first intra-species B. hyodysenteriae comparative genomic analysis.     

Differential Exoproteome Analysis of Two Corynebacterium pseudotuberculosis Biovar Ovis Strains Isolated from Goat (1002) and Sheep (C231)

Corynebacterium pseudotuberculosis is the etiologic agent of caseous lymphadenitis a chronic infectious disease affecting small ruminants. The 2D-DIGE technique was used to compare the exoproteomes of two C. pseudotuberculosis biovar ovis strains isolated from goat (strain 1002) and sheep (strain C231). Seventeen proteins differentially produced were identified here. Nine proteins appeared over-produced in the exoproteome of 1002 goat strain and 8 in that of C231 sheep strain. These proteins were related to various biological functions, such as the cell envelope, respiratory metabolism and proteolysis. This proteomic analysis revealed strain-specific exoproteins although each of the corresponding genes was found in both strain genomes. Such differential expression pattern may reflect inter-strain differences in adaptation to a specific host, in pathogenicity and or in antigenicity of this pathogenic bacterium.     

Horizontally Acquired Glycosyltransferase Operons Drive Salmonellae Lipopolysaccharide Diversity

The immunodominant lipopolysaccharide is a key antigenic factor for Gram-negative pathogens such as salmonellae where it plays key roles in host adaptation, virulence, immune evasion, and persistence. Variation in the lipopolysaccharide is also the major differentiating factor that is used to classify Salmonella into over 2600 serovars as part of the Kaufmann-White scheme. While lipopolysaccharide diversity is generally associated with sequence variation in the lipopolysaccharide biosynthesis operon, extraneous genetic factors such as those encoded by the glucosyltransferase (gtr) operons provide further structural heterogeneity by adding additional sugars onto the O-antigen component of the lipopolysaccharide. Here we identify and examine the O-antigen modifying glucosyltransferase genes from the genomes of Salmonella enterica and Salmonella bongori serovars. We show that Salmonella generally carries between 1 and 4 gtr operons that we have classified into 10 families on the basis of gtrC sequence with apparent O-antigen modification detected for five of these families. The gtr operons localize to bacteriophage-associated genomic regions and exhibit a dynamic evolutionary history driven by recombination and gene shuffling events leading to new gene combinations. Furthermore, evidence of Dam- and OxyR-dependent phase variation of gtr gene expression was identified within eight gtr families. Thus, as O-antigen modification generates significant intra- and inter-strain phenotypic diversity, gtr-mediated modification is fundamental in assessing Salmonella strain variability. This will inform appropriate vaccine and diagnostic approaches, in addition to contributing to our understanding of host-pathogen interactions.