Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.     

Sequence variation patterns along a latitudinal cline in Scots pine (Pinus sylvestris): signs of clinal adaptation?

Adaptive clinal variation is abundant in nature, and its genetic basis is of great interest. The polygenic nature of this variation poses a challenge for identifying the causative loci underlying these adaptations. Here, we have examined the patterns of sequence variation in ten candidate genes for timing of bud set and cold tolerance, traits that form strong latitudinal clines in Scots pine (Pinus sylvestris). A set of reference genes was used for comparison and to infer a simplified demographic background model with approximate Bayesian computation methods. Against the resulting bottleneck model, the neutrality tests show little signs for local adaptation, but species-wide selection is suggested in some of the studied genes. In line with the theoretical expectations, we see little evidence for adaptive differentiation with F _ST methods. Instead, allele frequency clines were found in three genes (dhn1, ftl2, and prr1). Our results add to the discussion on which molecular signals best characterize a gene subject to clinally varying selection. This will be especially relevant when these kinds of observations can be examined in parallel with association study results.     

THE GENETICS OF ADAPTATION: THE GENETIC BASIS OF RESISTANCE TO WASP PARASITISM IN DROSOPHILA MELANOGASTER

There have been very few genetic analyses of “natural” adaptations, that is, those not involving artificial selection or responses to human disturbance. Here we analyze the genetic basis of geographic variation in Drosophila melanogaster's resistance to parasitism by a wasp, Asobara tabida. Our results suggest that population differences in ability to encapsulate parasitoid eggs have a fairly simple genetic basis: 60% of the D. melanogaster genome plays no role in differences between resistant and susceptible populations. Instead, resistance gene(s) are restricted to chromosome two, and may be further restricted to the centromeric region of this chromosome. This finding suggests that natural adaptations—like many responses to artificial selection and human disturbance—sometimes have a simple genetic basis.     

Joint effects of genes underlying a temperature specialization tradeoff in yeast

A central goal of evolutionary genetics is to understand, at the molecular level, how organisms adapt to their environments. For a given trait, the answer often involves the acquisition of variants at unlinked sites across the genome. Genomic methods have achieved landmark successes in pinpointing these adaptive loci. To figure out how a suite of adaptive alleles work together, and to what extent they can reconstitute the phenotype of interest, requires their transfer into an exogenous background. We studied the joint effect of adaptive, gain-of-function thermotolerance alleles at eight unlinked genes from Saccharomyces cerevisiae, when introduced into a thermosensitive sister species, S. paradoxus. Although the loci damped each other’s beneficial impact (that is, they were subject to negative epistasis), most boosted high-temperature growth alone and in combination, and none was deleterious. The complete set of eight genes was sufficient to confer ~15% of the S. cerevisiae thermotolerance phenotype in the S. paradoxus background. The same loci also contributed to a heretofore unknown advantage in cold growth by S. paradoxus. Together, our data establish temperature resistance in yeasts as a model case of a genetically complex evolutionary tradeoff, which can be partly reconstituted from the sequential assembly of unlinked underlying loci.     

Genome-Wide Delineation of Natural Variation for Pod Shatter Resistance in Brassica napus

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.     

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.     

Production of congenic mouse strains carrying genomic intervals containing SLE-susceptibility genes derived from the SLE-prone NZM2410 strain

Systemic lupus erythematosus is inherited as a complex polygenic trait. Four genomic intervals containing major SLE-susceptibility loci were previously identified by interval mapping in the NZM2410 mouse model. In this paper, we utilized a marker-assisted selection protocol to produce four congenic mouse strains, each carrying an NZM2410-derived SLE-susceptibility interval on a C57BL/6-resistant background. Each strain carries only one susceptibility allele derived from this polygenic model and consequently can be used to characterize the specific component phenotypes contributed by individual SLE-susceptibility genes. We illustrate the efficacy of this approach with phenotypic data for one of our congenic strains, B6.NZMH2 ^z . Our results indicate that this single genomic interval from Chromosome (Chr) 17 of NZM2410 can mediate increased levels of IgG autoantibodies specific for chromatin and that, similar to results obtained in our original genetic cross, B6.NZMH2 ^z/b heterozygotes are more prone than B6.NZMH2 ^z homozygotes to the development of humoral autoimmunity to nuclear antigens. These results illustrate the feasibility of using congenic strains to dissect the complex pathogenic mechanisms that mediate polygenic SLE. These congenic strains will be valuable tools in the genetic analysis of SLE susceptibility. In future studies, these congenic strains will be interbred to produce bi- and tri-congenic strains in order to assess the role of genetic interactions in the expression of specific components of SLE pathogenesis. They will also be instrumental to the positional cloning and identification of the genes responsible for SLE susceptibility, via the production of congenic recombinants. Received: 1 September 1995 / Accepted: 20 December 1995     

Species and Population Level Molecular Profiling Reveals Cryptic Recombination and Emergent Asymmetry in the Dimorphic Mating Locus of C. reinhardtii

Heteromorphic sex-determining regions or mating-type loci can contain large regions of non-recombining sequence where selection operates under different constraints than in freely recombining autosomal regions. Detailed studies of these non-recombining regions can provide insights into how genes are gained and lost, and how genetic isolation is maintained between mating haplotypes or sex chromosomes. The Chlamydomonas reinhardtii mating-type locus (MT) is a complex polygenic region characterized by sequence rearrangements and suppressed recombination between its two haplotypes, MT+ and MT−. We used new sequence information to redefine the genetic contents of MT and found repeated translocations from autosomes as well as sexually controlled expression patterns for several newly identified genes. We examined sequence diversity of MT genes from wild isolates of C. reinhardtii to investigate the impacts of recombination suppression. Our population data revealed two previously unreported types of genetic exchange in Chlamydomonas MT—gene conversion in the rearranged domains, and crossover exchanges in flanking domains—both of which contribute to maintenance of genetic homogeneity between haplotypes. To investigate the cause of blocked recombination in MT we assessed recombination rates in crosses where the parents were homozygous at MT. While normal recombination was restored in MT+×MT+ crosses, it was still suppressed in MT−×MT− crosses. These data revealed an underlying asymmetry in the two MT haplotypes and suggest that sequence rearrangements are insufficient to fully account for recombination suppression. Together our findings reveal new evolutionary dynamics for mating loci and have implications for the evolution of heteromorphic sex chromosomes and other non-recombining genomic regions.     

Mitochondrial OXPHOS genes provides insights into genetics basis of hypoxia adaptation in anchialine cave shrimps

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (TST complex) comprises twenty cave-adapted taxa, which mainly occur in the anchialine environment. Anchialine habitats may undergo drastic environmental fluctuations, including spatial and temporal changes in salinity, temperature, and dissolved oxygen content. Previous studies of crustaceans from anchialine caves suggest that they have possessed morphological, behavioral, and physiological adaptations to cope with the extreme conditions, similar to other cave-dwelling crustaceans. However, the genetic basis has not been thoroughly explored in crustaceans from anchialine habitats, which can experience hypoxic regimes. To test whether the TST shrimp-complex hypoxia adaptations matched adaptive evolution of mitochondrial OXPHOS genes. The 13 OXPHOS genes from mitochondrial genomes of 98 shrimps and 1 outgroup were examined. For each of these genes was investigated and compared to orthologous sequences using both gene (i.e. branch-site and Datamonkey) and protein (i.e. TreeSAAP) level approaches. Positive selection was detected in 11 of the 13 candidate genes, and the radical amino acid changes sites scattered throughout the entire TST complex phylogeny. Additionally, a series of parallel/convergent amino acid substitutions were identified in mitochondrial OXPHOS genes of TST complex shrimps, which reflect functional convergence or similar genetic mechanisms of cave adaptation. The extensive occurrence of positive selection is suggestive of their essential role in adaptation to hypoxic anchialine environment, and further implying that TST complex shrimps might have acquired a finely capacity for energy metabolism. These results provided some new insights into the genetic basis of anchialine hypoxia adaptation.     

Genetics of Soybean-Heterodera glycines Interactions

The soybean cyst nematode, Heterodera glycines, is one of the most economically important pathogens of soybean. Effective management of the nematode is often dependent on the planting of resistant soybean cultivars. During the past 40 years, more than 60 soybean genotypes and plant introductions (PI) have been reported as resistant to H. glycines. About 130 modern soybean cultivars registered in the United States are resistant to certain races of H. glycines. Several resistance genes have been identified and genetically mapped; however, resistance levels in many soybean cultivars are not durable. Some older cultivars are no longer resistant to certain H. glycines populations in many production areas, especially if a soybean monoculture has been practiced. Past soybean registration reports show that all resistant cultivars developed in public institutions from the mid-1960s to the present have been derived from five PIs. This narrow genetic background is fragile. To further complicate the issue, soybean-H. glycines genetic interactions are complex and poorly understood. Studies to identify soybean resistance genes sometimes have overlapped, and the same genes may have been reported several times and designated by different names. Nevertheless, many potential resistance genes in existing germplasm resources have not yet been characterized. Clearly, it is necessary to identify new resistance genes, develop more precise selection methods, and integrate these resistance genes into new cultivars. Rational deployment of resistant cultivars is critical to future sustained soybean production.     

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent‐wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV‐driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome‐wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre‐RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV‐infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus‐driven selection in response to a suddenly emerging infectious disease.     

The impact of cis-acting polymorphisms on the human phenotype

Cis-acting polymorphisms that affect gene expression are now known to be frequent, although the extent and mechanisms by which such variation affects the human phenotype are, as yet, only poorly understood. Key signatures of cis-acting variation are differences in gene expression that are tightly associated with regulatory SNPs or expression Quantitative Trait Loci (eQTL) and an imbalance of allelic expression (AEI) in heterozygous samples. Such cis-acting sequence differences appear often to have been under selection within and between populations and are also thought to be important in speciation. Here we describe the example of lactase persistence. In medical research, variants that affect regulation in cis have been implicated in both monogenic and polygenic disorders, and in the metabolism of drugs. In this review we suggest that by further understanding common regulatory variations and how they interact with other genetic and environmental variables it will be possible to gain insight into important mechanisms behind complex disease, with the potential to lead to new methods of diagnosis and treatments.     

The [URE3] prion is not conserved among Saccharomyces species

The [URE3] prion of Saccharomyces cerevisiae is a self-propagating inactive form of the nitrogen catabolism regulator Ure2p. To determine whether the [URE3] prion is conserved in S. cerevisiae-related yeast species, we have developed genetic tools allowing the detection of [URE3] in Saccharomyces paradoxus and Saccharomyces uvarum. We found that [URE3] is conserved in S. uvarum. In contrast, [URE3] was not detected in S. paradoxus. The inability of S. paradoxus Ure2p to switch to a prion isoform results from the primary sequence of the protein and not from the lack of cellular cofactors as heterologous Ure2p can propagate [URE3] in this species. Our data therefore demonstrate that [URE3] is conserved only in a subset of Saccharomyces species. Implications of our finding on the physiological and evolutionary meaning of the yeast [URE3] prion are discussed.     

Genetic Basis of Haloperidol Resistance in Saccharomyces cerevisiae Is Complex and Dose Dependent

The genetic basis of most heritable traits is complex. Inhibitory compounds and their effects in model organisms have been used in many studies to gain insights into the genetic architecture underlying quantitative traits. However, the differential effect of compound concentration has not been studied in detail. In this study, we used a large segregant panel from a cross between two genetically divergent yeast strains, BY4724 (a laboratory strain) and RM11_1a (a vineyard strain), to study the genetic basis of variation in response to different doses of a drug. Linkage analysis revealed that the genetic architecture of resistance to the small-molecule therapeutic drug haloperidol is highly dose-dependent. Some of the loci identified had effects only at low doses of haloperidol, while other loci had effects primarily at higher concentrations of the drug. We show that a major QTL affecting resistance across all concentrations of haloperidol is caused by polymorphisms in SWH1, a homologue of human oxysterol binding protein. We identify a complex set of interactions among the alleles of the genes SWH1, MKT1, and IRA2 that are most pronounced at a haloperidol dose of 200 µM and are only observed when the remainder of the genome is of the RM background. Our results provide further insight into the genetic basis of drug resistance.     

Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.