Reconstructing sibling relationships in wild populations

Reconstruction of sibling relationships from genetic data is an important component of many biological applications. In particular, the growing application of molecular markers (microsatellites) to study wild populations of plant and animals has created the need for new computational methods of establishing pedigree relationships, such as sibgroups, among individuals in these populations. Most current methods for sibship reconstruction from microsatellite data use statistical and heuristic techniques that rely on a priori knowledge about various parameter distributions. Moreover, these methods are designed for data with large number of sampled loci and small family groups, both of which typically do not hold for wild populations. We present a deterministic technique that parsimoniously reconstructs sibling groups using only Mendelian laws of inheritance. We validate our approach using both simulated and real biological data and compare it to other methods. Our method is highly accurate on real data and compares favorably with other methods on simulated data with few loci and large family groups. It is the only method that does not rely on a priori knowledge about the population under study. Thus, our method is particularly appropriate for reconstructing sibling groups in wild populations.     

A new version of PRT software for sibling groups reconstruction with comments regarding several issues in the sibling reconstruction problem

Pedigree reconstruction using genotypic markers has become an important tool for the study of natural populations. The nonstandard nature of the underlying statistical problems has led to the necessity of developing specialized statistical and computational methods. In this article, a new version of pedigree reconstruction tools (PRT 2.0) is presented. The software implements algorithms proposed in Almudevar & Field (Journal of Agricultural Biological and Environmental Statistics, 4, 1999, 136) and Almudevar (Biometrics, 57, 2001a, 757) for the reconstruction of single generation sibling groups (SG). A wider range of enumeration algorithms is included, permitting improved computational performance. In particular, an iterative version of the algorithm designed for larger samples is included in a fully automated form. The new version also includes expanded simulation utilities, as well as extensive reporting, including half-sibling compatibility, parental genotype estimates and flagging of potential genotype errors. A number of alternative algorithms are described and demonstrated. A comparative discussion of the underlying methodologies is presented. Although important aspects of this problem remain open, we argue that a number of methodologies including maximum likelihood estimation (COLONY 1.2 and 2.0) and the set cover formulation (KINALYZER) exhibit undesirable properties in the sibling reconstruction problem. There is considerable evidence that large sets of individuals not genetically excluded as siblings can be inferred to be a true sibling group, but it is also true that unrelated individuals may be genetically compatible with a true sibling group by chance. Such individuals may be identified on a statistical basis. PRT 2.0, based on these sound statistical principles, is able to efficiently match or exceed the highest reported accuracy rates, particularly for larger SG. The new version is available at http://www.urmc.rochester.edu/biostat/people/faculty/almudevar.cfm.     

Boosting Probabilistic Graphical Model Inference by Incorporating Prior Knowledge from Multiple Sources

Inferring regulatory networks from experimental data via probabilistic graphical models is a popular framework to gain insights into biological systems. However, the inherent noise in experimental data coupled with a limited sample size reduces the performance of network reverse engineering. Prior knowledge from existing sources of biological information can address this low signal to noise problem by biasing the network inference towards biologically plausible network structures. Although integrating various sources of information is desirable, their heterogeneous nature makes this task challenging. We propose two computational methods to incorporate various information sources into a probabilistic consensus structure prior to be used in graphical model inference. Our first model, called Latent Factor Model (LFM), assumes a high degree of correlation among external information sources and reconstructs a hidden variable as a common source in a Bayesian manner. The second model, a Noisy-OR, picks up the strongest support for an interaction among information sources in a probabilistic fashion. Our extensive computational studies on KEGG signaling pathways as well as on gene expression data from breast cancer and yeast heat shock response reveal that both approaches can significantly enhance the reconstruction accuracy of Bayesian Networks compared to other competing methods as well as to the situation without any prior. Our framework allows for using diverse information sources, like pathway databases, GO terms and protein domain data, etc. and is flexible enough to integrate new sources, if available.     

Limited usefulness of microsatellite markers from the malaria vector Anopheles gambiae when applied to the closely related species Anopheles melas

Anopheles melas is a brackish water mosquito found in coastal West Africa where it is a dominant malaria vector locally. In order to facilitate genetic studies of this species, 45 microsatellite loci originally developed for Anopheles gambiae were sequenced in An. melas. Those that were suitable based on repeat number and flanking regions were examined in 2 natural populations from Equatorial Guinea. Only 15 loci were eventually deemed suitable as polymorphic markers in An. melas populations. These loci were screened in 4 populations from a wider geographic range. Heterozygosity estimates ranged from 0.18 to 0.79, and 2.5-15 average alleles were observed per locus, yielding 13 highly polymorphic markers and 2 loci with lower variability. To examine the usefulness of microsatellite markers when applied in a sibling species, the original An. gambiae specific markers were used to amplify 5 loci in An. melas. Null alleles were found for 1 An. gambiae marker. We discuss the pitfalls of using microsatellite loci across closely related species and conclude that in addition to the problem of null alleles associated with this practice, many loci may prove to be of very limited use as polymorphic markers even when used in a sibling species.     

Bayesian sparse‐based reconstruction in bioluminescence tomography improves localization accuracy and reduces computational time

Bioluminescence tomography (BLT) provides fundamental insight into biological processes in vivo. To fully realize its potential, it is important to develop image reconstruction algorithms that accurately visualize and quantify the bioluminescence signals taking advantage of limited boundary measurements. In this study, a new 2‐step reconstruction method for BLT is developed by taking advantage of the sparse a priori information of the light emission using multispectral measurements. The first step infers a wavelength‐dependent prior by using all multi‐wavelength measurements. The second step reconstructs the source distribution based on this developed prior. Simulation, phantom and in vivo results were performed to assess and compare the accuracy and the computational efficiency of this algorithm with conventional sparsity‐promoting BLT reconstruction algorithms, and results indicate that the position errors are reduced from a few millimeters down to submillimeter, and reconstruction time is reduced by 3 orders of magnitude in most cases, to just under a few seconds. The recovery of single objects and multiple (2 and 3) small objects is simulated, and the recovery of images of a mouse phantom and an experimental animal with an existing luminescent source in the abdomen is demonstrated. Matlab code is available at https://github.com/jinchaofeng/code/tree/master .      

Molecular genetic insights on cheetah (Acinonyx jubatus) ecology and conservation in Namibia

The extent and geographic patterns of molecular genetic diversity of the largest remaining free-ranging cheetah population were described in a survey of 313 individuals from throughout Namibia. Levels of relatedness, including paternity/maternity (parentage), were assessed across all individuals using 19 polymorphic microsatellite loci, and unrelated cheetahs (n = 89) from 7 regions were genotyped at 38 loci to document broad geographical patterns. There was limited differentiation among regions, evidence that this is a generally panmictic population. Measures of genetic variation were similar among all regions and were comparable with Eastern African cheetah populations. Parentage analyses confirmed several observations based on field studies, including 21 of 23 previously hypothesized family groups, 40 probable parent/offspring pairs, and 8 sibling groups. These results also verified the successful integration and reproduction of several cheetahs following natural dispersal or translocation. Animals within social groups (family groups, male coalitions, or sibling groups) were generally related. Within the main study area, radio-collared female cheetahs were more closely interrelated than similarly compared males, a pattern consistent with greater male dispersal. The long-term maintenance of current patterns of genetic variation in Namibia depends on retaining habitat characteristics that promote natural dispersal and gene flow of cheetahs.     

Seascape and life-history traits do not predict self-recruitment in a coral reef fish

The persistence and resilience of many coral reef species are dependent on rates of connectivity among sub-populations. However, despite increasing research efforts, the spatial scale of larval dispersal remains unpredictable for most marine metapopulations. Here, we assess patterns of larval dispersal in the angelfish Centropyge bicolor in Kimbe Bay, Papua New Guinea, using parentage and sibling reconstruction analyses based on 23 microsatellite DNA loci. We found that, contrary to previous findings in this system, self-recruitment (SR) was virtually absent at both the reef (0.4–0.5% at 0.15 km²) and the lagoon scale (0.6–0.8% at approx. 700 km²). While approximately 25% of the collected juveniles were identified as potential siblings, the majority of sibling pairs were sampled from separate reefs. Integrating our findings with earlier research from the same system suggests that geographical setting and life-history traits alone are not suitable predictors of SR and that high levels of localized recruitment are not universal in coral reef fishes.     

Genetic structure and heterozygosity in the kuru region,Eastern Highlands of New Guinea

Genetic structure of 20 of the 37 linguistic groups in the Eastern Highlands of New Guinea including the kuru region is analyzed using information on blood groups and serum protein polymorphisms. The average individual is heterozygous at 28.6% of loci and the average number of alleles per locus is 1.234. Coefficients of kinship for linguistic groups range from 0.005 for the sweet potato cultivating North Fore to 0.075 for the isolated Pawaians whose dietary staple is sago and who depend more on hunting and gathering. As one selects linguistic groups with smaller and smaller population size and increasing isolation, one finds that kinship coefficients rise as much as tenfold, but there is no concomitant loss of heterozygosity or trend toward fixation of alleles. Genetic relationships established by genetic distance trees and by principal components analysis are comparable and are consistent with other anthropological observations.     

Combining Harvest and Genetics to Estimate Reproduction in Wolves

Parameters of reproductive success are important to the management of wildlife populations. Genetic monitoring can be an effective approach for acquiring this important demographic information when traditional methods are unsuccessful, inefficient, or too expensive. This study demonstrates a novel application of genetic data opportunistically collected from harvested game to estimate a minimum annual count of breeding packs of gray wolves (Canis lupus) and to provide a coarse index of harvest vulnerability of young of the year (YOY) across packs. We used 18 microsatellite loci to genotype 98 gray wolf YOY from 2014 and 105 from 2015 harvested in Idaho, USA. Using this genotype data, we reconstructed sibling groups for each cohort using the program COLONY and treated full-sibling litters as proxies for unique packs. In addition to evaluating our marker panel using simulations, we assessed the accuracy of empirical relationship assignments by adding YOY of known relationship from long-term study packs to the dataset (27 individuals from 2014 and 61 from 2015) and tracking correctly reconstructed relationships. We varied COLONY input parameters to evaluate the power of relationship assignments under conditions that may be encountered when working with empirical data. We also compared COLONY's estimates of effective number of breeders based on sibship frequency to estimates based on a commonly used linkage-disequilibrium method. All COLONY runs for both cohorts correctly identified the known sibling relationships. Among the other individuals, changes in the geographic clustering of putative siblings, probabilities of inclusion and exclusion for reconstructed sibling groups, and consistency of relationship assignments across COLONY runs suggested that marker number had a larger effect on accuracy than access to population-level genetic data. Our estimates of breeding packs subjected to harvest within the state (52 for 2014 and 63 for 2015) differed from estimates reported by Idaho Department of Fish and Game by ≤6 for both years. Among packs that had pups harvested, most packs had 1–2 YOY harvested, whereas other packs had as many as 5 YOY harvested. All estimates of the number of effective breeders had overlapping confidence intervals regardless of method, though sibship frequency-based estimates had larger confidence intervals than estimates using the linkage disequilibrium method. Our study shows that sibling relationships can be accurately and reliably reconstructed from harvested gray wolves, and demonstrates a valuable new use of samples collected through harvest. © 2020 The Wildlife Society.     

Multipoint analysis of human quantitative genetic variation.

A unique method of partitioning human quantitative genetic variation into effects due to specific chromosomal regions is presented. This method is based on estimating the proportion of genetic material, R, shared identical by descent (IBD) by sibling pairs in a specified chromosomal region, on the basis of their marker genotypes at a set of marker loci spanning the region. The mean and variance of the distribution of R conditional on IBD status and recombination pattern between two marker loci are derived as a function of the distance between the two loci. The distribution of the estimates of R is exemplified using data on 22 loci on chromosome 7. A method of using the estimated R values and observed values of a quantitative trait in a set of sibships to estimate the proportion of total genetic variance explained by loci in the region of interest is presented. Monte Carlo simulation techniques are used to show that this method is more powerful than existing methods of quantitative linkage analysis based on sib pairs. It is also shown through simulation studies that the proposed method is sensitive to genetic variation arising from both a single locus of large effect as well as from several loosely linked loci of moderate phenotypic effect.     

Indicative value ofPediastrum and other coccal green algae in palaeoecology

Sporopollenin layers in the cell wall of coccal green algae are responsible for the resistance of cell walls to destructive processes during fossilization as well as during chemical preparation of samples for pollen-analysis. Pollen slides of samples from limnic sediments thus also contain some algal cell walls. Although some pollen-analysts tried to stress this fact, the finds of algae in pollen slides have not been paid systematic attention yet, despite their potential use for a more accurate palaeoecological reconstruction. The article summarizes the results of palaeoecological studies showing how the algae can be used in palaeoecological reconstruction of past environments. The possibility of utilizing the indicative value of algal finds is demonstrated on examples of algal communities from fossil, subrecent and recent sediments from different longitudes, latitudes, and altitudes. The identification and indicative values of species and varieties ofPediastrum are included in a special review (Komárek &; Jankovská, Biblioth. Phycol., in press). The contemporary knowledge of ecological requirements of the given taxa, completed by information from their fossil finds, makes possible the reconstruction of trophic and temperature conditions and of the purity of the water environment of the past water biotopes.     

Natural skip oviposition of the mosquito Aedes aegypti indicated by codominant genetic markers

This study examines the use of codominant restriction fragment length polymorphism (RFLP) markers to estimate the number of sibling families found within and among oviposition sites used by the mosquito Aedes aegypti (L) (Diptera: Culicidae). Estimates were made using pairwise relatedness (rxy) calculations based on alleles shared between individuals. Genotypes for eight laboratory mosquito families were determined at six RFLP loci and the observed allele frequencies were used to generate simulated distributions of rxy from full-sibling and unrelated pairs of individuals. The midpoint (mp) between the means of the pairwise rxy distributions was used to discriminate full-sibling families from unrelated families. Clusters of individuals with rxy values higher than the mp value were grouped as putative sibling families. This method was tested by calculating actual rxy for all pairwise comparisons of the known laboratory full-sibling and paternal half-sibling families, followed by upgma cluster analysis to group sibling families. The technique was then used for sibling estimations on wild caught mosquitoes collected at three locations in Trinidad, West Indies. From field populations, 35 families were estimated among 122 individuals tested with an average of 6.2 families per container. Members of 19 predicted families clustered as groups across multiple containers, providing molecular evidence for skip-oviposition behaviour in Ae. aegypti females, whereby individual females oviposit in more than one container.     

gerud 2.0: a computer program for the reconstruction of parental genotypes from half‐sib progeny arrays with known or unknown parents

Large groups of related progeny, as can be collected from discrete egg masses or fruits, present excellent opportunities for parentage analysis by allowing the reconstruction of parental genotypes. Current techniques of parental genotypic reconstruction require the knowledge of at least one parental genotype. Here, I present a new computer program that reconstructs parental genotypes when no parents are known, provided that the progeny array contains only full and half siblings. The prospects for successfully reconstructing parental genotypes for such progeny arrays with unknown parents are nearly as good as those for the case when one parent is known with certainty.     

Reconstructing species phylogeny of the carabid beetles Ohomopterus using multiple nuclear DNA sequences: heterogeneous information content and the performance of simultaneous analyses

We attempted a phylogenetic reconstruction for the carabid subgenus Ohomopterus (genus Carabus), a notable case of radiation with mitochondrial introgression across species. Sequence data from five nuclear single copy loci were used, including wingless (Wg), phosphoenolpyruvate carboxykinase (PepCK), cytochrome c (Cytc), elongation factor-1alpha (EF-1alpha), and an anonymous single copy locus (Carab1). Sequences of Cytc, EF-1alpha, and Carab1 included intron or intron-like parts with length variation. The analysis of individual loci resulted in low resolution of the phylogenetic relationships, and the monophyly of several morphologically recognized species for which multiple specimens were analyzed was not revealed. Several specimens were heterozygous, with non-monophyletic alleles observed in three of the five loci at which alleles in heterozygotes were separated. In a simultaneous analysis of the five loci with ambiguously aligned parts eliminated and heterozygotic sites treated as missing, the resulting tree was well resolved, but the branch support was generally weak because of conflicting phylogenetic signals from different loci. We also attempted to incorporate allelic sequence data plus the ambiguously aligned parts in the analysis, by using all possible combinations of alleles from different loci in heterozygotic individuals, but the resultant tree was not supported more strongly. Nonetheless, these simultaneous analyses provided support for the monophyly of several species and species groups, and revealed the basic evolutionary trend of OHOMOPTERUS: initial widespread groups with simpler genitalia and the origination of exaggerated genitalia in a derived clade. This study exemplifies problems inherent in the phylogenetic reconstruction of closely related organisms where low levels of variation limit the information content from each locus, while heterozygosity, different phylogenetic history of multiple loci, and alignment ambiguity further hamper phylogenetic reconstruction unless several loci converge on a uniform signal.     

Accurate partition of individuals into full-sib families from genetic data without parental information.

Two Markov chain Monte Carlo algorithms are proposed that allow the partitioning of individuals into full-sib groups using single-locus genetic marker data when no parental information is available. These algorithms present a method of moving through the sibship configuration space and locating the configuration that maximizes an overall score on the basis of pairwise likelihood ratios of being full-sib or unrelated or maximizes the full joint likelihood of the proposed family structure. Using these methods, up to 757 out of 759 Atlantic salmon were correctly classified into 12 full-sib families of unequal size using four microsatellite markers. Large-scale simulations were performed to assess the sensitivity of the procedures to the number of loci and number of alleles per locus, the allelic distribution type, the distribution of families, and the independent knowledge of population allelic frequencies. The number of loci and the number of alleles per locus had the most impact on accuracy. Very good accuracy can be obtained with as few as four loci when they have at least eight alleles. Accuracy decreases when using allelic frequencies estimated in small target samples with skewed family distributions with the pairwise likelihood approach. We present an iterative approach that partly corrects that problem. The full likelihood approach is less sensitive to the precision of allelic frequencies estimates but did not perform as well with the large data set or when little information was available (e.g., four loci with four alleles).     

A Novel Gaussian Extrapolation Approach for 2D Gel Electrophoresis Saturated Protein Spots

Analysis of images obtained from two-dimensional gel electrophoresis (2D-GE) is a topic of utmost importance in bioinformatics research, since commercial and academic software available currently has proven to be neither completely effective nor fully automatic, often requiring manual revision and refinement of computer generated matches. In this work, we present an effective technique for the detection and the reconstruction of over-saturated protein spots. Firstly, the algorithm reveals overexposed areas, where spots may be truncated, and plateau regions caused by smeared and overlapping spots. Next, it reconstructs the correct distribution of pixel values in these overexposed areas and plateau regions, using a two-dimensional least-squares fitting based on a generalized Gaussian distribution. Pixel correction in saturated and smeared spots allows more accurate quantification, providing more reliable image analysis results. The method is validated for processing highly exposed 2D-GE images, comparing reconstructed spots with the corresponding non-saturated image, demonstrating that the algorithm enables correct spot quantification.     

Assignment of paternity groups without access to parental genotypes: multiple mating and developmental plasticity in squid

We present a novel approach to investigating sibling relationships and reconstructing parental genotypes from a progeny array. The Bayesian method we have employed is flexible and may be applicable to a variety of situations in addition to the one presented here. While mutation rates and breeding population allele frequencies can be taken into account, the model requires relatively few loci and makes few assumptions. Paternity of 270 veined squid (Loligo forbesi) hatchlings from three egg strings collected from one location was assigned using five microsatellite loci. Paternal and maternal genotypes reconstructed for each of the three strings were identical, strongly indicating that a single female produced the strings that were fertilized by the same four males. The proportion of eggs fertilized was not equal between males in all three strings, with male 1 siring most offspring (up to 68% in string 1), through to male 4 siring the least (as low as 2.4% in string 1). Although temperature had a profound effect on incubation time, paternity did not affect this trait at 12 degrees C or 8 degrees C.     

Three-dimensional reconstruction of the actin cytoskeleton from stereo images

In eucaryotic cells, actin filaments are abundant components in the cytoskeleton where they form a complex three dimensional (3D) structural network that provides the cell with its shape and mechanical properties. However, understanding the structural and mechanical properties of actin filaments composing the cell cytoskeleton is often hampered by the inability to faithfully reconstruct the three-dimensional geometric relationships. This paper presents a vision-based reconstruction approach that automatically reconstitutes the three-dimensional structures of cytoskeletal polymers from stereo image pairs taken at the different tilt angles. The approach finds corresponding points between two images and recovers the depth information about the structures. The computational process consists of three major procedures: feature representation, stereo matching, and disparity refinement, implemented in a multi-resolution manner based on a coarse-to-fine strategy. The reconstruction depicts the three-dimensional structure of cytoskeletal polymers and their geometric relationships. New and useful information becomes available and allows quantitative analysis of the structure. Measurement of the cytoskeleton geometrical properties and the filament concentration in a defined volume are obtained by direct calculation.     

Disambiguating Multi–Modal Scene Representations Using Perceptual Grouping Constraints

In its early stages, the visual system suffers from a lot of ambiguity and noise that severely limits the performance of early vision algorithms. This article presents feedback mechanisms between early visual processes, such as perceptual grouping, stereopsis and depth reconstruction, that allow the system to reduce this ambiguity and improve early representation of visual information. In the first part, the article proposes a local perceptual grouping algorithm that — in addition to commonly used geometric information — makes use of a novel multi–modal measure between local edge/line features. The grouping information is then used to: 1) disambiguate stereopsis by enforcing that stereo matches preserve groups; and 2) correct the reconstruction error due to the image pixel sampling using a linear interpolation over the groups. The integration of mutual feedback between early vision processes is shown to reduce considerably ambiguity and noise without the need for global constraints.