Comparison of Methods for Handling Censored Records in Beef Fertility Data: Field Data

The purpose of this study was to compare methods for handling censored days to calving records in beef cattle data, and verify results of an earlier simulation study. Data were records from naturalservice matings of 33,176 first-calf females in Australian Angus herds.Three methods for handling censored records were evaluated. Censored records (records on noncalving females) were assigned penalty values on a within-contemporary group basis under the first method (DCPEN). Under the second method (DCSIM), censored records were drawn from their respective predictive truncated normal distributions, whereas censored records were deleted under the third method (DCMISS). Data were analyzed using a mixed linear model that included the fixed effects of contemporary group and sex of calf, linear and quadratic covariates for age at mating, and random effects of animal andresidual error. A Bayesian approach via Gibbs sampling was used to estimate variance components and predict breeding values.Posterior means (PM) (SD) of additive genetic variance for DCPEN, DCSIM, and DCMISS were 22.6d2 (4.2d2), 26.1d2(3.6d2), and 13.5d2(2.9d2),respectively. The PM (SD) of residual variance forDCPEN, DCSIM, and DCMISS were 431.4d2(5.0d2),371.4d2 (4.5d2), and 262.2d2(3.4d2), respectively. ThePM (SD) of heritability for DCPEN, DCSIM, andDCMISS were 0.05 (0.01), 0.07 (0.01), and 0.05 (0.01),respectively. Simulating trait records for noncalvingfemales resulted in similar heritability to the penaltymethod but lower residual variance. Pearson correlationsbetween posterior means of animal effects for sireswith more than 20 daughters with records were 0.99between DCPEN and DCSIM, 0.77 between DCPENand DCMISS, and 0.81 between DCSIM and DCMISS.Of the 424 sires ranked in the top 10% and bottom 10%of sires in DCPEN, 91% and 89%, respectively, werealso ranked in the top 10% and bottom 10% in DCSIM.Little difference was observed between DCPEN andDCSIM for correlations between posterior means of animaleffects for sires, indicating that no major rerankingof sires would be expected. This finding suggests littledifference between these two censored data handlingtechniques for use in genetic evaluation of days to calving

Studies on the value of incorporating the effect of dominance in genetic evaluations of dairy cattle, beef cattle and swine

Nonadditive genetic effects are currently ignored in national genetic evaluations of farm animals because of ignorance of thelevel of dominance variance for traits of interest and the difficult computational problems involved. Potential gains fromincluding the effects of dominance in genetic evaluations include “purification” of additive values and availability ofpredictions of specific combining abilities for each pair of prospective parents. This study focused on making evaluation withdominance effects feasible computationally and on ascertaining benefits of such an evaluation for dairy cattle, beef cattle,and swine. Using iteration on data, computing costs for evaluation with dominance effects included costs could be less thantwice expensive as with only an additive model. With Method Â, variance components could be estimated for problemsinvolving up to 10 millions equations. Dominance effects accounted for up to 10% of phenotypic variance; estimates werelarger for growth traits. As a percentage of additive variance, the estimate of dominance variance reached 78% for 21-d litterweight of swine and 47% for post weaning weight of beef cattle. When dominance effects are ignored, additive evaluationsare “contaminated”; effects are greatest for evaluations of dams in a single large family. These changes in ranking wereimportant for dairy cattle, especially for dams of full-sibs, but were less important for swine. Specific combining abilitiescannot be included in sire evaluations and need to be computed separately for each set of parents. The predictions of specificcombining abilities could be used in computerized mating programs via the Internet. Gains from including the dominanceeffect in genetic evaluations would be moderate but would outweigh expenditures to produce those evaluations

Sequencing and Analysis of the Mediterranean Amphioxus (Branchiostoma lanceolatum) Transcriptome

BACKGROUND: The basally divergent phylogenetic position of amphioxus (Cephalochordata), as well as its conserved morphology, development and genetics, make it the best proxy for the chordate ancestor. Particularly, studies using the amphioxus model help our understanding of vertebrate evolution and development. Thus, interest for the amphioxus model led to the characterization of both the transcriptome and complete genome sequence of the American species, Branchiostoma floridae. However, recent technical improvements allowing induction of spawning in the laboratory during the breeding season on a daily basis with the Mediterranean species Branchiostoma lanceolatum have encouraged European Evo-Devo researchers to adopt this species as a model even though no genomic or transcriptomic data have been available. To fill this need we used the pyrosequencing method to characterize the B. lanceolatum transcriptome and then compared our results with the published transcriptome of B. floridae. RESULTS: Starting with total RNA from nine different developmental stages of B. lanceolatum, a normalized cDNA library was constructed and sequenced on Roche GS FLX (Titanium mode). Around 1.4 million of reads were produced and assembled into 70,530 contigs (average length of 490 bp). Overall 37% of the assembled sequences were annotated by BlastX and their Gene Ontology terms were determined. These results were then compared to genomic and transcriptomic data of B. floridae to assess similarities and specificities of each species. CONCLUSION: We obtained a high-quality amphioxus (B. lanceolatum) reference transcriptome using a high throughput sequencing approach. We found that 83% of the predicted genes in the B. floridae complete genome sequence are also found in the B. lanceolatum transcriptome, while only 41% were found in the B. floridae transcriptome obtained with traditional Sanger based sequencing. Therefore, given the high degree of sequence conservation between different amphioxus species, this set of ESTs may now be used as the reference transcriptome for the Branchiostoma genus

Complete mitochondrial DNA sequences provide new insights into the Polynesian motif and the peopling of Madagascar

More than a decade of mitochondrial DNA (mtDNA) studies have given the 'Polynesian motif' renowned status as a marker for tracing the late-Holocene expansion of Austronesian speaking populations. Despite considerable research on the Polynesian motif in Oceania, there has been little equivalent work on the western edge of its expansion - leaving major issues unresolved regarding the motif's evolutionary history. This has also led to considerable uncertainty regarding the settlement of Madagascar. In this study, we assess mtDNA variation in 266 individuals from three Malagasy ethnic groups: the Mikea, Vezo, and Merina. Complete mtDNA genome sequencing reveals a new variant of the Polynesian motif in Madagascar; two coding region mutations define a Malagasy-specific sub-branch. This newly defined 'Malagasy motif' occurs at high frequency in all three ethnic groups (13-50%), and its phylogenetic position, geographic distribution, and estimated age all support a recent origin, but without conclusively identifying a specific source region. Nevertheless, the haplotype's limited diversity, similar to those of other mtDNA haplogroups found in our Malagasy groups, best supports a small number of initial settlers arriving to Madagascar through the same migratory process. Finally, the discovery of this lineage provides a set of new polymorphic positions to help localize the Austronesian ancestors of the Malagasy, as well as uncover the origin and evolution of the Polynesian motif itself

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

A chemical survey of exoplanets with ARIEL

Thousands of exoplanets have now been discovered with a huge range of masses, sizes and orbits: from rocky Earth-like planets to large gas giants grazing the surface of their host star. However, the essential nature of these exoplanets remains largely mysterious: there is no known, discernible pattern linking the presence, size, or orbital parameters of a planet to the nature of its parent star. We have little idea whether the chemistry of a planet is linked to its formation environment, or whether the type of host star drives the physics and chemistry of the planet’s birth, and evolution. ARIEL was conceived to observe a large number (~1000) of transiting planets for statistical understanding, including gas giants, Neptunes, super-Earths and Earth-size planets around a range of host star types using transit spectroscopy in the 1.25–7.8 μm spectral range and multiple narrow-band photometry in the optical. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres which should show minimal condensation and sequestration of high-Z materials compared to their colder Solar System siblings. Said warm and hot atmospheres are expected to be more representative of the planetary bulk composition. Observations of these warm/hot exoplanets, and in particular of their elemental composition (especially C, O, N, S, Si), will allow the understanding of the early stages of planetary and atmospheric formation during the nebular phase and the following few million years. ARIEL will thus provide a representative picture of the chemical nature of the exoplanets and relate this directly to the type and chemical environment of the host star. ARIEL is designed as a dedicated survey mission for combined-light spectroscopy, capable of observing a large and well-defined planet sample within its 4-year mission lifetime. Transit, eclipse and phase-curve spectroscopy methods, whereby the signal from the star and planet are differentiated using knowledge of the planetary ephemerides, allow us to measure atmospheric signals from the planet at levels of 10–100 part per million (ppm) relative to the star and, given the bright nature of targets, also allows more sophisticated techniques, such as eclipse mapping, to give a deeper insight into the nature of the atmosphere. These types of observations require a stable payload and satellite platform with broad, instantaneous wavelength coverage to detect many molecular species, probe the thermal structure, identify clouds and monitor the stellar activity. The wavelength range proposed covers all the expected major atmospheric gases from e.g. H2O, CO2, CH4 NH3, HCN, H2S through to the more exotic metallic compounds, such as TiO, VO, and condensed species. Simulations of ARIEL performance in conducting exoplanet surveys have been performed – using conservative estimates of mission performance and a full model of all significant noise sources in the measurement – using a list of potential ARIEL targets that incorporates the latest available exoplanet statistics. The conclusion at the end of the Phase A study, is that ARIEL – in line with the stated mission objectives – will be able to observe about 1000 exoplanets depending on the details of the adopted survey strategy, thus confirming the feasibility of the main science objectives.Peer reviewedFinal Published versio

A Genome-Wide RNAi Screen for Factors Involved in Neuronal Specification in Caenorhabditis elegans

One of the central goals of developmental neurobiology is to describe and understand the multi-tiered molecular events that control the progression of a fertilized egg to a terminally differentiated neuron. In the nematode Caenorhabditis elegans, the progression from egg to terminally differentiated neuron has been visually traced by lineage analysis. For example, the two gustatory neurons ASEL and ASER, a bilaterally symmetric neuron pair that is functionally lateralized, are generated from a fertilized egg through an invariant sequence of 11 cellular cleavages that occur stereotypically along specific cleavage planes. Molecular events that occur along this developmental pathway are only superficially understood. We take here an unbiased, genome-wide approach to identify genes that may act at any stage to ensure the correct differentiation of ASEL. Screening a genome-wide RNAi library that knocks-down 18,179 genes (94% of the genome), we identified 245 genes that affect the development of the ASEL neuron, such that the neuron is either not generated, its fate is converted to that of another cell, or cells from other lineage branches now adopt ASEL fate. We analyze in detail two factors that we identify from this screen: (1) the proneural gene hlh-14, which we find to be bilaterally expressed in the ASEL/R lineages despite their asymmetric lineage origins and which we find is required to generate neurons from several lineage branches including the ASE neurons, and (2) the COMPASS histone methyltransferase complex, which we find to be a critical embryonic inducer of ASEL/R asymmetry, acting upstream of the previously identified miRNA lsy-6. Our study represents the first comprehensive, genome-wide analysis of a single neuronal cell fate decision. The results of this analysis provide a starting point for future studies that will eventually lead to a more complete understanding of how individual neuronal cell types are generated from a single-cell embryo

Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P