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.

De Novo Assembly And Transcriptome Analysis Of The Rubber Tree (hevea Brasiliensis) And Snp Markers Development For Rubber Biosynthesis Pathways.

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

Microsatellite Markers For Urochloa Humidicola (poaceae) And Their Transferability To Other Urochloa Species.

Urochloa humidicola is a warm-season grass commonly used as forage in the tropics and is recognized for its tolerance to seasonal flooding. This grass is an important forage species for the Cerrado and Amazon regions of Brazil. U. humidicola is a polyploid species with variable ploidy (6X-9X) and facultative apomixis with high phenotypic plasticity. However, this apomixis and ploidy, as well as the limited knowledge of the genetic basis of the germplasm collection, have constrained genetic breeding activities, yet microsatellite markers may enable a better understanding of the species' genetic composition. This study aimed to develop and characterize new polymorphic microsatellite molecular markers in U. humidicola and to evaluate their transferability to other Urochloa species. A set of microsatellite markers for U. humidicola was identified from two new enriched genomic DNA libraries: the first library was constructed from a single sexual genotype and the second from a pool of eight apomictic genotypes selected on the basis of previous results. Of the 114 loci developed, 72 primer pairs presented a good amplification product, and 64 were polymorphic among the 34 genotypes tested. The number of bands per simple sequence repeat (SSR) locus ranged from 1 to 29, with a mean of 9.6 bands per locus. The mean polymorphism information content (PIC) of all loci was 0.77, and the mean discrimination power (DP) was 0.87. STRUCTURE analysis revealed differences among U. humidicola accessions, hybrids, and other Urochloa accessions. The transferability of these microsatellites was evaluated in four species of the genus, U. brizantha, U. decumbens, U. ruziziensis, and U. dictyoneura, and the percentage of transferability ranged from 58.33% to 69.44% depending on the species. This work reports new polymorphic microsatellite markers for U. humidicola that can be used for breeding programs of this and other Urochloa species, including genetic linkage mapping, quantitative trait loci identification, and marker-assisted selection.

QTL mapping for yield components in a tropical maize population using microsatellite markers

QTL mapping provides usefull information for breeding programs since it allows the estimation of genomic locations and genetic effects of chromossomal regions related to the expression of quantitative traits. The objective of this study was to map QTL related to several agronomic important traits associated with grain yield: ear weight (EW), prolificacy (PROL), ear number (NE), ear length (EL) and diameter (ED), number of rows on the ear (NRE) and number of kernels per row on the ear (NKPR). Four hundred F-2:3 tropical maize progenies were evaluated in five environments in Piracicaba, Sao Paulo, Brazil. The genetic map was previously estimated and had 117 microssatelite loci with average distance of 14 cM. Data was analysed using Composite Interval Mapping for each trait. Thirty six QTL were mapped and related to the expression of EW (2), PROL (3), NE (2), EL (5), ED (5), NRE (10), NKPR (5). Few QTL were mapped since there was high GxE interaction. Traits EW, PROL and EN showed high genetic correlation with grain yield and several QTL mapped to similar genomic regions, which could cause the observed correlation. However, further analysis using apropriate statistical models are required to separate linked versus pleiotropic QTL. Five QTL (named Ew1, Ne1, Ed3, Nre3 and Nre10) had high genetic effects, explaining from 10.8% (Nre3) to 16.9% (Nre10) of the phenotypic variance, and could be considered in further studies.

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Some factors complicate comparisons between linkage maps from different studies. This problem can be resolved if measures of precision, such as confidence intervals and frequency distributions, are associated with markers. We examined the precision of distances and ordering of microsatellite markers in the consensus linkage maps of chromosomes 1, 3 and 4 from two F 2 reciprocal Brazilian chicken populations, using bootstrap sampling. Single and consensus maps were constructed. The consensus map was compared with the International Consensus Linkage Map and with the whole genome sequence. Some loci showed segregation distortion and missing data, but this did not affect the analyses negatively. Several inversions and position shifts were detected, based on 95% confidence intervals and frequency distributions of loci. Some discrepancies in distances between loci and in ordering were due to chance, whereas others could be attributed to other effects, including reciprocal crosses, sampling error of the founder animals from the two populations, F(2) population structure, number of and distance between microsatellite markers, number of informative meioses, loci segregation patterns, and sex. In the Brazilian consensus GGA1, locus LEI1038 was in a position closer to the true genome sequence than in the International Consensus Map, whereas for GGA3 and GGA4, no such differences were found. Extending these analyses to the remaining chromosomes should facilitate comparisons and the integration of several available genetic maps, allowing meta-analyses for map construction and quantitative trait loci (QTL) mapping. The precision of the estimates of QTL positions and their effects would be increased with such information.

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production

Bioethanol is a biofuel produced mainly from the fermentation of carbohydrates derived from agricultural feedstocks by the yeast Saccharomyces cerevisiae. One of the most widely adopted strains is PE-2, a heterothallic diploid naturally adapted to the sugar cane fermentation process used in Brazil. Here we report the molecular genetic analysis of a PE-2 derived diploid (JAY270), and the complete genome sequence of a haploid derivative (JAY291). The JAY270 genome is highly heterozygous (similar to 2 SNPs/kb) and has several structural polymorphisms between homologous chromosomes. These chromosomal rearrangements are confined to the peripheral regions of the chromosomes, with breakpoints within repetitive DNA sequences. Despite its complex karyotype, this diploid, when sporulated, had a high frequency of viable spores. Hybrid diploids formed by outcrossing with the laboratory strain S288c also displayed good spore viability. Thus, the rearrangements that exist near the ends of chromosomes do not impair meiosis, as they do not span regions that contain essential genes. This observation is consistent with a model in which the peripheral regions of chromosomes represent plastic domains of the genome that are free to recombine ectopically and experiment with alternative structures. We also explored features of the JAY270 and JAY291 genomes that help explain their high adaptation to industrial environments, exhibiting desirable phenotypes such as high ethanol and cell mass production and high temperature and oxidative stress tolerance. The genomic manipulation of such strains could enable the creation of a new generation of industrial organisms, ideally suited for use as delivery vehicles for future bioenergy technologies.

Establishing a soybean germplasm core collection

Core collections are of strategic importance as they allow the use of a small part of a germplasm collection that is representative of the total collection. The objective of this study was to develop a soybean core collection of the USDA Soybean Germplasm Collection by comparing the results of random, proportional, logarithmic, multivariate proportional and multivariate logarithmic sampling strategies. All but the random sampling strategy used stratification of the entire collection based on passport data and maturity group classification. The multivariate proportional and multivariate logarithmic strategies made further use of qualitative and quantitative trait data to select diverse accessions within each stratum. The 18 quantitative trait data distribution parameters were calculated for each core and for the entire collection for pairwise comparison to validate the sampling strategies. All strategies were adequate for assembling a core collection. The random core collection best represented the entire collection in statistical terms. Proportional and logarithmic strategies did not maximize statistical representation but were better in selecting maximum variability. Multivariate proportional and multivariate logarithmic strategies produced the best core collections as measured by maximum variability conservation. The soybean core collection was established using the multivariate proportional selection strategy. (C) 2010 Elsevier B.V. All rights reserved.

Analysis of genomic and functional RFLP derived markers associated with sucrose content, fiber and yield QTLs in a sugarcane (Saccharum spp.) commercial cross

Expressed sequence tags derived markers have a great potential to be used in functional map construction and QTL tagging. In the present work, sugarcane genomic probes and expressed sequence tags having homology to genes, mostly involved in carbohydrate metabolism were used in RFLP assays to identify putative QTLs as well as their epistatic interactions for fiber content, cane yield, pol and tones of sugar per hectare, at two crop cycles in a progeny derived from a bi-parental cross of sugarcane elite materials. A hundred and twenty marker trait associations were found, of which 26 at both crop cycle and 32 only at first ratoon cane. A sucrose synthase derived marker was associated with a putative QTL having a high negative effect on cane yield and also with a QTL having a positive effect on Pol at both crop cycles. Fifty digenic epistatic marker interactions were identified for the four traits evaluated. Of these, only two were observed at both crop cycles.

Evaluation of algorithms used to order markers on genetic maps

When building genetic maps, it is necessary to choose from several marker ordering algorithms and criteria, and the choice is not always simple. In this study, we evaluate the efficiency of algorithms try (TRY), seriation (SER), rapid chain delineation (RCD), recombination counting and ordering (RECORD) and unidirectional growth (UG), as well as the criteria PARF (product of adjacent recombination fractions), SARF (sum of adjacent recombination fractions), SALOD (sum of adjacent LOD scores) and LHMC (likelihood through hidden Markov chains), used with the RIPPLE algorithm for error verification, in the construction of genetic linkage maps. A linkage map of a hypothetical diploid and monoecious plant species was simulated containing one linkage group and 21 markers with fixed distance of 3 cM between them. In all, 700 F(2) populations were randomly simulated with and 400 individuals with different combinations of dominant and co-dominant markers, as well as 10 and 20% of missing data. The simulations showed that, in the presence of co-dominant markers only, any combination of algorithm and criteria may be used, even for a reduced population size. In the case of a smaller proportion of dominant markers, any of the algorithms and criteria (except SALOD) investigated may be used. In the presence of high proportions of dominant markers and smaller samples (around 100), the probability of repulsion linkage increases between them and, in this case, use of the algorithms TRY and SER associated to RIPPLE with criterion LHMC would provide better results. Heredity (2009) 103, 494-502; doi:10.1038/hdy.2009.96; published online 29 July 2009

QTL mapping for reaction to Phaeosphaeria leaf spot in a tropical maize population

Phaeosphaeria leaf spot (PLS) is an important disease in tropical and subtropical maize (Zea mays, L.) growing areas, but there is limited information on its inheritance. Thus, this research was conducted to study the inheritance of the PLS disease in tropical maize by using QTL mapping and to assess the feasibility of using marker-assisted selection aimed to develop genotypes resistance to this disease. Highly susceptible L14-04B and highly resistant L08-05F inbred lines were crossed to develop an F(2) population. Two-hundred and fifty six F(2) plants were genotyped with 143 microsatellite markers and their F(2:3) progenies were evaluated at seven environments. Ten plants per plot were evaluated 30 days after silk emergence following a rating scale, and the plot means were used for analyses. The heritability coefficient on a progeny mean basis was high (91.37%), and six QTL were mapped, with one QTL on chromosomes 1, 3, 4, and 6, and two QTL on chromosome 8. The gene action of the QTL ranged from additive to partial dominance, and the average level of dominance was partial dominance; also a dominance x dominance epistatic effect was detected between the QTL mapped on chromosome 8. The phenotypic variance explained by each QTL ranged from 2.91 to 11.86%, and the joint QTL effects explained 41.62% of the phenotypic variance. The alleles conditioning resistance to PLS disease of all mapped QTL were in the resistant parental inbred L08-05F. Thus, these alleles could be transferred to other elite maize inbreds by marker-assisted backcross selection to develop hybrids resistant to PLS disease.

Is there an association between normal cognitive functioning and trinucleotide repeat expansion at loci involved in neurodegenerative disorders?

Recent reports have shown neurodegenerative disorders to be associated with abnormal expansions of a CAG trinucleotide repeat allele at various autosomal loci. While normal chromosomes have 14 to 44 repeats, disease chromosomes may have 60 to 84 repeats. The number of CAG repeats on mutant chromosomes correlates with increasing severity of disease or decreasing age at onset of symptoms. Since we are interested in identifying the many quantitative trait loci (QTL) influencing brain functioning, we examined the possibility that the number of CAG repeats in the normal size range at these loci are relevant to "normal" neural functioning. We have used 150 pairs of adolescent (aged 16 years) twins and their parents to examine allele size at the MJD, SCA1, and DRPLA loci in heterozygous normal individuals. These are part of a large ongoing project using cognitive and physiological measures to investigate the genetie influences on cognition, and an extensive protocol of tests is employed to assess some of the key components of intellectual functioning. This study selected to examine full-scale psychometric IQ (FSIQ) and a measure of information processing (choice reaction time) and working memory (slow wave amplitude). CAG repeat size was determined on an ABI Genescan system following multiplex PCR amplification. Quantitative genetic analyses were performed to determine QTL effects of MJD, SCA1, and DRPLA on cognitive functioning. Analyses are in progress and will be discussed.

Functional and Evolutionary Insights from the Genomes of Three Parasitoid Nasonia Species

We report here genome sequences and comparative analyses of three closely related parasitoid wasps: Nasonia vitripennis, N. giraulti, and N. longicornis. Parasitoids are important regulators of arthropod populations, including major agricultural pests and disease vectors, and Nasonia is an emerging genetic model, particularly for evolutionary and developmental genetics. Key findings include the identification of a functional DNA methylation tool kit; hymenopteran-specific genes including diverse venoms; lateral gene transfers among Pox viruses, Wolbachia, and Nasonia; and the rapid evolution of genes involved in nuclear-mitochondrial interactions that are implicated in speciation. Newly developed genome resources advance Nasonia for genetic research, accelerate mapping and cloning of quantitative trait loci, and will ultimately provide tools and knowledge for further increasing the utility of parasitoids as pest insect-control agents.

Congenic strains provide evidence that four mapped loci in chromosomes 2, 4, and 16 influence hypertension in the SHR

Aneas I, Rodrigues MV, Pauletti BA, Silva GJ, Carmona R, Cardoso L, Kwitek AE, Jacob HJ, Soler JM, Krieger JE. Congenic strains provide evidence that four mapped loci in chromosomes 2, 4, and 16 influence hypertension in the SHR. Physiol Genomics 37: 52-57, 2009. First published January 6, 2009; doi: 10.1152/physiolgenomics.90299.2008. - To dissect the genetic architecture controlling blood pressure (BP) regulation in the spontaneously hypertensive rat (SHR) we derived congenic rat strains for four previously mapped BP quantitative trait loci (QTLs) in chromosomes 2, 4, and 16. Target chromosomal regions from the Brown Norway rat (BN) averaging 13 - 29 cM were introgressed by marker-assisted breeding onto the SHR genome in 12 or 13 generations. Under normal salt intake, QTLs on chromosomes 2a, 2c, and 4 were associated with significant changes in systolic BP (13, 20, and 15 mmHg, respectively), whereas the QTL on chromosome 16 had no measurable effect. On high salt intake (1% NaCl in drinking water for 2 wk), the chromosome 16 QTL had a marked impact on SBP, as did the QTLs on chromosome 2a and 2c (18, 17, and 19 mmHg, respectively), but not the QTL on chromosome 4. Thus these four QTLs affected BP phenotypes differently: 1) in the presence of high salt intake (chromosome 16), 2) only associated with normal salt intake (chromosome 4), and 3) regardless of salt intake (chromosome 2c and 2a). Moreover, salt sensitivity was abrogated in congenics SHR. BN2a and SHR. BN16. Finally, we provide evidence for the influence of genetic background on the expression of the mapped QTLs individually or as a group. Collectively, these data reveal previously unsuspected nuances of the physiological roles of each of the four mapped BP QTLs in the SHR under basal and/or salt loading conditions unforeseen by the analysis of the F2 cross.

Genetic influence on the variance in P3 amplitude and latency

The P3(00) event-related potential (ERP) component is widely used as a measure of cognitive functioning and provides a sensitive electrophysiological index of the attentional and working memory demands of a task. This study investigated what proportion of the variance in the amplitude and latency of the P3, elicited in a delayed response working memory task, could be attributed to genetic factors. In 335 adolescent twin pairs and 48 siblings, the amplitude and latency of the P3 were examined at frontal, central, and parietal sites. Additive genetic factors accounted for 48% to 61% of the variance in P3 amplitude. Approximately one-third of the genetic variation at frontal sites was mediated by a common genetic factor that also influenced the genetic variation at parietal and central sites. Familial resemblance in P3 latency was due to genetic influence that accounted for 44% to 50% of the variance. Genetic covariance in P3 latency across sites was substantial, with a large part of the variance found at parietal, central, and frontal sites attributed to a common genetic factor. The findings provide further evidence that the P3 is a promising phenotype of neural activity of the brain and has the potential to be used in linkage and association analysis in the search for quantitative trait loci (QTLs) influencing cognition.