Runs of homozygosity and population history in cattle

ABSTRACT: BACKGROUND: Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring. The extent and frequency of ROHs may inform on the ancestry of an individual and its population. Here we use high density (n = 777,962) bi-allelic SNPs in a range of cattle breed samples to correlate ROH with the pedigree-based inbreeding coefficients and to validate subsequent analyses using 54,001 SNP genotypes. This study provides a first testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for several breeds. RESULTS: All animals genotyped on the HD panel displayed at least one ROH that was between 1--5 Mb in length with certain regions of the genome more likely to be involved in a ROH than others. Strong correlations (r = 0.75, p < 0.0001) existed between the pedigree-based inbreeding coefficient and a statistic based on sum of ROH of length > 0.5 KB and suggests that in the absence of an animal's pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples. CONCLUSIONS: Our findings suggest that ROH are frequent across all breeds but differing patterns of ROH length and burden illustrate variations in breed origins and recent management.     

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome‐wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.     

The evolutionary conservation of the mammalian peroxidase genes

Myeloperoxidase (MPO), eosinophil peroxidase (EPX) and lactoperoxidase (LPO) are mammalian peroxidase enzymes possessing similar structures and functions. Here, we demonstrate that the genes encoding these molecules form a cluster on mouse chromosome 11. Genomic sequence analysis revealed that the mouse LPO gene has similar genomic organization to the corresponding human gene. Our data strongly suggest the evolutionary conservation of mammalian peroxidase genes.     

Runs of homozygosity and distribution of functional variants in the cattle genome

Background

Recent developments in sequencing technology have facilitated widespread investigations of genomic variants, including continuous stretches of homozygous genomic regions. For cattle, a large proportion of these runs of homozygosity (ROH) are likely the result of inbreeding due to the accumulation of elite alleles from long-term selective breeding programs. In the present study, ROH were characterized in four cattle breeds with whole genome sequence data and the distribution of predicted functional variants was detected in ROH regions and across different ROH length classes.

Results

On average, 19.5 % of the genome was located in ROH across four cattle breeds. There were an average of 715.5 ROH per genome with an average size of ~750 kbp, ranging from 10 (minimum size considered) to 49,290 kbp. There was a significant correlation between shared short ROH regions and regions putatively under selection (p < 0.001). By investigating the relationship between ROH and the predicted deleterious and non-deleterious variants, we gained insight into the distribution of functional variation in inbred (ROH) regions. Predicted deleterious variants were more enriched in ROH regions than predicted non-deleterious variants, which is consistent with observations in the human genome. We also found that increased enrichment of deleterious variants was significantly higher in short (<100 kbp) and medium (0.1 to 3 Mbp) ROH regions compared with long (>3 Mbp) ROH regions (P < 0.001), which is different than what has been observed in the human genome.

Conclusions

This study illustrates the distribution of ROH and functional variants within ROH in cattle populations. These patterns are different from those in the human genome but consistent with the natural history of cattle populations, which is confirmed by the significant correlation between shared short ROH regions and regions putatively under selection. These findings contribute to understanding the effects of inbreeding and probably selection in shaping the distribution of functional variants in the cattle genome.

Electronic supplementary material

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

Runs of homozygosity do not influence survival to old age

Runs of homozygosity (ROH) are extended tracts of adjacent homozygous single nucleotide polymorphisms (SNPs) that are more common in unrelated individuals than previously thought. It has been proposed that estimating ROH on a genome-wide level, by making use of the genome-wide single nucleotide polymorphism (SNP) data, will enable to indentify recessive variants underlying complex traits. Here, we examined ROH larger than 1.5 Mb individually and in combination for association with survival in 5974 participants of the Rotterdam Study. In addition, we assessed the role of overall homozygosity, expressed as a percentage of the autosomal genome that is in ROH longer than 1.5 Mb, on survival during a mean follow-up period of 12 years. None of these measures of homozygosity was associated with survival to old age.     

Runs of homozygosity implicate autozygosity as a schizophrenia risk factor

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ~17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time.     

Runs of homozygosity reveal genome‐wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium‐density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH = 0.099), whereas Comisana showed the lowest (F_ROH = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Cavity types and microclimate: implications for ecological,evolutionary, and conservation studies

The abiotic conditions of the immediate environment of organisms are key factors for a better understanding of ecological and evolutionary processes. Yet, information in this regard is biased towards some habitat types, landscapes, and organisms. Here, we present a 2-year comparative study of the microclimatic properties (temperature, relative humidity, and their fluctuation) of three cavity types (nest boxes, cavities in bridges, and burrows in sandy cliffs) in an arid environment. We found marked and consistent months-long differences in microclimate among the three cavity types. Nest boxes were colder than the other cavity types, with temperature oscillations being an order of magnitude higher than in other cavity types. In contrast, microclimate was very stable in burrows and cavities in bridges, the former being generally warmer and drier than the latter. We also discuss the biological implications of microclimatic conditions and its variation in different cavity types by presenting two case studies, namely the temperature–humidity index and water vapor pressure during the hatching period of an endotherm and the chilling period during the diapause of an ectotherm ectoparasite. We stress the need for comparative studies of the same organisms subjected to different microclimates given the important ecological, evolutionary, and conservation implications.     

Considering evolutionary processes in conservation biology

Conservation biologists assign population distinctiveness by classifying populations as evolutionarily significant units (ESUs). Historically, this classification has included ecological and genetic data. However, recent ESU concepts, coupled with increasing availability of data on neutral genetic variation, have led to criteria based exclusively on molecular phylogenies. We argue that the earlier definitions of ESUs, which incorporated ecological data and genetic variation of adaptive significance, are more relevant for conservation. Furthermore, this dichotomous summary (ESU or not) of a continuum of population differentiation is not adequate for determining appropriate management actions. We argue for a broader categorization of population distinctiveness based on concepts of ecological and genetic exchangeability (sensu Templeton).     

Lineage-specific evolutionary rate in mammalian mtDNA

The existence of a lineage-specific nucleotide substitution rate in mammalian mtDNA has been investigated by analyzing the mtDNA of all available species, that is, 35 complete mitochondrial genomes from 14 mammalian orders. A detailed study of their evolutionary dynamics has been carried out on both ribosomal RNA and first and second codon positions (P12) of H-strand protein-coding genes by using two different types of relative-rate tests. Results are quite congruent between ribosomal and P12 sites. Significant rate variations have been observed among orders and among species of the same order. However, rate variation does not exceed 1.8-fold between the fastest (Proboscidea and Primates) and the slowest (Perissodactyla) evolving orders. Thus, the observed mitochondrial rate variations among taxa do not invalidate the suitability of mtDNA for drawing mammalian phylogeny. Dependence of evolutionary rate differences on variations in mutation and/or fixation rates was examined. Body size, generation time, and metabolic rate were tested, and no significant correlation was observed between them and the taxon-specific evolutionary rates, most likely because the latter might be influenced by multiple overlapping variable constraints.     

The conservation and evolutionary modularity of metabolism

Background  

Cellular metabolism is a fundamental biological system consisting of myriads of enzymatic reactions that together fulfill the basic requirements of life. The recent availability of vast amounts of sequence data from diverse sets of organisms provides an opportunity to systematically examine metabolism from a comparative perspective. Here we supplement existing genome and protein resources with partial genome datasets derived from 193 eukaryotes to present a comprehensive survey of the conservation of metabolism across 26 taxa representing the three domains of life.     

Transmembrane topography and evolutionary conservation of synaptophysin

Synaptophysin is the major integral membrane protein of small synaptic vesicles. Its primary structure deduced from rat and human complementary DNA sequences predicts that synaptophysin contains four transmembrane regions and a carboxyl-terminal domain having a novel repetitive structure. To elucidate the transmembrane organization of this protein in the synaptic vesicle, five antipeptide antibodies were raised. The site-specific antibodies were used to map the cognate sequences to the cytoplasmic or intravesicular side of the synaptic vesicle membrane by determining the susceptibility of the epitopes to proteolysis. The results confirm a topographic model for synaptophysin in which the protein spans the vesicle membrane four times, with both the amino and carboxyl terminus being cytoplasmic. In addition, the evolutionary conservation of the synaptophysin domains was addressed as a function of their membrane localization. To this end the primary structure of bovine synaptophysin was determined. Sequence comparisons between bovine, rat, and human synaptophysin revealed that only the intravesicular loops showed a significant number of amino acid substitutions (22%), while the transmembrane regions and cytoplasmic sequences were highly conserved (3% substitutions). These results depict synaptophysin as a protein with multiple membrane spanning regions whose functional site is likely to reside in highly conserved intramembranous and cytoplasmic sequences.     

Expression and evolutionary conservation of nanos-related genes in Hydra

The Drosophila gene nanos encodes two particular zinc finger motifs which are also found in germline-associated factors from nematodes to vertebrates. We cloned two nanos (nos)-related genes, Cnnos1 and Cnnos2 from Hydra magnipapillata. Using whole-mount in situ hybridization, the expression of Cnnos1 and Cnnos2 was examined. Cnnos1 was specifically expressed in multipotent stem cells and germline cells, but not in somatic cells. Cnnos2 was weakly expressed in germline cells and more specifically in the endoderm of the hypostome where it appears to be involved in head morphogenesis. In addition to structural conservation in the zinc finger domain of nanos-related genes, functional conservation of Cnnos1 was also demonstrated by the finding that a Cnnos1 transgene can partially rescue the nos ^RC phenotype that is defective in the egg production of Drosophila. Thus, the function of nanos-related genes in the germline appears to be well conserved from primitive to highly evolved metazoans. Received: 28 April 2000 / Accepted: 1 July 2000     

Genome-wide evolutionary conservation of N-glycosylation sites

Although posttranslational protein modifications are generally thought to perform important cellular functions, recent studies showed that a large fraction of phosphorylation sites are not evolutionarily conserved. Whether the same is true for other protein modifications, such as N-glycosylation is an open question. N-glycosylation is a form of cotranslational and posttranslational modification that occurs by enzymatic addition of a polysaccharide, or glycan, to an asparagine (N) residue of a protein. Examining a large set of experimentally determined mouse N-glycosylation sites, we find that the evolutionary rate of glycosylated asparagines is significantly lower than that of nonglycosylated asparagines of the same proteins. We further confirm that the conservation of glycosylated asparagines is accompanied by the conservation of the canonical motif sequence for glycosylation, suggesting that the above substitution rate difference is related to glycosylation. Interestingly, when solvent accessibility is considered, the substitution rate disparity between glycosylated and nonglycosylated asparagines is highly significant at solvent accessible sites but not at solvent inaccessible sites. Thus, although the solvent inaccessible glycosylation sites were experimentally identified, they are unlikely to be genuine or physiologically important. For solvent accessible asparagines, our analysis reveals a widespread and strong functional constraint on glycosylation, unlike what has been observed for phosphorylation sites in most studies, including our own analysis. Because the majority of N-glycosylation occurs at solvent accessible sites, our results show an overall functional importance for N-glycosylation.     

The evolutionary impacts of conservation actions

Conservation management for environmental sustainability is now ubiquitous. The ecological effects of these actions are well-intentioned and well-known. Although conservation biologists and managers increasingly incorporate evolutionary considerations into management plans, the evolutionary consequences of management strategies have remained relatively unexplored and unconsidered. But what are the evolutionary consequences? Here, we advocate a new research agenda focused on identifying, predicting, and countering the evolutionary consequences of conservation management. We showcase the examples of park creation and invasive species management, and speculate further on five other major methods of management. Park creation may cause selection for altered dispersal and behavior that utilizes human foods and structures. Management of invasive species may favor the evolution of resistance to or tolerance of control methods. In these and other cases, evolution may cause deviations from the predicted consequences of management strategies optimized without considering evolution, particularly when management results in or coincides with major environmental change, if population size change strongly, or if life histories are short enough to allow more rapid evolution. We call for research focused on: (1) experimental predictions and tests of evolution under particular management strategies, (2) widespread monitoring of managed populations and communities, and (3) meta-analysis and theoretical study aimed at simplifying the process of evolutionary prediction, particularly at systematizing a means of identifying traits likely to evolve due to likely existing genetic variance or high mutation rates. Ultimately, conservation biologists should incorporate evolutionary prediction into management planning to prevent the evolutionary domestication of the species that they are trying to protect.     

Assessing the utility of conserving evolutionary history

It is often claimed that conserving evolutionary history is more efficient than species‐based approaches for capturing the attributes of biodiversity that benefit people. This claim underpins academic analyses and recommendations about the distribution and prioritization of species and areas for conservation, but evolutionary history is rarely considered in practical conservation activities. One impediment to implementation is that arguments related to the human‐centric benefits of evolutionary history are often vague and the underlying mechanisms poorly explored. Herein we identify the arguments linking the prioritization of evolutionary history with benefits to people, and for each we explicate the purported mechanism, and evaluate its theoretical and empirical support. We find that, even after 25 years of academic research, the strength of evidence linking evolutionary history to human benefits is still fragile. Most – but not all – arguments rely on the assumption that evolutionary history is a useful surrogate for phenotypic diversity. This surrogacy relationship in turn underlies additional arguments, particularly that, by capturing more phenotypic diversity, evolutionary history will preserve greater ecosystem functioning, capture more of the natural variety that humans prefer, and allow the maintenance of future benefits to humans. A surrogate relationship between evolutionary history and phenotypic diversity appears reasonable given theoretical and empirical results, but the strength of this relationship varies greatly. To the extent that evolutionary history captures unmeasured phenotypic diversity, maximizing the representation of evolutionary history should capture variation in species characteristics that are otherwise unknown, supporting some of the existing arguments. However, there is great variation in the strength and availability of evidence for benefits associated with protecting phenotypic diversity. There are many studies finding positive biodiversity–ecosystem functioning relationships, but little work exists on the maintenance of future benefits or the degree to which humans prefer sets of species with high phenotypic diversity or evolutionary history. Although several arguments link the protection of evolutionary history directly with the reduction of extinction rates, and with the production of relatively greater future biodiversity via increased adaptation or diversification, there are few direct tests. Several of these putative benefits have mismatches between the relevant spatial scales for conservation actions and the spatial scales at which benefits to humans are realized. It will be important for future work to fill in some of these gaps through direct tests of the arguments we define here.     

Arabidopsis ovule development and its evolutionary conservation

Ovules have an important role during the life cycle of the plant, and they provide an excellent model for studying organogenesis in plants. As such, the molecular control of ovule development has been studied for many years. Recent studies in Arabidopsis have revealed important new data concerning ovule primordia formation, ovule identity determination, and patterning. Furthermore, interesting results about ovule development in other species, such as Petunia and rice, have been published recently. In this review, we discuss these recent findings in reference to ovule development in Arabidopsis. We compare available data with those of other species to investigate the evolutionary conservation of the regulatory pathways.