Genotyping strategy matters when analyzing hypervariable major histocompatibility complex‐Experience from a passerine bird

Genotyping of classical major histocompatibility complex (MHC) genes is challenging when they are hypervariable and occur in multiple copies. In this study, we used several different approaches to genotype the moderately variable MHC class I exon 3 (MHCIe3) and the highly polymorphic MHC class II exon 2 (MHCIIβe2) in the bluethroat (Luscinia svecica). Two family groups (eight individuals) were sequenced in replicates at both markers using Ion Torrent technology with both a single‐ and a dual‐indexed primer structure. Additionally, MHCIIβe2 was sequenced on Illumina MiSeq. Allele calling was conducted by modifications of the pipeline developed by Sommer et al. (BMC Genomics, 14, 2013, 542) and the software AmpliSAS. While the different genotyping strategies gave largely consistent results for MHCIe3, with a maximum of eight alleles per individual, MHCIIβe2 was remarkably complex with a maximum of 56 MHCIIβe2 alleles called for one individual. Each genotyping strategy detected on average 50%–82% of all MHCIIβe2 alleles per individual, but dropouts were largely allele‐specific and consistent within families for each strategy. The discrepancies among approaches indicate PCR biases caused by the platform‐specific primer tails. Further, AmpliSAS called fewer alleles than the modified Sommer pipeline. Our results demonstrate that allelic dropout is a significant problem when genotyping the hypervariable MHCIIβe2. As these genotyping errors are largely nonrandom and method‐specific, we caution against comparing genotypes across different genotyping strategies. Nevertheless, we conclude that high‐throughput approaches provide a major advance in the challenging task of genotyping hypervariable MHC loci, even though they may not reveal the complete allelic repertoire.     

Inference of genetic architecture from chromosome partitioning analyses is sensitive to genome variation,sample size,heritability and effect size distribution

Genomewide association studies have contributed immensely to our understanding of the genetic basis of complex traits. One major conclusion arising from these studies is that most traits are controlled by many loci of small effect, confirming the infinitesimal model of quantitative genetics. A popular approach to test for polygenic architecture involves so‐called “chromosome partitioning” where phenotypic variance explained by each chromosome is regressed on the size of the chromosome. First developed for humans, this has now been repeatedly used in other species, but there has been no evaluation of the suitability of this method in species that can differ in their genome characteristics such as number and size of chromosomes. Nor has the influence of sample size, heritability of the trait, effect size distribution of loci controlling the trait or the physical distribution of the causal loci in the genome been examined. Using simulated data, we show that these characteristics have major influence on the inferences of the genetic architecture of traits we can infer using chromosome partitioning analyses. In particular, small variation in chromosome size, small sample size, low heritability, a skewed effect size distribution and clustering of loci can lead to a loss of power and consequently altered inference from chromosome partitioning analyses. Future studies employing this approach need to consider and derive an appropriate null model for their study system, taking these parameters into consideration. Our simulation results can provide some guidelines on these matters, but further studies examining a broader parameter space are needed.     

What shapes Eurasian lynx distribution in human dominated landscapes: selecting prey or avoiding people?

In the multi-use landscape of southern Norway, the distribution of lynx is likely to be determined both by the abundance of their favoured prey – the roe deer – and the risk associated with the presence of humans because most lynx mortalities are caused by humans (recreational harvest, poaching, vehicle collisions). We described the distribution of the reproductive portion of the lynx population based on snow-track observations of females with dependent kittens collected over 10  yr (1997–2006) in southern Norway. We used the ecological-niche factor analysis to examine how lynx distribution was influenced by roe deer, human activity, habitat type, environmental productivity and elevation. Our first prediction that lynx should be found in areas of relatively high roe deer abundance was supported. However, our second prediction that lynx should avoid human activity was rejected, and lynx instead occupied areas more disturbed in average than those available (with the exception of the most densely occupied areas). Lynx, however, avoided the most disturbed areas and our third prediction of a trade-off between abundance of prey and avoidance of human activity was supported. On the one hand, roe deer in the most disturbed areas benefit to a large extent from current human land use practices, potentially allowing them to escape predation from lynx. On the other hand, the situation is not so favourable for the predators who are restricted in competition refuges with medium to low prey densities. The consequence is that lynx conservation will have to be achieved in a human modifed environment where the potential for a range of conflicts and high human-caused mortality will remain a constant threat.     

Molecular and phenotypic divergence in the bluethroat (Luscinia svecica) subspecies complex

Subspecies complexes may provide valuable insights into the early stages of the speciation process. The bluethroat (Luscinia svecica) consists of many morphologically distinct subspecies that differ most strikingly in the ornamental colour pattern of the male throat. We investigated the genetic and phenotypic differentiation in this subspecies complex, using (i) microsatellite genotyping (11 loci) of a sample of 364 individuals from bluethroat populations in Europe and Asia, and (ii) spectrometric and morphological measurements of a sample of 131 museum skin specimens. Population genetic analyses, based on microsatellite allele frequency variation, revealed a slight but significant overall population differentiation (F(ST) = 0.042). There was a well-differentiated southern group of subspecies with white or no throat spots and a less-differentiated northern group of chestnut-spotted populations. Phylogenetic analyses indicated that the southern all-blue and white-spotted forms are ancestral to the chestnut-spotted subspecies. In addition to the qualitative variation in throat plumage pattern already described in the literature, we found significant quantitative variation among subspecies in hue, chroma and brightness of the ultraviolet (UV)/blue throat coloration, and this variation seemed to be unrelated to the phylogenetic distance between subspecies.     

Genetic mapping of DArT markers in the Festuca�CLolium complex and their use in freezing tolerance association analysis

Species belonging to the Festuca-Lolium complex are important forage and turf species and as such, have been studied intensively. However, their out-crossing nature and limited availability of molecular markers make genetic studies difficult. Here, we report on saturation of F. pratensis and L. multiflorum genetic maps using Diversity Array Technology (DArT) markers and the DArTFest array.The 530 and 149 DArT markers were placed on genetic maps of L. multiflorum and F. pratensis, respectively, with overlap of 20 markers, which mapped in both species. The markers were sequenced and comparative sequence analysis was performed between L. multiflorum, rice and Brachypodium. The utility of the DArTFest array was then tested on a Festulolium population FuRs0357 in an integrated analysis using the DArT marker map positions to study associations between markers and freezing tolerance. Ninety six markers were significantly associated with freezing tolerance and five of these markers were genetically mapped to chromosomes 2, 4 and 7. Three genomic loci associated with freezing tolerance in the FuRs0357 population co-localized with chromosome segments and QTLs previously identified to be associated with freezing tolerance. The present work clearly confirms the potential of the DArTFest array in genetic studies of the Festuca-Lolium complex. The annotated DArTFest array resources could accelerate further studies and improvement of desired traits in Festuca-Lolium species.     

Proteomic identification of secreted proteins from human skeletal muscle cells and expression in response to strength training

Regular physical activity protects against several types of diseases. This may involve altered secretion of signaling proteins from skeletal muscle. Our aim was to identify the most abundantly secreted proteins in cultures of human skeletal muscle cells and to monitor their expression in muscles of strength-training individuals. A total of 236 proteins were detected by proteome analysis in medium conditioned by cultured human myotubes, which was narrowed down to identification of 18 classically secreted proteins expressed in skeletal muscle, using the SignalP 3.0 and Human Genome Expression Profile databases together with a published mRNA-based reconstruction of the human skeletal muscle secretome. For 17 of the secreted proteins, expression was confirmed at the mRNA level in cultured human myotubes as well as in biopsies of human skeletal muscles. RT-PCR analyses showed that 15 of the secreted muscle proteins had significantly enhanced mRNA expression in m. vastus lateralis and/or m. trapezius after 11 wk of strength training among healthy volunteers. For example, secreted protein acidic and rich in cysteine, a secretory protein in the membrane fraction of skeletal muscle fibers, was increased 3- and 10-fold in m. vastus lateralis and m. trapezius, respectively. Identification of proteins secreted by skeletal muscle cells in vitro facilitated the discovery of novel responses in skeletal muscles of strength-training individuals.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

Species‐specific communication bars interspecific mating between syntopic species of Zwicknia stoneflies (Plecoptera: Capniidae)

Most Northern Hemisphere stoneflies have species‐specific mating signals that are generally thought to constitute a barrier against interspecific mating. We tested this hypothesis in two species of the genus Zwicknia that have only very recently been recognised as distinct species, and that were found to occur together in a stream in Lower Saxony, Germany. Analyses of molecular markers COI and 28S in combination with wing length (distinguishing males of both species) and mating signals revealed no instance of hybridisation among 23 studied specimens. In addition, eleven further males identified on the basis of morphology alone all produced the expected species‐specific signal. Females and males of both species were presented with played back conspecific and heterospecific signals and duetting sequences, and responded only to conspecific stimuli. This lends support to the hypothesis that the intersexual communication system functions as an important pre‐mating barrier against gene flow, although post‐mating isolation cannot be excluded. Interspecific mating did occur when a mixed pair was confined together in a small container. Males of both species were found to call in response to played back duetting sequences with stereotypic latencies that are clearly longer than the latencies in male‐female duets. We interpret this as an indication of eavesdropping behaviour coupled with attempts to take over the perceived duet. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 969–980.     

Moose body mass variation revisited: disentangling effects of environmental conditions and genetics

Large-scale geographical variation in phenotypic traits within species is often correlated to local environmental conditions and population density. Such phenotypic variation has recently been shown to also be influenced by genetic structuring of populations. In ungulates, large-scale geographical variation in phenotypic traits, such as body mass, has been related to environmental conditions and population density, but little is known about the genetic influences. Research on the genetic structure of moose suggests two distinct genetic lineages in Norway, structured along a north-south gradient. This corresponds with many environmental gradients, thus genetic structuring provides an additional factor affecting geographical phenotypic variation in Norwegian moose. We investigated if genetic structure explained geographical variation in body mass in Norwegian moose while accounting for environmental conditions, age and sex, and if it captured some of the variance in body mass that previously was attributed to environmental factors. Genetic structuring of moose was the most important variable in explaining the geographic variation in body mass within age and sex classes. Several environmental variables also had strong explanatory power, related to habitat diversity, environmental seasonality and winter harshness. The results suggest that environmental conditions, landscape characteristics, and genetic structure should be evaluated together when explaining large-scale patterns in phenotypic characters or life history traits. However, to better understand the role of genetic and environmental effects on phenotypic traits in moose, an extended individual-based study of variation in fitness-related characters is needed, preferably in an area of convergence between different genetic lineages.     

Synthesis, molecular modeling studies and biological evaluation of fluorine substituted analogs of GW 501516

(±)-2-Fluoro-2-(2-methyl-4-(((4-methyl-2-(4-(trifluoromethyl)phenyl)thiazol-5-yl)methyl)thio)phenoxy)acetic acid (2a) has been prepared and subjected to biological testing against all three subtypes of the PPARs. This compound exhibited agonist effects with EC(50) values of 560 and 55 nM against PPARα and PPARδ, respectively, in a luciferase assay. Moreover, compound (±)-2a also exhibited potent ability to induce oleic acid oxidation in a human myotube cell assay with EC(50)=3.7 nM. Compound (±)-2a can be classified as a dual PPARα/δ agonist with a 10-fold higher potency against the PPARδ receptor than against the PPARα receptor. Molecular modeling studies revealed that both enantiomers of 2a bind to the PPARδ receptor with similar binding energies.