The unconventional secretion of stress-inducible protein 1 by a heterogeneous population of extracellular vesicles

The co-chaperone stress-inducible protein 1 (STI1) is released by astrocytes, and has important neurotrophic properties upon binding to prion protein (PrPC). However, STI1 lacks a signal peptide and pharmacological approaches pointed that it does not follow a classical secretion mechanism. Ultracentrifugation, size exclusion chromatography, electron microscopy, vesicle labeling, and particle tracking analysis were used to identify three major types of extracellular vesicles (EVs) released from astrocytes with sizes ranging from 20–50, 100–200, and 300–400 nm. These EVs carry STI1 and present many exosomal markers, even though only a subpopulation had the typical exosomal morphology. The only protein, from those evaluated here, present exclusively in vesicles that have exosomal morphology was PrPC. STI1 partially co-localized with Rab5 and Rab7 in endosomal compartments, and a dominant-negative for vacuolar protein sorting 4A (VPS4A), required for formation of multivesicular bodies (MVBs), impaired EV and STI1 release. Flow cytometry and PK digestion demonstrated that STI1 localized to the outer leaflet of EVs, and its association with EVs greatly increased STI1 activity upon PrPC-dependent neuronal signaling. These results indicate that astrocytes secrete a diverse population of EVs derived from MVBs that contain STI1 and suggest that the interaction between EVs and neuronal surface components enhances STI1–PrPC signaling

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

A multiparental cross population for mapping QTL for relevant agronomic traits in durum wheat (Triticum durum Desf.)

Multiparental cross designs for mapping quantitative trait loci (QTL) in crops are efficient alternatives to conventional biparental experimental populations because they exploit a broader genetic basis and higher mapping resolution. We describe the development and deployment of a multiparental recombinant inbred line (RIL) population in durum wheat (Triticum durum Desf.) obtained by crossing four elite cultivars characterized by different traits of agronomic value. A linkage map spanning 2,663 cM and including 7,594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs with a wheat-dedicated 90k SNP chip. A cluster file was developed for correct allele calling in the framework of the tetraploid durum wheat genome. Based on phenotypic data collected over four field experiments, a multi-trait quantitative trait loci (QTL) analysis was carried out for 18 traits of agronomic relevance (including yield, yield-components, morpho-physiological and seed quality traits). Across environments, a total of 63 QTL were identified and characterized in terms of the four founder haplotypes. We mapped two QTL for grain yield across environments and 23 QTL for grain yield components. A novel major QTL for number of grain per spikelet/ear was mapped on chr 2A and shown to control up to 39% of phenotypic variance in this cross. Functionally different QTL alleles, in terms of direction and size of genetic effect, were distributed among the four parents. Based on the occurrence of QTL-clusters, we characterized the breeding values (in terms of effects on yield) of most of QTL for heading and maturity as well as yield component and quality QTL. This multiparental RIL population provides the wheat community with a highly informative QTL mapping resource enabling the dissection of the genetic architecture of multiple agronomic relevant traits in durum wheat.

Misleading Population Estimates: Biases and Consistency of Visual Surveys and Matrix Modelling in the Endangered Bearded Vulture

Conservation strategies for long-lived vertebrates require accurate estimates of parameters relative to the populations' size, numbers of non-breeding individuals (the “cryptic” fraction of the population) and the age structure. Frequently, visual survey techniques are used to make these estimates but the accuracy of these approaches is questionable, mainly because of the existence of numerous potential biases. Here we compare data on population trends and age structure in a bearded vulture (Gypaetus barbatus) population from visual surveys performed at supplementary feeding stations with data derived from population matrix-modelling approximations. Our results suggest that visual surveys overestimate the number of immature (<2 years old) birds, whereas subadults (3–5 y.o.) and adults (>6 y.o.) were underestimated in comparison with the predictions of a population model using a stable-age distribution. In addition, we found that visual surveys did not provide conclusive information on true variations in the size of the focal population. Our results suggest that although long-term studies (i.e. population matrix modelling based on capture-recapture procedures) are a more time-consuming method, they provide more reliable and robust estimates of population parameters needed in designing and applying conservation strategies. The findings shown here are likely transferable to the management and conservation of other long-lived vertebrate populations that share similar life-history traits and ecological requirements.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Heritability, determinants and reference values of renal length: a family-based population study

OBJECTIVES: In this population-based study, reference values were generated for renal length, and the heritability and factors associated with kidney length were assessed. METHODS: Anthropometric parameters and renal ultrasound measurements were assessed in randomly selected nuclear families of European ancestry (Switzerland). The adjusted narrow sense heritability of kidney size parameters was estimated by maximum likelihood assuming multivariate normality after power transformation. Gender-specific reference centiles were generated for renal length according to body height in the subset of non-diabetic non-obese participants with normal renal function. RESULTS: We included 374 men and 419 women (mean ± SD, age 47 ± 18 and 48 ± 17 years, BMI 26.2 ± 4 and 24.5 ± 5 kg/m(2), respectively) from 205 families. Renal length was 11.4 ± 0.8 cm in men and 10.7 ± 0.8 cm in women; there was no difference between right and left renal length. Body height, weight and estimated glomerular filtration rate (eGFR) were positively associated with renal length, kidney function negatively, age quadratically, whereas gender and hypertension were not. The adjusted heritability estimates of renal length and volume were 47.3 ± 8.5 % and 45.5 ± 8.8 %, respectively (P < 0.001). CONCLUSION: The significant heritability of renal length and volume highlights the familial aggregation of this trait, independently of age and body size. Population-based references for renal length provide a useful guide for clinicians. KEY POINTS: • Renal length and volume are heritable traits, independent of age and size. • Based on a European population, gender-specific reference values/percentiles are provided for renal length. • Renal length correlates positively with body length and weight. • There was no difference between right and left renal lengths in this study. • This negates general teaching that the left kidney is larger and longer.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

DEVELOPMENT OF AN OPTIMIZATION MODEL FOR BIOFUEL FACILITY SIZE AND LOCATION AND A SIMULATION MODEL FOR DESIGN OF A BIOFUEL SUPPLY CHAIN

A range of societal issues have been caused by fossil fuel consumption in the transportation sector in the United States (U.S.), including health related air pollution, climate change, the dependence on imported oil, and other oil related national security concerns. Biofuels production from various lignocellulosic biomass types such as wood, forest residues, and agriculture residues have the potential to replace a substantial portion of the total fossil fuel consumption. This research focuses on locating biofuel facilities and designing the biofuel supply chain to minimize the overall cost. For this purpose an integrated methodology was proposed by combining the GIS technology with simulation and optimization modeling methods. The GIS based methodology was used as a precursor for selecting biofuel facility locations by employing a series of decision factors. The resulted candidate sites for biofuel production served as inputs for simulation and optimization modeling. As a precursor to simulation or optimization modeling, the GIS-based methodology was used to preselect potential biofuel facility locations for biofuel production from forest biomass. Candidate locations were selected based on a set of evaluation criteria, including: county boundaries, a railroad transportation network, a state/federal road transportation network, water body (rivers, lakes, etc.) dispersion, city and village dispersion, a population census, biomass production, and no co-location with co-fired power plants. The simulation and optimization models were built around key supply activities including biomass harvesting/forwarding, transportation and storage. The built onsite storage served for spring breakup period where road restrictions were in place and truck transportation on certain roads was limited. Both models were evaluated using multiple performance indicators, including cost (consisting of the delivered feedstock cost, and inventory holding cost), energy consumption, and GHG emissions. The impact of energy consumption and GHG emissions were expressed in monetary terms to keep consistent with cost. Compared with the optimization model, the simulation model represents a more dynamic look at a 20-year operation by considering the impacts associated with building inventory at the biorefinery to address the limited availability of biomass feedstock during the spring breakup period. The number of trucks required per day was estimated and the inventory level all year around was tracked. Through the exchange of information across different procedures (harvesting, transportation, and biomass feedstock processing procedures), a smooth flow of biomass from harvesting areas to a biofuel facility was implemented. The optimization model was developed to address issues related to locating multiple biofuel facilities simultaneously. The size of the potential biofuel facility is set up with an upper bound of 50 MGY and a lower bound of 30 MGY. The optimization model is a static, Mathematical Programming Language (MPL)-based application which allows for sensitivity analysis by changing inputs to evaluate different scenarios. It was found that annual biofuel demand and biomass availability impacts the optimal results of biofuel facility locations and sizes.

Impact of vessel size on angiographic and clinical outcomes of revascularization with biolimus-eluting stent with biodegradable polymer and sirolimus-eluting stent with durable polymer the LEADERS trial substudy

OBJECTIVES: We assessed the impact of vessel size on outcomes of stenting with biolimus-eluting degradable polymer stent (BES) and sirolimus-eluting permanent polymer stent (SES) within a randomized multicenter trial (LEADERS). BACKGROUND: Stenting of small vessels might be associated with higher rates of adverse events. METHODS: "All-comer" patients (n = 1,707) were randomized to BES and SES. Post-hoc-stratified analysis of angiographic and clinical outcomes at 9 months and 1 year, respectively, was performed for vessels with reference diameter 2.75 mm. RESULTS: Of 1,707 patients, 429 patients in the BES group with 576 lesions and 434 patients in the SES group with 557 lesions had only small vessels treated (50.6% of the patient cohort). In patients with small vessels there was no significant difference in overall major adverse cardiac events (MACE) rate (12.1% vs. 11.8%; p = 0.89) or target lesion revascularization (TLR) rate (9.6% vs. 7.4%; p = 0.26) between BES and SES. The MACE and TLR rates in the small-vessel patient population were higher than in the large-vessel population. The TLR rate was 9.6% versus 2.6%, and MACE rate was 12.1% versus 7.1% for small versus large vessels in the BES arm (TLR: hazard ratio [HR] = 3.724, p = 0.0013; MACE: HR = 1.720, p = 0.0412). In the SES arm, TLR was 7.4% versus 5.1%, and MACE was 11.8% versus 10.3% in small versus large vessels (TLR: HR = 1.435, p = 0.2594; MACE: HR = 1.149, p = 0.5546). CONCLUSIONS: Prevalence of small vessel disease is high in an "all-comer" population with higher TLR and MACE rates. The BES and SES seem equivalent in treatment outcomes of small vessels in this "all-comer" patient population.

Clinical and intravascular imaging outcomes at 1 and 2 years after implantation of absorb everolimus eluting bioresorbable vascular scaffolds in small vessels. Late lumen enlargement: does bioresorption matter with small vessel size? Insight from the ABSORB cohort B trial

BACKGROUND The long-term results after second generation everolimus eluting bioresorbable vascular scaffold (Absorb BVS) placement in small vessels are unknown. Therefore, we investigated the impact of vessel size on long-term outcomes, after Absorb BVS implantation. METHODS In ABSORB Cohort B Trial, out of the total study population (101 patients), 45 patients were assigned to undergo 6-month and 2-year angiographic follow-up (Cohort B1) and 56 patients to have angiographic follow-up at 1-year (Cohort B2). The pre-reference vessel diameter (RVD) was <2.5 mm (small-vessel group) in 41 patients (41 lesions) and ≥2.5 mm (large-vessel group) in 60 patients (61 lesions). Outcomes were compared according to pre-RVD. RESULTS At 2-year angiographic follow-up no differences in late lumen loss (0.29±0.16 mm vs 0.25±0.22 mm, p=0.4391), and in-segment binary restenosis (5.3% vs 5.3% p=1.0000) were demonstrated between groups. In the small-vessel group, intravascular ultrasound analysis showed a significant increase in vessel area (12.25±3.47 mm(2) vs 13.09±3.38 mm(2) p=0.0015), scaffold area (5.76±0.96 mm(2) vs 6.41±1.30 mm(2) p=0.0008) and lumen area (5.71±0.98 mm(2) vs 6.20±1.27 mm(2) p=0.0155) between 6-months and 2-year follow-up. No differences in plaque composition were reported between groups at either time point. At 2-year clinical follow-up, no differences in ischaemia-driven major adverse cardiac events (7.3% vs 10.2%, p=0.7335), myocardial infarction (4.9% vs 1.7%, p=0.5662) or ischaemia-driven target lesion revascularisation (2.4% vs 8.5%, p=0.3962) were reported between small and large vessels. No deaths or scaffold thrombosis were observed. CONCLUSIONS Similar clinical and angiographic outcomes at 2-year follow-up were reported in small and large vessel groups. A significant late lumen enlargement and positive vessel remodelling were observed in small vessels.

A spatial-temporal approach to surveillance of prostate cancer disparities in population subgroups

BACKGROUND: Prostate cancer mortality disparities exist among racial/ethnic groups in the United States, yet few studies have explored the spatiotemporal trend of the disease burden. To better understand mortality disparities by geographic regions over time, the present study analyzed the geographic variations of prostate cancer mortality by three Texas racial/ethnic groups over a 22-year period. METHODS: The Spatial Scan Statistic developed by Kulldorff et al was used. Excess mortality was detected using scan windows of 50% and 90% of the study period and a spatial cluster size of 50% of the population at risk. Time trend was analyzed to examine the potential temporal effects of clustering. Spatial queries were used to identify regions with multiple racial/ethnic groups having excess mortality. RESULTS: The most likely area of excess mortality for blacks occurred in Dallas-Metroplex and upper east Texas areas between 1990 and 1999; for Hispanics, in central Texas between 1992 and 1996: and for non-Hispanic whites, in the upper south and west to central Texas areas between 1990 and 1996. Excess mortality persisted among all racial/ethnic groups in the identified counties. The second scan revealed that three counties in west Texas presented an excess mortality for Hispanics from 1980-2001. Many counties bore an excess mortality burden for multiple groups. There is no time trend decline in prostate cancer mortality for blacks and non-Hispanic whites in Texas. CONCLUSION: Disparities in prostate cancer mortality among racial/ethnic groups existed in Texas. Central Texas counties with excess mortality in multiple subgroups warrant further investigation.

The Role of Cell Sterilization in Population Based Studies of Radiogenic Second Cancers Following Radiation Therapy

Advances in radiotherapy have generated increased interest in comparative studies of treatment techniques and their effectiveness. In this respect, pediatric patients are of specific interest because of their sensitivity to radiation induced second cancers. However, due to the rarity of childhood cancers and the long latency of second cancers, large sample sizes are unavailable for the epidemiological study of contemporary radiotherapy treatments. Additionally, when specific treatments are considered, such as proton therapy, sample sizes are further reduced due to the rareness of such treatments. We propose a method to improve statistical power in micro clinical trials. Specifically, we use a more biologically relevant quantity, cancer equivalent dose (DCE), to estimate risk instead of mean absorbed dose (DMA). Our objective was to demonstrate that when DCE is used fewer subjects are needed for clinical trials. Thus, we compared the impact of DCE vs. DMA on sample size in a virtual clinical trial that estimated risk for second cancer (SC) in the thyroid following craniospinal irradiation (CSI) of pediatric patients using protons vs. photons. Dose reconstruction, risk models, and statistical analysis were used to evaluate SC risk from therapeutic and stray radiation from CSI for 18 patients. Absorbed dose was calculated in two ways: with (1) traditional DMA and (2) with DCE. DCE and DMA values were used to estimate relative risk of SC incidence (RRCE and RRMA, respectively) after proton vs. photon CSI. Ratios of RR for proton vs. photon CSI (RRRCE and RRRMA) were then used in comparative estimations of sample size to determine the minimal number of patients needed to maintain 80% statistical power when using DCE vs. DMA. For all patients, we found that protons substantially reduced the risk of developing a second thyroid cancer when compared to photon therapy. Mean RRR values were 0.052±0.014 and 0.087±0.021 for RRRMA and RRRCE, respectively. However, we did not find that use of DCE reduced the number of patents needed for acceptable statistical power (i.e, 80%). In fact, when considerations were made for RRR values that met equipoise requirements and the need for descriptive statistics, the minimum number of patients needed for a micro-clinical trial increased from 17 using DMA to 37 using DCE. Subsequent analyses revealed that for our sample, the most influential factor in determining variations in sample size was the experimental standard deviation of estimates for RRR across the patient sample. Additionally, because the relative uncertainty in dose from proton CSI was so much larger (on the order of 2000 times larger) than the other uncertainty terms, it dominated the uncertainty in RRR. Thus, we found that use of corrections for cell sterilization, in the form of DCE, may be an important and underappreciated consideration in the design of clinical trials and radio-epidemiological studies. In addition, the accurate application of cell sterilization to thyroid dose was sensitive to variations in absorbed dose, especially for proton CSI, which may stem from errors in patient positioning, range calculation, and other aspects of treatment planning and delivery.

Pharmacological and non-pharmacological therapy for arrhythmias in the pediatric population: EHRA and AEPC-Arrhythmia Working Group joint consensus statement

In children with structurally normal hearts, the mechanisms of arrhythmias are usually the same as in the adult patient. Some arrhythmias are particularly associated with young age and very rarely seen in adult patients. Arrhythmias in structural heart disease may be associated either with the underlying abnormality or result from surgical intervention. Chronic haemodynamic stress of congenital heart disease (CHD) might create an electrophysiological and anatomic substrate highly favourable for re-entrant arrhythmias. As a general rule, prescription of antiarrhythmic drugs requires a clear diagnosis with electrocardiographic documentation of a given arrhythmia. Risk-benefit analysis of drug therapy should be considered when facing an arrhythmia in a child. Prophylactic antiarrhythmic drug therapy is given only to protect the child from recurrent supraventricular tachycardia during this time span until the disease will eventually cease spontaneously. In the last decades, radiofrequency catheter ablation is progressively used as curative therapy for tachyarrhythmias in children and patients with or without CHD. Even in young children, procedures can be performed with high success rates and low complication rates as shown by several retrospective and prospective paediatric multi-centre studies. Three-dimensional mapping and non-fluoroscopic navigation techniques and enhanced catheter technology have further improved safety and efficacy even in CHD patients with complex arrhythmias. During last decades, cardiac devices (pacemakers and implantable cardiac defibrillator) have developed rapidly. The pacing generator size has diminished and the pacing leads have become progressively thinner. These developments have made application of cardiac pacing in children easier although no dedicated paediatric pacing systems exist.

Male Dominance in the New Zealand Longfin Eel Population of a New Zealand River: Probable Causes and Implications for Management

The endemic New Zealand longfin eel Anguilla dieffenbachi (hereafter, longfin eel), is overfished, and in southern South Island, New Zealand, rivers have recently become predominated by males. This study examined length and age at sexual differentiation in male eels in the Aparima River catchment (area, 1,375 km(2); mean flow, 20 m(3.)s(-1)) and the sex ratio and distribution of eels throughout the catchment. Longfin eels differentiated into males mostly at lengths from 300 to 460 mm and ages from 10 to 25+ years. Females were rare: Of 738 eels examined for sexual differentiation, 466 were males and 5 were females, and a few others, not examined, were large enough to be female. These counts suggest a male : female ratio among differentiated longfin eels of 68:1. Of 31 differentiated shortfin eels A. australis, less common in the Aparima River, 26 were females. Male longfin eels were distributed throughout the main stern and tributaries; undifferentiated eels were more prevalent in lower and middle reaches and in the main stem than in upper reaches and tributaries. In other studies, male longfin eels predominated commercial catches in the Aparima and four other southernmost rivers, by 2.4:1 to 13.6:1 males to females. The Aparima River had the most skewed sex ratio. Longfin eel catches from the Aparima River will become more male predominated because few sublegal-size females were present. The length-frequency distributions of eels in the present samples and in the commercial catches were truncated just above minimum legal size (about 460 mm), showing that few females escape the fishery. Historically, females predominated these rivers. The recent change in sex ratio is attributable partly to selective harvest of females, and partly to changes in the structure of the population from fishing, such that differentiation into males has been favored. Longevity, delayed sexual maturity, semel-parity, and endemism with restricted range make the longfin eel particularly vulnerable to overfishing.

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.