Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

Parentage determination of an isolated Yangtze finless porpoise population Neophocaena phocaenoides asiaeorientalis in the Yangtze Tian-e-Zhou Baiji National Natural Reserve based on molecular data

Reproductive behaviors are poorly known for the Yangtze finless porpoise Neophocaena phocaenoides asiaeorientalis. In this study, the parentage of an isolated Yangtze finless porpoise population inhabiting the Yangtze Tian-e-Zhou Baiji National Natural Reserve is determined by analysis of microsatellite loci and mitochondrial DNA (mtDNA) control region sequences, and the porpoise's reproductive behaviors are studied. Overall 4 full parentage assignments and additional 3 single parentage assignments were determined for the population of 23 individuals. The analysis shows that their estimated reproductive cycle is shorter than that reported previously and there probably exists an overlapping between gestation and lactation period. The Study also shows that the female does not show fidelity to a particular male for breeding and vice versa, the oldest males did not monopolize mating and the dominance rank could not be so strict for the porpoise society. Moreover, the porpoise's mating pattern and relatedness among candidate parents are discussed here. These results provide important information for making guidelines of management and conservation for this protected population.

Genetic assessment of parentage in the caridean rock shrimp Rhynchocinetes typus based on microsatellite markers

Over the past decade, the common rock shrimp, Rhynchocinetes typus H. Milne Edwards, 1837, has been the focus of extensive investigations on mating behaviour. The species is now perceived as a model system for the study of reproductive strategies and sexual conflict in crustaceans displaying external fertilization. Using molecular markers, the current study assesses whether social mating behaviour in common rock shrimp translates into true genetic parentage. In a large mesocosm tank with >200 individuals of both sexes, the analysis of 15 families (22 eggs per female) for three informative microsatellites unambiguously confirmed multiple paternity in 11 instances (73%) involving, in each case, two to four males. Where more than one male was identified siring a particular brood, reproductive skew was apparent towards a single individual. Results suggest that multiple paternity in this species results from subordinate male coercive behaviour, female solicitation of multiple male matings or a combination of both.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Analysis of disputed single-parent/child and sibling relationships using 16 STR loci

This study describes the validation of short tandem repeat (STR) systems for the resolution of cases of disputed parentage where only a single parent is available for testing or where the claimed relationship of both parents is in doubt and also cases where sibship must be tested. Three separate multiplex systems the Second Generation Multiplex, Powerplex 1.2 and FFFL have been employed, giving a total of 16 STR loci. Both empirical and theoretical approaches to the validation have been adopted. Appropriate equations have been derived to calculate likelihood ratios for different relationships, incorporating a correction for subpopulation effects. An F(ST) point estimate of 1% has been applied throughout. Empirically, 101 cases of alleged father, alleged mother and child where analysed using six SLP systems and also using the three multiplex STR systems. Of the 202 relationships tested, 197 were independently resolved by both systems, providing either clear evidence of non-parentage or strong support for the relationship.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level

Abstract Background: The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North+East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe. Results: A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST=0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST=0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars. Conclusions: The variation found at group and sub-group levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

Rapid and accurate G M1-gangliosidosis diagnosis using a parentage testing microsatellite

The G M1-gangliosidosis is an autosomal recessive lysosomal storage disease caused by structural defects of the beta-galactosidase gene (GLB1) which lead to a severe phenotypical impairment in homozygous individuals, whereas heterozygous carriers remain clinically normal. Currently employed DNA parentage tests include the analysis of microsatellites, which also have a diagnostic predictive value. The aim of this study was to provide a reliable tool for genotyping the canine GLB1 which can be effectively integrated in parentage testing investigations. For this purpose the association between the GLB1 gene and the AHT K253 microsatellite was analyzed in 30 Alaskan huskies (11 GLB1+/+, 17 GLB1+/- and 2 GLB1-/- dogs). The 143 bp AHT K253 microsatellite allele was identified only in GLB1+/- and GLB1-/- animals and was in strong linkage disequilibrium with the causative mutation for G M1-gangliosidosis, a 19 bp duplication within exon 15 of the GLB1 gene. The results of the present study revealed a 100% concordance between the previous established genotypes and those obtained after the analysis of the AHT K253 microsatellite. Thus, the genotype of the AHT K253 microsatellite, which is routinely determined during dog parentage testing, has a high predictive value for the G M1-gangliosidosis carrier status.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Expressed sequence tag analysis of genes expressed during development of the tropical abalone Haliotis asinina

The tropical abalone. Haliotis asinina. is,in ideal species to investigate the molecular mechanisms that control development. growth, reproduction and shell formation in all cultured haliotids. Here we describe the analysis of 232 expressed sequence tags (EST) obtained front a developmental H. asinina cDNA library intended for future microarray studies. From this data set we identified 183 unique gene Clusters. Of these, 90 clusters showed significant homology with sequences lodged in GenBank, ranging in function from general housekeeping to signal transduction, gene regulation and cell-cell communication. Seventy-one clusters possessed completely novel ORFs greater than 50 codons in length, highlighting the paucity of sequence data from molluscs and other lophotrochozoans. This study of developmental gene expression in H. asinina provides the foundation for further detailed analyses of abalone growth, development and reproduction.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

Determination of preventive maintenance lead time using hybrid analysis

The time for conducting Preventive Maintenance (PM) on an asset is often determined using a predefined alarm limit based on trends of a hazard function. In this paper, the authors propose using both hazard and reliability functions to improve the accuracy of the prediction particularly when the failure characteristic of the asset whole life is modelled using different failure distributions for the different stages of the life of the asset. The proposed method is validated using simulations and case studies.