Accuracy and precision of methods to estimate the number of parents contributing to a half-sib progeny array

Molecular technologies have made feasible large-scale studies of genetic parentage in nature by permitting the genotypic examination of hundreds or thousands of progeny. One common goal of such studies is to estimate the true number of unshared parents who contributed to a large half-sib progeny array. Here we introduce computer programs designed to count the number of gametotypes contributed by unshared parents to each such progeny array, as well as assess the accuracy and precision of various estimators for the true number of unshared parents via computer simulation. These simulations indicate that under most biological conditions (1) a traditional approach (the multilocus MINIMUM METHOD) that merely counts the number of distinct haplotypes in offspring and divides by 2L, where L is the number of loci assayed, often vastly underestimates the true number of unshared parents who contributed to a half-sib progeny array; (2) a recently developed HAPLOTYPES estimator is a considerable improvement over the MINIMUM METHOD when parental numbers are high; and (3) the accuracy and precision of the HAPLOTYPES estimator increase as marker polymorphism and sample size increase, or as reproductive skew and the number of parents contributing to the progeny array decrease. Generally, HAPLOTYPES-based estimates of parental numbers in large half-sib cohorts should improve the characterization of organismal reproductive strategies and mating systems from genetic data.     

Worker policing limits the number of reproductives in a ponerine ant

Reproductive division of labour is an essential feature of insect sociality, but the regulation of sterility among colony members remains incompletely understood. Ant workers and queens are morphologically divergent and workers are only capable of producing males in a colony, although they usually do not do so. Worker policing is one mechanism proposed for their infertility and it can be expressed as either aggressive inhibition of ovarian activity among workers or destruction of worker-laid eggs. A few studies have shown that workers with developed ovaries are preferentially attacked by nest-mates, but adequate demonstration of worker policing also requires evidence that these attacks result in the suppression of ovarian activity or death. We investigated worker policing in the ponerine ant Harpegnathos saltator in which workers are able to mate and replace the founding queen. Five colonies were each divided into two groups, one of which consisted exclusively of infertile workers. Some individuals in the orphaned groups began laying eggs during the three-week separation and upon reunification these were vigorously attacked by infertile workers of the other groups. The ovarian activity of these new egg layers became inhibited, as revealed by subsequent dissection of marked individuals. Worker policing in H. saltator appears to function primarily in preventing an excess of reproductive workers.     

Variation in chromosome number in the seedling progeny of a somaclone of Paspalum dilatatum

The somaclone,C39,derived by tissue culture from the obligate apomict Paspalum dilatatum cv Raki(2n=50),had 50 chromosomes and a karyotype apparently identical to Raki.SC2 seedlings of C39 showed a high degree of phenotypic variation which was often associated with increased chromosome numbers,but some of the variant seedlings were karyotypically indistinguishable from Raki or C39.Plants with increased chromosome numbers exhibited a high degree of intraplant chromosome variation(aneusomaty).In one of the SC2 seedlings,the chromosome number of root tip cells varied from 58 to 82 and in several other seedlings the range was more than 10.The results suggested that the ability to form seed apomictically was much reduced in C39 and that this plant showed some capacity for sexual reproduction and the resulting seedlings,with a chromosome number of about 70,were genetically unstable.Of 11 SC2 seedlings examined cytologically,6 did not produce any viable seed.Seedlings grown from seed of the remaining 5 plants showed that aneusomaty persisted in the SC3 generation.SC3 seedlings which were phenotypically similar to their maternal parent showed a similar range of chromosome numbers to that parent.Some of the SC3 seedlings exhibited an even wider range of chromosome numbers(e.g.56-136),and these plants were all dwarfs.     

Estimating the effective number of breeders from heterozygote excess in progeny

The heterozygote-excess method is a recently published method for estimating the effective population size (Ne). It is based on the following principle: When the effective number of breeders (Neb) in a population is small, the allele frequencies will (by chance) be different in males and females, which causes an excess of heterozygotes in the progeny with respect to Hardy-Weinberg equilibrium expectations. We evaluate the accuracy and precision of the heterozygote-excess method using empirical and simulated data sets from polygamous, polygynous, and monogamous mating systems and by using realistic sample sizes of individuals (15-120) and loci (5-30) with varying levels of polymorphism. The method gave nearly unbiased estimates of Neb under all three mating systems. However, the confidence intervals on the point estimates of Neb were sufficiently small (and hence the heterozygote-excess method useful) only in polygamous and polygynous populations that were produced by <10 effective breeders, unless samples included > approximately 60 individuals and 20 multiallelic loci.     

金叶银杏半同胞子代无性系的叶色和色素含量变化及呈色机制分析

对4月份至11月份金叶银杏‘万年金'( Ginkgo biloba ‘Wannianjin')32个半同胞子代无性系与亲本的叶色差异进行比较;比较了不同色系叶片的色素含量和比值及叶色参数(L*、a*和b*)的变化,分析了叶色参数与叶片色素含量的相关性;并观察了不同色系的叶绿体超微结构。结果表明：32个半同胞子代无性系可被分为金黄、浅黄、草绿和蓝绿4个色系。随时间推移,草绿和蓝绿色系叶片的总叶绿素( Chl)、叶绿素a( Chla)、叶绿素b( Chlb)和类胡萝卜素( Car)含量均呈“双峰型”变化趋势,Car/Chl和Car/Chla比值的变幅均较小;而金黄和浅黄色系叶片的上述色素含量呈“升高—降低—升高”变化趋势,Car/Chl和Car/Chla比值总体呈“迅速下降—相对稳定—缓慢升高”的变化趋势。各色系叶片的上述色素含量在夏季均不同程度下降,Car/Chlb比值变化差异较大,且金黄和浅黄色系的各色素含量均低于草绿和蓝绿色系。随时间推移,金黄和浅黄色系叶片的L*、a*和b*值以及草绿和蓝绿色系叶片的L*和b*值均先降低后升高,后2个色系的a*值则先升高后下降;并且,前2个色系的L*和b*值总体上显著高于后2个色系,而a*值则总体上低于后2个色系。金黄色系的Chla和Chl含量与L*和a*值显著负相关,而其Car/Chl和Car/Chla比值则与L*、a*和b*值显著或极显著正相关;浅黄色系的Chlb含量与a*值显著负相关,其Car/Chla比值与L*和b*值以及Car/Chlb比值与a*值均显著正相关;草绿色系的Chla含量与L*值显著负相关,其Car/Chla比值与L*和b*值以及Car/Chlb比值与a*值均显著正相关;这3个色系叶片的其余指标间以及蓝绿色系叶片的各指标间均无显著相关性。观察结果显示：金黄和浅黄色系的叶绿体基粒片层发育不健全,基粒片层可见但排列较疏松,且无明显垛叠,分布范围小而稀疏;蓝绿和草绿色系叶绿体的基粒类囊体垛叠层数均较多,基粒片层发达且排列紧致、整齐,分布范围大而稠密。综合分析结果表明：‘万年金'4个色系半同胞子代无性系叶片的呈色差异和叶色变化由多种因素控制,其中,Car/Chl和Car/Chla比值高且叶绿体基粒片层发育不健全是叶片呈黄色的主要原因。     

Integrating copy number polymorphisms into array CGH analysis using a robust HMM

MOTIVATION: Array comparative genomic hybridization (aCGH) is a pervasive technique used to identify chromosomal aberrations in human diseases, including cancer. Aberrations are defined as regions of increased or decreased DNA copy number, relative to a normal sample. Accurately identifying the locations of these aberrations has many important medical applications. Unfortunately, the observed copy number changes are often corrupted by various sources of noise, making the boundaries hard to detect. One popular current technique uses hidden Markov models (HMMs) to divide the signal into regions of constant copy number called segments; a subsequent classification phase labels each segment as a gain, a loss or neutral. Unfortunately, standard HMMs are sensitive to outliers, causing over-segmentation, where segments erroneously span very short regions. RESULTS: We propose a simple modification that makes the HMM robust to such outliers. More importantly, this modification allows us to exploit prior knowledge about the likely location of "outliers", which are often due to copy number polymorphisms (CNPs). By "explaining away" these outliers with prior knowledge about the locations of CNPs, we can focus attention on the more clinically relevant aberrated regions. We show significant improvements over the current state of the art technique (DNAcopy with MergeLevels) on previously published data from mantle cell lymphoma cell lines, and on published benchmark synthetic data augmented with outliers. AVAILABILITY: Source code written in Matlab is available from http://www.cs.ubc.ca/~sshah/acgh.     

Evolution of polyandry by reduction in progeny number variance in structured populations

When there is a variation in the quality of males in a population, multiple mating can lead to an increase in the genetic fitness of a female by reducing the variance of the progeny number. The extent of selective advantage obtainable by this process is investigated for a population subdivided into structured demes. It is seen that for a wide range of model parameters (deme size, distribution of male quality, local resource level), multiple mating leads to a considerable increase in the fitness. Frequency-dependent selection or a stable coexistence between polyandry and monandry can also result when the possible costs involved in multiple mating are taken into account.     

Exploiting noise in array CGH data to improve detection of DNA copy number change

Developing effective methods for analyzing array-CGH data to detect chromosomal aberrations is very important for the diagnosis of pathogenesis of cancer and other diseases. Current analysis methods, being largely based on smoothing and/or segmentation, are not quite capable of detecting both the aberration regions and the boundary break points very accurately. Furthermore, when evaluating the accuracy of an algorithm for analyzing array-CGH data, it is commonly assumed that noise in the data follows normal distribution. A fundamental question is whether noise in array-CGH is indeed Gaussian, and if not, can one exploit the characteristics of noise to develop novel analysis methods that are capable of detecting accurately the aberration regions as well as the boundary break points simultaneously? By analyzing bacterial artificial chromosomes (BACs) arrays with an average 1 mb resolution, 19 k oligo arrays with the average probe spacing <100 kb and 385 k oligo arrays with the average probe spacing of about 6 kb, we show that when there are aberrations, noise in all three types of arrays is highly non-Gaussian and possesses long-range spatial correlations, and that such noise leads to worse performance of existing methods for detecting aberrations in array-CGH than the Gaussian noise case. We further develop a novel method, which has optimally exploited the character of the noise, and is capable of identifying both aberration regions as well as the boundary break points very accurately. Finally, we propose a new concept, posteriori signal-to-noise ratio (p-SNR), to assign certain confidence level to an aberration region and boundaries detected.     

On the optimum structure of progeny tests for recessive alleles

Summary A method and an alternative for determining the optimum structure of a progeny test program designed to reduce the frequency of recessive alleles are described. Some of the limitations and possible applications are discussed.     

Continuous-index hidden Markov modelling of array CGH copy number data

MOTIVATION: In recent years, a range of techniques for analysis and segmentation of array comparative genomic hybridization (aCGH) data have been proposed. For array designs in which clones are of unequal lengths, are unevenly spaced or overlap, the discrete-index view typically adopted by such methods may be questionable or improved. RESULTS: We describe a continuous-index hidden Markov model for aCGH data as well as a Monte Carlo EM algorithm to estimate its parameters. It is shown that for a dataset from the BT-474 cell line analysed on 32K BAC tiling microarrays, this model yields considerably better model fit in terms of lag-1 residual autocorrelations compared to a discrete-index HMM, and it is also shown how to use the model for e.g. estimation of change points on the base-pair scale and for estimation of conditional state probabilities across the genome. In addition, the model is applied to the Glioblastoma Multiforme data used in the comparative study by Lai et al. (Lai,W.R. et al. (2005) Comparative analysis of algorithms for identifying amplifications and deletions in array CGH data. Bioinformatics, 21, 3763-3370.) giving result similar to theirs but with certain features highlighted in the continuous-index setting.     

On the basic reproduction number in a random environment

The concept of basic reproduction number $R_0$ in population dynamics is studied in the case of random environments. For simplicity the dependence between successive environments is supposed to follow a Markov chain. $R_0$ is the spectral radius of a next-generation operator. Its position with respect to 1 always determines population growth or decay in simulations, unlike another parameter suggested in a recent article (Hernandez-Suarez et al., Theor Popul Biol, doi:10.1016/j.tpb.2012.05.004, 2012). The position of the latter with respect to 1 determines growth or decay of the population’s expectation. $R_0$ is easily computed in the case of scalar population models without any structure. The main emphasis is on discrete-time models but continuous-time models are also considered.     

Analysis of copy number variations in Holstein cows identify potential mechanisms contributing to differences in residual feed intake

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. In this study, we performed an initial analysis of copy number variations (CNVs) using BovineHD SNP genotyping data from 147 Holstein cows identified as having high or low feed efficiency as estimated by residual feed intake (RFI). We detected 443 candidate CNV regions (CNVRs) that represent 18.4?Mb (0.6?%) of the genome. To investigate the functional impacts of CNVs, we created two groups of 30 individual animals with extremely low or high estimated breeding values (EBVs) for RFI, and referred to these groups as low intake (LI; more efficient) or high intake (HI; less efficient), respectively. We identified 240 (~9.0?Mb) and 274 (~10.2?Mb) CNVRs from LI and HI groups, respectively. Approximately 30–40?% of the CNVRs were specific to the LI group or HI group of animals. The 240 LI CNVRs overlapped with 137 Ensembl genes. Network analyses indicated that the LI-specific genes were predominantly enriched for those functioning in the inflammatory response and immunity. By contrast, the 274 HI CNVRs contained 177 Ensembl genes. Network analyses indicated that the HI-specific genes were particularly involved in the cell cycle, and organ and bone development. These results relate CNVs to two key variables, namely immune response and organ and bone development. The data indicate that greater feed efficiency relates more closely to immune response, whereas cattle with reduced feed efficiency may have a greater capacity for organ and bone development.     

DsbA: a protein-folding catalyst contributing to bacterial virulence

DsbA, a periplasmic thiol:disulphide oxidoreductase, catalyses the folding of various factors, among which are virulence determinants or the components of type III secretory machinery. It is also necessary for intracellular survival and cell-to-cell spread of the intracellular pathogen Shigella flexneri.     

Microsatellite analysis of maternal half-sib families of Quercus robur, pedunculate oak: II. inferring the number of pollen donors from the offspring

We present an approach to infer the number of pollen donors directly from genotype data of open- pollinated progeny of Quercus robur (pedunculate oak), a highly outcrossing tree species. The approach is based on closely linked, highly polymorphic codominant microsatellite markers. Initially the close linkages between three previously mapped microsatellite loci were confirmed by studies of linkage disequilibrium (LD). Then an approach to track the pollen donors contributing to maternal half-sib families (open-pollinated families) was developed by analysing haplotype arrays of closely linked microsatellite markers transmitted from the fathers to the progeny. Simulated data of five linked microsatellite loci segregating in eight open-pollinated families were used to study the relationship between the number of paternal chromosomes detected by this ”haplotype approach” and the number of diploid fathers contributing to the families. The results showed that the number of diploid pollen donors can be expressed as an exponential function of the number of paternal chromosomes inferred from the progeny. The 95% confidence interval of this regression function is used to determine the minimum number of fathers contributing to a genotyped open-pollinated family of Quercus robur. Finally this open-pollinated family is used to demonstrate the resolution obtained with the ”haplotype approach”. Six independent microsatellite loci were used to study relatedness among all pairs of pollen gametes that share a haplotype of three linked markers. The results suggest that the majority of such gametes are identical by descent from the same father. The ”haplotype approach” presented here can be used to monitor the number of contributing pollen donors in commercial seedlot samples from oak or any other outcrossing tree species for which closely linked, highly polymorphic, codominant genetic markers are available. Received: 5 June 1999 / Accepted: 30 July 1999     

Control of the change from group to pair behaviour by Hodotermes mossambicus reproductives

Groups of Hodotermes mossambicus reproductives exhibit social behaviour. Maintenance of this behaviour by an individual depends on frequent contact of its body by the antennae of other termites. The antennae of the recipient do not play a rôle. Pairs on the other hand exhibit antisocial behaviour. Development of pair behaviour occurs when antenna to body contact is reduced or absent. Furthermore, at least the proximal parts of the antennae are required for full development of pair behaviour. Pheromones are not implicated in the change from group to pair behaviour.     

Predicting the size of the progeny mapping population required to positionally clone a gene

A key frustration during positional gene cloning (map-based cloning) is that the size of the progeny mapping population is difficult to predict, because the meiotic recombination frequency varies along chromosomes. We describe a detailed methodology to improve this prediction using rice (Oryza sativa L.) as a model system. We derived and/or validated, then fine-tuned, equations that estimate the mapping population size by comparing these theoretical estimates to 41 successful positional cloning attempts. We then used each validated equation to test whether neighborhood meiotic recombination frequencies extracted from a reference RFLP map can help researchers predict the mapping population size. We developed a meiotic recombination frequency map (MRFM) for approximately 1400 marker intervals in rice and anchored each published allele onto an interval on this map. We show that neighborhood recombination frequencies (R-map, >280-kb segments) extracted from the MRFM, in conjunction with the validated formulas, better predicted the mapping population size than the genome-wide average recombination frequency (R-avg), with improved results whether the recombination frequency was calculated as genes/cM or kb/cM. Our results offer a detailed road map for better predicting mapping population size in diverse eukaryotes, but useful predictions will require robust recombination frequency maps based on sampling more progeny.