How do SNP ascertainment schemes and population demographics affect inferences about population history?

Background

The selection of variable sites for inclusion in genomic analyses can influence results, especially when exemplar populations are used to determine polymorphic sites. We tested the impact of ascertainment bias on the inference of population genetic parameters using empirical and simulated data representing the three major continental groups of cattle: European, African, and Indian. We simulated data under three demographic models. Each simulated data set was subjected to three ascertainment schemes: (I) random selection; (II) geographically biased selection; and (III) selection biased toward loci polymorphic in multiple groups. Empirical data comprised samples of 25 individuals representing each continental group. These cattle were genotyped for 47,506 loci from the bovine 50 K SNP panel. We compared the inference of population histories for the empirical and simulated data sets across different ascertainment conditions using F_ST and principal components analysis (PCA).

Results

Bias toward shared polymorphism across continental groups is apparent in the empirical SNP data. Bias toward uneven levels of within-group polymorphism decreases estimates of F_ST between groups. Subpopulation-biased selection of SNPs changes the weighting of principal component axes and can affect inferences about proportions of admixture and population histories using PCA. PCA-based inferences of population relationships are largely congruent across types of ascertainment bias, even when ascertainment bias is strong.

Conclusions

Analyses of ascertainment bias in genomic data have largely been conducted on human data. As genomic analyses are being applied to non-model organisms, and across taxa with deeper divergences, care must be taken to consider the potential for bias in ascertainment of variation to affect inferences. Estimates of F_ST, time of separation, and population divergence as estimated by principal components analysis can be misleading if this bias is not taken into account.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1469-5) contains supplementary material, which is available to authorized users.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1595-0) contains supplementary material, which is available to authorized users.     

Geographic distribution of haplotype diversity at the bovine casein locus

The genetic diversity of the casein locus in cattle was studied on the basis of haplotype analysis. Consideration of recently described genetic variants of the casein genes which to date have not been the subject of diversity studies, allowed the identification of new haplotypes. Genotyping of 30 cattle breeds from four continents revealed a geographically associated distribution of haplotypes, mainly defined by frequencies of alleles at CSN1S1 and CSN3. The genetic diversity within taurine breeds in Europe was found to decrease significantly from the south to the north and from the east to the west. Such geographic patterns of cattle genetic variation at the casein locus may be a result of the domestication process of modern cattle as well as geographically differentiated natural or artificial selection. The comparison of African Bos taurus and Bos indicus breeds allowed the identification of several Bos indicus specific haplotypes (CSN1S1*C-CSN2*A²-CSN3*A_I/CSN3*H) that are not found in pure taurine breeds. The occurrence of such haplotypes in southern European breeds also suggests that an introgression of indicine genes into taurine breeds could have contributed to the distribution of the genetic variation observed.     

Genomic copy number variation in Mus musculus

Background

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

Results

We found 9634 putative autosomal CNVs across the samples affecting 6.87 % of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

Conclusion

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1713-z) contains supplementary material, which is available to authorized users.     

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-934) contains supplementary material, which is available to authorized users.     

On the association of common and rare genetic variation influencing body mass index: a combined SNP and CNV analysis

Background

As the architecture of complex traits incorporates a widening spectrum of genetic variation, analyses integrating common and rare variation are needed. Body mass index (BMI) represents a model trait, since common variation shows robust association but accounts for a fraction of the heritability. A combined analysis of single nucleotide polymorphisms (SNP) and copy number variation (CNV) was performed using 1850 European and 498 African-Americans from the Study of Addiction: Genetics and Environment. Genetic risk sum scores (GRSS) were constructed using 32 BMI-validated SNPs and aggregate-risk methods were compared: count versus weighted and proxy versus imputation.

Results

The weighted SNP-GRSS constructed from imputed probabilities of risk alleles performed best and was highly associated with BMI (p = 4.3×10⁻¹⁶) accounting for 3% of the phenotypic variance. In addition to BMI-validated SNPs, common and rare BMI/obesity-associated CNVs were identified from the literature. Of the 84 CNVs previously reported, only 21-kilobase deletions on 16p12.3 showed evidence for association with BMI (p = 0.003, frequency = 16.9%), with two CNVs nominally associated with class II obesity, 1p36.1 duplications (OR = 3.1, p = 0.009, frequency 1.2%) and 5q13.2 deletions (OR = 1.5, p = 0.048, frequency 7.7%). All other CNVs, individually and in aggregate, were not associated with BMI or obesity. The combined model, including covariates, SNP-GRSS, and 16p12.3 deletion accounted for 11.5% of phenotypic variance in BMI (3.2% from genetic effects). Models significantly predicted obesity classification with maximum discriminative ability for morbid-obesity (p = 3.15×10⁻¹⁸).

Conclusion

Results show that incorporating validated effect sizes and allelic probabilities improve prediction algorithms. Although rare-CNVs did not account for significant phenotypic variation, results provide a framework for integrated analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-368) contains supplementary material, which is available to authorized users.     

Genomic comparative analysis of the environmental Enterococcus mundtii against enterococcal representative species

Background

Enterococcus mundtii is a yellow-pigmented microorganism rarely found in human infections. The draft genome sequence of E. mundtii was recently announced. Its genome encodes at least 2,589 genes and 57 RNAs, and 4 putative genomic islands have been detected. The objective of this study was to compare the genetic content of E. mundtii with respect to other enterococcal species and, more specifically, to identify genes coding for putative virulence traits present in enterococcal opportunistic pathogens.

Results

An in-depth mining of the annotated genome was performed in order to uncover the unique properties of this microorganism, which allowed us to detect a gene encoding the antimicrobial peptide mundticin among other relevant features. Moreover, in this study a comparative genomic analysis against commensal and pathogenic enterococcal species, for which genomic sequences have been released, was conducted for the first time. Furthermore, our study reveals significant similarities in gene content between this environmental isolate and the selected enterococci strains (sharing an “enterococcal gene core” of 805 CDS), which contributes to understand the persistence of this genus in different niches and also improves our knowledge about the genetics of this diverse group of microorganisms that includes environmental, commensal and opportunistic pathogens.

Conclusion

Although E. mundtii CRL1656 is phylogenetically closer to E. faecium, frequently responsible of nosocomial infections, this strain does not encode the most relevant relevant virulence factors found in the enterococcal clinical isolates and bioinformatic predictions indicate that it possesses the lowest number of putative pathogenic genes among the most representative enterococcal species. Accordingly, infection assays using the Galleria mellonella model confirmed its low virulence.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-489) contains supplementary material, which is available to authorized users.     

Genomic variation in tomato,from wild ancestors to contemporary breeding accessions

Background

Domestication modifies the genomic variation of species. Quantifying this variation provides insights into the domestication process, facilitates the management of resources used by breeders and germplasm centers, and enables the design of experiments to associate traits with genes. We described and analyzed the genetic diversity of 1,008 tomato accessions including Solanum lycopersicum var. lycopersicum (SLL), S. lycopersicum var. cerasiforme (SLC), and S. pimpinellifolium (SP) that were genotyped using 7,720 SNPs. Additionally, we explored the allelic frequency of six loci affecting fruit weight and shape to infer patterns of selection.

Results

Our results revealed a pattern of variation that strongly supported a two-step domestication process, occasional hybridization in the wild, and differentiation through human selection. These interpretations were consistent with the observed allele frequencies for the six loci affecting fruit weight and shape. Fruit weight was strongly selected in SLC in the Andean region of Ecuador and Northern Peru prior to the domestication of tomato in Mesoamerica. Alleles affecting fruit shape were differentially selected among SLL genetic subgroups. Our results also clarified the biological status of SLC. True SLC was phylogenetically positioned between SP and SLL and its fruit morphology was diverse. SLC and “cherry tomato” are not synonymous terms. The morphologically-based term “cherry tomato” included some SLC, contemporary varieties, as well as many admixtures between SP and SLL. Contemporary SLL showed a moderate increase in nucleotide diversity, when compared with vintage groups.

Conclusions

This study presents a broad and detailed representation of the genomic variation in tomato. Tomato domestication seems to have followed a two step-process; a first domestication in South America and a second step in Mesoamerica. The distribution of fruit weight and shape alleles supports that domestication of SLC occurred in the Andean region. Our results also clarify the biological status of SLC as true phylogenetic group within tomato. We detect Ecuadorian and Peruvian accessions that may represent a pool of unexplored variation that could be of interest for crop improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1444-1) contains supplementary material, which is available to authorized users.     

Nitrogen and water addition reduce leaf longevity of steppe species

Background and aims

Changes in supplies of resources will modify plant functional traits. However, few experimental studies have addressed the effects of nitrogen and water variations, either singly or in combination, on functional traits.

Methods

A 2-year field experiment was conducted to test the effects of nitrogen and water addition on leaf longevity and other functional traits of the two dominant (Agropyron cristatum and Stipa krylovii) and three most common species (Cleistogenes squarrosa, Melilotoides ruthenica and Potentilla tanacetifolia) in a temperate steppe in northern China.

Key Results

Additional nitrogen and water increased leaf nitrogen content and net photosynthetic rate, and changed other measured functional traits. Leaf longevity decreased significantly with both nitrogen addition (–6 days in 2007 and –5·4 days in 2008; both P < 0·001) and watering (–13 days in 2007 and –9·9 days in 2008; both P < 0·001), and significant differences in leaf longevity were also found among species. Nitrogen and water interacted to affect leaf longevity and other functional traits. Soil water content explained approx. 70 % of the shifts in leaf longevity. Biomass at both species and community level increased under water and nitrogen addition because of the increase in leaf biomass production per individual plant.

Conclusions

The results suggest that additional nitrogen and water supplies reduce plant leaf longevity. Soil water availability might play a fundamental role in determining leaf longevity and other leaf functional traits, and its effects can be modified by soil nitrogen availability in semi-arid areas. The different responses of species to resource alterations may cause different global change ramifications under future climate change scenarios.     

How slug herbivory of juvenile hybrid willows alters chemistry,growth and subsequent susceptibility to diverse plant enemies

Background and Aims

Selective feeding by herbivores, especially at the seedling or juvenile phase, has the potential to change plant traits and ultimately the susceptibility of surviving plants to other enemies. Moreover, since hybridization is important to speciation and can lead to introgression of traits between plant species, differential feeding (herbivore-induced mortality) can influence the expression of resistance traits of hybrids and ultimately determine the consequences of hybridization. While it would be expected that herbivore-induced mortality would lead to greater resistance, there may be trade-offs whereby resistance to one herbivore increases susceptibility to others. The hypothesis was tested that the exotic slug, Arion subfuscus, causes non-random survival of hybrid willows and alters plant: (1) susceptibility to slugs; (2) secondary and nutritional chemistry, and growth; and (3) susceptibility to other phytophages.

Methods

Two populations of plants, control and selected, were created by placing trays of juvenile willows in the field and allowing slugs access to only some. When ≤10 individuals/tray remained (approx. 85 % mortality), ‘selected’ and undamaged ‘control’ trays were returned to a common area. Traits of these populations were then examined in year 1 and in subsequent years.

Key Results

The selected population was less palatable to slugs. Surprisingly, foliar concentrations of putative defence traits (phenolic glycosides and tannins) did not differ between treatments, but the selected population had higher foliar nitrogen and protein, lower carbon to nitrogen ratio and greater above-ground biomass, indicating that vigorously growing plants were inherently more resistant to slugs. Interestingly, selected plants were more susceptible to three phytophages: an indigenous pathogen (Melampsora epitea), a native herbivorous beetle (Chrysomela knabi) and an exotic willow leaf beetle (Plagiodera versicolora).

Conclusions

This exotic slug changed the population structure of F₂ hybrid willows in unanticipated ways. Defence expression remained unchanged, while nutritional and growth traits changed. These changes caused plants to be more susceptible to other plant enemies. Other exotic herbivore species are anticipated to have similar direct and indirect effects on native plant populations.     

On the Origin of Indonesian Cattle

Background

Two bovine species contribute to the Indonesian livestock, zebu (Bos indicus) and banteng (Bos javanicus), respectively. Although male hybrid offspring of these species is not fertile, Indonesian cattle breeds are supposed to be of mixed species origin. However, this has not been documented and is so far only supported by preliminary molecular analysis.

Methods and Findings

Analysis of mitochondrial, Y-chromosomal and microsatellite DNA showed a banteng introgression of 10–16% in Indonesian zebu breeds. East-Javanese Madura and Galekan cattle have higher levels of autosomal banteng introgression (20–30%) and combine a zebu paternal lineage with a predominant (Madura) or even complete (Galekan) maternal banteng origin. Two Madura bulls carried taurine Y-chromosomal haplotypes, presumably of French Limousin origin. In contrast, we did not find evidence for zebu introgression in five populations of the Bali cattle, a domestic form of the banteng.

Conclusions

Because of their unique species composition Indonesian cattle represent a valuable genetic resource, which potentially may also be exploited in other tropical regions.     

Virulence related sequences; insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background

Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains.

Results

Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection.

Conclusion

The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1512-6) contains supplementary material, which is available to authorized users.     

Novel adaptation to hawkmoth pollinators in Clarkia reduces efficiency,not attraction of diurnal visitors

Background and Aims

Plant populations experiencing divergent pollination environments may be under selection to modify floral traits in ways that increase both attractiveness to and efficiency of novel pollinators. These changes may come at the cost of reducing overall effectiveness of other pollinators. The goal of this study was to examine differences in attractiveness and efficiency between Clarkia concinna and C. breweri, sister species of annual plants with parapatric distributions.

Methods

An assessment was made as to whether observed differences in visitors between natural populations are driven by differences in floral traits or differences in the local pollination environment. Differences in floral attractiveness were quantified by setting out arrays of both species in the geographical range of each species and exposing both species to nocturnal hawkmoths (Hyles lineata) in flight cages. Differences in visitor efficiency were estimated by measuring stigma–visitor contact frequency and pollen loads for diurnal visitors, and pollen deposition on stigmas for hawkmoths.

Key Results

The composition of visitors to arrayed plants was similar between plant species at any particular site, but highly divergent among sites, and reflected differences in visitors to natural populations. Diurnal insects visited both species, but were more common at C. concinna populations. Hummingbirds and hawkmoths were only observed visiting within the range of C. breweri. Despite attracting similar species when artificially presented together, C. concinna and C. breweri showed large differences in pollinator efficiency. All visitors except hawkmoths pollinated C. concinna more efficiently.

Conclusions

Differences in the available pollinator community may play a larger role than differences in floral traits in determining visitors to natural populations of C. concinna and C. breweri. However, floral traits mediate differences in pollinator efficiency. Increased effectiveness of the novel hawkmoth pollinator on C. breweri comes at relatively little cost in attractiveness to other visitors, but at large cost in their efficiency as pollinators.     

Quantitative imaging of rhizosphere pH and CO2 dynamics with planar optodes

Background and Aims

Live imaging methods have become extremely important for the exploration of biological processes. In particular, non-invasive measurement techniques are key to unravelling organism–environment interactions in close-to-natural set-ups, e.g. in the highly heterogeneous and difficult-to-probe environment of plant roots: the rhizosphere. pH and CO₂ concentration are the main drivers of rhizosphere processes. Being able to monitor these parameters at high spatio-temporal resolution is of utmost importance for relevant interpretation of the underlying processes, especially in the complex environment of non-sterile plant–soil systems. This study introduces the application of easy-to-use planar optode systems in different set-ups to quantify plant root–soil interactions.

Methods

pH- and recently developed CO₂-sensors were applied to rhizobox systems to investigate roots with different functional traits, highlighting the potential of these tools. Continuous and highly resolved real-time measurements were made of the pH dynamics around Triticum turgidum durum (durum wheat) roots, Cicer arietinum (chickpea) roots and nodules, and CO₂ dynamics in the rhizosphere of Viminaria juncea.

Key Results

Wheat root tips acidified slightly, while their root hair zone alkalized their rhizosphere by more than 1 pH unit and the effect of irrigation on soil pH could be visualized as well. Chickpea roots and nodules acidified the surrounding soil during N₂ fixation and showed diurnal changes in acidification activity. A growing root of V. juncea exhibited a large zone of influence (mm) on soil CO₂ content and therefore on its biogeochemical surrounding, all contributing to the extreme complexity of the root–soil interactions.

Conclusions

This technique provides a unique tool for future root research applications and overcomes limitations of previous systems by creating quantitative maps without, for example, interpolation and time delays between single data points.