Genetic correlations and heritability estimates of reproductive traits of Nellore bulls

Estimates of phenotypic, genetics and residual variances for reproductive traits in 5903 Nellore bulls were obtained. The experimental model used was multiple trait derivative-free restricted maximum likelihood. The values obtained for heritability were 0.24 +/- 0.05 for scrotal circumference at 450 days of age and 0.37 +/- 0.05 at 21 months for age at the time of the breeding soundness evaluation; 0.24 +/- 0.05 and 0.26 +/- 0.05 for left and right testicle length; 0.29 +/- 0.05 and 0.31 +/- 0.05 for left and right testicle width; 0.12 +/- 0.04 for testicle format; 0.33 +/- 0.06 for testicle volume; 0.11 +/- 0.03 for gross motility; 0.08 +/- 0.03 for individual motility and 0.05 +/- 0.02 for spermatic vigor; 0.20 +/- 0.04, 0.03 +/- 0.02 and 0.19 +/- 0.04 for larger defects, smaller defects and total defects, respectively. The values for heritability for testicular biometric characteristics were moderate to high while the seminal characteristics, presented low values. Genetic correlations between scrotal circumference with all the reproductive traits were favorable, suggesting the scrotal circumference as a feature of choice in the selection of bulls.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Genetic diversity in an indigenous horse breed: implications for mating strategies and the control of future inbreeding

The Franches-Montagnes is an indigenous Swiss horse breed, with approximately 2500 foalings per year. The stud book is closed, and no introgression from other horse breeds was conducted since 1998. Since 2006, breeding values for 43 different traits (conformation, performance and coat colour) are estimated with a best linear unbiased prediction (BLUP) multiple trait animal model. In this study, we evaluated the genetic diversity for the breeding population, considering the years from 2003 to 2008. Only horses with at least one progeny during that time span were included. Results were obtained based on pedigree information as well as from molecular markers. A series of software packages were screened to combine best the best linear unbiased prediction (BLUP) methodology with optimal genetic contribution theory. We looked for stallions with highest breeding values and lowest average relationship to the dam population. Breeding with such stallions is expected to lead to a selection gain, while lowering the future increase in inbreeding within the breed.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Barnacle recurit density, biomass, and reproductive output measured in rocky intertidal habitats in Glasgow Head, on the Atlantic coast of Nova Scotia, Canada

Through a field experiment, we show that a predator has negative nonconsumptive effects (NCEs) on different life-history stages of the same prey species. Shortly before the recruitment season of the barnacle Semibalanus balanoides (May-June), we established experimental cages in rocky intertidal habitats in Nova Scotia, Canada. The cages were used to manipulate the presence and absence of dogwhelks, Nucella lapillus, the main predators of barnacles. At the centre of each cage, we installed a tile where barnacle pelagic larvae could settle and the resulting recruits grow. Mesh prevented caged dogwhelks from accessing the tiles, but allowed waterborne dogwhelk cues to reach the tiles. Results in May indicated that barnacle larvae settled preferentially on tiles from cages without dogwhelks. In November, at the end of the dogwhelk activity period and once the barnacle recruits had grown to adult size, barnacle body mass was lower in the presence of dogwhelks. This limitation may have resulted from a lower barnacle feeding activity with nearby dogwhelks, as found by a previous study. The observed larval and adult responses in barnacles are consistent with attempts to decrease predation risk. November data also indicated that dogwhelk cues limited barnacle reproductive output, a possible consequence of the limited growth of barnacles. Overall, this study suggests that a predator species might influence trait evolution in a prey species through NCEs on different life-history stages.

A simple method to localise pleiotropic susceptibility loci using univariate linkage analyses of correlated traits

Univariate linkage analysis is used routinely to localise genes for human complex traits. Often, many traits are analysed but the significance of linkage for each trait is not corrected for multiple trait testing, which increases the experiment-wise type-I error rate. In addition, univariate analyses do not realise the full power provided by multivariate data sets. Multivariate linkage is the ideal solution but it is computationally intensive, so genome-wide analysis and evaluation of empirical significance are often prohibitive. We describe two simple methods that efficiently alleviate these caveats by combining P-values from multiple univariate linkage analyses. The first method estimates empirical pointwise and genome-wide significance between one trait and one marker when multiple traits have been tested. It is as robust as an appropriate Bonferroni adjustment, with the advantage that no assumptions are required about the number of independent tests performed. The second method estimates the significance of linkage between multiple traits and one marker and, therefore, it can be used to localise regions that harbour pleiotropic quantitative trait loci (QTL). We show that this method has greater power than individual univariate analyses to detect a pleiotropic QTL across different situations. In addition, when traits are moderately correlated and the QTL influences all traits, it can outperform formal multivariate VC analysis. This approach is computationally feasible for any number of traits and was not affected by the residual correlation between traits. We illustrate the utility of our approach with a genome scan of three asthma traits measured in families with a twin proband.

Independent sources of condition dependency and multiple pathways determine a composite trait: lessons from carotenoid-based plumage colouration.

Many colour ornaments are composite traits consisting of at least four components, which themselves may be more complex, determined by independent evolutionary pathways, and potentially being under different environmental control. To date, little evidence exists that several different components of colour elaboration are condition dependent and no direct evidence exists that different ornamental components are affected by different sources of variation. For example, in carotenoid-based plumage colouration, one of the best-known condition-dependent ornaments, colour elaboration stems from both condition-dependent pigment concentration and structural components. Some environmental flexibility of these components has been suggested, but specifically which and how they are affected remains unknown. Here, we tested whether multiple colour components may be condition dependent, by using a comprehensive 3 × 2 experimental design, in which we carotenoid supplemented and immune challenged great tit nestlings (Parus major) and quantified effects on different components of colouration. Plumage colouration was affected by an interaction between carotenoid availability and immune challenge. Path analyses showed that carotenoid supplementation increased plumage saturation via feather carotenoid concentration and via mechanisms unrelated to carotenoid deposition, while immune challenge affected feather length, but not carotenoid concentration. Thus, independent condition-dependent pathways, affected by different sources of variation, determine colour elaboration. This provides opportunities for the evolution of multiple signals within components of ornamental traits. This finding indicates that the selective forces shaping the evolution of different components of a composite trait and the trait's signal content may be more complex than believed so far, and that holistic approaches are required for drawing comprehensive evolutionary conclusions.

Radiation of the Australian flora: What can comparisons of molecular phylogenies across multiple taxa tell us about the evolution of diversity in present-day communities?

The Australian fossil record shows that from ca. 25 Myr ago, the aseasonal-wet biome (rainforest and wet heath) gave way to the unique Australian sclerophyll biomes dominated by eucalypts, acacias and casuarinas. This transition coincided with tectonic isolation of Australia, leading to cooler, drier, more seasonal climates. From 3 Myr ago, aridification caused rapid opening of the central Australian and zone. Molecular phylogenies with dated nodes have provided new perspectives on how these events could have affected the evolution of the Australian flora. During the Mid-Cenozoic (25-10 Myr ago) period of climatic change, there were rapid radiations in sclerophyll taxa, such as Banksia, eucalypts, pea-flowered legumes and Allocasuarina. At the same time, taxa restricted to the aseasonal-wet biome (Nothofagus, Podocarpaceae and Araucariaceae) did not radiate or were depleted by extinction. During the Pliocene aridification, two Eremean biome taxa (Lepidium and Chenopodiaceae) radiated rapidly after dispersing into Australia from overseas. It is clear that the biomes have different histories. Lineages in the aseasonal-wet biome are species poor, with sister taxa that are species rich, either outside Australia or in the sclerophyll biomes. In conjunction with the fossil record, this indicates depletion of the Australian aseasonal-wet biome from the Mid-Cenozoic. In the sclerophyll biomes, there have been multiple exchanges between the southwest and southeast, rather than single large endemic radiations after a vicariance event. There is need for rigorous molecular phylogenetic studies so that additional questions can be addressed, such as how interactions between biomes may have driven the speciation process during radiations. New studies should include the hither-to neglected monsoonal tropics.

The gall-inducing habit has evolved multiple times among the eriococcid scale insects (Sternorrhyncha : Coccoidea : Eriococcidae)

The habit of inducing plant galls has evolved multiple times among insects but most species diversity occurs in only a few groups, such as gall midges and gall wasps. This phylogenetic clustering may reflect adaptive radiations in insect groups in which the trait has evolved. Alternatively, multiple independent origins of galling may suggest a selective advantage to the habit. We use DNA sequence data to examine the origins of galling among the most speciose group of gall-inducing scale insects, the eriococcids. We determine that the galling habit has evolved multiple times, including four times in Australian taxa, suggesting that there has been a selective advantage to galling in Australia. Additionally, although most gall-inducing eriococcid species occur on Myrtaceae, we found that lineages feeding on Myrtaceae are no more likely to have evolved the galling habit than those feeding on other plant groups. However, most gall-inducing species-richness is clustered in only two clades (Apiomorpha and Lachnodius + Opisthoscelis), all of which occur exclusively on Eucalyptus s.s. The Eriococcidae and the large genus Eriococcus were determined to be non-monophyletic and each will require revision. (C) 2004 The Linnean Society of London.

Experimental autoimmune encephalomyelitis evolution was not modified by multiple infections with Strongyloides venezuelensis

P>According to the hygiene hypothesis, the increased incidence of allergic and autoimmune diseases in developed countries is mainly explained by the decreased contact between the human population and certain environmental agents as lactobacillus, mycobacteria and helminths. In this study, we evaluated the effect of multiple infections with Strongyloides venezuelensis on the development of experimental autoimmune encephalomyelitis (EAE) in Lewis rats. Multiple infections before EAE induction were not able to change the evolution of the disease. No alterations were observed in weight loss, clinical score and inflammation intensity at the central nervous system. The presence of significant levels of parasite-specific IgG1 but not IgG2b suggested a Th2 polarization. However, the percentage and absolute number of CD4+CD25+Foxp3+ T cells were not changed, being their levels in the spleen and lymph nodes of infected rats comparable to the ones found in normal animals. These results suggest that a Th2-polarized response without concomitant expansion of Foxp3+ regulatory T cells was not able to modify EAE progression. Even though these results do not threaten the hygiene hypothesis, they suggest that this paradigm might be an oversimplification. They also emphasize the need of a study to compare the immunoregulatory ability associated with different helminth spp.

Chromosomal structure: a selectable trait for evolution

Dissertation presented to obtain the Ph.D degree in Evolutionary Biology

Evolution of the gene regulatory network controlling flower dorsoventral asymmetry

Tese de Doutoramento em Biologia de Plantas.

Multiple causes of the evolution of dispersal in metapopulations

Abstract This work investigates the outcome of the interaction of the multiple causes of selection acting on dispersal in metapopulations. Dispersal, defined here as the ability of individuals to move out of their natal population to reproduce in an other one, has three main causes. First, population variability, as caused by random population extinctions, induces high incentives to disperse through the probability to recolonize an empty population and thus to escape competition for space. This adds to the second cause, kin competition avoidance where individuals in a crowded patch will benefit from the release of competition with relatives caused by dispersal. Dispersal may thus be viewed as an altruistic act. Third, dispersal might evolve as a strategy of avoiding inbred matings which are expected to bear fitness costs due to the presence of a mutation load. The interaction of inbreeding avoidance and kin competition is explored in chapter 2. Conditions conducive to the establishment of a high relatedness within population are expected to induce high dispersal through both kin competition avoidance and inbreeding avoidance. However, the dynamics of inbreeding depression is bound to depend on the level of gene flow as well as on the deleterious mutation parameters. Mutations more prone to settle a high level of inbreeding depression will select for increased dispersal. Chapter 3 investigates the effect of the mating system on the joint dynamics of dispersal and inbreeding depression. Higher inbreeding rates as those found in various mating systems lead to a more efficient purge of the deleterious mutations. However, this decrease in the costs of inbreeding are usually accompanied by a higher within deme relatedness which balances the decreased effect of inbreeding avoidance on the evolution of dispersal. Finally, population turnover, as found in most natural populations has a dual effect on dispersal. Indeed, it increases dispersal by the increased probability of winning a breeding slot in extinct demes it creates but, on the other hand, it counter-selects for dispersal through the slow establishment of unsaturated demic conditions which contribute to lower the local competition for space. Résumé Ce travail se propose d'étudier les effets conjoints des multiples causes de l'évolution de la dispersion en métapopulation. La dispersion, définie ici comme étant la capacité de quitter sa population d'origine pour se reproduire dans une antre population, possède trois principales causes. Premièrement, l'extinction aléatoire de populations sélectionne pour plus de dispersion car elle augmente la Probabilité de recoloniser un patch éteint et donc d'échapper à la compétition locale. La seconde cause, l'évitement de la compétition de parentèle, sélectionne pour plus de dispersion par les bénéfices qu'elle apporte par diminution de la compétition entre individus apparentés. Troisièmement, la dispersion évolue "comme stratégie d'évitement de la dépression de consanguinité présente dans des petites populations isolées. L'interaction entre l'évitement de la consanguinité et de la compétition de parentèle est étudiée dans le chapitre 2. Les conditions conduisant à l'établissement d'un fort apparentement à l'intérieur des populations sont celles qui génèrent le plus de sélection pour la dispersion. Cependant, la dynamique de la dépression de consanguinité est dépendante de la dispersion entre populations ainsi que des paramètres des mutations délétères. Les mutations créant le plus de dépression de consanguinité sont celles qui sélectionneront le plus pour de la dispersion. Le chapitre 3 s'intéresse aux effets du système de reproduction sur la dynamique conjointe du fardeau de mutation et de la dispersion. La purge des mutations délétère étant plus sévère dans des conditions de forte consanguinité, elle diminue les coûts de la consanguinité mais est habituellement accompagné par une augmentation de l'apparentement et donc l'effet peut être neutre sur la dispersion. Finalement, le turnover de populations a un effet dual sur la dispersion. La dispersion est sélectionnée par l'augmentation de la probabilité de gagner une place de reproduction dans des patchs éteints mais elle est également contre sélectionnée par la désaturation des patchs causée par l'extinction et la diminution de la compétition pour l'espace qui intervient dans ce cas.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.