Evaluation of breeding strategies for polledness in dairy cattle using a newly developed simulation framework for quantitative and Mendelian traits

Background: Intensified selection of polled individuals has recently gained importance in predominantly horned dairy cattle breeds as an alternative to routine dehorning. The status quo of the current polled breeding pool of genetically-closely related artificial insemination sires with lower breeding values for performance traits raises questions regarding the effects of intensified selection based on this founder pool. Methods: We developed a stochastic simulation framework that combines the stochastic simulation software QMSim and a self-designed R program named QUALsim that acts as an external extension. Two traits were simulated in a dairy cattle population for 25 generations: one quantitative (QMSim) and one qualitative trait with Mendelian inheritance (i.e. polledness, QUALsim). The assignment scheme for qualitative trait genotypes initiated realistic initial breeding situations regarding allele frequencies, true breeding values for the quantitative trait and genetic relatedness. Intensified selection for polled cattle was achieved using an approach that weights estimated breeding values in the animal best linear unbiased prediction model for the quantitative trait depending on genotypes or phenotypes for the polled trait with a user-defined weighting factor. Results: Selection response for the polled trait was highest in the selection scheme based on genotypes. Selection based on phenotypes led to significantly lower allele frequencies for polled. The male selection path played a significantly greater role for a fast dissemination of polled alleles compared to female selection strategies. Fixation of the polled allele implies selection based on polled genotypes among males. In comparison to a base breeding scenario that does not take polledness into account, intensive selection for polled substantially reduced genetic gain for this quantitative trait after 25 generations. Reducing selection intensity for polled males while maintaining strong selection intensity among females, simultaneously decreased losses in genetic gain and achieved a final allele frequency of 0.93 for polled. Conclusions: A fast transition to a completely polled population through intensified selection for polled was in contradiction to the preservation of high genetic gain for the quantitative trait. Selection on male polled genotypes with moderate weighting, and selection on female polled phenotypes with high weighting, could be a suitable compromise regarding all important breeding aspects.

Effect of family relatedness on characteristics of estimated IBD probabilities in relation to precision of QTL estimates

Background: A random QTL effects model uses a function of probabilities that two alleles in the same or in different animals at a particular genomic position are identical by descent (IBD). Estimates of such IBD probabilities and therefore, modeling and estimating QTL variances, depend on marker polymorphism, strength of linkage and linkage disequilibrium of markers and QTL, and the relatedness of animals in the pedigree. The effect of relatedness of animals in a pedigree on IBD probabilities and their characteristics was examined in a simulation study. Results: The study based on nine multi-generational family structures, similar to a pedigree structure of a real dairy population, distinguished by an increased level of inbreeding from zero to 28 % across the studied population. Highest inbreeding level in the pedigree, connected with highest relatedness, was accompanied by highest IBD probabilities of two alleles at the same locus, and by lower relative variation coefficients. Profiles of correlation coefficients of IBD probabilities along the marked chromosomal segment with those at the true QTL position were steepest when the inbreeding coefficient in the pedigree was highest. Precision of estimated QTL location increased with increasing inbreeding and pedigree relatedness. A method to assess the optimum level of inbreeding for QTL detection is proposed, depending on population parameters. Conclusions: An increased overall relationship in a QTL mapping design has positive effects on precision of QTL position estimates. But the relationship of inbreeding level and the capacity for QTL detection depending on the recombination rate of QTL and adjacent informative marker is not linear. © 2010 Freyer et al., licensee BioMed Central Ltd.

The future for genetic studies in reproduction

Genetic factors contribute to risk of many common diseases affecting reproduction and fertility. In recent years, methods for genome-wide association studies(GWAS) have revolutionized gene discovery forcommontraits and diseases. Results of GWAS are documented in the Catalog of Published Genome-Wide Association Studies at the National Human Genome Research Institute and report over 70 publications for 32 traits and diseases associated with reproduction. These include endometriosis, uterine fibroids, age at menarche and age at menopause. Results that pass appropriate stringent levels of significance are generally well replicated in independent studies. Examples of genetic variation affecting twinning rate, infertility, endometriosis and age at menarche demonstrate that the spectrum of disease-related variants for reproductive traits is similar to most other common diseases.GWAS 'hits' provide novel insights into biological pathways and the translational value of these studies lies in discovery of novel gene targets for biomarkers, drug development and greater understanding of environmental factors contributing to disease risk. Results also show that genetic data can help define sub-types of disease and co-morbidity with other traits and diseases. To date, many studies on reproductive traits have used relatively small samples. Future genetic marker studies in large samples with detailed phenotypic and clinical information will yield new insights into disease risk, disease classification and co-morbidity for many diseases associated with reproduction and infertility.

The contribution of genes to cortical thickness and volume

We analyzed brain MRI data from 372 young adult twins toidentify cortical regions in which gray matter thickness and volume are influenced by genetics. This was achieved using an A/C/E structural equation model that divides the variance of these traits, at each point on the cortex, into additive genetic (A), shared (C), and unique environmental (E) components. A strong genetic influencewas found in frontal and parietal regions. Inaddition, we correlated cortical thickness with full-scale intelligence quotient for comparison with the A/C/E maps, and several regions where cortical structure was correlated with intelligence quotient are under genetic control. These cortical measures may be useful phenotypes to narrow the searchfor quantitative trait lociinfluencing brain structure.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

Genetic analyses in a sample of individuals with high or low BMD shows association with multiple Wnt pathway genes

INTRODUCTION Although the high heritability of BMD variation has long been established, few genes have been conclusively shown to affect the variation of BMD in the general population. Extreme truncate selection has been proposed as a more powerful alternative to unselected cohort designs in quantitative trait association studies. We sought to test these theoretical predictions in studies of the bone densitometry measures BMD, BMC, and femoral neck area, by investigating their association with members of the Wnt pathway, some of which have previously been shown to be associated with BMD in much larger cohorts, in a moderate-sized extreme truncate selected cohort (absolute value BMD Z-scores = 1.5-4.0; n = 344). MATERIALS AND METHODS Ninety-six tag-single nucleotide polymorphism (SNPs) lying in 13 Wnt signaling pathway genes were selected to tag common genetic variation (minor allele frequency [MAF] > 5% with an r(2) > 0.8) within 5 kb of all exons of 13 Wnt signaling pathway genes. The genes studied included LRP1, LRP5, LRP6, Wnt3a, Wnt7b, Wnt10b, SFRP1, SFRP2, DKK1, DKK2, FZD7, WISP3, and SOST. Three hundred forty-four cases with either high or low BMD were genotyped by Illumina Goldengate microarray SNP genotyping methods. Association was tested either by Cochrane-Armitage test for dichotomous variables or by linear regression for quantitative traits. RESULTS Strong association was shown with LRP5, polymorphisms of which have previously been shown to influence total hip BMD (minimum p = 0.0006). In addition, polymorphisms of the Wnt antagonist, SFRP1, were significantly associated with BMD and BMC (minimum p = 0.00042). Previously reported associations of LRP1, LRP6, and SOST with BMD were confirmed. Two other Wnt pathway genes, Wnt3a and DKK2, also showed nominal association with BMD. CONCLUSIONS This study shows that polymorphisms of multiple members of the Wnt pathway are associated with BMD variation. Furthermore, this study shows in a practical trial that study designs involving extreme truncate selection and moderate sample sizes can robustly identify genes of relevant effect sizes involved in BMD variation in the general population. This has implications for the design of future genome-wide studies of quantitative bone phenotypes relevant to osteoporosis.

The effect of transforming growth factor β1 gene polymorphisms in ankylosing spondylitis

Objectives. To determine whether genetic polymorphisms in or near the transforming growth factor β1 (TGFB1) locus were associated d with susceptibility to or severity of ankylosing spondylitis (AS). Methods. Five intragenic single-nucleotide polymorphisms (SNP) and three microsatellite markers flanking the TGFB1 locus were genotyped. Seven hundred and sixty-two individuals from 184 multiplex families were genotyped for the microsatellite markers and two of the promoter SNPs. One thousand and two individuals from 212 English and 170 Finnish families with AS were genotyped for all five intragenic SNPs. A structured questionnaire was used to assess the age of symptom onset, disease duration and disease severity scores, including the BASDAI (Bath Ankylosing Spondylitis Disease Activity Index) and BASFI (Bath Ankylosing Spondylitis Functional Index). Results. A weak association was noted between the rare TGFB1 + 1632 T allele and AS in the Finnish population (P = 0.04) and in the combined data set (P = 0.03). No association was noted between any other SNPs or SNP haplotype and AS, even among those families with positive non-parametric linkage scores. The TGFB1 +1632 polymorphism was also associated with a younger age of symptom onset (English population, allele 2 associated with age of onset greater by 4.2 yr, P = 0.05; combined data set, allele 2 associated with age of onset greater by 3.2 yr, P = 0.02). A haplotype of coding region SNPs (TGFB1 +869/ +915+1632 alleles 2/1/2) was associated with age of symptom onset in both the English parent-case trios and the combined data set (English data set, haplotype 2/1/2 associated with age of onset greater by 4.9 yr, P = 0.03; combined data set, haplotype 2/1/2 associated with greater age of onset by 4.2 yr, P = 0.006). Weak linkage with AS susceptibility was noted and the peak LOD score was 1.3 at distance 2 cM centromeric to the TGFB1 gene. No other linkage or association was found between quantitative traits and the markers. Conclusion. This study suggests that the polymorphisms within the TGFB1 gene play at most a small role in AS and that other genes encoded on chromosome 19 are involved in susceptibility to the disease.

Genome-wide association and genetic functional studies identify autism susceptibility candidate 2 gene (AUTS2) in the regulation of alcohol consumption

Alcohol consumption is a moderately heritable trait, but the genetic basis in humans is largely unknown, despite its clinical and societal importance. We report a genome-wide association study meta-analysis of approximately 2.5 million directly genotyped or imputed SNPs with alcohol consumption (gram per day per kilogram body weight) among 12 population-based samples of European ancestry, comprising 26,316 individuals, with replication genotyping in an additional 21,185 individuals. SNP rs6943555 in autism susceptibility candidate 2 gene (AUTS2) was associated with alcohol consumption at genome-wide significance (P = 4 x 10(-8) to P = 4 x 10(-9)). We found a genotype-specific expression of AUTS2 in 96 human prefrontal cortex samples (P = 0.026) and significant (P < 0.017) differences in expression of AUTS2 in whole-brain extracts of mice selected for differences in voluntary alcohol consumption. Down-regulation of an AUTS2 homolog caused reduced alcohol sensitivity in Drosophila (P < 0.001). Our finding of a regulator of alcohol consumption adds knowledge to our understanding of genetic mechanisms influencing alcohol drinking behavior.

Pre-emptive breeding to combine superior eating quality in tropical super sweet corn with resistance to major diseases

This report presents the process and outcomes of a five year project, which employed genetics and breeding approach for integrating disease resistance,agronomy and quality traits that enhances sustainable productivity improvement in sweet corn production. The report outlines a molecular markers based approach to introgress quantitative traits loci that are believed to contribute to resistance to downy mildew, a potentially devastating disease that threatens sweet corn and other similar crops. It also details the approach followed to integrate resistances for other major diseases such as southern rust (caused by Puccinia polysora Underw), Northern Corn Leaf Blight (Exserohilum turcicum) with improved agronomy and eating quality. The report explains the importance of heterosis (hybrid vigour) and combining ability in the development of useful sweet corn hybrids. It also explains the relevance of parental performance to predict its breeding value and the performance of its hybrids.

Compound heterozygous mutations in RIPPLY2 associated with vertebral segmentation defects

Segmentation defects of the vertebrae (SDV) are caused by aberrant somite formation during embryogenesis and result in irregular formation of the vertebrae and ribs. The Notch signal transduction pathway plays a critical role in somite formation and patterning in model vertebrates. In humans, mutations in several genes involved in the Notch pathway are associated with SDV, with both autosomal recessive (MESP2, DLL3, LFNG, HES7) and autosomal dominant (TBX6) inheritance. However, many individuals with SDV do not carry mutations in these genes. Using whole-exome capture and massive parallel sequencing, we identified compound heterozygous mutations in RIPPLY2 in two brothers with multiple regional SDV, with appropriate familial segregation. One novel mutation (c.A238T:p.Arg80*) introduces a premature stop codon. In transiently transfected C2C12 mouse myoblasts, the RIPPLY2 mutant protein demonstrated impaired transcriptional repression activity compared with wild-type RIPPLY2 despite similar levels of expression. The other mutation (c.240-4T>G), with minor allele frequency <0.002, lies in the highly conserved splice site consensus sequence 5' to the terminal exon. Ripply2 has a well-established role in somitogenesis and vertebral column formation, interacting at both gene and protein levels with SDV-associated Mesp2 and Tbx6. We conclude that compound heterozygous mutations in RIPPLY2 are associated with SDV, a new gene for this condition. © The Author 2014.

Phenotypic characters of rice landraces reveal independent lineages of short-grain aromatic indica rice

Rice landraces are lineages developed by farmers through artificial selection during the long-term domestication process. Despite huge potential for crop improvement, they are largely understudied in India. Here, we analyse a suite of phenotypic characters from large numbers of Indian landraces comprised of both aromatic and non-aromatic varieties. Our primary aim was to investigate the major determinants of diversity, the strength of segregation among aromatic and non-aromatic landraces as well as that within aromatic landraces. Using principal component analysis, we found that grain length, width and weight, panicle weight and leaf length have the most substantial contribution. Discriminant analysis can effectively distinguish the majority of aromatic from non-aromatic landraces. More interestingly, within aromatic landraces long-grain traditional Basmati and short-grain non-Basmati aromatics remain morphologically well differentiated. The present research emphasizes the general patterns of phenotypic diversity and finds out the most important characters. It also confirms the existence of very unique short-grain aromatic landraces, perhaps carrying signatures of independent origin of an additional aroma quantitative trait locus in the indica group, unlike introgression of specific alleles of the BADH2 gene from the japonica group as in Basmati. We presume that this parallel origin and evolution of aroma in short-grain indica landraces are linked to the long history of rice domestication that involved inheritance of several traits from Oryza nivara, in addition to O. rufipogon. We conclude with a note that the insights from the phenotypic analysis essentially comprise the first part, which will likely be validated with subsequent molecular analysis.

Transgenerational Effects of Early Life Starvation on Growth, Reproduction, and Stress Resistance in Caenorhabditis elegans.

Starvation during early development can have lasting effects that influence organismal fitness and disease risk. We characterized the long-term phenotypic consequences of starvation during early larval development in Caenorhabditis elegans to determine potential fitness effects and develop it as a model for mechanistic studies. We varied the amount of time that larvae were developmentally arrested by starvation after hatching ("L1 arrest"). Worms recovering from extended starvation grew slowly, taking longer to become reproductive, and were smaller as adults. Fecundity was also reduced, with the smallest individuals most severely affected. Feeding behavior was impaired, possibly contributing to deficits in growth and reproduction. Previously starved larvae were more sensitive to subsequent starvation, suggesting decreased fitness even in poor conditions. We discovered that smaller larvae are more resistant to heat, but this correlation does not require passage through L1 arrest. The progeny of starved animals were also adversely affected: Embryo quality was diminished, incidence of males was increased, progeny were smaller, and their brood size was reduced. However, the progeny and grandprogeny of starved larvae were more resistant to starvation. In addition, the progeny, grandprogeny, and great-grandprogeny were more resistant to heat, suggesting epigenetic inheritance of acquired resistance to starvation and heat. Notably, such resistance was inherited exclusively from individuals most severely affected by starvation in the first generation, suggesting an evolutionary bet-hedging strategy. In summary, our results demonstrate that starvation affects a variety of life-history traits in the exposed animals and their descendants, some presumably reflecting fitness costs but others potentially adaptive.

Genetic and genomic analyses of musculoskeletal differences between BEH and BEL strains

Berlin high (BEH) and Berlin low (BEL) strains selected for divergent growth differ 3-fold in body weight. We aimed at examining muscle mass, which is a major contributor to body weight, by exploring anatomical characteristics of the soleus muscle, its fiber numbers and their cross sectional area (CSA), by analysing transcriptome of the gastrocnemius and by initiating quantitative trait locus (QTL) mapping. BEH muscles were 4-to-8 times larger compared to BEL strain. In sub-strain BEH+/+, mutant myostatin was replaced with a wild type allele, however, BEH+/+muscles still were 2-to-4 times larger compared to the BEL strain. BEH soleus contained 2-times more (P<0.0001) and 2-times larger in CSA (P<0.0001) fibers compared to BEL strain. In addition, soleus femoral attachment anomaly (SFAA) was observed in all BEL mice. One significant (chromosome 1) and four suggestive (chromosomes 3, 4, 6 and 9) muscle weight QTLs were mapped in 21-day old F2 intercross (n=296) between BEH and BEL strains. The frequency of SFAA incidence in the F2 and in the backcross to BEL strain (BCL) suggested the presence of more than one causative gene. Two suggestive SFAA QTLs were mapped in BCL, however, their peak markers were not associated with the phenotype in F2. RNA-Seq analysis revealed 2,148 differentially expressed (P<0.1) genes and 45,673 SNPs and >2,000 indels between BEH+/+ and BEL males. In conclusion, contrasting muscle traits, genomic and gene expression differences between BEH and BEL strains provide a promising model for the search of genes involved in muscle growth and musculoskeletal morphogenesis.

The role of adiposity in cardiometabolic traits:a Mendelian randomization analysis

The association between adiposity and cardiometabolic traits is well known from epidemiological studies. Whilst the causal relationship is clear for some of these traits, for others it is not. We aimed to determine whether adiposity is causally related to various cardiometabolic traits using the Mendelian randomization approach.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.