Fitting mathematical functions to extended lactation curves and forecasting late-lactation milk yields of dairy cows.

A 305-d lactation followed by a 60-d dry period has traditionally been considered economically optimal, yet dairy cows in modern intensive dairy systems are frequently dried off while still producing significant quantities of milk. Managing cows for an extended lactation has reported production, welfare, and economic benefits, but not all cows are suitable for an extended lactation. Implementation of an extended lactation strategy on-farm could benefit from use of a decision support system, based on a mathematical lactation model, that can identify suitable cows during early lactation that have a high likelihood of producing above a target milk yield (MY) at 305 d in milk (DIM). Therefore, our objectives were (1) to compare the suitability of 3 commonly used lactation models for modeling extended lactations (Dijkstra, Wood, and Wilmink) in primiparous and multiparous cows under a variety of lactation lengths, and (2) to determine the amount of early-lactation daily MY data needed to accurately forecast MY at d 305 by using the most suitable model and determine whether this is sufficient for identifying cows suitable for an extended lactation before the end of a typical voluntary waiting period (50-90 d). Daily MY data from 467 individual Holstein-Friesian lactations (DIM >305 d; 379 ± 65-d lactation length [mean ± SD]) were fitted by the 3 lactation models using a nonlinear regression procedure. The parameter estimates of these models, lactation characteristics (peak yield, time to peak yield, and persistency), and goodness-of-fit were compared between parity and different lactation lengths. The models had similar performance, and differences between parity groups were consistent with previous literature. Then, data from only the first i DIM for each individual lactation, where i was incremented by 30 d from 30 to 150 DIM and by 50 d from 150 to 300 DIM, were fitted by each model to forecast MY at d 305. The Dijkstra model was selected for further analysis, as it had superior goodness-of-fit statistics for i= 30 and 60. The data set was fit twice by the Dijkstra model, with parameter bounds either unconstrained or constrained. The quality of predictions of MY at d 305 improved with increasing data availability for both models and assisting the model fitting procedure with more biologically relevant constraints on parameters improved the predictions, but neither was reliable enough for practical use on-farm due to the high uncertainty of forecasted predictions. Using 90 d of data, the constrained model correctly classified 66% of lactations as being above or below a target MY at d 305 of 25 kg/d, with a probability threshold of 0.95. The proportion of correct classifications became smaller at lower targets of MY at d 305 and became greater when using more lactation days. Overall, further work is required to develop a model that can forecast late-lactation MY with sufficient accuracy for practical use. We envisage that a hybridized machine learning and mechanistic model that incorporates additional historical and genetic information with early-lactation MY could produce meaningful lactation curve forecasts.

Looking beyond the surface: Muir Torre syndrome.

Muir-Torre Syndrome (MTS) is associated with multiple visceral malignancies. Initial presentation may be a benign skin tumor mimicking a sebaceous cyst. This case report highlights the importance of early diagnosis, genetic testing, and multidisciplinary screening. A 67-year-old man was diagnosed with MTS following excision of a skin lesion (sebaceoma). He was declined both screening colonoscopy and genetic testing. Subsequently, advanced colon cancer was found following presentation with iron deficiency anemia, which ultimately led to palliation despite successful surgery. MTS can present insidiously with skin lesions clinically diagnosed as sebaceous cysts. Once MTS is suspected on histology, genetic testing and screening for MTS-related cancers is warranted. Better understanding of the genetic variants for MTS can aid in earlier diagnosis thus not dismissing the need for screening for MTS-related cancers.

Reduction in mucosa thickness is associated with changes in immune function in the colon mucosa during the weaning transition in Holstein bull dairy calves.

This study aims to characterize changes in the structure and the molecules related to immune function in the colon mucosa in dairy calves during the weaning transition (weaned at week 6 of age). Colon mucosa thickness, measured at week 5 to 8 and 12 of age, decreased for 2 weeks after weaning, but then recovered. Colon mucosa's transcriptome profiling at week 5, 7, and 12 of age was obtained using RNA-sequencing. Functional analysis showed that pathways related to immune function were up-regulated postweaning. A weighted gene co-expression network analysis identified 17 immune function related genes, expressed higher postweaning, which were negatively correlated with colon mucosa thickness, suggesting that these genes may be involved in colon mucosa inflammation and recovery from mucosa thickness decrement during the weaning transition. As such, it is important to determine the function of immune cells in the colon mucosa during the weaning transition in dairy calves.

Transcriptome profiling revealed that key rumen epithelium functions change in relation to short-chain fatty acids and rumen epithelium-attached microbiota during the weaning transition.

This study aims to characterize the functional changes of the rumen epithelium associated with ruminal short-chain fatty acid (SCFA) concentration and epithelium-attached microbes during the weaning transition in dairy calves. Ruminal SCFA concentrations were determined, and transcriptome and microbiota profiling in biopsied rumen papillae were obtained from Holstein calves before and after weaning using RNA- and amplicon sequencing. Metabolic pathway analysis showed that pathways related to SCFA metabolism and cell apoptosis were up- and down-regulated postweaning, respectively. Functional analysis showed that genes related to SCFA absorption, metabolism, and protective roles against oxidative stress were positively correlated with ruminal SCFA concentrations. The relative abundance of epithelium-attached Rikenellaceae RC9 gut group and Campylobacter was positively correlated with genes involved in SCFA absorption and metabolism, suggesting that these microbes can cooperatively affect host functions. Future research should examine the contribution of attenuated apoptosis on rumen epithelial functional shifts during the weaning transition.

Feeding an unsalable carrot total-mixed ration altered bacterial amino acid degradation in the rumen of lambs.

The objective of this study was to determine the influence of a total-mixed ration including unsalable carrots at 45% DM on the rumen microbiome; and the plasma, rumen and liver metabolomes. Carrots discarded at processing were investigated as an energy-dense substitute for barley grain in a conventional feedlot diet, and improved feed conversion efficiency by 25%. Here, rumen fluid was collected from 34 Merino lambs at slaughter (n = 16 control; n = 18 carrot) after a feeding period of 11-weeks. The V4 region of the 16S rRNA gene was sequenced to profile archaeal and bacterial microbe communities. Further, a comprehensive, targeted profile of known metabolites was constructed for blood plasma, rumen fluid and biopsied liver metabolites using a gas chromatography mass spectrometry (GC-MS) metabolomics approach. An in vitro batch culture was used to characterise ruminal fermentation including gas and methane (CH4) production. In vivo rumen microbial community structure of carrot fed lambs was dissimilar (P < 0.01; PERMANOVA), and all measures of alpha diversity were greater (P < 0.01), compared to those fed the control diet. Unclassified genera in Bacteroidales (15.9 ± 6.74% relative abundance; RA) were more abundant (P < 0.01) in the rumen fluid of carrot-fed lambs, while unclassified taxa in the Succinivibrionaceae family (11.1 ± 3.85% RA) were greater (P < 0.01) in the control. The carrot diet improved in vitro ruminal fermentation evidenced as an 8% increase (P < 0.01) in DM digestibility and a 13.8% reduction (P = 0.01) in CH4 on a mg/ g DM basis, while the control diet increased (P = 0.04) percentage of propionate within total VFA by 20%. Fourteen rumen fluid metabolites and 27 liver metabolites were influenced (P ≤ 0.05) by diet, while no effect (P ≥ 0.05) was observed in plasma metabolites. The carrot diet enriched (impact value = 0.13; P = 0.01) the tyrosine metabolism pathway (acetoacetic acid, dopamine and pyruvate), while the control diet enriched (impact value = 0.42; P ≤ 0.02) starch and sucrose metabolism (trehalose and glucose) in rumen fluid. This study demonstrated that feeding 45% DM unsalable carrots diversified bacterial communities in the rumen. These dietary changes influenced pathways of tyrosine degradation, such that previous improvements in feed conversion efficiency in lambs could be explained.

The 'Tommy Atkins' mango genome reveals candidate genes for fruit quality.

BACKGROUND: Mango, Mangifera indica L., an important tropical fruit crop, is grown for its sweet and aromatic fruits. Past improvement of this species has predominantly relied on chance seedlings derived from over 1000 cultivars in the Indian sub-continent with a large variation for fruit size, yield, biotic and abiotic stress resistance, and fruit quality among other traits. Historically, mango has been an orphan crop with very limited molecular information. Only recently have molecular and genomics-based analyses enabled the creation of linkage maps, transcriptomes, and diversity analysis of large collections. Additionally, the combined analysis of genomic and phenotypic information is poised to improve mango breeding efficiency. RESULTS: This study sequenced, de novo assembled, analyzed, and annotated the genome of the monoembryonic mango cultivar 'Tommy Atkins'. The draft genome sequence was generated using NRGene de-novo Magic on high molecular weight DNA of 'Tommy Atkins', supplemented by 10X Genomics long read sequencing to improve the initial assembly. A hybrid population between 'Tommy Atkins' x 'Kensington Pride' was used to generate phased haplotype chromosomes and a highly resolved phased SNP map. The final 'Tommy Atkins' genome assembly was a consensus sequence that included 20 pseudomolecules representing the 20 chromosomes of mango and included ~ 86% of the ~ 439 Mb haploid mango genome. Skim sequencing identified ~ 3.3 M SNPs using the 'Tommy Atkins' x 'Kensington Pride' mapping population. Repeat masking identified 26,616 genes with a median length of 3348 bp. A whole genome duplication analysis revealed an ancestral 65 MYA polyploidization event shared with Anacardium occidentale. Two regions, one on LG4 and one on LG7 containing 28 candidate genes, were associated with the commercially important fruit size characteristic in the mapping population. CONCLUSIONS: The availability of the complete 'Tommy Atkins' mango genome will aid global initiatives to study mango genetics.

Ocean currents and the population genetic signature of fish migrations.

Animal migrations are a fascinating and global phenomenon, yet they are often difficult to study and sometimes poorly understood. Here, we build on classic ecological theory by hypothesizing that some enigmatic spawning migrations across coastal marine habitats can be inferred from the population genetic signature of larval dispersal by ocean currents. We test this assumption by integrating spatially realistic simulations of alternative spawning migration routes, associated patterns of larval dispersal, and associated variation in the population genetic structure of eastern Australian sea mullet (Mugil cephalus). We then use simulation results to assess the implications of alternative spawning destinations for larval replenishment, and we contrast simulated against measured population genetic variation. Both analyses suggest that the spawning migrations of M. cephalus in eastern Australia are likely to be localized (approximately 100 km along the shore), and that spawning is likely to occur in inshore waters. Our conclusions are supported by multiple lines of evidence available through independent studies, but they challenge the more traditional assumption of a single, long-distance migration event with subsequent offshore spawning in the East Australian Current. More generally, our study operationalizes classic theory on the relationship between fish migrations, ocean currents, and reproductive success. However, rather than confirming the traditionally assumed adaptation of migratory behavior to dominant ocean current flow, our findings support the concept of a genetically measurable link between fish migrations and local oceanographic conditions, specifically water temperature and coastal retention of larvae. We believe that future studies using similar approaches for high resolution and spatially realistic ecological-genetic scenario testing can help rapidly advance our understanding of key ecological processes in many other marine species.

The Sorghum QTL Atlas: a powerful tool for trait dissection, comparative genomics and crop improvement.

KEY MESSAGE: We describe the development and application of the Sorghum QTL Atlas, a high-resolution, open-access research platform to facilitate candidate gene identification across three cereal species, sorghum, maize and rice. The mechanisms governing the genetic control of many quantitative traits are only poorly understood and have yet to be fully exploited. Over the last two decades, over a thousand QTL and GWAS studies have been published in the major cereal crops including sorghum, maize and rice. A large body of information has been generated on the genetic basis of quantitative traits, their genomic location, allelic effects and epistatic interactions. However, such QTL information has not been widely applied by cereal improvement programs and genetic researchers worldwide. In part this is due to the heterogeneous nature of QTL studies which leads QTL reliability variation from study to study. Using approaches to adjust the QTL confidence interval, this platform provides access to the most updated sorghum QTL information than any database available, spanning 23 years of research since 1995. The QTL database provides information on the predicted gene models underlying the QTL CI, across all sorghum genome assembly gene sets and maize and rice genome assemblies and also provides information on the diversity of the underlying genes and information on signatures of selection in sorghum. The resulting high-resolution, open-access research platform facilitates candidate gene identification across 3 cereal species, sorghum, maize and rice. Using a number of trait examples, we demonstrate the power and resolution of the resource to facilitate comparative genomics approaches to provide a bridge between genomics and applied breeding.

Genetic Map of Mango: A Tool for Mango Breeding.

Mango (Mangifera indica) is an economically and nutritionally important tropical/subtropical tree fruit crop. Most of the current commercial cultivars are selections rather than the products of breeding programs. To improve the efficiency of mango breeding, molecular markers have been used to create a consensus genetic map that identifies all 20 linkage groups in seven mapping populations. Polyembryony is an important mango trait, used for clonal propagation of cultivars and rootstocks. In polyembryonic mango cultivars, in addition to a zygotic embryo, several apomictic embryos develop from maternal tissue surrounding the fertilized egg cell. This trait has been associated with linkage group 8 in our consensus genetic map and has been validated in two of the seven mapping populations. In addition, we have observed a significant association between trait and single nucleotide polymorphism (SNP) markers for the vegetative trait of branch habit and the fruit traits of bloom, ground skin color, blush intensity, beak shape, and pulp color.

Sequence diversity and differential expression of major phenylpropanoid-flavonoid biosynthetic genes among three mango varieties.

BACKGROUND: Mango fruits contain a broad spectrum of phenolic compounds which impart potential health benefits; their biosynthesis is catalysed by enzymes in the phenylpropanoid-flavonoid (PF) pathway. The aim of this study was to reveal the variability in genes involved in the PF pathway in three different mango varieties Mangifera indica L., a member of the family Anacardiaceae: Kensington Pride (KP), Irwin (IW) and Nam Doc Mai (NDM) and to determine associations with gene expression and mango flavonoid profiles. RESULTS: A close evolutionary relationship between mango genes and those from the woody species poplar of the Salicaceae family (Populus trichocarpa) and grape of the Vitaceae family (Vitis vinifera), was revealed through phylogenetic analysis of PF pathway genes. We discovered 145 SNPs in total within coding sequences with an average frequency of one SNP every 316 bp. Variety IW had the highest SNP frequency (one SNP every 258 bp) while KP and NDM had similar frequencies (one SNP every 369 bp and 360 bp, respectively). The position in the PF pathway appeared to influence the extent of genetic diversity of the encoded enzymes. The entry point enzymes phenylalanine lyase (PAL), cinnamate 4-mono-oxygenase (C4H) and chalcone synthase (CHS) had low levels of SNP diversity in their coding sequences, whereas anthocyanidin reductase (ANR) showed the highest SNP frequency followed by flavonoid 3'-hydroxylase (F3'H). Quantitative PCR revealed characteristic patterns of gene expression that differed between mango peel and flesh, and between varieties. CONCLUSIONS: The combination of mango expressed sequence tags and availability of well-established reference PF biosynthetic genes from other plant species allowed the identification of coding sequences of genes that may lead to the formation of important flavonoid compounds in mango fruits and facilitated characterisation of single nucleotide polymorphisms between varieties. We discovered an association between the extent of sequence variation and position in the pathway for up-stream genes. The high expression of PAL, C4H and CHS genes in mango peel compared to flesh is associated with high amounts of total phenolic contents in peels, which suggest that these genes have an influence on total flavonoid levels in mango fruit peel and flesh. In addition, the particularly high expression levels of ANR in KP and NDM peels compared to IW peel and the significant accumulation of its product epicatechin gallate (ECG) in those extracts reflects the rate-limiting role of ANR on ECG biosynthesis in mango.

Complete genome sequence and intracellular protein localization of Datura yellow vein nucleorhabdovirus.

A limited number of plant rhabdovirus genomes have been fully sequenced, making taxonomic classification, evolutionary analysis and molecular characterization of this virus group difficult. We have for the first time determined the complete genome sequence of 13,188 nucleotides of Datura yellow vein nucleorhabdovirus (DYVV). DYVV genome organization resembles that of its closest relative, Sonchus yellow net virus (SYNV), with six ORFs in antigenomic orientation, separated by highly conserved intergenic regions and flanked by complementary 3' leader and 5' trailer sequences. As is typical for nucleorhabdoviruses, all viral proteins, except the glycoprotein, which is targeted to the endoplasmic reticulum, are localized to the nucleus. Nucleocapsid (N) protein, matrix (M) protein and polymerase, as components of nuclear viroplasms during replication, have predicted strong canonical nuclear localization signals, and N and M proteins exclusively localize to the nucleus when transiently expressed as GFP fusions. As in all nucleorhabdoviruses studied so far, N and phosphoprotein P interact when co-expressed, significantly increasing P nuclear localization in the presence of N protein. This research adds to the list of complete genomes of plant-infecting rhabdoviruses, provides molecular tools for further characterization and supports classification of DYVV as a nucleorhabdovirus closely related to but with some distinct differences from SYNV.

Grain sorghum proteomics: integrated approach toward characterization of endosperm storage proteins in kafirin allelic variants.

Grain protein composition determines quality traits, such as value for food, feedstock, and biomaterials uses. The major storage proteins in sorghum are the prolamins, known as kafirins. Located primarily on the periphery of the protein bodies surrounding starch, cysteine-rich ß- and γ-kafirins may limit enzymatic access to internally positioned α-kafirins and starch. An integrated approach was used to characterize sorghum with allelic variation at the kafirin loci to determine the effects of this genetic diversity on protein expression. Reversed-phase high performance liquid chromatography and lab-on-a-chip analysis showed reductions in alcohol-soluble protein in ß-kafirin null lines. Gel-based separation and liquid chromatography-tandem mass spectrometry identified a range of redox active proteins affecting storage protein biochemistry. Thioredoxin, involved in the processing of proteins at germination, has reported impacts on grain digestibility and was differentially expressed across genotypes. Thus, redox states of endosperm proteins, of which kafirins are a subset, could affect quality traits in addition to the expression of proteins.

A population of sexual Daphnia pulex resists invasion by asexual clones.

Asexual reproduction avoids the costs associated with sex, predicting that invading asexual clones can quickly replace sexual populations. Daphnia pulex populations in the Great Lakes area are predominately asexual, but the elimination of sexual populations by invading clones is poorly understood. Asexual clones were detected at low frequency in one rare sexual population in 1995, with some increase in frequency during 2003 and 2004. However, these clones remained at low frequency during further yearly sampling (2005-2013) with no evidence that the resident sexual population was in danger of elimination. There was evidence for hybridization between rare males produced by asexual clones and sexual females with the potential to produce new asexual genotypes and spread the genetic factors for asexuality. In a short-term laboratory competition experiment, the two most common asexual clones did not increase in frequency relative to a genetically diverse sexual population due in part to a greater investment in diapausing eggs that trades-off current population growth for increased contribution to the egg bank. Our results suggest that a successful invasion can be prolonged, requiring a combination of clonal genotypes with high fitness, persistence of clones in the egg bank and negative factors affecting the sexual population such as inbreeding depression resulting from population bottlenecks.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes.

BACKGROUND: Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL. RESULTS: In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL. CONCLUSIONS: NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the "arms race" with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.

Whole-genome sequencing reveals untapped genetic potential in Africa's indigenous cereal crop sorghum.

Sorghum is a food and feed cereal crop adapted to heat and drought and a staple for 500 million of the world's poorest people. Its small diploid genome and phenotypic diversity make it an ideal C4 grass model as a complement to C3 rice. Here we present high coverage (16-45 × ) resequenced genomes of 44 sorghum lines representing the primary gene pool and spanning dimensions of geographic origin, end-use and taxonomic group. We also report the first resequenced genome of S. propinquum, identifying 8 M high-quality SNPs, 1.9 M indels and specific gene loss and gain events in S. bicolor. We observe strong racial structure and a complex domestication history involving at least two distinct domestication events. These assembled genomes enable the leveraging of existing cereal functional genomics data against the novel diversity available in sorghum, providing an unmatched resource for the genetic improvement of sorghum and other grass species.

The role of hybridization in the origin and spread of asexuality in Daphnia.

The molecular mechanisms leading to asexuality remain little understood despite their substantial bearing on why sexual reproduction is dominant in nature. Here, we examine the role of hybridization in the origin and spread of obligate asexuality in Daphnia pulex, arguably the best-documented case of contagious asexuality. Obligately parthenogenetic (OP) clones of D. pulex have traditionally been separated into 'hybrid' (Ldh SF) and 'nonhybrid' (Ldh SS) forms because the lactase dehydrogenase (Ldh) locus distinguishes the cyclically parthenogenetic (CP) lake dwelling Daphnia pulicaria (Ldh FF) from its ephemeral pond dwelling sister species D. pulex (Ldh SS). The results of our population genetic analyses based on microsatellite loci suggest that both Ldh SS and SF OP individuals can originate from the crossing of CP female F1 (D. pulex × D. pulicaria) and backcross with males from OP lineages carrying genes that suppress meiosis specifically in female offspring. In previous studies, a suite of diagnostic markers was found to be associated with OP in Ldh SS D. pulex lineages. Our association mapping supports a similar genetic mechanism for the spread of obligate parthenogenesis in Ldh SF OP individuals. Interestingly, our study shows that CP D. pulicaria carry many of the diagnostic microsatellite alleles associated with obligate parthenogenesis. We argue that the assemblage of mutations that suppress meiosis and underlie obligate parthenogenesis in D. pulex originated due to a unique historical hybridization and introgression event between D. pulex and D. pulicaria.

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus).

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional-scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS-FLX Titanium platform. A total of 172,837 raw reads (17.4 Mbp) were retrieved, 95,500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. FST = 0.0025) can be resolved with a likely power ≥67% by using standard sample sizes of 50 and validated subsets of available markers.

Evolutionary factors affecting Lactate dehydrogenase A and B variation in the Daphnia pulex species complex.

BACKGROUND: Evidence for historical, demographic and selective factors affecting enzyme evolution can be obtained by examining nucleotide sequence variation in candidate genes such as Lactate dehydrogenase (Ldh). Two closely related Daphnia species can be distinguished by their electrophoretic Ldh genotype and habitat. Daphnia pulex populations are fixed for the S allele and inhabit temporary ponds, while D. pulicaria populations are fixed for the F allele and inhabit large stratified lakes. One locus is detected in most allozyme surveys, but genome sequencing has revealed two genes, LdhA and LdhB. RESULTS: We sequenced both Ldh genes from 70 isolates of these two species from North America to determine if the association between Ldh genotype and habitat shows evidence for selection, and to elucidate the evolutionary history of the two genes. We found that alleles in the pond-dwelling D. pulex and in the lake-dwelling D. pulicaria form distinct groups at both loci, and the substitution of Glutamine (S) for Glutamic acid (F) at amino acid 229 likely causes the electrophoretic mobility shift in the LDHA protein. Nucleotide diversity in both Ldh genes is much lower in D. pulicaria than in D. pulex. Moreover, the lack of spatial structuring of the variation in both genes over a wide geographic area is consistent with a recent demographic expansion of lake populations. Neutrality tests indicate that both genes are under purifying selection, but the intensity is much stronger on LdhA. CONCLUSIONS: Although lake-dwelling D. pulicaria hybridizes with the other lineages in the pulex species complex, it remains distinct ecologically and genetically. This ecological divergence, coupled with the intensity of purifying selection on LdhA and the strong association between its genotype and habitat, suggests that experimental studies would be useful to determine if variation in molecular function provides evidence that LDHA variants are adaptive.

Quantitative analysis of sperm plane circular movement in the blue mussels Mytilus edulis, M. trossulus and their hybrids.

Fertilization success of free spawning organisms such as Mytilus species depends on gamete interactions. Therefore, gamete traits such as sperm movement are important for determining fertilization success in free spawning organisms. Since little is known about sperm movement pattern in Mytilus species, the purpose of this study was to investigate sperm movement pattern of blue mussel M. edulis, M. trossulus and their hybrids using computer-assisted sperm movement video analysis. Sperm of all genotypes were found to conduct circular movement in a two-dimensional plane. Furthermore, new sperm movement parameters, real time radius (R), angle change rate (θ) and the center of circular track (O(t)) were developed to verify and quantitatively describe the plane circular movement pattern using software (Image-J) that may be widely applied to sperm movement study in other organisms. Angle change rate was positively correlated to fertilization success. However, no correlation between fertilization and real time radius was detected. Although no interspecific differences were found in the radius, the F¹ (first generation) hybrid sperm had a lower angle change rate than M. edulis and M. trossulus. Published studies have shown that sperm circular movement is more prevalent in aquatic broadcast spawning species than in species with mating behavior or internal fertilization. Therefore, a two-dimensional circular movement pattern in sperm may represent a trait that increases fertilization success for broadcast spawning species by either increasing gamete interaction rate at a small scale and/or avoiding swimming further away from the eggs before sperm detects the chemoattractant gradient.

Incomplete reproductive isolation in the blue mussel (Mytilus edulis and M. trossulus) hybrid zone in the Northwest Atlantic: role of gamete interactions and larval viability.

Reproductive isolation in free-spawning organisms may involve only small changes in the gamete surface molecules that control fertilization, linking gamete incompatibility and speciation. Most studies have focused on species in which natural hybrids are absent and reproductive isolation is complete, but how gamete incompatibility evolves remains unclear. Reproductive isolation is incomplete between two sympatric mussel species (Mytilus edulis, M. trossulus) that hybridize in nature. In this study prezygotic and postzygotic components of reproductive incompatibility were examined in laboratory crosses. Conspecific crosses showed significantly greater rates of fertilization than heterospecific crosses, although some females of both species showed heterospecific gamete compatibility. The proportion of fertilized eggs developing into normal larvae was not significantly different between conspecific and heterospecific crosses, but survival of normal larvae was greater for conspecific crosses. Mixed-species sperm experiments suggested that conspecific sperm preference may further limit hybridization. The different components of reproductive incompatibility and total incompatibility varied among females of both species. Although our study has shown that partial reproductive isolation between M. edulis and M. trossulus involves both prezygotic gamete interactions and postzygotic larval survival, further research is required to determine the potential role of gamete incompatibility in the evolution of complete reproductive isolation between these species.