Ancestral haplotype reconstruction in endogamous populations using identity-by-descent

In this work we develop a novel algorithm for reconstructing the genomes of ancestral individuals, given genotype or sequence data from contemporary individuals and an extended pedigree of family relationships. A pedigree with complete genomes for every individual enables the study of allele frequency dynamics and haplotype diversity across generations, including deviations from neutrality such as transmission distortion. When studying heritable diseases, ancestral haplotypes can be used to augment genome-wide association studies and track disease inheritance patterns. The building blocks of our reconstruction algorithm are segments of Identity-By-Descent (IBD) shared between two or more genotyped individuals. The method alternates between identifying a source for each IBD segment and assembling IBD segments placed within each ancestral individual. Unlike previous approaches, our method is able to accommodate complex pedigree structures with hundreds of individuals genotyped at millions of SNPs.We apply our method to an Old Order Amish pedigree from Lancaster, Pennsylvania, whose founders came to North America from Europe during the early 18th century. The pedigree includes 1338 individuals from the past 12 generations, 394 with genotype data. The motivation for reconstruction is to understand the genetic basis of diseases segregating in the family through tracking haplotype transmission over time. Using our algorithm thread, we are able to reconstruct an average of 224 ancestral individuals per chromosome. For these ancestral individuals, on average we reconstruct 79% of their haplotypes. We also identify a region on chromosome 16 that is difficult to reconstruct—we find that this region harbors a short Amish-specific copy number variation and the gene HYDIN. thread was developed for endogamous populations, but can be applied to any extensive pedigree with the recent generations genotyped. We anticipate that this type of practical ancestral reconstruction will become more common and necessary to understand rare and complex heritable diseases in extended families.     

lra: A long read aligner for sequences and contigs

It is computationally challenging to detect variation by aligning single-molecule sequencing (SMS) reads, or contigs from SMS assemblies. One approach to efficiently align SMS reads is sparse dynamic programming (SDP), where optimal chains of exact matches are found between the sequence and the genome. While straightforward implementations of SDP penalize gaps with a cost that is a linear function of gap length, biological variation is more accurately represented when gap cost is a concave function of gap length. We have developed a method, lra, that uses SDP with a concave-cost gap penalty, and used lra to align long-read sequences from PacBio and Oxford Nanopore (ONT) instruments as well as de novo assembly contigs. This alignment approach increases sensitivity and specificity for SV discovery, particularly for variants above 1kb and when discovering variation from ONT reads, while having runtime that are comparable (1.05-3.76×) to current methods. When applied to calling variation from de novo assembly contigs, there is a 3.2% increase in Truvari F1 score compared to minimap2+htsbox. lra is available in bioconda (https://anaconda.org/bioconda/lra) and github (https://github.com/ChaissonLab/LRA).     

Correction to: Sex and race differences of cerebrospinal fluid metabolites in healthy individuals

Metabolomics - Metabolomics applications to the aquaculture research are increasing steadily. The use of standardized proton nuclear magnetic resonance (1H NMR) spectroscopy can provide the...     

Lichenochemical Screening and Antioxidant Capacity of Four Tunisian Lichen Species

The phenolic composition and antioxidant capacity of four Tunisian lichen species, Cladonia rangiformis, Flavoparmelia caperata, Squamarina cartilaginea and Xanthoria parietina, were determined in order to provide a better understanding of their lichenochemical composition. Powdered material of F. caperata was the richest in total phenolic content (956.68 μg GAE g⁻¹ DW) and S. cartilaginea in proanthocyanidin content (77.31 μg CE g⁻¹ DW), while the acetone extract of X. parietina showed the highest flavonoid content (9.56 μg CE g⁻¹ DW). The antioxidant capacity of all lichen extracts and crude material was evaluated by DPPH^. scavenging, iron-chelating, and iron-reducing powers. Results showed that methanol extracts of S. cartilaginea had the highest DPPH^. antioxidant capacity (IC₅₀=0.9 μg mL⁻¹) and the highest iron-reducing power was attributed to the acetone extract of this species. All extracts of all species were further screened by Fourier-transform infrared spectroscopy (FT-IR) and nuclear resonance spectroscopy (NMR); results showed an abundance of phenols, aromatic compounds, and fatty acids. Overall, our results showed that the investigated species are a rich source of potentially bioactive compounds with valuable properties.     

Specialized soil types affect host acceptability and performance of weed biocontrol candidates: implications for host specificity assessments

The Eurasian gall-forming weevil Ceutorhynchus cardariae Korotyeav (Coleoptera: Curculionidae) is a biological control candidate for the invasive Eurasian Lepidium draba L. (Brassicaceae) in the western USA. Among 157 nontarget plant species that have been tested, some North American Caulanthus and Streptanthus species, confamilial with Lepidium, were found to be at potential risk of attack by C. cardariae. Many Caulanthus and Streptanthus species grow on serpentine soils, which are characterized by low nutrient content and high concentrations of various combinations of heavy metals. Some of these species accumulate heavy metals, which have been shown to act as deterrents against insect herbivory. Standard pre-release host specificity tests with C. cardariae used plants propagated on horticultural soils, which could have inflated performance by C. cardariae on Caulanthus and Streptanthus species. To examine this possibility, we assessed the performance of C. cardariae on three Caulanthus species, the federally listed threatened and endangered Streptanthus glandulosus ssp. albidus, and Lepidium draba, on plants propagated in horticultural soil or in native serpentine soil. Our study showed that native serpentine soil influenced C. cardariae attack. All plant species, including L. draba, received less feeding damage and gall formation when grown in serpentine soil. In addition, feeding by C. cardariae was much less and fewer galls were formed on the confamilial species than on L. draba, regardless of soil type. Our data show that native confamilial species restricted to specialized soil types may be at less risk of herbivore attack than predicted based on tests conducted in horticultural soil.

     

Lake morphometry determines Dreissena invasion dynamics

Predicting the ecosystem effects of invasive species and the best control strategies requires understanding population dynamics and population regulation. Invasive bivalves zebra and quagga mussels (Dreissena spp.) are considered the most aggressive invaders in freshwaters and have become major drivers of ecosystem processes in the Laurentian Great Lakes. Combining all lake-wide studies of Dreissena spp. conducted in the Great Lakes, we found that invasion dynamics are largely governed by lake morphometry. Where both species are present, quagga mussels generally become dominant in 8–13 years. Thereafter, zebra mussels remain common in shallow lakes and embayments and lake-wide Dreissena density may remain similar, while in deep lakes quagga led to a near-complete displacement of zebra mussels and an ensuing dramatic increase in overall dreissenid density. In deep lakes, overall Dreissena biomass peaked later and achieved?~?threefold higher levels than in shallow lakes. Comparison with 21 waterbodies in North America and Europe colonized by both dreissenids confirmed that patterns of invasion dynamics found in the Great Lakes are very consistent with other waterbodies, and thus can be generalized to other lakes. Our biophysical model predicted that the long-term reduction in primary producers by mussel grazing may be fourfold less in deep compared to shallow lakes due to thermal stratification and a smaller proportion of the epilimnion in contact with the bottom. While this impact remains greatest in shallow areas, we show that when lakes are vertically well-mixed, dreissenid grazing impact may be greatest offshore, revealing a potentially strong offshore carbon and phosphorus sink.

     

Aquatic models of human ciliary diseases

Cilia are microtubule‐based structures that either transmit information into the cell or move fluid outside of the cell. There are many human diseases that arise from malfunctioning cilia. Although mammalian models provide vital insights into the underlying pathology of these diseases, aquatic organisms such as Xenopus and zebrafish provide valuable tools to help screen and dissect out the underlying causes of these diseases. In this review we focus on recent studies that identify or describe different types of human ciliopathies and outline how aquatic organisms have aided our understanding of these diseases.     

Punishment is strongly motivated by revenge and weakly motivated by inequity aversion

There are two broad functional explanations for second-party punishment: fitness-leveling and deterrence. The former suggests that people punish to reduce fitness differences, while the latter suggests that people punish in order to reciprocate losses and deter others from inflicting losses on them in the future. We explore the relative roles of these motivations using a pre-registered, two-player experiment with 2426 US participants from Amazon Mechanical Turk. Participants played as the “responder” and were assigned to either a Take or Augment condition. In the Take condition, the “partner” could steal money from the responder's bonus or do nothing. In the Augment condition, the partner could augment the responder's bonus by giving them money at no cost to themselves or do nothing. We also manipulated the responders' starting endowments, such that after the partner's decision, responders experienced different payoff outcomes: advantageous inequity, equality, or varying degrees of disadvantageous inequity. Responders then decided whether to pay a cost to punish the partner. Punishment was clearly influenced by theft and was most frequent when theft resulted in disadvantageous inequity. However, people also punished in the absence of theft, particularly when confronted with disadvantageous inequity. While the effect of inequity on punishment was small, our results suggest that punishment is motivated by more than just the desire to reciprocate losses. These findings highlight the multiple motivations undergirding punishment and bear directly on functional explanations for the existence of punishment in human societies.     

Diversity of dimethylsulfide-degrading methanogens and sulfate-reducing bacteria in anoxic sediments along the Medway Estuary,UK

Methane is a powerful greenhouse gas but the microbial diversity mediating methylotrophic methanogenesis is not well-characterized. One overlooked route to methane is via the degradation of dimethylsulfide (DMS), an abundant organosulfur compound in the environment. Methanogens and sulfate-reducing bacteria (SRB) can degrade DMS in anoxic sediments depending on sulfate availability. However, we know little about the underlying microbial community and how sulfate availability affects DMS degradation in anoxic sediments. We studied DMS-dependent methane production along the salinity gradient of the Medway Estuary (UK) and characterized, for the first time, the DMS-degrading methanogens and SRB using cultivation-independent tools. DMS metabolism resulted in high methane yield (39%–42% of the theoretical methane yield) in anoxic sediments regardless of their sulfate content. Methanomethylovorans, Methanolobus and Methanococcoides were dominant methanogens in freshwater, brackish and marine incubations respectively, suggesting niche-partitioning of the methanogens likely driven by DMS amendment and sulfate concentrations. Adding DMS also led to significant changes in SRB composition and abundance in the sediments. Increases in the abundance of Sulfurimonas and SRB suggest cryptic sulfur cycling coupled to DMS degradation. Our study highlights a potentially important pathway to methane production in sediments with contrasting sulfate content and sheds light on the diversity of DMS degraders.     

Natural variation in colony inbreeding does not influence susceptibility to a fungal pathogen in a termite

Reduced genetic diversity through inbreeding can negatively affect pathogen resistance. This relationship becomes more complicated in social species, such as social insects, since the chance of disease transmission increases with the frequency of interactions among individuals. However, social insects may benefit from social immunity, whereby individual physiological defenses may be bolstered by collective‐level immune responses, such as grooming or sharing of antimicrobial substance through trophallaxis. We set out to determine whether differences in genetic diversity between colonies of the subterranean termite, Reticulitermes flavipes, accounts for colony survival against pathogens. We sampled colonies throughout the United States (Texas, North Carolina, Maryland, and Massachusetts) and determined the level of inbreeding of each colony. To assess whether genetically diverse colonies were better able to survive exposure to diverse pathogens, we challenged groups of termite workers with two strains of a pathogenic fungus, one local strain present in the soil surrounding sampled colonies and another naïve strain, collected outside the range of this species. We found natural variation in the level of inbreeding between colonies, but this variation did not explain differences in susceptibility to either pathogen. Although the naïve strain was found to be more hazardous than the local strain, colony resistance was correlated between two strains, meaning that colonies had either relatively high or low susceptibility to both strains regardless of their inbreeding coefficient. Overall, our findings may reflect differential virulence between the strains, immune priming of the colonies via prior exposure to the local strain, or a coevolved resistance toward this strain. They also suggest that colony survival may rely more upon additional factors, such as different behavioral response thresholds or the influence of a specific genetic background, rather than the overall genetic diversity of the colony.