POPREP: a generic report for population management

Genetic variation provides a basis upon which populations can be genetically improved. Management of animal genetic resources in order to minimize loss of genetic diversity both within and across breeds has recently received attention at different levels, e. g., breed, national and international levels. A major need for sustainable improvement and conservation programs is accurate estimates of population parameters, such as rate of inbreeding and effective population size. A software system (POPREP) is presented that automatically generates a typeset report. Key parameters for population management, such as age structure, generation interval, variance in family size, rate of inbreeding, and effective population size form the core part of this report. The report includes a default text that describes definition, computation and meaning of the various parameters. The report is summarized in two pdf files, named Population Structure and Pedigree Analysis Reports. In addition, results (e. g., individual inbreeding coefficients, rate of inbreeding and effective population size) are stored in comma-separate-values files that are available for further processing. Pedigree data from eight livestock breeds from different species and countries were used to describe the potential of POPREP and to highlight areas for further research.

LIFE CYCLE AND POPULATION STRUCTURE OF AEGLA PAULENSIS (DECAPODA: ANOMURA: AEGLIDAE)

We describe growth, longevity, sex ratio, reproductive period, and recruitment of Aegla paulensis from Jaragua Stale Park, Sao Paulo, Brazil (23 degrees 27'27.9 '' S; 46 degrees 45'32.3 '' W). The population was sampled monthly (September 2007 through August 2009) with the aid of traps. Over five thousand individuals were captured, sexed, measured (carapace length = CL) and inspected for reproductive traits (females only), and then released back to the sampling site. The pattern of the reproductive cycle was strongly seasonal (austral mid autumn through late winter), with a single recruitment pulse per year. The obtained von Bertalanffy growth equations were CL = 21.25[1-e(-0.041(t + 1.250))] and CL = 16.52[1-e(-0.049(t + 1.823))] for males and females, respectively. Males (mean CL +/- SD = 11.86 +/- 2.79 mm) attain larger sizes than females (mean CL +/- SD = 10.84 +/- 2.36 mm). Aegla paulensis reproduces twice during an estimated life span of 40.2 months for females and 33.9 months for males. Temporal variation of sex ratio showed a distinctive pattern characterized by a sequence of three distinct periods that repeated from one year to another, and which suggested that a behavioral component influence the proportion of sex in adult specimens sampled with traps during reproductive and non-reproductive periods.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Successful maintenance of a stingless bee population despite a severe genetic bottleneck

Stingless bees play an important ecological role as pollinators of many wild plant species in the tropics and have significant potential for the pollination of agricultural crops. Nevertheless, conservation efforts as well as commercial breeding programmes require better guidelines on the amount of genetic variation that is needed to maintain viable populations. In this context, we carried out a long-term genetic study on the stingless bee Melipona scutellaris to evaluate the population viability consequences of prolonged breeding from a small number of founder colonies. In particular, it was artificially imposed a genetic bottleneck by setting up a population starting from only two founder colonies, and continued breeding from it for a period of over 10 years in a location outside its natural area of occurrence. We show that despite a great reduction in the number of alleles present at both neutral microsatellite loci and the sex-determining locus relative to its natural source population, and an increased frequency in the production of sterile diploid males, the genetically impoverished population could be successfully bred and maintained for at least 10 years. This shows that in stingless bees, breeding from a small stock of colonies may have less severe consequences than previously suspected. In addition, we provide a simulation model to determine the number of colonies that are needed to maintain a certain number of sex alleles in a population, thereby providing useful guidelines for stingless bee breeding and conservation efforts.

Modeling associations between genetic markers using Bayesian networks

Motivation: Understanding the patterns of association between polymorphisms at different loci in a population ( linkage disequilibrium, LD) is of fundamental importance in various genetic studies. Many coefficients were proposed for measuring the degree of LD, but they provide only a static view of the current LD structure. Generative models (GMs) were proposed to go beyond these measures, giving not only a description of the actual LD structure but also a tool to help understanding the process that generated such structure. GMs based in coalescent theory have been the most appealing because they link LD to evolutionary factors. Nevertheless, the inference and parameter estimation of such models is still computationally challenging. Results: We present a more practical method to build GM that describe LD. The method is based on learning weighted Bayesian network structures from haplotype data, extracting equivalence structure classes and using them to model LD. The results obtained in public data from the HapMap database showed that the method is a promising tool for modeling LD. The associations represented by the learned models are correlated with the traditional measure of LD D`. The method was able to represent LD blocks found by standard tools. The granularity of the association blocks and the readability of the models can be controlled in the method. The results suggest that the causality information gained by our method can be useful to tell about the conservability of the genetic markers and to guide the selection of subset of representative markers.

Reconfiguration of distribution networks to minimize loss and disruption costs using genetic algorithms

In this paper a computational implementation of an evolutionary algorithm (EA) is shown in order to tackle the problem of reconfiguring radial distribution systems. The developed module considers power quality indices such as long duration interruptions and customer process disruptions due to voltage sags, by using the Monte Carlo simulation method. Power quality costs are modeled into the mathematical problem formulation, which are added to the cost of network losses. As for the EA codification proposed, a decimal representation is used. The EA operators, namely selection, recombination and mutation, which are considered for the reconfiguration algorithm, are herein analyzed. A number of selection procedures are analyzed, namely tournament, elitism and a mixed technique using both elitism and tournament. The recombination operator was developed by considering a chromosome structure representation that maps the network branches and system radiality, and another structure that takes into account the network topology and feasibility of network operation to exchange genetic material. The topologies regarding the initial population are randomly produced so as radial configurations are produced through the Prim and Kruskal algorithms that rapidly build minimum spanning trees. (C) 2009 Elsevier B.V. All rights reserved.

Structure of a Population of Hydromedusa maximiliani (Testudines, Chelidae) from Parque Estadual da Serra do Mar, an Atlantic Rainforest Preserve in Southeastern Brazil

We studied the structure of a population of Hydromedusa maximiliani associated with a stream in Parque Estadual da Serra do Mar, Nucleo Itutinga-Piloes, southeastern Brazil, between October 2004 and October 2005. Twenty-five individuals were captured, and a population size of 43.72 +/- 23.7 individuals was estimated. This value is similar to that of the population of Parque Estadual Carlos Botelho, another Atlantic forest reserve of southeastern Brazil. Males were recaptured more frequently than females, suggesting higher activity and/or greater movement of males.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Genetic variation for anatomic traits and volume shrinkage and their correlation with wood basic density in an Eucalyptus camaldulensis Dehn base population

Eucalyptus camaldulensis has great importance in Brazil because of their phenotypic plasticity for different environmental conditions, as soils, altitudes and rainfall. This study is an investigation of a base population of E. camaldulensis from Australia through a progeny test implanted in Selviria, MS. The trial was established in a randomized block design, with 25 families and 60 replications of single tree plots. Genetic parameters for anatomic traits and volume shrinkage were estimated, as well as their correlations with wood basic density. No significant differences among progenies were observed for the traits studied. The additive genetic variation coefficient at individual and among progeny levels ranged from low (0.26%) to high (16.98%). The narrow sense heritability at individual and family means levels also ranged from low (0.01) to high (0.87). This indicates that some traits are under strong genetic control and can be improved by selection. In the present situation, in order to attain the highest genetic gains, the sequential selection among and within progeny would be recommended.

Genetic population data of 12 STR loci of the PowerPlex (R) Y system in the state of Sao Paulo population (Southeast of Brazil)

Allele frequency distributions and population data for 12 Y-chromosomal short tandem repeats (STRs) included in the PowerPlex (R) Y Systems (Promega) were obtained for a sample of 200 healthy unrelated males living in S (a) over tildeo Paulo State (Southeast of Brazil). A total of 192 haplotypes were identified, of which 184 were unique and 8 were found in 2 individuals. The average gene diversity of the 12 Y-STR was 0.6746 and the haplotype diversity was 0.9996. Pairwise analysis confirmed that our population is more similar with the Italy, North Portugal and Spain, being more distant of the Japan. (c) 2007 Elsevier Ireland Ltd. All rights reserved.

The Genetic Speciation of Archaeological Fish Bone: A Feasibility Study from Southeast Queensland

Current genetic methods enable highly specific identification of DNA from modern fish bone. The applicability of these methods to the identification of archaeological fish bone was investigated through a study of a sample from late Holocene southeast Queensland sites. The resultant overall success rate of 2% indicates that DNA analysis is, as yet, not feasible for identifying fish bone from any given site. Taphonomic issues influencing the potential of genetic identification methods are raised and discussed in light of this result.

An update on genetic, structural and functional studies of arylamine N-acetyltransferases in eucaryotes and procaryotes

Arylamine N-acetyltransferase (NAT) was first identified as the inactivator of the anti-tubercular drug isoniazid, The enzyme was shown to catalyse the transfer of an acetyl group from acetyl-CoA to the terminal nitrogen of the hydrazine drug. The rate of inactivation of isoniazid was polymorphically distributed in the population and was one of the first examples of pharmacogenetic variation, NAT was identified recently in Mycobacterium tuberculosis and is a candidate for; modulating the response to isoniazid, Genome sequences have revealed many homologous members of this unique family of enzymes. The first three-dimensional structure of a member of the NAT family identifies a catalytic triad consisting of aspartate, histidine and cysteine proposed to form the activation mechanism. So far, all procaryotic NATs resemble the human enzyme which acetylates isoniazid (NAT2), Human NAT2 is characteristic of drug-metabolizing enzymes: it is found in liver and intestine, In humans and other mammals, there are up to three different isoenzymes. If only one isoenzyme is present, it is like human NAT1. Human NAT1 and its murine equivalent specifically acetylate the folate catabolite p-amino-benzoylglutamate. NAT1 and its murine homologue each have a ubiquitous tissue distribution and are expressed early in development at the blastocyst stage, During murine embryonic development, NAT is expressed in the developing neural tube. The proposed endogenous role of NAT in folate metabolism, and its multi-allelic nature, indicate that its role in development should be assessed further.

Comparative analysis of population structure of endemic Western Atlantic hermit crab Loxopagurus loxochelis (Decapoda, Anomura) in two regions of state of Sao Paulo, Brazil

This work characterized the population structure of the hermit crab Loxopagurus loxochelis (Moreira, 1901) in terms of size frequency distribution and sex ratio. Specimens were collected monthly, over a period of one year (from July 2002 to June 2003), in seven transects (from 5 to 35 m of depth) using fishing boat equipped with two double-rig trawl nets, in Caraguatatuba and Ubatuba regions (state of Sao Paulo, Brazil). A total of 366 hermit crabs were collected in Caraguatatuba [222 males (60.65%), 114 non-ovigerous females (31.15%) and 30 ovigerous females (8.20%)] and 126 hermit crabs in Ubatuba [81 males (64.28%), 38 non-ovigerous females (30.16%) and seven ovigerous females (5.56%)]. In Caraguatatuba the highest incidence of ovigerous females occurred during winter (July 2002), whereas in Ubatuba, the number was incipient. The cephalothoracic shield length ranged from 2.0 to 7.9mm (5.29 +/- 0.96mm) in Caraguatatuba, and from 2.7 to 7.5mm (5.32 +/- 0.95mm) in Ubatuba. The mean size of males was significantly larger than the mean size of females in both regions. Overall sex ratio was in favor of males (1.54:1 in Caraguatatuba and 1.9:1 in Ubatuba). Sexual dimorphism was recorded to L. loxochelis by the presence of males in the largest size classes, following the standard pattern observed in Decapoda. There was an unimodal size distribution for both sexes, with normal distributions in both regions. The higher number of males in relation to females may indicate the existence of different growth and mortality rates between the sexes. Despite of the different geomorphologic characteristics between Caraguatatuba and Ubatuba regions, the dynamics of development was similar for both populations.