Comparison of breeding value prediction for two traits in a Nellore-Angus crossbred population using different Bayesian modeling methodologies

The objectives of this study were to 1) compare four models for breeding value prediction using genomic or pedigree information and 2) evaluate the impact of fixed effects that account for family structure. Comparisons were made in a Nellore-Angus population comprising F₂, F₃ and half-siblings to embryo transfer F₂ calves with records for overall temperament at weaning (TEMP; n = 769) and Warner-Bratzler shear force (WBSF; n = 387). After quality control, there were 34,913 whole genome SNP markers remaining. Bayesian methods employed were BayesB ($\tilde{\pi }$ = 0.995 or 0.997 for WBSF or TEMP, respectively) and BayesC (π = 0 and $\tilde{\pi }$), where $\tilde{\pi }$ is the ideal proportion of markers not included. Direct genomic values (DGV) from single trait Bayesian analyses were compared to conventional pedigree-based animal model breeding values. Numerically, BayesC procedures (using $\tilde{\pi }$) had the highest accuracy of all models for WBSF and TEMP ($\widehat{\rho }$_gĝ = 0.843 and 0.923, respectively), but BayesB had the least bias (regression of performance on prediction closest to 1, $\widehat{\beta }$_y,x = 2.886 and 1.755, respectively). Accounting for family structure decreased accuracy and increased bias in prediction of DGV indicating a detrimental impact when used in these prediction methods that simultaneously fit many markers.     

Mixed Venous and Arterial Pco2

The rebreathing method of measuring oxygenated mixed venous Pco₂ (P$\bar{v}$co₂) was originally introduced as a bloodless way to estimate arterial Pco₂ (Paco₂). It has become common practice to subtract 6 mm Hg from the P$\bar{v}$co₂ to obtain the Paco₂ but there are many circumstances in which this leads to an overestimate of the Paco₂. Measurements of P$\bar{v}$co₂ and Paco₂ in 19 patients have shown that a better approximation to Paco₂ under normal conditions of cardiac output and arterial O₂ saturation is Paco₂ = 0·8 P$\bar{v}$co₂. These studies also showed that the P$\bar{v}$co₂ — Paco₂ difference may be much wider, particularly in the presence of arterial unsaturation and a low cardiac output.The factors governing the venoarterial Pco₂ difference are reviewed and their magnitude is calculated to emphasize the complementary roles of measurements of P$\bar{v}$co₂ and Paco₂ in the assessment of patients with cardiorespiratory disease.     

A Variable Sampling Interval Synthetic Xbar Chart for the Process Mean

The usual practice of using a control chart to monitor a process is to take samples from the process with fixed sampling interval (FSI). In this paper, a synthetic X¯ control chart with the variable sampling interval (VSI) feature is proposed for monitoring changes in the process mean. The VSI synthetic X¯ chart integrates the VSI X¯ chart and the VSI conforming run length (CRL) chart. The proposed VSI synthetic X¯ chart is evaluated using the average time to signal (ATS) criterion. The optimal charting parameters of the proposed chart are obtained by minimizing the out-of-control ATS for a desired shift. Comparisons between the VSI synthetic X¯ chart and the existing X¯, synthetic X¯, VSI X¯ and EWMA X¯ charts, in terms of ATS, are made. The ATS results show that the VSI synthetic X¯ chart outperforms the other X¯ type charts for detecting moderate and large shifts. An illustrative example is also presented to explain the application of the VSI synthetic X¯ chart.     

Understanding the genetic diversity, spatial genetic structure and mating system at the hierarchical levels of fruits and individuals of a continuous Theobroma cacao population from the Brazilian Amazon

Understanding the mating patterns of populations of tree species is a key component of ex situ genetic conservation. In this study, we analysed the genetic diversity, spatial genetic structure (SGS) and mating system at the hierarchical levels of fruits and individuals as well as pollen dispersal patterns in a continuous population of Theobroma cacao in Pará State, Brazil. A total of 156 individuals in a 0.56 ha plot were mapped and genotyped for nine microsatellite loci. For the mating system analyses, 50 seeds were collected from nine seed trees by sampling five fruits per tree (10 seeds per fruit). Among the 156 individuals, 127 had unique multilocus genotypes, and the remaining were clones. The population was spatially aggregated; it demonstrated a significant SGS up to 15 m that could be attributed primarily to the presence of clones. However, the short seed dispersal distance also contributed to this pattern. Population matings occurred mainly via outcrossing, but selfing was observed in some seed trees, which indicated the presence of individual variation for self-incompatibility. The matings were also correlated, especially within ($\widehat{r}$_p(m)=0.607) rather than among the fruits ($\widehat{r}$_p(m)=0.099), which suggested that a small number of pollen donors fertilised each fruit. The paternity analysis suggested a high proportion of pollen migration (61.3%), although within the plot, most of the pollen dispersal encompassed short distances (28 m). The determination of these novel parameters provides the fundamental information required to establish long-term ex situ conservation strategies for this important tropical species.     

Relationship between the Absorption and Emission Spectra and the "Red Drop" in the Action Spectra of Fluorescence In Vivo

The Stepanov equation, relating the intensity of emission, f_e($\bar{v}$), at a given frequency, and that of absorption, k($\bar{v}$), at the same frequency, is applied, in its modified form (see equation 3 in text) to suspensions of Chlorella, Porphyridium, and Anacystis and to chlorophyll solutions. This application can reveal whether the yield of fluorescence, Φ($\bar{v}$), is constant, or changes with frequency. In Chlorella (green alga) a sharp drop of Φ($\bar{v}$) is indicated towards the lower frequencies (longer waves), beginning around $\bar{v}$ = 1.48 × 10⁴cm^-1 (680 mμ); the Φ($\bar{v}$) function calculated from the Stepanov equation is in fair agreement with the directly determined action spectrum for the excitation of chlorophyll fluorescence in this organism. In Porphyridium (red alga) and Anacystis (blue-green alga) application of the Stepanov equation supports the conclusions, derived from direct measurements, of a much earlier “red drop” of the fluorescence excitation spectra. Direct measurements suggest that the drop in Porphyridium may begin at about 1.53 × 10⁴cm^-1 (654 mμ); in Anacystis, it may begin already above 1.57 × 10⁴cm^-1 (<637mμ). These results confirm the relation, postulated earlier by Duysens and others, between the action spectra of photosynthesis and of chlorophyll a fluorescence in algal cells. The relation of these findings to spectroscopic evidence, suggesting the existence of two main chlorophyll a components in vivo, in green as well as in red and blue-green algae, is discussed.     

Modeling Caribou Movements: Seasonal Ranges and Migration Routes of the Central Arctic Herd

Migration is an important component of the life history of many animals, but persistence of large-scale terrestrial migrations is being challenged by environmental changes that fragment habitats and create obstacles to animal movements. In northern Alaska, the Central Arctic herd (CAH) of barren-ground caribou (Rangifer tarandus granti) is known to migrate over large distances, but the herd’s seasonal distributions and migratory movements are not well documented. From 2003–2007, we used GPS radio-collars to determine seasonal ranges and migration routes of 54 female caribou from the CAH. We calculated Brownian bridges to model fall and spring migrations for each year and used the mean of these over all 4 years to identify areas that were used repeatedly. Annual estimates of sizes of seasonal ranges determined by 90% fixed kernel utilization distributions were similar between summer and winter ($\bar{X}$ = 27,929 SE = 1,064 and $\bar{X}$ = 26,585 SE = 4912 km², respectively). Overlap between consecutive summer and winter ranges varied from 3.3–18.3%. Percent overlap between summer ranges used during consecutive years ($\bar{X}$ = 62.4% SE = 3.7%) was higher than for winter ranges ($\bar{X}$ = 42.8% SE = 5.9%). Caribou used multiple migration routes each year, but some areas were used by caribou during all years, suggesting that these areas should be managed to allow for continued utilization by caribou. Restoring migration routes after they have been disturbed or fragmented is challenging. However, prior knowledge of movements and threats may facilitate maintenance of migratory paths and seasonal ranges necessary for long-term persistence of migratory species.     

Determination of sedimentation coefficients of subcellular particles of rat liver homogenates. A theoretical and experimental approach

A simple and convenient method is described for the determination of the average sedimentation coefficients ($\overline{S}?values$) of subcellular particles in a homogenous solution by velocity sedimentation in the preparative ultracentrifuge (swinging-bucket rotor). Theoretically the method is based on the distribution of intact organelles between sediment and supernatant and the convergence of their sedimentation at high centrifugal effects. The data have been treated according to the rate equation of parallel sedimentation in horizontal cylindrical tubes which has been converted to a linear function. By measuring activities of marker enzymes in the sediments, the theory has been confirmed experimentally in the determination of $\overline{S}?values$ and the distribution between supernatant and sediment of intact mitochondria, lysosomes and peroxisomes of rat liver homogenates. The effect of particle concentration on the $\overline{S}?values$ of mitochondria has been determined.     

Rhinos in the Parks: An Island-Wide Survey of the Last Wild Population of the Sumatran Rhinoceros

In the 200 years since the Sumatran rhinoceros was first scientifically described (Fisher 1814), the range of the species has contracted from a broad region in Southeast Asia to three areas on the island of Sumatra and one in Kalimantan, Indonesia. Assessing population and spatial distribution of this very rare species is challenging because of their elusiveness and very low population number. Using an occupancy model with spatial dependency, we assessed the fraction of the total landscape occupied by Sumatran rhinos over a 30,345-km² survey area and the effects of covariates in the areas where they are known to occur. In the Leuser Landscape (surveyed in 2007), the model averaging result of conditional occupancy estimate was ψ^(SE[ψ^])=0.151(0.109) or 2,371.47 km², and the model averaging result of replicated level detection probability p^(SE[p^])=0.252(0.267); in Way Kambas National Park—2008: ψ^(SE[ψ^])=0.468(0.165) or 634.18 km², and p^(SE[p^])=0.138(0.571); and in Bukit Barisan Selatan National Park—2010: ψ^(SE[ψ^])=0.322(0.049) or 819.67 km², and p^(SE[p^])=0.365(0.42). In the Leuser Landscape, rhino occurrence was positively associated with primary dry land forest and rivers, and negatively associated with the presence of a road. In Way Kambas, occurrence was negatively associated with the presence of a road. In Bukit Barisan Selatan, occurrence was negatively associated with presence of primary dryland forest and rivers. Using the probabilities of site occupancy, we developed spatially explicit maps that can be used to outline intensive protection zones for in-situ conservation efforts, and provide a detailed assessment of conserving Sumatran rhinos in the wild. We summarize our core recommendation in four points: consolidate small population, strong protection, determine the percentage of breeding females, and recognize the cost of doing nothing. To reduce the probability of poaching, here we present only the randomized location of site level occupancy in our result while retaining the overall estimation of occupancy for a given area.     

The effects of fasting on swimming performance in juvenile qingbo (Spinibarbus sinensis) at two temperatures

We measured the following variables to investigate the effects of fasting and temperature on swimming performance in juvenile qingbo (Spinibarbus sinensis): the critical swimming speed (U_crit), resting metabolic rate ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$ and active metabolic rate ${(\dot{\mathrm{M}}O}_{2active})$ of fish fasting for 0 (control), 1, 2 and 4 weeks at low and high acclimation temperatures (15 and 25 °C). Both fasting treatment and temperature acclimation had significant effects on all parameters measured (P<0.05). Fasting at the higher temperature had a negative effect on all measured parameters after 1 week (P<0.05). However, when acclimated to the lower temperature, fasting had a negative effect on U_crit until week 2 and on ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$, ${(\dot{\mathrm{M}}O}_{\mathrm{2active}})$ and metabolic scope (MS, ${(\dot{\mathrm{M}}O}_{\mathrm{2active}})$–${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$) until week 4 (P<0.05). The values of all parameters at the lower temperature were significantly lower than those at the higher temperature in the identical fasting period groups except for ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$ of the fish that fasted for 2 weeks. The relationship between fasting time (T) and U_crit was described as U_crit(15)=−0.302T²−0.800T+35.877 (r=0.781, n=32, P<0.001) and U_crit(25)=0.471T²−3.781T+50.097 (r=0.766, n=32, P<0.001₎ at 15 and 25 °C, respectively. The swimming performance showed less decrease in the early stage of fasting but more decrease in the later stage at the low temperature compared to the high temperature, which might be related to thermal acclimation time, resting metabolism, respiratory capacity, energy stores, enzyme activity in muscle tissue and energy substrate utilization changes with fasting between low and high temperatures. The divergent response of the swimming performance to fasting in qingbo at different temperatures might be an adaptive strategy to seasonal temperature and food resource variation in their habitat.     

A Global Analysis of Deforestation in Moist Tropical Forest Protected Areas

Protected areas (PAs) have been established to conserve tropical forests, but their effectiveness at reducing deforestation is uncertain. To explore this issue, we combined high resolution data of global forest loss over the period 2000–2012 with data on PAs. For each PA we quantified forest loss within the PA, in buffer zones 1, 5, 10 and 15 km outside the PA boundary as well as a 1 km buffer within the PA boundary. We analysed 3376 tropical and subtropical moist forest PAs in 56 countries over 4 continents. We found that 73% of PAs experienced substantial deforestation pressure, with >0.1% a⁻¹ forest loss in the outer 1 km buffer. Forest loss within PAs was greatest in Asia (0.25% a⁻¹) compared to Africa (0.1% a⁻¹), the Neotropics (0.1% a⁻¹) and Australasia (Australia and Papua New Guinea; 0.03% a⁻¹). We defined performance (P) of a PA as the ratio of forest loss in the inner 1 km buffer compared to the loss that would have occurred in the absence of the PA, calculated as the loss in the outer 1 km buffer corrected for any difference in deforestation pressure between the two buffers. To remove the potential bias due to terrain, we analysed a subset of PAs (n = 1804) where slope and elevation in inner and outer 1 km buffers were similar (within 1° and 100 m, respectively). We found 41% of PAs in this subset reduced forest loss in the inner buffer by at least 25% compared to the expected inner buffer forest loss (P<0.75). Median performance (P˜) of subset reserves was 0.87, meaning a reduction in forest loss within the PA of 13%. We found PAs were most effective in Australasia (P˜=0.16), moderately successful in the Neotropics (P˜=0.72) and Africa (P˜=0.83), but ineffective in Asia (P˜=1). We found many countries have PAs that give little or no protection to forest loss, particularly in parts of Asia, west Africa and central America. Across the tropics, the median effectiveness of PAs at the national level improved with gross domestic product per capita. Whilst tropical and subtropical moist forest PAs do reduce forest loss, widely varying performance suggests substantial opportunities for improved protection, particularly in Asia.     

Angle of Inclination of Tank-Treading Red Cells: Dependence on Shear Rate and Suspending Medium

Red cells suspended in solutions much more viscous than blood plasma assume an almost steady-state orientation when sheared above a threshold value of shear rate. This orientation is a consequence of the motion of the membrane around the red cell called tank-treading. Observed along the undisturbed vorticity of the shear flow, tank-treading red cells appear as slender bodies. Their orientation can be quantified as an angle of inclination (θ) of the major axis with respect to the undisturbed flow direction. We measured θ using solution viscosities (η₀) and shear rates (γ˙) covering one and three orders of magnitude, respectively. At the lower values of η₀, θ was almost independent of γ˙. At the higher values of η₀, θ displayed a maximum at intermediate shear rates. The respective maximal values of θ increased by ∼10° from 10.7 to 104 mPas. After accounting for the absent membrane viscosity in models by using an increased cytoplasmic viscosity, their predictions of θ agree qualitatively with our data. Comparison of the observed variation of θ at constant γ˙ with model results suggests a change in the reference configuration of the shear stiffness of the membrane.     

Hunting,Exotic Carnivores,and Habitat Loss: Anthropogenic Effects on a Native Carnivore Community,Madagascar

The wide-ranging, cumulative, negative effects of anthropogenic disturbance, including habitat degradation, exotic species, and hunting, on native wildlife has been well documented across a range of habitats worldwide with carnivores potentially being the most vulnerable due to their more extinction prone characteristics. Investigating the effects of anthropogenic pressures on sympatric carnivores is needed to improve our ability to develop targeted, effective management plans for carnivore conservation worldwide. Utilizing photographic, line-transect, and habitat sampling, as well as landscape analyses and village-based bushmeat hunting surveys, we provide the first investigation of how multiple forms of habitat degradation (fragmentation, exotic carnivores, human encroachment, and hunting) affect carnivore occupancy across Madagascar’s largest protected area: the Masoala-Makira landscape. We found that as degradation increased, native carnivore occupancy and encounter rates decreased while exotic carnivore occupancy and encounter rates increased. Feral cats (Felis species) and domestic dogs (Canis familiaris) had higher occupancy than half of the native carnivore species across Madagascar’s largest protected landscape. Bird and small mammal encounter rates were negatively associated with exotic carnivore occupancy, but positively associated with the occupancy of four native carnivore species. Spotted fanaloka (Fossa fossana) occupancy was constrained by the presence of exotic feral cats and exotic small Indian civet (Viverricula indica). Hunting was intense across the four study sites where hunting was studied, with the highest rates for the small Indian civet (x¯ = 90 individuals consumed/year), the ring-tailed vontsira (Galidia elegans) (x¯ = 58 consumed/year), and the fosa (Cryptoprocta ferox) (x¯ = 31 consumed/year). Our modeling results suggest hunters target intact forest where carnivore occupancy, abundance, and species richness, are highest. These various anthropogenic pressures and their effects on carnivore populations, especially increases in exotic carnivores and hunting, have wide-ranging, global implications and demand effective management plans to target the influx of exotic carnivores and unsustainable hunting that is affecting carnivore populations across Madagascar and worldwide.     

The Projection of a Test Genome onto a Reference Population and Applications to Humans and Archaic Hominins

We introduce a method for comparing a test genome with numerous genomes from a reference population. Sites in the test genome are given a weight, w, that depends on the allele frequency, x, in the reference population. The projection of the test genome onto the reference population is the average weight for each x, w¯(x). The weight is assigned in such a way that, if the test genome is a random sample from the reference population, then w¯(x)=1. Using analytic theory, numerical analysis, and simulations, we show how the projection depends on the time of population splitting, the history of admixture, and changes in past population size. The projection is sensitive to small amounts of past admixture, the direction of admixture, and admixture from a population not sampled (a ghost population). We compute the projections of several human and two archaic genomes onto three reference populations from the 1000 Genomes project—Europeans, Han Chinese, and Yoruba—and discuss the consistency of our analysis with previously published results for European and Yoruba demographic history. Including higher amounts of admixture between Europeans and Yoruba soon after their separation and low amounts of admixture more recently can resolve discrepancies between the projections and demographic inferences from some previous studies.