Brent W. Murray

Nonlinkage of major histocompatibility complex class I and class II loci in bony fishes.

Abstract In tetrapods, the functional (classical) class I and class II B loci of the major histocompatibility complex (Mhc) are tightly linked in a single chromosomal region. In an earlier study, we demonstrated that in the zebrafish, Danio rerio, order Cypriniformes, the two classes are present on different chromosomes. Here, we show that the situation is similar in the stickleback, Gasterosteus aculeatus, order Gasterosteiformes, the common guppy, Poecilia reticulata, order Cyprinodontiformes, and the cichlid fish Oreochromis niloticus, order Perciformes. These data, together with unpublished results from other laboratories suggest that in all Euteleostei, the classical class I and class II B loci are in separate linkage groups, and that in at least some of these taxa, the class II loci are in two different groups. Since Euteleostei are at least as numerous as tetrapods, in approximately one-half of jawed vertebrates, the class I and class II regions are not linked.

A Contig Map of the Mhc Class I Genomic Region in the Zebrafish Reveals Ancient Synteny

In contrast to the human and mouse Mhc, in which the clusters of class I and class II loci reside in close vicinity to one another, in the zebrafish, Danio rerio, they are found in different linkage groups. Chromosome walking using BAC (bacterial artificial chromosome) and PAC (P1 artificial chromosome) clones reveals the zebrafish class I region to occupy a segment of ∼450 kb and to encompass at least 19 loci. These include three class I (Dare-UDA, -UEA, -UFA), five proteasome subunit β (PSMB8, -9A, -9C, -11, -12), two TAPs (TAP2A, TAP2B), and one TAP binding protein (TAPBP). This arrangement contrasts with the arrangements found in human and mouse Mhc, in which the orthologues of the PSMB, TAP, and TAPBP loci reside within the class II region. In addition to this main zebrafish class I contig, a shorter contig of about 150 kb contains two additional class I (UBA, UCA) and at least five other loci. It probably represents a different haplotype of part of the class I region. The previously identified UAA gene shares an identical 5′ part with UEA, but the two genes differ in their 3′ parts. One of them is probably the result of an unequal crossing over. The described organization has implications for the persistence of syntenic relationships, coevolution of loci, and interpretation of the origin of the human/mouse Mhc organization.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 April 2010-31 May 2010

This article documents the addition of 396 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Anthocidaris crassispina, Aphis glycines, Argyrosomus regius, Astrocaryum sciophilum, Dasypus novemcinctus, Delomys sublineatus, Dermatemys mawii, Fundulus heteroclitus, Homalaspis plana, Jumellea rossii, Khaya senegalensis, Mugil cephalus, Neoceratitis cyanescens, Phalacrocorax aristotelis, Phytophthora infestans, Piper cordulatum, Pterocarpus indicus, Rana dalmatina, Rosa pulverulenta, Saxifraga oppositifolia, Scomber colias, Semecarpus kathalekanensis, Stichopus monotuberculatus, Striga hermonthica, Tarentola boettgeri and Thermophis baileyi. These loci were cross‐tested on the following species: Aphis gossypii, Sooretamys angouya, Euryoryzomys russatus, Fundulus notatus, Fundulus olivaceus, Fundulus catenatus, Fundulus majalis, Jumellea fragrans, Jumellea triquetra Jumellea recta, Jumellea stenophylla, Liza richardsonii, Piper marginatum, Piper aequale, Piper darienensis, Piper dilatatum, Rana temporaria, Rana iberica, Rana pyrenaica, Semecarpus anacardium, Semecarpus auriculata, Semecarpus travancorica, Spondias acuminata, Holigarna grahamii, Holigarna beddomii, Mangifera indica, Anacardium occidentale, Tarentola delalandii, Tarentola caboverdianus and Thermophis zhaoermii.

Sequence variation at the major histocompatibility complex DRB loci in beluga (Delphinapterus leucas) and narwhal (Monodon monoceros).

Abstract The variation at loci with similarity to DRB class II major histocompatibility complex loci was assessed in 313 beluga collected from 13 sampling locations across North America, and 11 narwhal collected in the Canadian high Arctic. Variation was assessed by amplification of exon 2, which codes for the peptide binding region, via the polymerase chain reaction, followed by either cloning and DNA sequencing or single-stranded conformation polymorphism analysis. Two DRB loci were identified in beluga: DRB1, a polymorphic locus, and, DRB2, a monomorphic locus. Eight alleles representing five distinct lineages (based on sequence similarity) were found at the beluga DRB1 locus. Although the relative number of alleles is low when compared with terrestrial mammals, the amino acid variation found among the lineages is moderate. At the DRB1 locus, the average number of nonsynonymous substitutions per site is greater than the average number of synonymous substitutions per site (0.0806 : 0.0207, respectively;P<0.01). Most of the 31 amino acid substitutions do not conserve the physiochemical properties of the residue, and 21 of these are located at positions implicated as forming pockets responsible for the selective binding of foreign peptide side chains. Only DRB1 variation was examined in 11 narwhal, revealing a low amount of variation. These data are consistent with an important role for the DRB1 locus in the cellular immune response of beluga. In addition, the ratio of nonsynonymous to synonymous substitutions is similar to that among primate alleles, arguing against a reduction in the balancing selection pressure in the marine environment. Two hypotheses may explain the modest amount of Mhc variation when compared with terrestrial mammals: small population sizes at speciation or a reduced neutral substitution rate in cetaceans.

Evolution of glutamine synthetase in vertebrates: multiple glutamine synthetase genes expressed in rainbow trout (Oncorhynchus mykiss).

SUMMARY Glutamine synthetase (GSase) is a key enzyme in nitrogen metabolism and encoded by a single gene in mammals. Using PCR cloning techniques, including RT-PCR from total RNA and PCR from a cDNA library, we find evidence of four expressed GSase mRNAs for the tetraploid rainbow trout. For two of these mRNAs (Onmy-GS01, -GS02) we characterize the full-length coding regions, and for two others (Onmy-GS03, -GS04), we describe partial sequences. Northern analysis of Onmy-GS01, -GS02, -GS03 and -GS04 indicates that (1) Onmy-GS02 is expressed at higher levels relative to the other transcripts in most adult tissues, with the exception of brain and gill, where Onmy-GS01 is at the highest level, and (2) the tissue with the highest level of expression of all four transcripts is the brain, with decreasing levels in the intestine, liver, red muscle, gill/kidney, white muscle and heart. Clearly, rainbow trout possess multiple GSase genes with differing levels of tissue expression, implying manifold potential routes of regulation for this octameric enzyme. Our data also indicate that caution should be taken when interpreting mRNA expression data of a single gene, unless multiple genes have been ruled out. Consistent with a southern blot, phylogenetic and intron sequence analyses imply that the trout genes are encoded by at least four separate loci, belonging to two distinct evolutionary branches. Our data on rainbow trout, together with those from two full-length zebrafish Danio rerio GSase genes compiled from GenBank ESTs, support the idea that fish GSases are polyphyletic and that gene duplications have occurred at multiple points and in independent lineages throughout the evolution of bony fishes.

Spatial Genetic Structure of a Symbiotic Beetle-Fungal System: Toward Multi-Taxa Integrated Landscape Genetics

Spatial patterns of genetic variation in interacting species can identify shared features that are important to gene flow and can elucidate co-evolutionary relationships. We assessed concordance in spatial genetic variation between the mountain pine beetle (Dendroctonus ponderosae) and one of its fungal symbionts, Grosmanniaclavigera, in western Canada using neutral genetic markers. We examined how spatial heterogeneity affects genetic variation within beetles and fungi and developed a novel integrated landscape genetics approach to assess reciprocal genetic influences between species using constrained ordination. We also compared landscape genetic models built using Euclidean distances based on allele frequencies to traditional pair-wise Fst. Both beetles and fungi exhibited moderate levels of genetic structure over the total study area, low levels of structure in the south, and more pronounced fungal structure in the north. Beetle genetic variation was associated with geographic location while that of the fungus was not. Pinevolume and climate explained beetle genetic variation in the northern region of recent outbreak expansion. Reciprocal genetic relationships were only detectedin the south where there has been alonger history of beetle infestations. The Euclidean distance and Fst-based analyses resulted in similar models in the north and over the entire study area, but differences between methods in the south suggest that genetic distances measures should be selected based on ecological and evolutionary contexts. The integrated landscape genetics framework we present is powerful, general, and can be applied to other systems to quantify the biotic and abiotic determinants of spatial genetic variation within and among taxa.

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north‐south population structure as evidenced by: (i) Bayesian‐based analyses, (ii) north‐south genetic relationships and diversity gradients; and (iii) a lack of isolation‐by‐distance in the northernmost cluster. The north‐south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate‐driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range‐wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 December 2010-31 January 2011: PERMANENT GENETIC RESOURCES NOTE

This article documents the addition of 238 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Alytes dickhilleni, Arapaima gigas, Austropotamobius italicus, Blumeria graminis f. sp. tritici, Cobitis lutheri, Dendroctonus ponderosae, Glossina morsitans morsitans, Haplophilus subterraneus, Kirengeshoma palmata, Lysimachia japonica, Macrolophus pygmaeus, Microtus cabrerae, Mytilus galloprovincialis, Pallisentis (Neosentis) celatus, Pulmonaria officinalis, Salminus franciscanus, Thais chocolata and Zootoca vivipara. These loci were cross‐tested on the following species: Acanthina monodon, Alytes cisternasii, Alytes maurus, Alytes muletensis, Alytes obstetricans almogavarii, Alytes obstetricans boscai, Alytes obstetricans obstetricans, Alytes obstetricans pertinax, Cambarellus montezumae, Cambarellus zempoalensis, Chorus giganteus, Cobitis tetralineata, Glossina fuscipes fuscipes, Glossina pallidipes, Lysimachia japonica var. japonica, Lysimachia japonica var. minutissima, Orconectes virilis, Pacifastacus leniusculus, Procambarus clarkii, Salminus brasiliensis and Salminus hilarii.

Major histocompatibility complex class II A genes in cichlid fishes: identification, expression, linkage relationships, and haplotype variation.

Abstract Two cichlid species, the haplochromine Aulonocara hansbaenschi and the tilapiine Oreochromis niloticus, were used to study the major histocompatibility complex (Mhc) class II A variation within this group. Multiple class II A sequences were recovered from A. hansbaenschi and O. niloticus cDNA libraries and three sequence families, DAA, DBA, and DCA, were identified. Sets of O. niloticus haploid embryo families were used to determine the linkage relationships of these genes. Two independently assorting linkage groups were detected, DAA and DBA/DCA, neither of which is linked to the previously described Mhc class I gene cluster. Three DCA genes and up to four DBA genes were found to segregate in different haplotypes, whereas DAA occurred as a single locus. Four DBA haplotypes, DBA*H1-H4, were identified and shown to co-segregate with the previously described class II B haplotypes. Four DCA haplotypes, DCA*H1-H4, were found at a distance of 37 cM from the DBA/class II B cluster; in one DCA haplotype, DCA*H5, the genes were tightly linked to the DBA/class II B clusters. Transcripts of DAA and DBA genes were found in O. niloticus hepatopancreas and spleen; transcripts of DCA genes were detected in the A. hansbaenschi cDNA library, but not in O. niloticus. These findings provide a basis for using class II haplotypes as markers in the study of adaptive radiation in the cichlid species flocks of the East African Great Lakes.

Sexual reproduction and gene flow in the pine pathogen Dothistroma septosporum in British Columbia.

Dothistroma septosporum has caused a serious needle blight epidemic in the lodgepole pine forests in northwest British Columbia over the past several years. Although ascocarps had been observed in British Columbia, nothing was known about the contribution of sexual reproduction, gene flow and long-distance dispersal to the epidemic. Amplified fragment length polymorphism and mating-type markers in 19 sites were used to generate population and reproductive data. Overall, evidence suggests a mixed mode of reproduction. Haplotypic diversity was high, with 79 unique and 56 shared haplotypes (possible clones) identified from 192 fungal isolates. Overall, mating-type segregation did not differ significantly from 1:1; however, random mating was rejected in most populations in the index of association and parsimony tree-length permutation analyses using the full data set and, when using clone-corrected data sets, more of the smaller populations showed random mating. Two of the smaller populations consistently showed random mating for both tests using both clone-corrected and noncorrected data. High gene flow is suggested by no differentiation between 14 of the 19 sites, several of which came from young plantations where the pathogen was not likely present prior to the current outbreak. The remaining five sites showed some level of divergence, possibly due to historic separation and endemic pathogen populations. Results indicate a high evolutionary potential and long-distance dispersal in this pathogen, important to consider in future forest management.