Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulae for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.

Growth data from a field trial of Quercus suber plants regenerated from selected trees and their half-sib progenies by somatic embryogenesis

The development of reliable clonal propagation technologies is a requisite for performing Multi-Varietal Forestry (MVF). Somatic embryogenesis is considered the tissue culture based method more suitable for operational breeding of forest trees. Vegetative propagation is very difficult when tissues are taken from mature donors, making clonal propagation of selected trees almost impossible. We have been able to induce somatic embryogenesis in leaves taken from mature oak trees, including cork oak (Quercus suber). This important species of the Mediterranean ecosystem produces cork regularly, conferring to this species a significant economic value. In a previous paper we reported the establishment of a field trial to compare the growth of plants of somatic origin vs zygotic origin, and somatic plants from mature trees vs somatic plants from juvenile seedlings. For that purpose somatic seedlings were regenerated from five selected cork oak trees and from young plants of their half-sib progenies by somatic embryogenesis. They were planted in the field together with acorn-derived plants of the same families. After the first growth period, seedlings of zygotic origin doubled the height of somatic seedlings, showing somatic plants of adult and juvenile origin similar growth. Here we provide data on height and diameter increases after two additional growth periods. In the second one, growth parameters of zygotic seedlings were also significantly higher than those of somatic ones, but there were not significant differences in height increase between seedlings and somatic plants of mature origin. In the third growth period, height and diameter increases of somatic seedlings cloned from the selected trees did not differ from those of zygotic seedlings, which were still higher than data from plants obtained from somatic embryos from the sexual progeny. Therefore, somatic seedlings from mature origin seem not to be influenced by a possible ageing effect, and plants from somatic embryos tend to minimize the initial advantage of plants from acorns

Half-sib progenies evaluation in velvet grass

O capim-lanudo (Holcus lanatus L.) possui grande potencial de utilização como pastagem de inverno nas regiões subtropicais, devido a seu bom estabelecimento, persistência, produção, resistência ao frio, palatabilidade e capacidade de afilhamento. O objetivo deste trabalho foi estimar a variabilidade genética e a herdabilidade de caracteres agronômicos e identificar progênies superiores. Foram avaliadas 60 progênies de meios-irmãos, para altura (NH) e diâmetro da planta na fase vegetativa, duração do ciclo vegetativo (ciclo), altura final (FH) e número de afilhos. Houve diferença entre progênies para todos os caracteres. As estimativas de herdabilidade foram de 38%, 32%, 92%, 57% e 64% para NH, diâmetro, ciclo, FH e afilhos, respectivamente. O maior ganho genético estimado foi de 30,77 % para o número de afilhos. Existe variabilidade para todos os caracteres e são esperados ganhos na seleção entre progênies.

Genetics of anthracnose panel canker disease resistance and its relationship with yield and growth characters in half-sib progenies of rubber tree (Hevea brasiliensis)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genotypic evaluation of Dent composite (Zea mays L) half-sib progenies in late harvest crop conditions

The present work was conducted at FCAV-UNESP, Campus de Jaboticabal, in late harvest crop time of 1991. One hundred interpopulational hybrids obtained from a top-cross between Dent Composite and Flint Composite populations were evaluated. The analyzed characteristics were: Spodoptera frugiperda damages, lodging, damaged ears percentage, and productivity. Estimates of heritability, variance between progenies and genetical gain for each character of available population were determined. Estimates of the progenies variance, like those from heritability for the different analysed characters, indicate that there is an adequate genetical variance for the utilization of that material on subsequent genetic breeding programs, though allowing genetical gains, on the following selection cycles for the characteristics: grain weight, lodging and damages caused by armyworm larvae, whereas for the other characteristics the obtained results were not satisfactory.

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y ps) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.

Genetic improvement of giant freshwater prawn in Vietnam

The giant freshwater prawn (Macrobrachium rosenbergii) or GFP is one of the most important freshwater crustacean species in the inland aquaculture sector of many tropical and subtropical countries. Since the 1990’s, there has been rapid global expansion of freshwater prawn farming, especially in Asian countries, with an average annual rate of increase of 48% between 1999 and 2001 (New, 2005). In Vietnam, GFP is cultured in a variety of culture systems, typically in integrated or rotational rice-prawn culture (Phuong et al., 2006) and has become one of the most common farmed aquatic species in the country, due to its ability to grow rapidly and to attract high market price and high demand. Despite potential for expanded production, sustainability of freshwater prawn farming in the region is currently threatened by low production efficiency and vulnerability of farmed stocks to disease. Commercial large scale and small scale GFP farms in Vietnam have experienced relatively low stock productivity, large size and weight variation, a low proportion of edible meat (large head to body ratio), scarcity of good quality seed stock. The current situation highlights the need for a systematic stock improvement program for GFP in Vietnam aimed at improving economically important traits in this species. This study reports on the breeding program for fast growth employing combined (between and within) family selection in giant freshwater prawn in Vietnam. The base population was synthesized using a complete diallel cross including 9 crosses from two local stocks (DN and MK strains) and a third exotic stock (Malaysian strain - MY). In the next three selection generations, matings were conducted between genetically unrelated brood stock to produce full-sib and (paternal) half-sib families. All families were produced and reared separately until juveniles in each family were tagged as a batch using visible implant elastomer (VIE) at a body size of approximately 2 g. After tags were verified, 60 to 120 juveniles chosen randomly from each family were released into two common earthen ponds of 3,500 m2 pond for a grow-out period of 16 to 18 weeks. Selection applied at harvest on body weight was a combined (between and within) family selection approach. 81, 89, 96 and 114 families were produced for the Selection line in the F0, F1, F2 and F3 generations, respectively. In addition to the Selection line, 17 to 42 families were produced for the Control group in each generation. Results reported here are based on a data set consisting of 18,387 body and 1,730 carcass records, as well as full pedigree information collected over four generations. Variance and covariance components were estimated by restricted maximum likelihood fitting a multi-trait animal model. Experiments assessed performance of VIE tags in juvenile GFP of different size classes and individuals tagged with different numbers of tags showed that juvenile GFP at 2 g were of suitable size for VIE tags with no negative effects evident on growth or survival. Tag retention rates were above 97.8% and tag readability rates were 100% with a correct assignment rate of 95% through to mature animal size of up to 170 g. Across generations, estimates of heritability for body traits (body weight, body length, cephalothorax length, abdominal length, cephalothorax width and abdominal width) and carcass weight traits (abdominal weight, skeleton-off weight and telson-off weight) were moderate and ranged from 0.14 to 0.19 and 0.17 to 0.21, respectively. Body trait heritabilities estimated for females were significantly higher than for males whereas carcass weight trait heritabilities estimated for females and males were not significantly different (P > 0.05). Maternal and common environmental effects for body traits accounted for 4 to 5% of the total variance and were greater in females (7 to 10%) than in males (4 to 5%). Genetic correlations among body traits were generally high in both sexes. Genetic correlations between body and carcass weight traits were also high in the mixed sexes. Average selection response (% per generation) for body weight (transformed to square root) estimated as the difference between the Selection and the Control group was 7.4% calculated from least squares means (LSMs), 7.0% from estimated breeding values (EBVs) and 4.4% calculated from EBVs between two consecutive generations. Favourable correlated selection responses (estimated from LSMs) were detected for other body traits (12.1%, 14.5%, 10.4%, 15.5% and 13.3% for body length, cephalothorax length, abdominal length, cephalothorax width and abdominal width, respectively) over three selection generations. Data in the second selection generation showed positive correlated responses for carcass weight traits (8.8%, 8.6% and 8.8% for abdominal weight, skeleton-off weight and telson-off weight, respectively). Data in the third selection generation showed that heritability for body traits were moderate and ranged from 0.06 to 0.11 and 0.11 to 0.22 at weeks 10 and 18, respectively. Body trait heritabilities estimated at week 10 were not significantly lower than at week 18. Genetic correlations between body traits within age and genetic correlations for body traits between ages were generally high. Overall our results suggest that growth rate responds well to the application of family selection and carcass weight traits can also be improved in parallel, using this approach. Moreover, selection for high growth rate in GFP can be undertaken successfully before full market size has been reached. The outcome of this study was production of an improved culture strain of GFP for the Vietnamese culture industry that will be trialed in real farm production environments to confirm the genetic gains identified in the experimental stock improvement program.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Detecting quantitative trait loci affecting beef tenderness on bovine chromosome 7 near calpastatin and lysyl oxidase

From a study of 3 large half-sib families of cattle, we describe linkage between DNA polymorphisms on bovine chromosome 7 and meat tenderness. Quantitative trait loci (QTL) for Longissimus lumborum peak force (LLPF) and Semitendonosis adhesion (STADH) were located to this map of DNA markers, which includes the calpastatin ( CAST) and lysyl oxidase (LOX) genes. The LLPF QTL has a maximum lodscore of 4.9 and allele substitution of approximately 0.80 of a phenotypic standard deviation, and the peak is located over the CAST gene. The STADH QTL has a maximum lodscore of 3.5 and an allele substitution of approximately 0.37 of a phenotypic standard deviation, and the peak is located over the LOX gene. This suggests 2 separate likelihood peaks on the chromosome. Further analyses of meat tenderness measures in the Longissimus lumborum, LLPF and LL compression (LLC), in which outlier individuals or kill groups are removed, demonstrate large shifts in the location of LLPF QTL, as well as confirming that there are indeed 2 QTL on bovine chromosome 7. We found that both QTL are reflected in both LLPF and LLC measurements, suggesting that both these components of tenderness, myofibrillar and connective tissue, are detected by both measurements in this muscle.

Quantitative trait loci detection and benefits from marker-assisted selection in dairy cattle

Knowing the chromosomal areas or actual genes affecting the traits under selection would add more information to be used in the selection decisions which would potentially lead to higher genetic response. The first objective of this study was to map quantitative trait loci (QTL) affecting economically important traits in the Finnish Ayrshire population. The second objective was to investigate the effects of using QTL information in marker-assisted selection (MAS) on the genetic response and the linkage disequilibrium between the different parts of the genome. Whole genome scans were carried out on a grand-daughter design with 12 half-sib families and a total of 493 sons. Twelve different traits were studied: milk yield, protein yield, protein content, fat yield, fat content, somatic cell score (SCS), mastitis treatments, other veterinary treatments, days open, fertility treatments, non-return rate, and calf mortality. The average spacing of the typed markers was 20 cM with 2 to 14 markers per chromosome. Associations between markers and traits were analyzed with multiple marker regression. Significance was determined by permutation and genome-wise P-values obtained by Bonferroni correction. The benefits from MAS were investigated by simulation: a conventional progeny testing scheme was compared to a scheme where QTL information was used within families to select among full-sibs in the male path. Two QTL on different chromosomes were modelled. The effects of different starting frequencies of the favourable alleles and different size of the QTL effects were evaluated. A large number of QTL, 48 in total, were detected at 5% or higher chromosome-wise significance. QTL for milk production were found on 8 chromosomes, for SCS on 6, for mastitis treatments on 1, for other veterinary treatments on 5, for days open on 7, for fertility treatments on 7, for calf mortality on 6, and for non-return rate on 2 chromosomes. In the simulation study the total genetic response was faster with MAS than with conventional selection and the advantage of MAS persisted over the studied generations. The rate of response and the difference between the selection schemes reflected clearly the changes in allele frequencies of the favourable QTL. The disequilibrium between the polygenes and QTL was always negative and it was larger with larger QTL size. The disequilibrium between the two QTL was larger with QTL of large effect and it was somewhat larger with MAS for scenarios with starting frequencies below 0.5 for QTL of moderate size and below 0.3 for large QTL. In conclusion, several QTL affecting economically important traits of dairy cattle were detected. Further studies are needed to verify these QTL, check their presence in the present breeding population, look for pleiotropy and fine map the most interesting QTL regions. The results of the simulation studies show that using MAS together with embryo transfer to pre-select young bulls within families is a useful approach to increase the genetic merit of the AI-bulls compared to conventional selection.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment. This fraction is result of undesirable genotype-by-environment interactions (GxE) and measured by the genetic correlation (rg) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of GxE over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels.

ScreeningCorymbiapopulations for resistance toPuccinia psidii

Disease screening to determine the threat Puccinia psidii poses to plantation and native eucalypts in Australia was undertaken in half-sib families of two contrasting eucalypt species, Eucalyptus cloeziana and E. argophloia. Artificial inoculation with a single-lesion isolate of P. psidii was used to screen these species for resistance to the biotype of P. psidii established in Australia. The objective was to characterize resistance to P. psidii within these two distinct species: E. argophloia, a vulnerable species with a narrow distribution, and E. cloeziana, a species with a broad and extensive distribution in Queensland. Results for E. cloeziana indicate that inland provenances are more resistant to P. psidii infection than provenances from coastal regions. Heritability estimates for the two assessment systems used (resistance on a 1-to-5 ordinal scale verses resistance on a 0-to-1 binomial scale) were low to high (0.24 to 0.63) for E. argophloia and moderate to high (0.4 to 0.91) for E. cloeziana, indicating a significant level of additive genetic variance for rust resistance within the populations. This study demonstrates the potential to select resistant families within the tested populations and indicates that P. psidii could detrimentally affect these species in native forests, nurseries, and plantations.

Screening Eucalyptus cloeziana and E. argophloia populations for resistance to Puccinia psidii

Disease screening to determine the threat Puccinia psidii poses to plantation and native eucalypts in Australia was undertaken in half-sib families of two contrasting eucalypt species, Eucalyptus cloeziana and E. argophloia. Artificial inoculation with a single-lesion isolate of P. psidii was used to screen these species for resistance to the biotype of P. psidii established in Australia. The objective was to characterize resistance to P. psidii within these two distinct species: E. argophloia, a vulnerable species with a narrow distribution, and E. cloeziana, a species with a broad and extensive distribution in Queensland. Results for E. cloeziana indicate that inland provenances are more resistant to P. psidii infection than provenances from coastal regions. Heritability estimates for the two assessment systems used (resistance on a 1-to-5 ordinal scale verses resistance on a 0-to-1 binomial scale) were low to high (0.24 to 0.63) for E. argophloia and moderate to high (0.4 to 0.91) for E. cloeziana, indicating a significant level of additive genetic variance for rust resistance within the populations. This study demonstrates the potential to select resistant families within the tested populations and indicates that P. psidii could detrimentally affect these species in native forests, nurseries, and plantations. Disease screening to determine the threat Puccinia psidii poses to plantation and native eucalypts in Australia was undertaken in half-sib families of two contrasting eucalypt species, Eucalyptus cloeziana and E. argophloia. Artificial inoculation with a single-lesion isolate of P. psidii was used to screen these species for resistance to the biotype of P. psidii established in Australia. The objective was to characterize resistance to P. psidii within these two distinct species: E. argophloia, a vulnerable species with a narrow distribution, and E. cloeziana, a species with a broad and extensive distribution in Queensland. Results for E. cloeziana indicate that inland provenances are more resistant to P. psidii infection than provenances from coastal regions. Heritability estimates for the two assessment systems used (resistance on a 1-to-5 ordinal scale verses resistance on a 0-to-1 binomial scale) were low to high (0.24 to 0.63) for E. argophloia and moderate to high (0.4 to 0.91) for E. cloeziana, indicating a significant level of additive genetic variance for rust resistance within the populations. This study demonstrates the potential to select resistant families within the tested populations and indicates that P. psidii could detrimentally affect these species in native forests, nurseries, and plantations.

Genetic parameter estimates for growth traits at early stage of Pacific abalone, Haliotis discus hannai Ino

Heritability and genetic and phenotypic correlations were estimated for juvenile growth traits of Pacific abalone Haliotis discus hannai Ino. The estimates were calculated from shell length and shell width measurements on progeny resulting from 12 half-sib families and 36 full-sib families obtained using artificial fertilization of mating three females to each male. The measurements were taken at 10, 20 and 30 d after fertilization. It was found that heritability estimates based on sire component ranged from 0.23 to 0.36 for shell length and 0.21 to 0.32 for shell width. Heritability estimates from dam component were larger than those from sire component at three ages, indicating presence of maternal effects, non-additive genetic effects and common environmental effects. Phenotypic correlations were significant at three ages (P < 0.05), with values of 0.92, 0.93 and 0.92, respectively. Genetic correlations from the paternal half-sib correlation analysis were highly positive at three ages, with values of 0.50, 0.78 and 0.81, respectively. The results suggest that selective breeding is an effective approach to improving growth traits of Pacific abalone stocks.