A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Effects of laboratory colonization on Bactrocera dorsalis (Diptera, Tephritidae) mating behaviour: ‘what a difference a year makes

Laboratory-reared insects are widely known to have significantly reduced genetic diversity in comparison to wild populations; however, subtle behavioural changes between laboratory-adapted and wild or ‘wildish’ (i.e., within one or very few generations of field collected material) populations are less well understood. Quantifying alterations in behaviour, particularly sexual, in laboratory-adapted insects is important for mass-reared insects for use in pest management strategies, especially those that have a sterile insect technique component. We report subtle changes in sexual behaviour between ‘wildish’ Bactrocera dorsalis flies (F1 and F2) from central and southern Thailand and the same colonies 12 months later when at six generations from wild. Mating compatibility tests were undertaken under standardised semi-natural conditions, with number of homo/heterotypic couples and mating location in field cages analysed via compatibility indices. Central and southern populations of B. dorsalis displayed positive assortative mating in the 2010 trials but mated randomly in the 2011 trials. ‘Wildish’ southern Thailand males mated significantly earlier than central Thailand males in 2010; this difference was considerably reduced in 2011, yet homotypic couples from southern Thailand still formed significantly earlier than all other couple combinations. There was no significant difference in couple location in 2010; however, couple location significantly differed among pair types in 2011 with those involving southern Thailand females occurring significantly more often on the tree relative to those with central Thailand females. Relative participation also changed with time, with more southern Thailand females forming couples relative to central Thailand females in 2010; this difference was considerably decreased by 2011. These results reveal how subtle changes in sexual behaviour, as driven by laboratory rearing conditions, may significantly influence mating behaviour between laboratory-adapted and recently colonised tephritid fruit flies over a relatively short period of time.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10−9 to P = 1.8 × 10−40) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10−3 to P = 1.2 × 10−13). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

Phenotyping novel stay-green traits to capture genetic variation in senescence dynamics

Stay-green plants retain green leaves longer after anthesis and can have improved yield, particularly under water limitation. As senescence is a dynamic process, genotypes with different senescence patterns may exhibit similar final normalised difference vegetative index (NDVI). By monitoring NDVI from as early as awn emergence to maturity, we demonstrate that analysing senescence dynamics improves insight into genotypic stay-green variation. A senescence evaluation tool was developed to fit a logistic function to NDVI data and used to analyse data from three environments for a wheat (Triticum aestivum L.) population whose lines contrast for stay-green. Key stay-green traits were estimated including, maximum NDVI, senescence rate and a trait integrating NDVI variation after anthesis, as well as the timing from anthesis to onset, midpoint and conclusion of senescence. The integrative trait and the timing to onset and mid-senescence exhibited high positive correlations with yield and a high heritability in the three studied environments. Senescence rate was correlated with yield in some environments, whereas maximum NDVI was associated with yield in a drought-stressed environment. Where resources preclude frequent measurements, we found that NDVI measurements may be restricted to the period of rapid senescence, but caution is required when dealing with lines of different phenology. In contrast, regular monitoring during the whole period after flowering allows the estimation of senescence dynamics traits that may be reliably compared across genotypes and environments. We anticipate that selection for stay-green traits will enhance genetic progress towards high-yielding, stay-green germplasm.

Genetic and epigenetic variants in the MTHFR gene are not associated with non-Hodgkin lymphoma

The methylenetetrahydrofolate reductase (MTHFR) gene codes for the MTHFR enzyme which plays a key role in the pathway of folate and methionine metabolism. Polymorphisms of genes in this pathway affect its regulation and have been linked to lymphoma. In this study we examined whether we could detect an association between two common non-synonomous MTHFR polymorphisms, 677C>T (rs1801133) and 1298A>C (rs1801131), and susceptibility to non-Hodgkin lymphoma (NHL) in an Australian case-control cohort. We found no significant differences between genotype or allele frequencies for either polymorphisms between lymphoma cases and controls. We also explored whether epigenetic modification of MTHFR, specifically DNA methylation of a CpG island in the MTHFR promoter region, is associated with NHL using blood samples from patients. No difference in methylation levels was detected between the case and control samples suggesting that although hypermethylation of MTHFR has been reported in tumour tissues, particularly in the diffuse large B-cell lymphoma subtype of NHL, methylation of this MTHFR promoter CpG island is not a suitable epigenetic biomarker for NHL diagnosis or prognosis in peripheral blood samples. Further studies into epigenetic variants could focus on genes that are robustly associated with NHL susceptibility.

SPECT-aided diagnostics and genetic risk factors in Parkinson's disease

Hypokinesia, rigidity, tremor, and postural instability are the cardinal symptoms of Parkinson s disease (PD). Since these symptoms are not specific to PD the diagnosis may be uncertain in early PD. Etiology and pathogenesis of PD remain unclear. There is no neuroprotective therapy. Genetic findings are expected to reveal metabolic routes in PD pathogenesis and thereby eventually lead to therapeutic innovations. In this thesis, we first aimed to study the usefulness and accuracy of 123I-b-CIT SPECT in the diagnosis of PD in a consecutive clinic-based material including various movement disorders. We subsequently a genetic project to identify genetic risk factors for sporadic PD using a candidate gene approach in a case-control setting including 147 sporadic PD patients and 137 spouse controls. Dopamine transporter imaging by 123I-b-CIT SPECT could distinguish PD from essential tremor, drug-induced parkinsonism, dystonia and psychogenic parkinsonism. However, b-CIT uptake in Parkinson plus syndromes (PSP and multiple system atrophy) and dementia with Lewy bodies was not significantly different from PD. 123I-b-CIT SPECT could not reliably differentiate PD from vascular parkinsonism. 123I-b-CIT SPECT was 100% sensitive and specific in the diagnosis of PD in patients younger than 55 years but less specific in older patients, due to differential distribution of the above conditions in the younger and older age groups. 123I-b-CIT SPECT correlated with symptoms and detected bilateral nigrostriatal defect in patients whose PD was still in unilateral stage. Thus, in addition to as a differential diagnostic aid, 123I-b-CIT SPECT may be used to detect PD early, even pre-symptomatically in at-risk individuals. 123I-b-CIT SPECT was used to aid in the collection of patients to the genetic studies. In the genetic part of this thesis we found an association between PD and a polymorphic CAG-repeat in POLG1 gene encoding the catalytic subunit of mitochondrial polymerase gamma. The CAG-repeat encodes a polyglutamine tract (polyQ), the two most common lengths of which are 10Q (86-90%) and 11Q. In our Finnish material, the rarer non-10Q or non-11Q length variants (6Q-9Q, 12Q-14Q, 4R+9Q) were more frequent in patients than in spouse controls (10% vs. 3.5 %, p=0.003), or population controls (p=0.001). Therefore, we performed a replication study in 652 North American PD patients and 292 controls. Non-10/11Q alleles were more common in the US PD patients compared to the controls but the difference did not reach statistical significance (p=0.07). This larger data suggested our original definition of variant length allele might need reconsideration. Most previous studies on phenotypic effects of POLG1 polyQ have defined 10Q as the only normal allele. Non-10Q alleles were significantly more common in patients compared to the controls (17.3% vs. 12.3 %, p= 0.005). This association between non-10Q length variants and PD remained significant when compared to a larger set of 1541 literature controls (p=0.00005). In conclusion, POLG1 polyQ alleles other than 10Q may predispose to PD. We did not find association between PD and parkin or DJ-1, genes underlying autosomal recessive parkinsonism. The functional Val158Met polymorphism, which affects the catalytic effect of COMT enzyme, and another coding polymorphism in COMT were not associated with PD in our patient material. The APOE e2/3/4 polymorphism modifies risk for Alzheimer s disease and prognosis of for example brain trauma. APOE promoter and enhancer polymorphisms 219G/T and +113G/C, and APOE e3 haplotypes, have also been shown to modify the risk of Alzheimer s disease but not reported in PD. No association was found between PD and APOE e2/3/4 polymorphism, the promoter or enhancer polymorphisms, or the e3 haplotypes.

Genetic background of bipolar disorder and related cognitive impairments

Bipolar disorder (BP) is a complex psychiatric disorder characterized by episodes of mania and depression. BP affects approximately 1% of the world’s population and shows no difference in lifetime prevalence between males and females. BP arises from complex interactions among genetic, developmental and environmental factors, and it is likely that several predisposing genes are involved in BP. The genetic background of BP is still poorly understood, although intensive and long-lasting research has identified several chromosomal regions and genes involved in susceptibility to BP. This thesis work aims to identify the genetic variants that influence bipolar disorder in the Finnish population by candidate gene and genome-wide linkage analyses in families with many BP cases. In addition to diagnosis-based phenotypes, neuropsychological traits that can be seen as potential endophenotypes or intermediate traits for BP were analyzed. In the first part of the thesis, we examined the role of the allelic variants of the TSNAX/DISC1 gene cluster to psychotic and bipolar spectrum disorders and found association of distinct allelic haplotypes with these two groups of disorders. The haplotype at the 5’ end of the Disrupted-in-Schizophrenia-1 gene (DISC1) was over-transmitted to males with psychotic disorder (p = 0.008; for an extended haplotype p = 0.0007 with both genders), whereas haplotypes at the 3’ end of DISC1 associated with bipolar spectrum disorder (p = 0.0002; for an extended haplotype p = 0.0001). The variants of these haplotypes also showed association with different cognitive traits. The haplotypes at the 5’ end associated with perseverations and auditory attention, while the variants at the 3’ end associated with several cognitive traits including verbal fluency and psychomotor processing speed. Second, in our complete set of BP families with 723 individuals we studied six functional candidate genes from three distinct signalling systems: serotonin-related genes (SLC6A4 and TPH2), BDNF -related genes (BDNF, CREB1 and NTRK2) and one gene related to the inflammation and cytokine system (P2RX7). We replicated association of the functional variant Val66Met of BDNF with BP and better performance in retention. The variants at the 5’ end of SLC6A4 also showed some evidence of association among males (p = 0.004), but the widely studied functional variants did not yield any significant results. A protective four-variant haplotype on P2RX7 showed evidence of association with BP and executive functions: semantic and phonemic fluency (p = 0.006 and p = 0.0003, respectively). Third, we analyzed 23 bipolar families originating from the North-Eastern region of Finland. A genome-wide scan was performed using the 6K single nucleotide polymorphism (SNP) array. We identified susceptibility loci at chromosomes 7q31 with a LOD score of 3.20 and at 9p13.1 with a LOD score of 4.02. We followed up both linkage findings in the complete set of 179 Finnish bipolar families. The finding on chromosome 9p13 was supported (maximum LOD score of 3.02), but the susceptibility gene itself remains unclarified. In the fourth part of the thesis, we wanted to test the role of the allelic variants that have associated with bipolar disorder in recent genome-wide association studies (GWAS). We could confirm findings for the DFNB31, SORCS2, SCL39A3, and DGKH genes. The best signal in this study comes from DFNB31, which remained significant after multiple testing corrections. Two variants of SORCS2 were allelic replications and presented the same signal as the haplotype analysis. However, no association was detected with the PALB2 gene, which was the most significantly associated region in the previous GWAS. Our results indicate that BP is heterogeneous and its genetic background may accordingly vary in different populations. In order to fully understand the allelic heterogeneity that underlies common diseases such as BP, complete genome sequencing for many individuals with and without the disease is required. Identification of the specific risk variants will help us better understand the pathophysiology underlying BP and will lead to the development of treatments with specific biochemical targets. In addition, it will further facilitate the identification of environmental factors that alter risk, which will potentially provide improved occupational, social and psychological advice for individuals with high risk of BP.

Origins of the sex differences in handedness and mental rotation ability : genetic, environmental and hormonal effects

In humans, well-replicated and robust sex differences in cognitive functions exist for handedness and mental rotation ability. A common characteristic in human cognitive functions is the lateralization of language functions. Handedness is a common measure of laterality and is related to language lateralization. The prevalence of left-handedness is higher in males than in females, the male to female ratio being about 1.2. Among cognitive abilities, the largest sex difference is evident in the Vandenberg and Kuse Mental Rotation Test (MRT), which requires the ability to rotate objects in mental space. On average, males achieve scores one standard deviation higher than females in the MRT. The present thesis investigated the origins of the sex differences in laterality and spatial ability as represented by handedness and mental rotation ability, respectively. Two population-based Finnish twin cohorts were utilized in this study. Handedness was studied in 25 810 twins and 4068 singletons born before 1958 from the Older Finnish Twin Cohort, and in 4736 twins born in 1983-87 from the FinnTwin12. MRT was studied in a sub-sample of 804 young adult participants from the FinnTwin12 sample. The main findings of this study were: 1) the prevalence of left-handedness was higher among males than among females in both singletons and in twins; 2) males had significantly higher scores than females in MRT; 3) about one quarter of the variance in handedness and about half of the variance in MRT was explained by genetic effects, whereas the remainder of the variance in these traits was explained by environmental effects unique to each individual. The magnitude of the genetic effects was similar in both sexes; 4) left-handedness was significantly less common in female co-twins of a male than in female co-twins of a female, and female co-twins of a male scored significantly higher than did female co-twins of a female in the Mental Rotation Test. This dissertation discusses whether these differences between females from opposite- and same-sex twin pairs are due to the prenatal transfer of testosterone from the male fetus in females with male co-twins or whether they arise from postnatal socialization effects.

Genetic Algorithm for Damage Assessment

Genetic Algorithms (GAs) are recognized as an alternative class of computational model, which mimic natural evolution to solve problems in a wide domain including machine learning, music generation, genetic synthesis etc. In the present study Genetic Algorithm has been employed to obtain damage assessment of composite structural elements. It is considered that a state of damage can be modeled as reduction in stiffness. The task is to determine the magnitude and location of damage. In a composite plate that is discretized into a set of finite elements, if a jth element is damaged, the GA based technique will predict the reduction in Ex and Ey and the location j. The fact that the natural frequency decreases with decrease in stiffness is made use of in the method. The natural frequency of any two modes of the damaged plates for the assumed damage parameters is facilitated by the use of Eigen sensitivity analysis. The Eigen value sensitivities are the derivatives of the Eigen values with respect to certain design parameters. If ωiu is the natural frequency of the ith mode of the undamaged plate and ωid is that of the damaged plate, with δωi as the difference between the two, while δωk is a similar difference in the kth mode, R is defined as the ratio of the two. For a random selection of Ex,Ey and j, a ratio Ri is obtained. A proper combination of Ex,Ey and j which makes Ri−R=0 is obtained by Genetic Algorithm.

Fast and accurate quantification using Genetic Algorithm optimized H-1-C-13 refocused constant-time INEPT

Using Genetic Algorithm, a global optimization method inspired by nature's evolutionary process, we have improved the quantitative refocused constant-time INEPT experiment (Q-INEPT-CT) of Makela et al. (JMR 204 (2010) 124-130) with various optimization constraints. The improved `average polarization transfer' and `min-max difference' of new delay sets effectively reduces the experimental time by a factor of two (compared with Q-INEPT-CT, Makela et al.) without compromising on accuracy. We also discuss a quantitative spectral editing technique based on average polarization transfer. (C) 2013 Elsevier Inc. All rights reserved.

Genetic differentiation between Atlantic salmon populations in the Windermere catchment

Genetic analysis, using single locus probes for genomic DNA, revealed that the juvenile Atlantic salmon populations in the Rivers Leven, Rothay and Troutbeck were related but genetically distinct. This genetic differentiation is greater than might be expected (by comparison with other salmon populations in the UK) and it is recommended that no action is taken which might promote genetic exchange between the three rivers. Thus, future fisheries management practices should treat the salmon from each site as separate genetic stocks. It is unlikely that any attempts to encourage fish currently spawning in the River Leven (downstream of Windermere) to utilize the upper catchment will be successful. The faster growth rate of juvenile salmon in the River Leven, compared with the River Rothay, probably results from a difference in temperature between the inflowing streams and the main outflow of Windermere. Precocious sexual maturation of some male parr was found in all three populations but the incidence (13-33%) is well within the range reported for other waters. Because of their enhanced growth rate, it is likely that some of the precocious males in the River Leven were 0+ fish. A very high incidence of hybridization (>18%) between Atlantic salmon and brown/sea trout was found in Troutbeck but not in the other rivers. Mitochondrial DNA analysis of these hybrids revealed them to be the product of several, independent cross-fertilizations involving both sexes of both species. The implications of this finding are discussed in relation to the availability of suitable spawning sites in Troutbeck.

Genetic variation of rougheye rockfish (Sebastes aleutianus) and shortraker rockfish (S. borealis) inferred from allozymes

Rougheye rockfish (Sebastes aleutianus) and shortraker rockfish (Sebastes borealis) were collected from the Washington coast, the Gulf of Alaska, the southern Bering Sea, and the eastern Kamchatka coast of Russia (areas encompassing most of their geographic distribution) for population genetic analyses. Using starch gel electrophoresis, we analyzed 1027 rougheye rockfish and 615 shortraker rockfish for variation at 29 proteincoding loci. No genetic heterogeneity was found among shortraker rockfish throughout the sampled regions, although shortraker in the Aleutian Islands region, captured at deeper depths, were found to be significantly smaller in size than the shortraker caught in shallower waters from Southeast Alaska. Genetic analysis of the rougheye rockfish revealed two evolutionary lineages that exist in sympatry with little or no gene f low between them. When analyzed as two distinct species, neither lineage exhibited heterogeneity among regions. Sebastes aleutianus seems to inhabit waters throughout the Gulf of Alaska and more southern waters, whereas S. sp. cf. aleutianus inhabits waters throughout the Gulf of Alaska, Aleutian Islands, and Asia. The distribution of the two rougheye rockfish lineages may be related to depth where they are sympatric. The paler color morph, S. aleutianus, is found more abundantly in shallower waters and the darker color morph, Sebastes sp. cf. aleutianus, inhabits deeper waters. Sebastes sp. cf. aleutianus, also exhibited a significantly higher prevalence of two parasites, N. robusta and T. trituba, than did Sebastes aleutianus, in the 2001 samples, indicating a possible difference in habitat and (or) resource use between the two lineages.

An environmentally based growth model that uses finite difference calculus with maximum likelihood method: its application to the brackish water bivalve Corbicula japonica in Lake Abashiri, Japan

We present a growth analysis model that combines large amounts of environmental data with limited amounts of biological data and apply it to Corbicula japonica. The model uses the maximum-likelihood method with the Akaike information criterion, which provides an objective criterion for model selection. An adequate distribution for describing a single cohort is selected from available probability density functions, which are expressed by location and scale parameters. Daily relative increase rates of the location parameter are expressed by a multivariate logistic function with environmental factors for each day and categorical variables indicating animal ages as independent variables. Daily relative increase rates of the scale parameter are expressed by an equation describing the relationship with the daily relative increase rate of the location parameter. Corbicula japonica grows to a modal shell length of 0.7 mm during the first year in Lake Abashiri. Compared with the attain-able maximum size of about 30 mm, the growth of juveniles is extremely slow because their growth is less susceptible to environmental factors until the second winter. The extremely slow growth in Lake Abashiri could be a geographical genetic variation within C. japonica.