Handling missingness when modeling the force of infection from clustered seroprevalence data

Modeling infectious diseases data is a relatively young research area in which clustering and stratification are key features. It is not unlikely for these data to have missing values. If values are missing completely at random, the analysis on the complete cases is valid. However, in practice this assumption is usually not fulfilled. This article shows the effect of ignoring missing data in modeling the force of infection of the bovine herpesvirus-1 in Belgian cattle and proposes the use of weighted generalized estimating equations with constrained fractional polynomials as a flexible modeling tool.     

In situ evaluation of fruit yield and estimation of repeatability coefficient for major fruit traits of umbu tree [Spondias tuberosa (Anacardiaceae)] in the semi-arid region of Brazil

The umbu tree (Spondias tuberosa) is a xerophytic woody species endemic to the northeast of Brazil. The present work reports a three-year in situ evaluation of major yield components in a random set of trees in Petrolina city, Pernambuco State, Brazil and the application of repeatability coefficient in the context of genetic resources utilization. Five traits were assessed: total fruit yield per plant (FY); total number of fruits per plant (NF); average fruit weight (FW); average pulp weight (PW) and average fruit diameter (FD). The values observed for FY ranged from 4.2–184 kg with mean of 61.5 kg. The values for NF ranged from 257 to 12,981 fruits/tree with mean 3,993 fruits. For FW, the range was from 10.82–23.36 g with mean of 16.03 g. The values for PW ranged from 7 to 17 g with mean of 11.2 g, while FD values varied from 2.5 to 3.5 cm with mean of 3.0 cm. Six methods were employed to estimate repeatability coefficient (r) for all parameters. The r ranged from 0.68–0.76 for NF; 0.87–0.89 for FW; 0.65–0.75 for FY; 0.64–0.78 for PW and 0.70–0.84 for FD. Narrower ranges across methods were observed to r-values greater than 0.80, as observed for FW and FD. Three to four years of measurements will be necessary for FY, NF, and PW and one year for FW to obtain a precision of 90% to provide a reliable identification of individuals for in situ or ex situ conservation.     

Effects of soil properties and trace metals on urease activities of calcareous soils

A study was conducted to investigate the relationship between urease activity and some physical, chemical, and microbiological properties of soils from central Iran. Inhibitory effects of Cr, Cd, and Pb on urease activity were also studied. Results indicated that no significant difference was observed between urease activity of field-moist and air-dried soils. Soil organic C and total N correlated highly significantly with urease activity, with r values of 0.899*** and 0.797***, respectively. There was also a significant correlation between urease activity and the number of bacteria grown on urea-agar media, with r value of 0.470*. A significant negative correlation (r =–0.492*) was observed between urease activity and electrical conductivity of saturation paste extracts. There were no significant correlations between urease activity and soil textural properties, pH, calcium carbonate equivalent, cation exchange capacity, and populations of soil bacteria on nutrient agar and population of soil fungi on potato dextrose agar. Both Cd and Pb inhibited urease activity to a similar extent and to a greater extent than did Cr.     

Incorporating parameter uncertainty into prediction intervals for spatial data modeled via a parametric variogram

In spatial predictions, researchers usually treat the estimated theoretical variogram parameters as known without error and ignore the variability of the parameter estimators. Although the prediction is still unbiased, the prediction error is usually underestimated. Therefore, the coverage probability of the prediction interval usually is lower than the nominal probability. A simulation study is performed to show how the coverage probability for prediction relates to the true range and sill of an exponential variogram. This article proposes two parametric bootstrap methods to incorporate the variability of the corresponding parameter estimators. A simulation study is performed to evaluate the coverage probability of these proposed methods. Finally, we apply the parametric bootstrap methods to a real dataset and compare the results with those from naive (i.e., treating estimated parameters as known) and Bayesian methods.     

Nitrogen fixation in biological soil crusts from southeast Utah,USA

Biological soil crusts can be the dominant source of N for arid land ecosystems. We measured potential N fixation rates biweekly for 2 years, using three types of soil crusts: (1) crusts whose directly counted cells were >98% Microcoleus vaginatus (light crusts); (2) crusts dominated by M. vaginatus, but with 20% or more of the directly counted cells represented by Nostoc commune and Scytonema myochrous (dark crusts); and (3) the soil lichen Collema sp. At all observation times, Collema had higher nitrogenase activity (NA) than dark crusts, which had higher NA than light crusts, indicating that species composition is critical when estimating N inputs. In addition, all three types of crusts generally responded in a similar fashion to climate conditions. Without precipitation within a week of collection, no NA was recorded, regardless of other conditions being favorable. Low (<1°C) and high (>26°C) temperatures precluded NA, even if soils were moist. If rain or snow melt had occurred 3 or less days before collection, NA levels were highly correlated with daily average temperatures of the previous 3 days (r²=0.93 for Collema crusts; r²=0.86 for dark crusts and r²=0.83 for light crusts) for temperatures between 1°C and 26°C. If a precipitation event followed a long dry period, NA levels were lower than if collection followed a time when soils were wet for extended periods (e.g., winter). Using a combination of data from a recording weather datalogger, time-domain reflectometry, manual dry-down curves, and N fixation rates at different temperatures, annual N input from the different crust types was estimated. Annual N input from dark crusts found at relatively undisturbed sites was estimated at 9 kg ha^-1 year^-1. With 20% cover of the N-fixing soil lichen Collema, inputs are estimated at 13 kg ha^-1 year^-1. N input from light crusts, generally indicating soil surface disturbance, was estimated at 1.4 kg ha^-1 year^-1. The rates in light crusts are expected to be highly variable, as disturbance history will determine cyanobacterial biomass and therefore N fixation rates.     

Comparison between phenetic characterisation using RAPD and ISSR markers and phenotypic data of cultivated chestnut (Castanea sativa Mill.)

Patterns of phenotypic and phenetic variability in sixPortuguese cultivars of chestnut (Castanea sativa Mill.) are evaluated. Morphological characterisationwas based on the quantification of seventeen traits. Varianceanalysis showed significant differences among cultivars, and cultivar× year for all the traits studied, and trees within cultivarsshowed also some significant differences for some of themorphological variables. A significant correlation was obtainedbetween length of the leaf blade and the percentage of unisexual andandrogynic inflorescence with the effective thermal index,accumulated rainfall from April to October and from July to October,or the accumulated temperature below seven during the dormant period.Principal Component and cluster analysis were performed to group thecultivars, according to their similarity coefficients. For molecularcharacterisation, 125 RAPD and 157 ISSR polymorphic markers wereamplified using 28 and 7 primers respectively. High level ofcongruence among the two marker systems (r =90.5%) was obtained from comparison of pheneticsimilarities based on the percentage of shared fragments. ISSRmarkers revealed important advantages over RAPDs, due to a higheffective multiplex ratio (12.5 for ISSR compared with 2.2 forRAPD analysis) and reproducibility. Although morphological andmolecular results are comparable, slight differences are showed incluster analysis UPGMA dendrograms. Molecular analysis explainedhomonym situations among 'Martainha' and'Longal' cultivars in Portugal.     

Variation and phylogeny of the triploid cultivated potatoSolanum chaucha Juz. et Buk. based on RAPD and isozyme markers

Ninety four accessions of the cultivated triploid potatoS. chaucha were analyzed and classified in genotypic groups using 9 isozyme loci and RAPD markers disclosed by 20 arbitrary 10-mer primers. Eight isozyme loci out of nine were polymorphic. A total of 22 allozymes were analyzed but none of them were specific for any genotypic group. About half (52%) of the 102 RAPD markers scored, were polymorphic, all of them showing polymorphism among groups and rarely within groups. Eighteen RAPD markers were specific for certain genotypes. The isozyme markers showed a certain amount of intra group variation which made classification less reliable than with RAPD markers. A total of 10 triploid genetic groups were discriminated using both techniques together. A single primer was found to be sufficient to distinguish all 10 groups. All varieties of a single group are considered to have been derived from the same cross and then clonally propagated, even though there is a high amount of morphological variation within a single genotypic group due probably to somatic mutations. RAPD markers have been shown to be more reliable in the classification of triploid potato varieties than other genetic markers like isozymes, proteins and morphological traits.     

QTL mapping of sheath blight resistance in a deep-water rice cultivar

Sheath blight, caused by the pathogen Rhizoctonia solani Kühn, is one of the most serious diseases of rice and leads to severe yield loss worldwide. A recombinant inbred line (RIL) population consisting of 121 lines was constructed from a cross between HH1B and RSB03, the latter of which is a deep-water rice variety. Five traits were used to evaluate sheath blight resistance, namely disease rating (DR), lesion length (LL), lesion height (LH), relative lesion length [RLL, the ratio of LL to plant height (PH)], and relative LH (RLH, the ratio of LH to PH). Using the RIL population and 123 molecular markers, we identified 28 quantitative trait loci (QTLs) for the five traits in two environments. These QTLs are located on nine chromosomes and most of them are environment specific. A major QTL for DR (qSBR1) on chromosome 1 was identified with contributions of 12.7% at Shanghai and 42.6% at Hainan, and it collocated with a QTL for PH. The allele at this locus from RSB03 enhances sheath blight resistance and increases PH. Another QTL for DR on chromosome 7 was adjacent to QTLs for heading date (HD) and four other disease traits. RSB03 also carries the resistant allele at this locus and shortens HD. The susceptible parent, HH1B, provides the resistance allele at the locus qSBR8, where QTLs for four other disease traits were identified. QTL mapping results showed that most QTLs for LL, LH, RLL, and RLH are collocated with QTLs for DR. Three QTLs for DR are independent from HD, PH, and four other disease traits, while four QTLs are closely related to HD and PH. Four QTLs for LL, LH, RLL, and RLH are independent from DR, HD, and PH, while there is only one region harboring QTLs for these four traits and HD. Correlation analysis and QTL mapping results indicated that LL, LH, RLL, and RLH might be important indices, like DR, for evaluating the level of resistance to rice sheath blight.     

Allele mining and haplotype discovery in barley candidate genes for drought tolerance

In the present study, allele mining was conducted on a panel of drought related candidate genes in a set of 96 barley genotypes using EcoTILLING, which is a variant of the targeting induced local lesions in genomes (TILLING) technology. Analyzing approximately 1.5 million basepairs in barley a total number of 94 verified unique haplotypes were identified in 18 amplicons designed for 9 genes. Overall, 185 single nucleotide polymorphisms (SNPs) and 46 insertions/deletions (INDELs) were detected with a mean of 1SNP/92 bp and 1INDEL/372 bp genomic sequence. Based on overlapping haplotype sequences, markers were developed for four candidate genes (HvARH1, HvSRG6, HvDRF1, HVA1), which allows distinguishing between the main haplotypes showing either differences in amino acid sequence or which have larger INDELs in the promoter region. As “proof of concept”, the HvARH1 and HvSRG6 haplotypes were tested for the level of abscisic acid-induced gene expression in subsets of genotypes belonging to different haplotype categories. An integrated database was developed to contain information about the genes, genotypes, and haplotypes analyzed in this study. The database supplies profound information about the natural variation in the tested drought related candidate genes providing a significant asset for further mapping studies dealing with this highly polygenic trait.     

Impact of epistasis and QTL × environmental interaction on the oil filling rate of soybean seed at different developmental stages

The oil accumulation in the developing soybean seed has been shown to be a dynamic process with different rates and activities at different phases affected by both genotype and environment. The objective of the present study was to investigate additive, epistatic and quantitative trait loci (QTL) × environment interaction (QE) effects of the QTL controlling oil filling rate in soybean seed. A total of 143 recombinant inbred lines (RILs) derived from the cross of Charleston and Dongnong 594 were used in this study to obtain 2 years of field data (2004 and 2005). A total of 26 QTL with significantly unconditional and conditional additive (a) effect and/or additive × environment interaction (ae) effect at different filling stages were identified on 14 linkage groups. Among the QTL with significant a effects, 18 QTL showed positive effects and 6 QTL had negative effects on seed filling rate of oil content during seed development. A total of 29 epistatic pairwise QTL underlying seed filling rate were identified at different filling stages. About 28 pairs of the QTL showed additive × additive epistatic (aa) effects and 14 pairs of the QTL showed aa × environment interaction (aae) effects at different filling stages. QTL with aa and aae (additive × additive × environment) effects appeared to vary at different filling stages. Our results demonstrated that oil filling rate in soybean seed were under genetic, developmental and environmental control.