THE GENOMIC TRAJECTORY OF HYBRID SWARMS: OUTCOMES OF REPEATED CROSSES BETWEEN POPULATIONS OF TIGRIOPUS CALIFORNICUS

Introgressive hybridization between genetically divergent populations is an important evolutionary process. The degree to which repeated hybridization events between the same parental taxa lead to similar genomic outcomes is unknown. This study addressed this question by following genomic trajectories of replicate hybrid swarms of the copepod Tigriopus californicus over many generations of free mating. Swarm composition was determined both by differential reproductive success of founder individuals and subsequent selection on hybrid genotypes. For one cross, between two populations showing differential fitness in the laboratory and no hybrid breakdown, the genetic trajectory was highly repeatable: replicates rapidly became dominated by alleles from the fitter parent. In a second cross, between two populations showing similar fitness and significant F2 hybrid breakdown, alleles from alternative populations dominated different replicates. Swarms exhibited a general temporal trend of decreasing cytonuclear mismatch. Some patterns of differential introgression across the genome were strikingly congruent amongst swarm replicates, both within and between cross types, and reflected patterns of segregation distortion previously observed within controlled crosses between the same parental populations. Differences in heterozygosity between the sexes, and evidence for a previously suspected sex‐distortion locus, suggest that complex interactions between sex and genotype influence hybrid swarm outcome.     

Genetic Mapping and Genomic Selection Using Recombination Breakpoint Data

The correct models for quantitative trait locus mapping are the ones that simultaneously include all significant genetic effects. Such models are difficult to handle for high marker density. Improving statistical methods for high-dimensional data appears to have reached a plateau. Alternative approaches must be explored to break the bottleneck of genomic data analysis. The fact that all markers are located in a few chromosomes of the genome leads to linkage disequilibrium among markers. This suggests that dimension reduction can also be achieved through data manipulation. High-density markers are used to infer recombination breakpoints, which then facilitate construction of bins. The bins are treated as new synthetic markers. The number of bins is always a manageable number, on the order of a few thousand. Using the bin data of a recombinant inbred line population of rice, we demonstrated genetic mapping, using all bins in a simultaneous manner. To facilitate genomic selection, we developed a method to create user-defined (artificial) bins, in which breakpoints are allowed within bins. Using eight traits of rice, we showed that artificial bin data analysis often improves the predictability compared with natural bin data analysis. Of the eight traits, three showed high predictability, two had intermediate predictability, and two had low predictability. A binary trait with a known gene had predictability near perfect. Genetic mapping using bin data points to a new direction of genomic data analysis.     

Analysis of Allele-Specific Expression in Mouse Liver by RNA-Seq: A Comparison With Cis-eQTL Identified Using Genetic Linkage

We report an analysis of allele-specific expression (ASE) and parent-of-origin expression in adult mouse liver using next generation sequencing (RNA-Seq) of reciprocal crosses of heterozygous F₁ mice from the parental strains C57BL/6J and DBA/2J. We found a 60% overlap between genes exhibiting ASE and putative cis-acting expression quantitative trait loci (cis-eQTL) identified in an intercross between the same strains. We discuss the various biological and technical factors that contribute to the differences. We also identify genes exhibiting parental imprinting and complex expression patterns. Our study demonstrates the importance of biological replicates to limit the number of false positives with RNA-Seq data.     

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

The number of markers and samples needed for detecting bottlenecks under realistic scenarios,with and without recovery: a simulation‐based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection methods exist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. ( 2012 ) analysed the power of two approaches, the M‐ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100), which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when postbottleneck recovery has occurred. For the realistic situations examined, we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals, (ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii) recovery after a moderate amount of time or gradual change in size reduces power, by up to one‐half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated, using simulation modelling, as a routine component of molecular ecology studies.     

Molecularly imprinted cryogels for human interferon‐alpha purification from human gingival fibroblast culture

Interferons are important proteins for the immune system because of their antiviral, anti‐proliferating and immunomodulatory activities. Therapeutic value of these proteins against certain types of tumors caused interest and investigations aimed to obtain highly purified interferons. Molecular imprinting is an efficient method for purification with high selectivity, specificity and good reproducibility. In this study, we utilized advantages of molecular imprinting technique for the purification of interferon from human gingival fibroblast culture. For this purpose, interferon α‐2b imprinted poly(hydroxyethyl methacrylate) cryogel (hIFN‐α‐MIP) was prepared. Optimum adsorption conditions were determined, and maximum adsorption capacity of hIFN‐α‐MIP cryogel was found as 254.8 × 10⁴ IU/g from aqueous solution. All interferon measurements are expressed as International Unit (IU), which is a unit measurement used to quantify biologically active substances like interferon based on their biological activity or effect. Selectivity experiments were performed using competitive proteins and repeated adsorption–desorption studies showed that the adsorption capacity maintained almost at a constant value after ten cycles. For the purification of interferon from human gingival fibroblast culture, fast protein liquid chromatography was used and the specific activity of the purified interferon α‐2b on HeLa cell line was found between the values 3.45 × 10⁸ IU/mg and 3.75 × 10⁸ IU/mg. The results are promising, and the molecular imprinting technique is effective for the purification of interferon α‐2b. Copyright © 2013 John Wiley & Sons, Ltd.