Genetic architecture of body size in mammals

Island demographics and trait associations in white-tailed deer

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

The distribution of fitness effects during adaptive walks using a simple genetic network

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.

The far-reaching effects of genetic process in a keystone predator species, grey wolves

Although detrimental genetic processes are known to adversely affect the viability of populations, little is known about how detrimental genetic processes in a keystone species can affect the functioning of ecosystems. Here, we assessed how changes in the genetic characteristics of a keystone predator, grey wolves, affected the ecosystem of Isle Royale National Park over two decades. Changes in the genetic characteristic of the wolf population associated with a genetic rescue event, followed by high levels of inbreeding, led to a rise and then fall in predation rates on moose, the primary prey of wolves and dominant mammalian herbivore in this system. Those changes in predation rate led to large fluctuations in moose abundance, which in turn affected browse rates on balsam fir, the dominant forage for moose during winter and an important boreal forest species. Thus, forest dynamics can be traced back to changes in the genetic characteristics of a predator population.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Large offspring have enhanced lifetime reproductive success: Long-term carry-over effects of weaning size in gray seals (Halichoerus grypus)

An individual's size in early stages of life may be an important source of individual variation in lifetime reproductive performance, as size effects on ontogenetic development can have cascading physiological and behavioral consequences throughout life. Here, we explored how size-at-young influences subsequent reproductive performance in gray seals (Halichoerus grypus) using repeated encounter and reproductive data on a marked sample of 363 females that were measured for length after weaning, at ~4 weeks of age, and eventually recruited to the Sable Island breeding colony. Two reproductive traits were considered: provisioning performance (mass of weaned offspring), modeled using linear mixed effects models; and reproductive frequency (rate at which a female returns to breed), modeled using mixed effects multistate mark-recapture models. Mothers with the longest weaning lengths produced pups 8 kg heavier and were 20% more likely to breed in a given year than mothers with the shortest lengths. Correlation in body lengths between weaning and adult life stages, however, is weak: Longer pups do not grow to be longer than average adults. Thus, covariation between weaning length and future reproductive performance appears to be a carry-over effect, where the size advantages afforded in early juvenile stages may allow enhanced long-term performance in adulthood.

Novel Haplotype in the HHEX Gene Promoter Associated with Body Length in Pigs

The screening of important candidate genes and the identification of genetic markers are important for molecular selection in the pig industry. The hematopoietically expressed homeobox (HHEX) gene plays an important role in embryonic development and organogenesis; however, the genetic variation and expression pattern of the porcine HHEX gene remains to be clarified. In this study, semiquantitative RT-PCR and immunohistochemistry results showed the specific expression of the HHEX gene in porcine cartilage tissues. A novel haplotype consisting of two SNPs rs80901185 (T > C) and rs80934526 (A > G) was detected in the promoter region of the HHEX gene. The expression of the HHEX gene was significantly higher in Yorkshire pigs (TA haplotype) than in Wuzhishan pigs (CG haplotype), and a population analysis showed that this haplotype was significantly associated with body length. An analysis subsequently revealed that the -586 to -1 bp region of the HHEX gene promoter showed the highest activity. Furthermore, we found that the activity of the TA haplotype was significantly higher than that of the CG haplotype by changing the potential binding of transcription factors YY1 and HDAC2. In summary, we conclude that the porcine HHEX gene may contribute to the breeding of pigs for body length traits.

Screening Indels from the whole genome to identify the candidates and their association with economic traits in several goat breeds

In the present study, the re-sequencing of our previous whole-genome sequencing (WGS) for selected individuals of Dazu-black goat (DBG) and Inner-Mongolia Cashmere goat (IMCG) breeds was used to detect and compare the differentiation in Indels depending on the reference genome of goat. Then, three selected candidate Indels rs668795676, rs657996810, and rs669452874 of the three genes SUFU, SYCP2L and GLIPR1L1, respectively, have been chosen, based on the results of prior GWAS across the genome, and examined for their association with body weight and dimensions (body height, body length, heart girth, chest width, cannon circumference, and chest depth) by kompetitive allele specific PCR assay for 342 goats from the three studied goat breeds (DBG, n = 203; ♂99, ♀104), IMCG (n = 65; 15♂, 50♀), and Hechuan white goat (HWG, n = 74; 34♂, 40♀) breeds. The analysis of 192.747 Gb WGS revealed an average 334,151 Indels in the whole genome of DBG and IMCG breeds. Chromosome 1 had a maximum number of mutations (Indels) of 58,497 and 55,527 for IMCG and DBG, respectively, while chromosome 25 had the least number of mutations of 15,680 and 16,103 for IMCG and DBG, respectively. The majority of Indels were either Ins or Del of short fragments of 1-5 bp, which covered 79.06 and 71.78% of the total number of Indels mutations in IMCG and DBG, respectively. Comparing the differences of Indels between the studied goat breeds revealed 100 and 110 unique Indels for IMCG and DBG, respectively. The Indels loci in the intron region were unique for both studied goat breeds which were related to 30 and 38 candidate genes in IMCG and DBG, respectively, including SUFU, SYCP2L, and GLIPR1L1 genes. Concerning rs669452874 locus, body height and body length of Del/Del genotype in DBG were significantly higher (P < 0.05) than that of Ins/Del genotype, while body height and body length of Del/Del genotype in IMCG were significantly higher (P < 0.01) than those of Ins/Ins genotype, whereas body height and body length and heart girth of Del/Del genotype in HWG were significantly higher (P < 0.01) than those of the Ins/Del and Ins/Ins genotypes. Thus, Del/Del genotype of rs669452874 locus can be used as a candidate molecular marker related to the body dimensions in the studied goat breeds.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

The genetic architecture of phenotypic diversity in the Betta fish (Betta splendens)

The Betta fish displays a remarkable variety of phenotypes selected during domestication. However, the genetic basis underlying these traits remains largely unexplored. Here, we report a high-quality genome assembly and resequencing of 727 individuals representing diverse morphotypes of the Betta fish. We show that current breeds have a complex domestication history with extensive introgression with wild species. Using a genome-wide association study, we identify the genetic basis of multiple traits, including coloration patterns, the "Dumbo" phenotype with pectoral fin outgrowth, extraordinary enlargement of body size that we map to a major locus on chromosome 8, the sex determination locus that we map to dmrt1, and the long-fin phenotype that maps to the locus containing kcnj15. We also identify a polygenic signal related to aggression, involving multiple neural system-related genes such as esyt2, apbb2, and pank2. Our study provides a resource for developing the Betta fish as a genetic model for morphological and behavioral research in vertebrates.

Maintenance of quantitative genetic variance in complex, multi-trait phenotypes: The contribution of rare, large effect variants in two Drosophila species

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.

Genomic diversity and signatures of selection in meat and fancy rabbit breeds based on high-density marker data

Background

Domestication of the rabbit (Oryctolagus cuniculus) has led to a multi-purpose species that includes many breeds and lines with a broad phenotypic diversity, mainly for external traits (e.g. coat colours and patterns, fur structure, and morphometric traits) that are valued by fancy rabbit breeders. As a consequence of this human-driven selection, distinct signatures are expected to be present in the rabbit genome, defined as signatures of selection or selective sweeps. Here, we investigated the genome of three Italian commercial meat rabbit breeds (Italian Silver, Italian Spotted and Italian White) and 12 fancy rabbit breeds (Belgian Hare, Burgundy Fawn, Champagne d’Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex, Rhinelander and Thuringian) by using high-density single nucleotide polymorphism data. Signatures of selection were identified based on the fixation index (F_ST) statistic with different approaches, including single-breed and group-based methods, the latter comparing breeds that are grouped based on external traits (different coat colours and body sizes) and types (i.e. meat vs. fancy breeds).

Results

We identified 309 genomic regions that contained signatures of selection and that included genes that are known to affect coat colour (ASIP, MC1R and TYR), coat structure (LIPH), and body size (LCORL/NCAPG, COL11A1 and HOXD) in rabbits and that characterize the investigated breeds. Their identification proves the suitability of the applied methodologies for capturing recent selection events. Other regions included novel candidate genes that might contribute to the phenotypic variation among the analyzed breeds, including genes for pigmentation-related traits (EDNRA, EDNRB, MITF and OCA2) and body size, with a strong candidate for dwarfism in rabbit (COL2A1).

Conclusions

We report a genome-wide view of genetic loci that underlie the main phenotypic differences in the analyzed rabbit breeds, which can be useful to understand the shift from the domestication process to the development of breeds in O. cuniculus. These results enhance our knowledge about the major genetic loci involved in rabbit external traits and add novel information to understand the complexity of the genetic architecture underlying body size in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00696-9.

Identification and validation of a regulatory mutation upstream of the BMP2 gene associated with carcass length in pigs

Background

Carcass length is very important for body size and meat production for swine, thus understanding the genetic mechanisms that underly this trait is of great significance in genetic improvement programs for pigs. Although many quantitative trait loci (QTL) have been detected in pigs, very few have been fine-mapped to the level of the causal mutations. The aim of this study was to identify potential causal single nucleotide polymorphisms (SNPs) for carcass length by integrating a genome-wide association study (GWAS) and functional assays.

Results

Here, we present a GWAS in a commercial Duroc × (Landrace × Yorkshire) (DLY) population that reveals a prominent association signal (P = 4.49E−07) on pig chromosome 17 for carcass length, which was further validated in two other DLY populations. Within the detected 1 Mb region, the BMP2 gene stood out as the most likely causal candidate because of its functions in bone growth and development. Whole-genome gene expression studies showed that the BMP2 gene was differentially expressed in the cartilage tissues of pigs with extreme carcass length. Then, we genotyped an additional 267 SNPs in 500 selected DLY pigs, followed by further whole-genome SNP imputation, combined with deep genome resequencing data on multiple pig breeds. Reassociation analyses using genotyped and imputed SNP data revealed that the rs320706814 SNP, located approximately 123 kb upstream of the BMP2 gene, was the strongest candidate causal mutation, with a large association with carcass length, with a ~ 4.2 cm difference in length across all three DLY populations (N = 1501; P = 3.66E−29). This SNP segregated in all parental lines of the DLY (Duroc, Large White and Landrace) and was also associated with a significant effect on body length in 299 pure Yorkshire pigs (P = 9.2E−4), which indicates that it has a major value for commercial breeding. Functional assays showed that this SNP is likely located within an enhancer and may affect the binding affinity of transcription factors, thereby regulating BMP2 gene expression.

Conclusions

Taken together, these results suggest that the rs320706814 SNP on pig chromosome 17 is a putative causal mutation for carcass length in the widely used DLY pigs and has great value in breeding for body size in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00689-0.

Effect of dietary protein on growth performance, and serum biochemical index in late pregnant Hu ewes and their offspring

The purpose of this study was to investigate the effects of different protein levels in late pregnancy on ewe and lamb growth performance, serum biochemical indexes. Thirty-three ewes (46.4 ± 1.38 kg initial weight) were randomly divided into 3 groups, with 11 ewes in each group. The protein levels of three diets formulated to provide components to meet 10.00 MJ/kg ME requirements diets were: 10.12%, 11.26%, 12.4%. Ewes were raised from the 90th day of pregnancy to the end of delivery, and the lambs were weaned at 60 days. Dietary protein levels had significant effects on blood urea nitrogen, glucose, ammonia nitrogen and triglyceride of ewes (p < 0.05). The height, chest depth, chest circumference, straight crown hip length and curved crown hip length of lambs decreased at first and then increased with the increase of protein. The body length, chest circumference, head width and head length of weaned lambs decreased at first and then increased with the increase of protein. The results showed that when the dietary protein level was increased to 12.4%, the amino acid, glucose and fat metabolism of ewes were affected. The body size development of lambs was better than 10.12% and 11.26% proteins.

The GWAS Analysis of Body Size and Population Verification of Related SNPs in Hu Sheep

Body size is an important indicator of growth and health in sheep. In the present study, we performed Genome-Wide Association Studies (GWAS) to detect significant single-nucleotide polymorphisms (SNPs) associated with Hu sheep's body size. After genotyping parental (G1) and offspring (G2) generation of the nucleus herd for meat production of Hu sheep and conducting GWAS on the body height, chest circumference, body length, tail length, and tail width of the two groups, 5 SNPs associated with body height and 4 SNPs correlated with chest circumference were identified at the chromosomal significance level. No SNPs were significantly correlated to body length, tail length, and width. Four out of the 9 SNPs were found to be located within the 4 genes. KITLG and CADM2 are considered as candidate functional genes related to body height; MCTP1 and COL4A6 are candidate functional genes related to chest circumference. The 9 SNPs found in GWAS were verified using the G3 generation of the nucleus herd for meat production. Nine products were amplified around the 9 sites, and 29 SNPs were found; 3 mutation sites, G > C mutation at 134 bp downstream of s554331, T > G mutation at 19 bp upstream of s26859.1, and A > G mutation at 81 bp downstream of s26859.1, were significantly correlated to the body height. Dual-luciferase reporter gene experiments showed that the 3 SNPs could significantly impact dual-luciferase and gene transcription activity.

Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population

Chicken growth traits are economically important, but the relevant genetic mechanisms have not yet been elucidated. Herein, we performed a genome-wide association study to identify the variants associated with growth traits. In total, 860 chickens from a Gushi-Anka F₂ resource population were phenotyped for 68 growth and carcass traits, and 768 samples were genotyped based on the genotyping-by-sequencing (GBS) method. Finally, 734 chickens and 321,314 SNPs remained after quality control and removal of the sex chromosomes, and these data were used to carry out a GWAS analysis. A total of 470 significant single-nucleotide polymorphisms (SNPs) for 43 of the 68 traits were detected and mapped on chromosomes (Chr) 1-6, -9, -10, -16, -18, -23, and -27. Of these, the significant SNPs in Chr1, -4, and -27 were found to be associated with more than 10 traits. Multiple traits shared significant SNPs, indicating that the same mutation in the region might have a large effect on multiple growth or carcass traits. Haplotype analysis revealed that SNPs within the candidate region of Chr1 presented a mosaic pattern. The significant SNPs and pathway enrichment analysis revealed that the MLNR, MED4, CAB39L, LDB2, and IGF2BP1 genes could be putative candidate genes for growth and carcass traits. The findings of this study improve our understanding of the genetic mechanisms regulating chicken growth and carcass traits and provide a theoretical basis for chicken breeding programs.

Using prior information from humans to prioritize genes and gene-associated variants for complex traits in livestock

Genome-Wide Association Studies (GWAS) in large human cohorts have identified thousands of loci associated with complex traits and diseases. For identifying the genes and gene-associated variants that underlie complex traits in livestock, especially where sample sizes are limiting, it may help to integrate the results of GWAS for equivalent traits in humans as prior information. In this study, we sought to investigate the usefulness of results from a GWAS on human height as prior information for identifying the genes and gene-associated variants that affect stature in cattle, using GWAS summary data on samples sizes of 700,000 and 58,265 for humans and cattle, respectively. Using Fisher's exact test, we observed a significant proportion of cattle stature-associated genes (30/77) that are also associated with human height (odds ratio = 5.1, p = 3.1e-10). Result of randomized sampling tests showed that cattle orthologs of human height-associated genes, hereafter referred to as candidate genes (C-genes), were more enriched for cattle stature GWAS signals than random samples of genes in the cattle genome (p = 0.01). Randomly sampled SNPs within the C-genes also tend to explain more genetic variance for cattle stature (up to 13.2%) than randomly sampled SNPs within random cattle genes (p = 0.09). The most significant SNPs from a cattle GWAS for stature within the C-genes did not explain more genetic variance for cattle stature than the most significant SNPs within random cattle genes (p = 0.87). Altogether, our findings support previous studies that suggest a similarity in the genetic regulation of height across mammalian species. However, with the availability of a powerful GWAS for stature that combined data from 8 cattle breeds, prior information from human-height GWAS does not seem to provide any additional benefit with respect to the identification of genes and gene-associated variants that affect stature in cattle.

Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds

Background

Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions.

Results

Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C).

Conclusions

We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Genetic control of temperament traits across species: association of autism spectrum disorder risk genes with cattle temperament

Background

Temperament traits are of high importance across species. In humans, temperament or personality traits correlate with psychological traits and psychiatric disorders. In cattle, they impact animal welfare, product quality and human safety, and are therefore of direct commercial importance. We hypothesized that genetic factors that contribute to variation in temperament among individuals within a species will be shared between humans and cattle. Using imputed whole-genome sequence data from 9223 beef cattle from three cohorts, a series of genome-wide association studies was undertaken on cattle flight time, a temperament phenotype measured as the time taken for an animal to cover a short-fixed distance after release from an enclosure. We also investigated the association of cattle temperament with polymorphisms in bovine orthologs of risk genes for neuroticism, schizophrenia, autism spectrum disorders (ASD), and developmental delay disorders in humans.

Results

Variants with the strongest associations were located in the bovine orthologous region that is involved in several behavioural and cognitive disorders in humans. These variants were also partially validated in independent cattle cohorts. Genes in these regions (BARHL2, NDN, SNRPN, MAGEL2, ABCA12, KIFAP3, TOPAZ1, FZD3, UBE3A, and GABRA5) were enriched for the GO term neuron migration and were differentially expressed in brain and pituitary tissues in humans. Moreover, variants within 100 kb of ASD susceptibility genes were associated with cattle temperament and explained 6.5% of the total additive genetic variance in the largest cattle cohort. The ASD genes with the most significant associations were GABRB3 and CUL3. Using the same 100 kb window, a weak association was found with polymorphisms in schizophrenia risk genes and no association with polymorphisms in neuroticism and developmental delay disorders risk genes.

Conclusions

Our analysis showed that genes identified in a meta-analysis of cattle temperament contribute to neuron development functions and are differentially expressed in human brain tissues. Furthermore, some ASD susceptibility genes are associated with cattle temperament. These findings provide evidence that genetic control of temperament might be shared between humans and cattle and highlight the potential for future analyses to leverage results between species.

Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color

Coat color is among the most distinctive phenotypes in cattle. Worldwide, several breeds share peculiar coat color features such as the presence of a fawn pigmentation of the calf at birth, turning over time to grey, and sexual dichromatism. The aim of this study was to search for polymorphisms under differential selection by contrasting grey cattle breeds displaying the above phenotype with non-grey cattle breeds, and to identify the underlying genes. Using medium-density SNP array genotype data, a multi-cohort F_ST-outlier approach was adopted for a total of 60 pair-wise comparisons of the 15 grey with 4 non-grey cattle breeds (Angus, Limousin, Charolais, and Holstein), with the latter selected as representative of solid and piebald phenotypes, respectively. Overall, more than 50 candidate genes were detected; almost all were either directly or indirectly involved in pigmentation, and some of them were already known for their role in phenotypes related with hair graying in mammals. Notably, 17 relevant genes, including SDR16C5, MOS, SDCBP, and NSMAF, were located in a signal on BTA14 convergently observed in all the four considered scenarios. Overall, the key stages of pigmentation (melanocyte development, melanogenesis, and pigment trafficking/transfer) were all represented among the pleiotropic functions of the candidate genes, suggesting the complex nature of the grey phenotype in cattle.

Adult survival in migratory caribou is negatively associated with MHC functional diversity

Genes of the major histocompatibility complex (MHC) are involved in acquired immunity in vertebrates. Only a few studies have investigated the fitness consequences of MHC gene diversity in wild populations. Here, we looked at the association between annual survival and body mass and MHC-DRB exon 2 (MHC-DRB) genetic diversity, obtained from high-throughput sequencing, in two declining migratory caribou (Rangifer tarandus) herds. To disentangle the potential direct and general effects of MHC-DRB genetic diversity, we compared different indices of diversity that were either based on DNA-sequence variation or on physicochemical divergence of the translated peptides, thereby covering a gradient of allelic-to-functional diversity. We found that (1) body mass was not related to MHC-DRB diversity or genotype, and (2) adult survival probability was negatively associated with point accepted mutation distance, a corrected distance that considers the likelihood of each amino acid substitution to be accepted by natural selection. In addition, we found no evidence of fluctuating selection over time on MHC-DRB diversity. We concluded that direct effects were involved in the negative relationship between MHC functional diversity and survival, although the mechanism underlying this result remains unclear. A possible explanation could be that individuals with higher MHC diversity suffer higher costs of immunity (immunopathology). Our results suggest that genetic diversity is not always beneficial even in genes that are likely to be strongly shaped by balancing selection.

Epigenomics and genotype-phenotype association analyses reveal conserved genetic architecture of complex traits in cattle and human

Background

Lack of comprehensive functional annotations across a wide range of tissues and cell types severely hinders the biological interpretations of phenotypic variation, adaptive evolution, and domestication in livestock. Here we used a combination of comparative epigenomics, genome-wide association study (GWAS), and selection signature analysis, to shed light on potential adaptive evolution in cattle.

Results

We cross-mapped 8 histone marks of 1300 samples from human to cattle, covering 178 unique tissues/cell types. By uniformly analyzing 723 RNA-seq and 40 whole genome bisulfite sequencing (WGBS) datasets in cattle, we validated that cross-mapped histone marks captured tissue-specific expression and methylation, reflecting tissue-relevant biology. Through integrating cross-mapped tissue-specific histone marks with large-scale GWAS and selection signature results, we for the first time detected relevant tissues and cell types for 45 economically important traits and artificial selection in cattle. For instance, immune tissues are significantly associated with health and reproduction traits, multiple tissues for milk production and body conformation traits (reflecting their highly polygenic architecture), and thyroid for the different selection between beef and dairy cattle. Similarly, we detected relevant tissues for 58 complex traits and diseases in humans and observed that immune and fertility traits in humans significantly correlated with those in cattle in terms of relevant tissues, which facilitated the identification of causal genes for such traits. For instance, PIK3CG, a gene highly specifically expressed in mononuclear cells, was significantly associated with both age-at-menopause in human and daughter-still-birth in cattle. ICAM, a T cell-specific gene, was significantly associated with both allergic diseases in human and metritis in cattle.

Conclusion

Collectively, our results highlighted that comparative epigenomics in conjunction with GWAS and selection signature analyses could provide biological insights into the phenotypic variation and adaptive evolution. Cattle may serve as a model for human complex traits, by providing additional information beyond laboratory model organisms, particularly when more novel phenotypes become available in the near future.

Whole genome sequence analysis reveals genetic structure and X-chromosome haplotype structure in indigenous Chinese pigs

Chinese indigenous pigs exhibit considerable phenotypic diversity, but their population structure and the genetic basis of agriculturally important traits need further exploration. Here, we sequenced the whole genomes of 24 individual pigs representing 22 breeds distributed throughout China. For comparison with European and commercial breeds (one pig per breed), we included seven published pig genomes with our new genomes for analyses. Our results showed that breeds grouped together based on morphological classifications are not necessarily more genetically similar to each other than to breeds from other groups. We found that genetic material from European pigs likely introgressed into five Chinese breeds. We have identified two new subpopulations of domestic pigs that encompass morphology-based criteria in China. The Southern Chinese subpopulation comprises the classical South Chinese Type and part of the Central China Type. In contrast, the Northern Chinese subpopulation comprises the North China Type, the Lower Yangtze River Basin Type, the Southwest Type, the Plateau Type, and the remainder of the Central China Type. Eight haplotypes and two recombination sites were identified within a conserved 40.09 Mb linkage-disequilibrium (LD) block on the X chromosome. Potential candidate genes (LEPR, FANCC, COL1A1, and PCCA) influencing body size were identified. Our findings provide insights into the phylogeny of Chinese indigenous pig breeds and benefit gene mining efforts to improve major economic traits.

The genetic architecture of the maize progenitor, teosinte, and how it was altered during maize domestication

The genetics of domestication has been extensively studied ever since the rediscovery of Mendel's law of inheritance and much has been learned about the genetic control of trait differences between crops and their ancestors. Here, we ask how domestication has altered genetic architecture by comparing the genetic architecture of 18 domestication traits in maize and its ancestor teosinte using matched populations. We observed a strongly reduced number of QTL for domestication traits in maize relative to teosinte, which is consistent with the previously reported depletion of additive variance by selection during domestication. We also observed more dominance in maize than teosinte, likely a consequence of selective removal of additive variants. We observed that large effect QTL have low minor allele frequency (MAF) in both maize and teosinte. Regions of the genome that are strongly differentiated between teosinte and maize (high FST) explain less quantitative variation in maize than teosinte, suggesting that, in these regions, allelic variants were brought to (or near) fixation during domestication. We also observed that genomic regions of high recombination explain a disproportionately large proportion of heritable variance both before and after domestication. Finally, we observed that about 75% of the additive variance in both teosinte and maize is "missing" in the sense that it cannot be ascribed to detectable QTL and only 25% of variance maps to specific QTL. This latter result suggests that morphological evolution during domestication is largely attributable to very large numbers of QTL of very small effect.

The LCORL Locus is under Selection in Large-Sized Pakistani Goat Breeds

Goat domestication and human selection for valued traits have formed diverse breeds with characteristic phenotypes. This process led to the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity-so-called "selection signatures". We previously reported an analysis of selection signatures based on pooled whole-genome sequencing data of 20 goat breeds and bezoar goats. In the present study, we reanalyzed the data and focused on a subset of eight Pakistani goat breeds (Angora, Barbari, Beetal, Dera Din Panah, Kamori, Nachi, Pahari, Teddy). We identified 749 selection signatures based on reduced heterozygosity in these breeds. A search for signatures that are shared across large-sized goat breeds revealed that five medium-to-large-sized Pakistani goat breeds had a common selection signature on chromosome 6 in a region harboring the LCORL gene, which has been shown to modulate height or body size in several mammalian species. Fine-mapping of the region confirmed that all five goat breeds with the selection signature were nearly fixed for the same haplotype in a ~191 kb region spanning positions 37,747,447-37,938,449. From the pool sequencing data, we identified a frame-shifting single base insertion into an isoform-specific exon of LCORL as a potential candidate causal variant mediating the size-increasing effect. If this preliminary result can be confirmed in independent replication studies, genotyping of this variant might be used to improve breeding programs and the selection for stature in goats.

Quantitative genetics of body size evolution on islands: an individual-based simulation approach

According to the island rule, small-bodied vertebrates will tend to evolve larger body size on islands, whereas the opposite happens to large-bodied species. This controversial pattern has been studied at the macroecological and biogeographical scales, but new developments in quantitative evolutionary genetics now allow studying the island rule from a mechanistic perspective. Here, we develop a simulation approach based on an individual-based model to model body size change on islands as a progressive adaptation to a moving optimum, determined by density-dependent population dynamics. We applied the model to evaluate body size differentiation in the pigmy extinct hominin Homo floresiensis, showing that dwarfing may have occurred in only about 360 generations (95% CI ranging from 150 to 675 generations). This result agrees with reports suggesting rapid dwarfing of large mammals on islands, as well as with the recent discovery that small-sized hominins lived in Flores as early as 700 kyr ago. Our simulations illustrate the power of analysing ecological and evolutionary patterns from an explicit quantitative genetics perspective.

Closing the 'phenotype gap' in precision medicine: improving what we measure to understand complex disease mechanisms

The central concept underlying precision medicine is a mechanistic understanding of each disease and its response to therapy sufficient to direct a specific intervention. To execute on this vision requires parsing incompletely defined disease syndromes into discrete mechanistic subsets and developing interventions to precisely address each of these etiologically distinct entities. This will require substantial adjustment of traditional paradigms which have tended to aggregate high-level phenotypes with very different etiologies. In the current environment, where diagnoses are not mechanistic, drug development has become so expensive that it is now impractical to imagine the cost-effective creation of new interventions for many prevalent chronic conditions. The vision of precision medicine also argues for a much more seamless integration of research and development with clinical care, where shared taxonomies will enable every clinical interaction to inform our collective understanding of disease mechanisms and drug responses. Ideally, this would be executed in ways that drive real-time and real-world discovery, innovation, translation, and implementation. Only in oncology, where at least some of the biology is accessible through surgical excision of the diseased tissue or liquid biopsy, has "co-clinical" modeling proven feasible. In most common germline disorders, while genetics often reveal the causal mutations, there still remain substantial barriers to efficient disease modeling. Aggregation of similar disorders under single diagnostic labels has directly contributed to the paucity of etiologic and mechanistic understanding by directly reducing the resolution of any subsequent studies. Existing clinical phenotypes are typically anatomic, physiologic, or histologic, and result in a substantial mismatch in information content between the phenomes in humans or in animal 'models' and the variation in the genome. This lack of one-to-one mapping of discrete mechanisms between disease and animal models causes a failure of translation and is one form of 'phenotype gap.' In this review, we will focus on the origins of the phenotyping deficit and approaches that may be considered to bridge the gap, creating shared taxonomies between human diseases and relevant models, using cardiovascular examples.

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.

Insights into body size variation in cetaceans from the evolution of body-size-related genes

BACKGROUND

Cetaceans exhibit an exceptionally wide range of body size, yet in this regard, their genetic basis remains poorly explored. In this study, 20 body-size-related genes for which duplication, mutation, or deficiency can cause body size change in mammals were chosen to preliminarily investigate the evolutionary mechanisms underlying the dramatic body size variation in cetaceans.

RESULTS

We successfully sequenced 20 body-size-related genes in six representative species of cetaceans. A total of 46 codons from 10 genes were detected and determined to be under strong positive selection, 32 (69.6%) of which were further found to be under radical physiochemical changes; moreover, some of these sites were localized in or near important functional regions. Interestingly, positively selected genes were well matched with body size evolution: for small cetaceans, strong evidence of positive selection was detected at ACAN, OBSL1, and GRB10, within which mutations or duplications could cause short stature; positive selection was found in large cetaceans at CBS and EIF2AK3, which could promote growth, and at the PLOD1 gene, within which mutations could cause tall stature. Importantly, relationship analyses revealed that the evolutionary rate of CBS was positively related to body length and body mass with statistical significance. Additionally, we identified 32 cetacean-specific amino acid changes in 10 genes.

CONCLUSIONS

This is the first study to investigate the molecular basis of dramatic body size variation in cetaceans. Our results provide evidence of the positive selection of several body-size-related genes in cetaceans, as well as divergent selection between large or small cetaceans, which suggest cetacean body size variation possibly associated with these genes. In addition, cetacean-specific amino acid changes might have played key roles in body size evolution after the divergence of cetaceans from their terrestrial relatives. Overall, the evolutionary pattern of these body-size-related genes could provide new insights into genetic mechanisms for the body size variation in cetaceans.

Common Disease Is More Complex Than Implied by the Core Gene Omnigenic Model

The evidence that most adult-onset common diseases have a polygenic genetic architecture fully consistent with robust biological systems supported by multiple back-up mechanisms is now overwhelming. In this context, we consider the recent "omnigenic" or "core genes" model. A key assumption of the model is that there is a relatively small number of core genes relevant to any disease. While intuitively appealing, this model may underestimate the biological complexity of common disease, and therefore, the goal to discover core genes should not guide experimental design. We consider other implications of polygenicity, concluding that a focus on patient stratification is needed to achieve the goals of precision medicine.

Hypomorphic mutation of the mouse Huntington's disease gene orthologue

Rare individuals with inactivating mutations in the Huntington's disease gene (HTT) exhibit variable abnormalities that imply essential HTT roles during organ development. Here we report phenotypes produced when increasingly severe hypomorphic mutations in the murine HTT orthologue Htt, (HdhneoQ20, HdhneoQ50, HdhneoQ111), were placed over a null allele (Hdhex4/5). The most severe hypomorphic allele failed to rescue null lethality at gastrulation, while the intermediate, though still severe, alleles yielded recessive perinatal lethality and a variety of fetal abnormalities affecting body size, skin, skeletal and ear formation, and transient defects in hematopoiesis. Comparative molecular analysis of wild-type and Htt-null retinoic acid-differentiated cells revealed gene network dysregulation associated with organ development that nominate polycomb repressive complexes and miRNAs as molecular mediators. Together these findings demonstrate that Htt is required both pre- and post-gastrulation to support normal development.

Population structure, genetic diversity and selection signatures within seven indigenous Pakistani goat populations

Goat farming in Pakistan depends on indigenous breeds that have adapted to specific agro-ecological conditions. Pakistan has a rich resource of goat breeds, and the genetic diversity of these goat breeds is largely unknown. In this study, genetic diversity and population structure were characterized from seven indigenous goat breeds using the goat 50K SNP chip. The genetic diversity analysis showed that Bugi toori goats have the highest inbreeding level, consistent with the highest linkage disequilibrium, lowest diversity and long run of heterozygosity segments. This indicates that this breed should be prioritized in future conservation activities. The population structure analysis revealed four fairly distinct clusters (including Bugi toori, Bari, Black Tapri and some Kamori) and three other breeds that are seemingly the results of admixture between these or related groups (some Kamori, Pateri, Tapri and White Tapri). The selection signatures were evaluated in each breed. A total of 2508 putative selection signals were reported. The 26 significant windows were identified in more than four breeds, and selection signatures spanned several genes that directly or indirectly influence traits included coat colour variation (KIT), reproduction (BMPR1B, GNRHR, INSL6, JAK2 and EGR4), body size (SOCS2), ear size (MSRB3) and milk composition (ABCG2, SPP1, CSN1S2, CSN2, CSN3 and PROLACTIN).

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals

Stature is affected by many polymorphisms of small effect in humans ¹ . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes^2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10^-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

Deciphering signature of selection affecting beef quality traits in Angus cattle

Artificial selection towards a desired phenotype/trait has modified the genomes of livestock dramatically that generated breeds that greatly differ in morphology, production and environmental adaptation traits. Angus cattle are among the famous cattle breeds developed for superior beef quality. This paper aimed at exploring genomic regions under selection in Angus cattle that are associated with meat quality traits and other associated phenotypes. The whole genome of 10 Angus cattle was compared with 11 Hanwoo (A-H) and 9 Jersey (A-J) cattle breeds using a cross-population composite likelihood ratio (XP-CLR) statistical method. The top 1% of the empirical distribution was taken as significant and annotated using UMD3.1. As a result, 255 and 210 genes were revealed under selection from A-H and A-J comparisons, respectively. The WebGestalt gene ontology analysis resulted in sixteen (A-H) and five (A-J) significantly enriched KEGG pathways. Several pathways associated with meat quality traits (insulin signaling, type II diabetes mellitus pathway, focal adhesion pathway, and ECM-receptor interaction), and feeding efficiency (olfactory transduction, tight junction, and metabolic pathways) were enriched. Genes affecting beef quality traits (e.g., FABP3, FTO, DGAT2, ACS, ACAA2, CPE, TNNI1), stature and body size (e.g., PLAG1, LYN, CHCHD7, RPS20), fertility and dystocia (e.g., ESR1, RPS20, PPP2R1A, GHRL, PLAG1), feeding efficiency (e.g., PIK3CD, DNAJC28, DNAJC3, GHRL, PLAG1), coat color (e.g., MC1-R) and genetic disorders (e.g., ITGB6, PLAG1) were found to be under positive selection in Angus cattle. The study identified genes and pathways that are related to meat quality traits and other phenotypes of Angus cattle. The findings in this study, after validation using additional or independent dataset, will provide useful information for the study of Angus cattle in particular and beef cattle in general.

Genomic regions and genes related to inter-population differences in body size in the ground beetle Carabus japonicus

Body size is a key trait in diversification among animal species, and revealing the gene regions responsible for body size diversification among populations or related species is important in evolutionary biology. We explored the genomic regions associated with body size differences in Carabus japonicus ground beetle populations by quantitative trait locus (QTL) mapping of F2 hybrids from differently sized parents from two populations using restriction site-associated DNA sequencing and de novo assembly of the beetle whole genome. The assembled genome had a total length of 191 Mb with a scaffold N50 of 0.73 Mb; 14,929 protein-coding genes were predicted. Three QTLs on different linkage groups had major effects on the overall size, which is composed chiefly of elytral length. In addition, we found QTLs on autosomal and X chromosomal linkage groups that affected head length and width, thoracic width, and elytral width. We determined the gene loci potentially related to control of body size in scaffolds of the genome sequence, which contained the QTL regions. The genetic basis of body size variation based on a small number of major loci would promote differentiation in body size in response to selection pressures related to variations in environmental conditions and inter-specific interactions.

Contribution of rare and low-frequency whole-genome sequence variants to complex traits variation in dairy cattle

Background

Whole-genome sequencing and imputation methodologies have enabled the study of the effects of genomic variants with low to very low minor allele frequency (MAF) on variation in complex traits. Our objective was to estimate the proportion of variance explained by imputed sequence variants classified according to their MAF compared with the variance explained by the pedigree-based additive genetic relationship matrix for 17 traits in Nordic Holstein dairy cattle.

Results

Imputed sequence variants were grouped into seven classes according to their MAF (0.001–0.01, 0.01–0.05, 0.05–0.1, 0.1–0.2, 0.2–0.3, 0.3–0.4 and 0.4–0.5). The total contribution of all imputed sequence variants to variance in deregressed estimated breeding values or proofs (DRP) for different traits ranged from 0.41 [standard error (SE) = 0.026] for temperament to 0.87 (SE = 0.011) for milk yield. The contribution of rare variants (MAF < 0.01) to the total DRP variance explained by all imputed sequence variants was relatively small (a maximum of 12.5% for the health index). Rare and low-frequency variants (MAF < 0.05) contributed a larger proportion of the explained DRP variances (>13%) for health-related traits than for production traits (<11%). However, a substantial proportion of these variance estimates across different MAF classes had large SE, especially when the variance explained by a MAF class was small. The proportion of DRP variance that was explained by all imputed whole-genome sequence variants improved slightly compared with variance explained by the 50 k Illumina markers, which are routinely used in bovine genomic prediction. However, the proportion of DRP variance explained by imputed sequence variants was lower than that explained by pedigree relationships, ranging from 1.5% for milk yield to 37.9% for the health index.

Conclusions

Imputed sequence variants explained more of the variance in DRP than the 50 k markers for most traits, but explained less variance than that captured by pedigree-based relationships. Although in humans partitioning variants into groups based on MAF and linkage disequilibrium was used to estimate heritability without bias, many of our bovine estimates had a high SE. For a reliable estimate of the explained DRP variance for different MAF classes, larger sample sizes are needed.

Electronic supplementary material

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

Study of male-mediated gene flow across a hybrid zone in the common shrew (Sorex araneus) using Y chromosome

Despite many studies, the impact of chromosome rearrangements on gene flow between chromosome races of the common shrew (Sorex araneus Linnaeus, 1758) remains unclear. Interracial hybrids form meiotic chromosome complexes that are associated with reduced fertility. Nevertheless comprehensive investigations of autosomal and mitochondrial markers revealed weak or no barrier to gene flow between chromosomally divergent populations. In a narrow zone of contact between the Novosibirsk and Tomsk races hybrids are produced with extraordinarily complex configurations at meiosis I. Microsatellite markers have not revealed any barrier to gene flow, but the phenotypic differentiation between races is greater than may be expected if gene flow was unrestricted. To explore this contradiction we analyzed the distribution of the Y chromosome SNP markers within this hybrid zone. The Y chromosome variants in combination with race specific autosome complements allow backcrosses to be distinguished and their proportion among individuals within the hybrid zone to be evaluated. The balanced ratio of the Y variants observed among the pure race individuals as well as backcrosses reveals no male mediated barrier to gene flow. The impact of reproductive unfitness of backcrosses on gene flow is discussed as a possible mechanism of the preservation of race-specific morphology within the hybrid zone.

Rensching cats and dogs: feeding ecology and fecundity trends explain variation in the allometry of sexual size dimorphism

The tendency for sexual size dimorphism (SSD) to increase with body mass in taxa where males are larger, and to decrease when females are larger, is known as Rensch's rule. In mammals, where the trend occurs, it is believed to be the result of a competitive advantage for larger males, while female mass is constrained by the energetics of reproduction. Here, we examine the allometry of SSD within the Felidae and Canidae, demonstrating distinctly different patterns: in felids, there is positive allometric scaling, while there is no trend in canids. We hypothesize that feeding ecology, via its effect on female spacing patterns, is responsible for the difference; larger male mass may be advantageous only where females are dispersed such that males can defend access to them. This is supported by the observation that felids are predominately solitary, and all are obligate carnivores. Similarly, carnivorous canids are more sexually dimorphic than insectivores and omnivores, but carnivory does not contribute to a Rensch effect as dietary variation occurs across the mass spectrum. The observed inter-familial differences are also consistent with reduced constraints on female mass in the canids, where litter size increases with body mass, versus no observable allometry in the felids.

Combined GWAS and 'guilt by association'-based prioritization analysis identifies functional candidate genes for body size in sheep

Background

Body size in sheep is an important indicator of productivity, growth and health as well as of environmental adaptation. It is a composite quantitative trait that has been studied with high-throughput genomic methods, i.e. genome-wide association studies (GWAS) in various mammalian species. Several genomic markers have been associated with body size traits and genes have been identified as causative candidates in humans, dog and cattle. A limited number of related GWAS have been performed in various sheep breeds and have identified genomic regions and candidate genes that partly account for body size variability. Here, we conducted a GWAS in Frizarta dairy sheep with phenotypic data from 10 body size measurements and genotypic data (from Illumina ovineSNP50 BeadChip) for 459 ewes.

Results

The 10 body size measurements were subjected to principal component analysis and three independent principal components (PC) were constructed, interpretable as width, height and length dimensions, respectively. The GWAS performed for each PC identified 11 significant SNPs, at the chromosome level, one on each of the chromosomes 3, 8, 9, 10, 11, 12, 19, 20, 23 and two on chromosome 25. Nine out of the 11 SNPs were located on previously identified quantitative trait loci for sheep meat, production or reproduction. One hundred and ninety-seven positional candidate genes within a 1-Mb distance from each significant SNP were found. A guilt-by-association-based (GBA) prioritization analysis (PA) was performed to identify the most plausible functional candidate genes. GBA-based PA identified 39 genes that were significantly associated with gene networks relevant to body size traits. Prioritized genes were identified in the vicinity of all significant SNPs except for those on chromosomes 10 and 12. The top five ranking genes were TP53, BMPR1A, PIK3R5, RPL26 and PRKDC.

Conclusions

The results of this GWAS provide evidence for 39 causative candidate genes across nine chromosomal regions for body size traits, some of which are novel and some are previously identified candidates from other studies (e.g. TP53, NTN1 and ZNF521). GBA-based PA has proved to be a useful tool to identify genes with increased biological relevance but it is subjected to certain limitations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0316-3) contains supplementary material, which is available to authorized users.

Functional confirmation of PLAG1 as the candidate causative gene underlying major pleiotropic effects on body weight and milk characteristics

A major pleiotropic quantitative trait locus (QTL) located at ~25 Mbp on bovine chromosome 14 affects a myriad of growth and developmental traits in Bos taurus and indicus breeds. These QTL have been attributed to two functional variants in the bidirectional promoter of PLAG1 and CHCHD7. Although PLAG1 is a good candidate for mediating these effects, its role remains uncertain given that these variants are also associated with expression of five additional genes at the broader locus. In the current study, we conducted expression QTL (eQTL) mapping of this region using a large, high depth mammary RNAseq dataset representing 375 lactating cows. Here we show that of the seven previously implicated genes, only PLAG1 and LYN are differentially expressed by QTL genotype, and only PLAG1 bears the same association signature of the growth and body weight QTLs. For the first time, we also report significant association of PLAG1 genotype with milk production traits, including milk fat, volume, and protein yield. Collectively, these data strongly suggest PLAG1 as the causative gene underlying this diverse range of traits, and demonstrate new effects for the locus on lactation phenotypes.

Imputation-Based Fine-Mapping Suggests That Most QTL in an Outbred Chicken Advanced Intercross Body Weight Line Are Due to Multiple, Linked Loci

The Virginia chicken lines have been divergently selected for juvenile body weight for more than 50 generations. Today, the high- and low-weight lines show a >12-fold difference for the selected trait, 56-d body weight. These lines provide unique opportunities to study the genetic architecture of long-term, single-trait selection. Previously, several quantitative trait loci (QTL) contributing to weight differences between the lines were mapped in an F₂-cross between them, and these were later replicated and fine-mapped in a nine-generation advanced intercross of them. Here, we explore the possibility to further increase the fine-mapping resolution of these QTL via a pedigree-based imputation strategy that aims to better capture the genetic diversity in the divergently selected, but outbred, founder lines. The founders of the intercross were high-density genotyped, and then pedigree-based imputation was used to assign genotypes throughout the pedigree. Imputation increased the marker density 20-fold in the selected QTL, providing 6911 markers for the subsequent analysis. Both single-marker association and multi-marker backward-elimination analyses were used to explore regions associated with 56-d body weight. The approach revealed several statistically and population structure independent associations and increased the mapping resolution. Further, most QTL were also found to contain multiple independent associations to markers that were not fixed in the founder populations, implying a complex underlying architecture due to the combined effects of multiple, linked loci perhaps located on independent haplotypes that still segregate in the selected lines.

Sexual size dimorphism in musteloids: An anomalous allometric pattern is explained by feeding ecology

Rensch's rule states that sexual size dimorphism (SSD) increases with body size in taxa where males are larger, and decreases when females are larger. The dominant explanation for the trend is currently that competitive advantage for males is greater in larger individuals, whereas female size is constrained by the energetics of rearing offspring. This rule holds for a variety of vertebrate taxa, and opposing trends are rare. We examine the allometry of SSD within the Musteloidea and demonstrate a hypo-allometry contrary to Rensch's rule, with lower SSD associated with larger body size. We provide evidence that feeding ecology is involved. Where diet promotes group-living, the optimal strategy for the males of larger species is often not to attempt to defend access to multiple females, obviating any competitive advantage of relatively greater size. We conclude that the effect of feeding ecology on mating systems may be a hitherto neglected factor explaining variation in SSD.

Comparison of gene-based rare variant association mapping methods for quantitative traits in a bovine population with complex familial relationships

BACKGROUND

There is growing interest in the role of rare variants in the variation of complex traits due to increasing evidence that rare variants are associated with quantitative traits. However, association methods that are commonly used for mapping common variants are not effective to map rare variants. Besides, livestock populations have large half-sib families and the occurrence of rare variants may be confounded with family structure, which makes it difficult to disentangle their effects from family mean effects. We compared the power of methods that are commonly applied in human genetics to map rare variants in cattle using whole-genome sequence data and simulated phenotypes. We also studied the power of mapping rare variants using linear mixed models (LMM), which are the method of choice to account for both family relationships and population structure in cattle.

RESULTS

We observed that the power of the LMM approach was low for mapping a rare variant (defined as those that have frequencies lower than 0.01) with a moderate effect (5 to 8 % of phenotypic variance explained by multiple rare variants that vary from 5 to 21 in number) contributing to a QTL with a sample size of 1000. In contrast, across the scenarios studied, statistical methods that are specialized for mapping rare variants increased power regardless of whether multiple rare variants or a single rare variant underlie a QTL. Different methods for combining rare variants in the test single nucleotide polymorphism set resulted in similar power irrespective of the proportion of total genetic variance explained by the QTL. However, when the QTL variance is very small (only 0.1 % of the total genetic variance), these specialized methods for mapping rare variants and LMM generally had no power to map the variants within a gene with sample sizes of 1000 or 5000.

CONCLUSIONS

We observed that the methods that combine multiple rare variants within a gene into a meta-variant generally had greater power to map rare variants compared to LMM. Therefore, it is recommended to use rare variant association mapping methods to map rare genetic variants that affect quantitative traits in livestock, such as bovine populations.

Natural history collections as windows on evolutionary processes

Natural history collections provide an immense record of biodiversity on Earth. These repositories have traditionally been used to address fundamental questions in biogeography, systematics and conservation. However, they also hold the potential for studying evolution directly. While some of the best direct observations of evolution have come from long-term field studies or from experimental studies in the laboratory, natural history collections are providing new insights into evolutionary change in natural populations. By comparing phenotypic and genotypic changes in populations through time, natural history collections provide a window into evolutionary processes. Recent studies utilizing this approach have revealed some dramatic instances of phenotypic change over short timescales in response to presumably strong selective pressures. In some instances, evolutionary change can be paired with environmental change, providing a context for potential selective forces. Moreover, in a few cases, the genetic basis of phenotypic change is well understood, allowing for insight into adaptive change at multiple levels. These kinds of studies open the door to a wide range of previously intractable questions by enabling the study of evolution through time, analogous to experimental studies in the laboratory, but amenable to a diversity of species over longer timescales in natural populations.

Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH ) with pedigree-based inbreeding coefficients (F_PED ). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.

Environmental temperature impact on bone and cartilage growth

Environmental temperature can have a surprising impact on extremity growth in homeotherms, but the underlying mechanisms have remained elusive for over a century. Limbs of animals raised at warm ambient temperature are significantly and permanently longer than those of littermates housed at cooler temperature. These remarkably consistent lab results closely resemble the ecogeographical tenet described by Allen's "extremity size rule," that appendage length correlates with temperature and latitude. This phenotypic growth plasticity could have adaptive significance for thermal physiology. Shortened extremities help retain body heat in cold environments by decreasing surface area for potential heat loss. Homeotherms have evolved complex mechanisms to maintain tightly regulated internal temperatures in challenging environments, including "facultative extremity heterothermy" in which limb temperatures can parallel ambient. Environmental modulation of tissue temperature can have direct and immediate consequences on cell proliferation, metabolism, matrix production, and mineralization in cartilage. Temperature can also indirectly influence cartilage growth by modulating circulating levels and delivery routes of essential hormones and paracrine regulators. Using an integrated approach, this article synthesizes classic studies with new data that shed light on the basis and significance of this enigmatic growth phenomenon and its relevance for treating human bone elongation disorders. Discussion centers on the vasculature as a gateway to understanding the complex interconnection between direct (local) and indirect (systemic) mechanisms of temperature-enhanced bone lengthening. Recent advances in imaging modalities that enable the dynamic study of cartilage growth plates in vivo will be key to elucidating fundamental physiological mechanisms of long bone growth regulation.

Alu-mediated recombination defect in IGF1R: haploinsufficiency in a patient with short stature

BACKGROUND

The insulin-like growth factor (IGF) receptor (IGF1R) is essential for normal development and growth. IGF1R mutations cause IGF-1 resistance resulting in intrauterine and postnatal growth failure. The phenotypic spectrum related to IGF1R mutations remains to be fully understood.

METHODS

Auxological and endocrinological data of a patient identified previously were assessed. The patient's fibroblasts were studied to characterize the IGF1R deletion, mRNA fate, protein expression and signalling capabilities.

RESULTS

The boy, who carries a heterozygous IGF1R exon 6 deletion caused by Alu element-mediated recombination and a heterozygous SHOX variant (p.Met240Ile), was born appropriate for gestational age but developed proportionate short stature postnatally. IGF-1 levels were low-normal. None of the stigmata associated with SHOX deficiency or sporadically observed in IGF1R mutation carriers were present. Nonsense-mediated mRNA decay led to a substantial decline of IGF1R dosage and IGF-1-dependent receptor autophosphorylation but not impaired downstream signalling.

CONCLUSION

We present the first detailed report of an intragenic IGF1R deletion identified in a patient who, apart from short stature, deviates from all established markers that qualify a growth-retarded child for IGF1R analysis. Although such children will usually escape routine clinical mutation screenings, they can contribute to the understanding of factors and mechanisms that cooperate with the IGF1R.

Describing the genetic architecture of epilepsy through heritability analysis

Epilepsy is highly heritable, but its genetic architecture is poorly understood. Speed et al. estimate the number of susceptibility loci, show that common variants account for the majority of heritability, and demonstrate that epilepsy consists of genetically distinct subtypes. They conclude that gene-based prediction models may have clinical utility in first-seizure settings.

Epilepsy is a disease with substantial missing heritability; despite its high genetic component, genetic association studies have had limited success detecting common variants which influence susceptibility. In this paper, we reassess the role of common variants on epilepsy using extensions of heritability analysis. Our data set consists of 1258 UK patients with epilepsy, of which 958 have focal epilepsy, and 5129 population control subjects, with genotypes recorded for over 4 million common single nucleotide polymorphisms. Firstly, we show that on the liability scale, common variants collectively explain at least 26% (standard deviation 5%) of phenotypic variation for all epilepsy and 27% (standard deviation 5%) for focal epilepsy. Secondly we provide a new method for estimating the number of causal variants for complex traits; when applied to epilepsy, our most optimistic estimate suggests that at least 400 variants influence disease susceptibility, with potentially many thousands. Thirdly, we use bivariate analysis to assess how similar the genetic architecture of focal epilepsy is to that of non-focal epilepsy; we demonstrate both significant differences (P = 0.004) and significant similarities (P = 0.01) between the two subtypes, indicating that although the clinical definition of focal epilepsy does identify a genetically distinct epilepsy subtype, there is also scope to improve the classification of epilepsy by incorporating genotypic information. Lastly, we investigate the potential value in using genetic data to diagnose epilepsy following a single epileptic seizure; we find that a prediction model explaining 10% of phenotypic variation could have clinical utility for deciding which single-seizure individuals are likely to benefit from immediate anti-epileptic drug therapy.

Heterozygous GHR gene mutation in a child with idiopathic short stature

Several monogenic defects have been reported to be associated with idiopathic short stature. Focusing on growth hormone receptor (GHR)-gene alterations, the heterozygosity of the same gene defect may be associated with a range of growth deficits. We found a heterozygous mutation (V144I) within exon 6 of the GHR gene in a patient with a low level of insulin-like growth factor I (IGF-I), normal level of GH, and severe short stature. Despite the lack of statistical difference, an overall tendency for reduced wt-GH-induction of GHR activation and Jak/Stat signalling in cells transiently expressing GHR-V144I alone or co-expressing wt-GHR compared to cells expressing only wt-GHR was found when GH doses were increased. Our results suggest that, although GHR sequence variants are responsible for some functional alterations commonly observed in children with idiopathic short stature, these changes may not explain all the height deficits observed in these subjects.