MEGASAT: automated inference of microsatellite genotypes from sequence data

MEGASAT is software that enables genotyping of microsatellite loci using next-generation sequencing data. Microsatellites are amplified in large multiplexes, and then sequenced in pooled amplicons. MEGASAT reads sequence files and automatically scores microsatellite genotypes. It uses fuzzy matches to allow for sequencing errors and applies decision rules to account for amplification artefacts, including nontarget amplification products, replication slippage during PCR (amplification stutter) and differential amplification of alleles. An important fea- ture of MEGASAT is the generation of histograms of the length–frequency distributions of amplification products for each locus and each individual. These histograms, analogous to electropherograms traditionally used to score microsatellite genotypes, enable rapid evaluation and editing of automatically scored genotypes. MEGASAT is written in Perl, runs on Windows, Mac OS X and Linux systems, and includes a simple graphical user interface. We demon- strate MEGASAT using data from guppy, Poecilia reticulata. We genotype 1024 guppies at 43 microsatellites per run on an Illumina MiSeq sequencer. We evaluated the accuracy of automatically called genotypes using two methods, based on pedigree and repeat genotyping data, and obtained estimates of mean genotyping error rates of 0.021 and 0.012. In both estimates, three loci accounted for a disproportionate fraction of genotyping errors; conversely, 26 loci were scored with 0–1 detected error (error rate ≤0.007). Our results show that with appropriate selection of loci, automated genotyping of microsatellite loci can be achieved with very high throughput, low genotyping error and very low genotyping costs

Genomics reveal population structure, evolutionary history, and signatures of selection in the northern bottlenose whale, Hyperoodon ampullatus

Funding: This work was supported by Fisheries and Oceans Canada (DFO) Maritimes and National Geographic emerging explorer grant to L.J.F, with support by and Natural Sciences and Engineering Research Council of Canada (NSERC) and Killam Nova Scotia Doctoral Scholarships. Work was also supported by US Office of Naval Research and US Strategic Environmental Research and Development Program (SERDP), DFO, University of Windsor, Crown-Indigenous Relations and Northern Affairs Canada, Nunavut Fisheries Association, Government of Nunavut, and NSERC. Funding and resources for sequencing the northern bottlenose whale genome was supported by the CanSeq150 program of Canada’s Genomics Enterprise.Information on wildlife population structure, demographic history, and adaptations are fundamental to understanding species evolution and informing conservation strategies. To study this ecological context for a cetacean of conservation concern, we conducted the first genomic assessment of the northern bottlenose whale, Hyperoodon ampullatus, using whole-genome resequencing data (n = 37) from five regions across the North Atlantic Ocean. We found a range-wide pattern of isolation-by-distance with a genetic subdivision distinguishing three subgroups: the Scotian Shelf, western North Atlantic, and Jan Mayen regions. Signals of elevated levels of inbreeding in the Endangered Scotian Shelf population indicate this population may be more vulnerable than the other two subgroups. In addition to signatures of inbreeding, evidence of local adaptation in the Scotian Shelf was detected across the genome. We found a long-term decline in effective population size for the species, which poses risks to their genetic diversity and may be exacerbated by the isolating effects of population subdivision. Protecting important habitat and migratory corridors should be prioritized to rebuild population sizes that were diminished by commercial whaling, strengthen gene flow, and ensure animals can move across regions in response to environmental changes.Publisher PDFPeer reviewe

Genome-wide changes in genetic diversity in a population of Myotis lucifugus affected by white-nose syndrome

Novel pathogens can cause massive declines in populations, and even extirpation of hosts. But disease can also act as a selective pressure on survivors, driving the evolution of resistance or tolerance. Bat white-nose syndrome (WNS) is a rapidly spreading wildlife disease in North America. The fungus causing the disease invades skin tissues of hibernating bats, resulting in disruption of hibernation behavior, premature energy depletion, and subsequent death. We used whole-genome sequencing to investigate changes in allele frequencies within a population of Myotis lucifugus in eastern North America to search for genetic resistance to WNS. Our results show low F-ST values within the population across time, i.e., prior to WNS (Pre-WNS) compared to the population that has survived WNS (Post-WNS). However, when dividing the population with a geographical cut-off between the states of Pennsylvania and New York, a sharp increase in values on scaffold GL429776 is evident in the Post-WNS samples. Genes present in the diverged area are associated with thermoregulation and promotion of brown fat production. Thus, although WNS may not have subjected the entire M. lucifugus population to selective pressure, it may have selected for specific alleles in Pennsylvania through decreased gene flow within the population. However, the persistence of remnant sub-populations in the aftermath of WNS is likely due to multiple factors in bat life history.Peer reviewe

Fragmentation and shear band formation by slow compression of brittle porous media

Localized fragmentation is an important phenomenon associated with the formation of shear bands and faults in granular media. It can be studied by empirical observation, by laboratory experiment, or by numerical simulation. Here we investigate the spatial structure and statistics of fragmentation using discrete element simulations of the strain-controlled uniaxial compression of cylindrical samples of different finite size. As the system approaches failure, damage localizes in a narrow shear band or synthetic fault "gouge" containing a large number of poorly sorted noncohesive fragments on a broad bandwidth of scales, with properties similar to those of natural and experimental faults. We determine the position and orientation of the central fault plane, the width of the shear band, and the spatial and mass distribution of fragments. The relative width of the shear band decreases as a power law of the system size, and the probability distribution of the angle of the central fault plane converges to around 30 degrees, representing an internal coefficient of friction of 0.7 or so. The mass of fragments is power law distributed, with an exponent that does not depend on scale, and is near that inferred for experimental and natural fault gouges. The fragments are in general angular, with a clear self-affine geometry. The consistency of this model with experimental and field results confirms the critical roles of preexisting heterogeneity, elastic interactions, and finite system size to grain size ratio on the development of shear bands and faults in porous media. © 2016 American Physical Society

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish

Funding Information Natural Sciences and Engineering Research Council of Canada (NSERC) Strategic Project Atlantic Canada Opportunities Agency and Department of Tourism, Culture, Industry and Innovation grants allocated to the Labrador Institute (MC) Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Genomics Research and Development Initiative (GRDI) Weston Family AwardPeer reviewedPublisher PD

The role of biophysical cohesion on subaqueous bed form size

Biologically active, fine-grained sediment forms abundant sedimentary deposits on Earth's surface, and mixed mud-sand dominates many coasts, deltas, and estuaries. Our predictions of sediment transport and bed roughness in these environments presently rely on empirically based bed form predictors that are based exclusively on biologically inactive cohesionless silt, sand, and gravel. This approach underpins many paleoenvironmental reconstructions of sedimentary successions, which rely on analysis of cross-stratification and bounding surfaces produced by migrating bed forms. Here we present controlled laboratory experiments that identify and quantify the influence of physical and biological cohesion on equilibrium bed form morphology. The results show the profound influence of biological cohesion on bed form size and identify how cohesive bonding mechanisms in different sediment mixtures govern the relationships. The findings highlight that existing bed form predictors require reformulation for combined biophysical cohesive effects in order to improve morphodynamic model predictions and to enhance the interpretations of these environments in the geological record

The Alberta Heart Failure Etiology and Analysis Research Team (HEART) study

Background Nationally, symptomatic heart failure affects 1.5-2% of Canadians, incurs $3 billion in hospital costs annually and the global burden is expected to double in the next 1–2 decades. The current one-year mortality rate after diagnosis of heart failure remains high at >25%. Consequently, new therapeutic strategies need to be developed for this debilitating condition. Methods/Design The objective of the Alberta HEART program (http://albertaheartresearch.ca) is to develop novel diagnostic, therapeutic and prognostic approaches to patients with heart failure with preserved ejection fraction. We hypothesize that novel imaging techniques and biomarkers will aid in describing heart failure with preserved ejection fraction. Furthermore, the development of new diagnostic criteria will allow us to: 1) better define risk factors associated with heart failure with preserved ejection fraction; 2) elucidate clinical, cellular and molecular mechanisms involved with the development and progression of heart failure with preserved ejection fraction; 3) design and test new therapeutic strategies for patients with heart failure with preserved ejection fraction. Additionally, Alberta HEART provides training and education for enhancing translational medicine, knowledge translation and clinical practice in heart failure. This is a prospective observational cohort study of patients with, or at risk for, heart failure. Patients will have sequential testing including quality of life and clinical outcomes over 12 months. After that time, study participants will be passively followed via linkage to external administrative databases. Clinical outcomes of interest include death, hospitalization, emergency department visits, physician resource use and/or heart transplant. Patients will be followed for a total of 5 years. Discussion Alberta HEART has the primary objective to define new diagnostic criteria for patients with heart failure with preserved ejection fraction. New criteria will allow for targeted therapies, diagnostic tests and further understanding of the patients, both at-risk for and with heart failure

Effect of cadence on locomotor–respiratory coupling during upper-body exercise

Introduction: Asynchronous arm-cranking performed at high cadences elicits greater cardiorespiratory responses compared to low cadences. This has been attributed to increased postural demand and locomotor–respiratory coupling (LRC), and yet, this has not been empirically tested. This study aimed to assess the effects of cadence on cardiorespiratory responses and LRC during upper-body exercise. Methods: Eight recreationally-active men performed arm-cranking exercise at moderate and severe intensities that were separated by 10 min of rest. At each intensity, participants exercised for 4 min at each of three cadences (50, 70, and 90 rev min−1) in a random order, with 4 min rest-periods applied in-between cadences. Exercise measures included LRC via whole- and half-integer ratios, cardiorespiratory function, perceptions of effort (RPE and dyspnoea), and diaphragm EMG using an oesophageal catheter. Results: The prevalence of LRC during moderate exercise was highest at 70 vs. 50 rev min−1 (27 ± 10 vs. 13 ± 9%, p = 0.000) and during severe exercise at 90 vs. 50 rev min−1 (24 ± 7 vs. 18 ± 5%, p = 0.034), with a shorter inspiratory time and higher mean inspiratory flow (p < 0.05) at higher cadences. During moderate exercise, (Formula presented.) and fC were higher at 90 rev min−1 (p < 0.05) relative to 70 and 50 rev min−1 ((Formula presented.) 1.19 ± 0.25 vs. 1.05 ± 0.21 vs. 0.97 ± 0.24 L min−1; fC 116 ± 11 vs. 101 ± 13 vs. 101 ± 12 b min−1), with concomitantly elevated dyspnoea. There were no discernible cadence-mediated effects on diaphragm EMG. Conclusion: Participants engage in LRC to a greater extent at moderate-high cadences which, in turn, increase respiratory airflow. Cadence rate should be carefully considered when designing aerobic training programmes involving the upper-limbs

Genome-Wide Changes in Genetic Diversity in a Population of Myotis lucifugus Affected by White-Nose Syndrome

Novel pathogens can cause massive declines in populations, and even extirpation of hosts. But disease can also act as a selective pressure on survivors, driving the evolution of resistance or tolerance. Bat white-nose syndrome (WNS) is a rapidly spreading wildlife disease in North America. The fungus causing the disease invades skin tissues of hibernating bats, resulting in disruption of hibernation behavior, premature energy depletion, and subsequent death. We used whole-genome sequencing to investigate changes in allele frequencies within a population of Myotis lucifugus in eastern North America to search for genetic resistance to WNS. Our results show low F-ST values within the population across time, i.e., prior to WNS (Pre-WNS) compared to the population that has survived WNS (Post-WNS). However, when dividing the population with a geographical cut-off between the states of Pennsylvania and New York, a sharp increase in values on scaffold GL429776 is evident in the Post-WNS samples. Genes present in the diverged area are associated with thermoregulation and promotion of brown fat production. Thus, although WNS may not have subjected the entire M. lucifugus population to selective pressure, it may have selected for specific alleles in Pennsylvania through decreased gene flow within the population. However, the persistence of remnant sub-populations in the aftermath of WNS is likely due to multiple factors in bat life history