Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (F_st) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). F_st was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of F_st values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.     

Leveraging whole genome sequencing data for demographic inference with approximate Bayesian computation

Accounting for historical demographic features, such as the strength and timing of gene flow and divergence times between closely related lineages, is vital for many inferences in evolutionary biology. Approximate Bayesian computation (ABC) is one method commonly used to estimate demographic parameters. However, the DNA sequences used as input for this method, often microsatellites or RADseq loci, usually represent a small fraction of the genome. Whole genome sequencing (WGS) data, on the other hand, have been used less often with ABC, and questions remain about the potential benefit of, and how to best implement, this type of data; we used pseudo‐observed data sets to explore such questions. Specifically, we addressed the potential improvements in parameter estimation accuracy that could be associated with WGS data in multiple contexts; namely, we quantified the effects of (a) more data, (b) haplotype‐based summary statistics, and (c) locus length. Compared with a hypothetical RADseq data set with 2.5 Mbp of data, using a 1 Gbp data set consisting of 100 Kbp sequences led to substantial gains in the accuracy of parameter estimates, which was mostly due to haplotype statistics and increased data. We also quantified the effects of including (a) locus‐specific recombination rates, and (b) background selection information in ABC analyses. Importantly, assuming uniform recombination or ignoring background selection had a negative effect on accuracy in many cases. Software and results from this method validation study should be useful for future demographic history analyses.     

Measuring dominance and diversity in ecological communities: choosing the right variables

Abstract. Although there are many studies in searching for the most useful community diversity index, the importance of choosing an appropriate parameter has been neglected. Here, we examined dominance and diversity in early post-fire chaparral communities using different variables, i.e. plant density, cover and biomass. Significant different results were produced by applying different parameters and the difference may be caused by the inconsistency in density, cover, and biomass allocated in each life form. Among the three parameters, biomass was most successful in detecting differences among communities because the apportionment of biomass among species was more variable than that of density. Although the three species variables represent different aspects of community properties and their relative performance may vary among communities, we recommend the use of biomass or productivity data as the most appropriate variable because it can best represent per capita resource use and resource partitioning among organisms in competitive situations.     

Mitogenomic analysis of Montipora cactus and Anacropora matthai (cnidaria; scleractinia; acroporidae) indicates an unequal rate of mitochondrial evolution among Acroporidae corals

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for specimens of the coral species Montipora cactus (Bernard 1897) and Anacropora matthai (Pillai 1973), representing two morphologically distinct genera of the family Acroporidae. These sequences were compared with the published mt genome sequence for the confamilial species, Acropora tenuis (Dana 1846). The size of the mt genome was 17,887 bp and 17,888 bp for M. cactus and A. matthai. Gene content and organization was found to be very similar among the three Acroporidae mt genomes with a group I intron occurring in the NADH dehyrogenase 5 (nad5) gene. The intergenic regions were also similar in length among the three corals. The control region located between the small ribosomal RNA (ms) and the cytochrome oxidase 3 (cox3) gene was significantly smaller in M. cactus and A. matthai (both 627 bp) than in A. tenuis (1086 bp). Only one set of repeated sequences was identified at the 3′-end of the control regions in M. cactus and A. matthai. A lack of the abundant repetitive elements which have been reported for A. tenuis, accounts for the relatively short control regions in M. cactus and A. matthai. Pairwise distances and relative rate analyses of 13 protein coding genes, the group I intron and the largest intergenic region, igr3, revealed significant differences in the rate of molecular evolution of the mt genome among the three species, with an extremely slow rate being seen between Montipora and Anacropora. It is concluded that rapid mt genome evolution is taking place in genus Acropora relative to the confamilial genera Montipora and Anacropora although all are within the relatively slow range thought to be typical of Anthozoa.     

Diet partitioning in sympatric Atlantic salmon and brown trout in streams with contrasting riparian vegetation

Prey intake by Atlantic salmon Salmo salar and brown trout Salmo trutta was measured across different riparian vegetation types: grassland, open canopy deciduous and closed canopy deciduous, in upland streams in County Mayo, Western Ireland. Fishes were collected by electrofishing while invertebrates were sampled from the benthos using a Surber sampler and drifting invertebrates collected in drift traps. Aquatic invertebrates dominated prey numbers in the diets of 0+ year Atlantic salmon and brown trout and 1+ year Atlantic salmon, whereas terrestrial invertebrates were of greater importance for diets of 1+ and 2+ year brown trout. Terrestrial prey biomass was generally greater than aquatic prey for 1+ and 2+ year brown trout across seasons and riparian types. Prey intake was greatest in spring and summer and least in autumn apart from 2+ year brown trout that sustained feeding into autumn. Total prey numbers captured tended to be greater for all age classes in streams with deciduous riparian canopy. Atlantic salmon consumed more aquatic prey and brown trout more terrestrial prey with an ontogenetic increase in prey species richness and diversity. Atlantic salmon and brown trout diets were most similar in summer. Terrestrial invertebrates provided an important energy subsidy particularly for brown trout. In grassland streams, each fish age class was strongly associated with aquatic, mainly benthic invertebrates. In streams with deciduous riparian canopy cover, diet composition partitioned between conspecifics with older brown trout associated with surface drifting terrestrial invertebrates and older Atlantic salmon associated with aquatic invertebrates with a high drift propensity in the water column and 0+ year fish feeding on benthic aquatic invertebrates. Deciduous riparian canopy cover may therefore facilitate vertical partitioning of feeding position within the water column between sympatric Atlantic salmon and brown trout. Implications for riparian management are discussed.     

Structure and expression of the overlapping ND4L and ND5 genes of Neurospora crassa mitochondria

Summary Genes homologous to the mammalian mitochondrial NADH dehydrogenase subunit genes ND4L and ND5 were identified in the mitochondrial genome of the filamentous fungus Neurospora crassa, and the structure and expression of these genes was examined. The ND4L gene (interrupted by one intervening sequence) potentially encodes an 89 residue long hydrophobic protein that shares about 26% homology (or 41% homology if conservative amino acid substitutions are allowed) with the analogous human mitochondrial protein. The ND5 gene (which contains two introns) encodes a 715 residue polypeptide that shares 23% homology with the human analogue; a 300 amino acid long region is highly conserved (50% homology) in the two ND5 proteins. The stop codon of the ND4L gene overlaps the initiation codon of the downstream ND5 gene, and the two genes are contranscribed and probably cotranslated. A presumed mature dicistronic (ND4L plus ND5) RNA was detected. The postulated mRNA (about 3.2 kb) contains 5 and 3 non-coding regions of about 86 and 730 nucleotides, respectively; this species is generated from very large precursor RNAs by a complex processing pathway. The ND4L and ND5 introns are all stable after their excision from the precursor species.Abbreviations bp base pairs - rRNA ribosomal RNA - ND NADH dehydrogenase - URF unidentified reading frame - kDal kilodaltons; a.a., amino acid     

Organization of Methanobacterium thermoautotrophicum bacteriophage psi M1 DNA

Summary M1 is a virulent bacteriophage of Methanobacterium thermoautotrophicum strain Marburg. Restriction enzyme analysis of the linear, 30.4 kb phage DNA led to a circular map of the 27.1 kb M1 genome. M1 is thus circularly permuted and exhibits terminal redundancy of approximately 3 kb. Packaging of M1 DNA from a concatemeric precursor initiates at the pac site which was identified at coordinate 4.6 kb on the circular genome map. It proceeds clockwise for at least five packaging rounds. Headful packaging was also shown for M2, a phage variant with a 0.7 kb deletion at coordinate 23.25 on the map.     

Effects of microbial floras on the distributions of five domestic Drosophila species across fruit resources

Summary The distributions of five Drosophila species and four components of the microflora have been compared across a total of 48 traps baited with four different fruit and vegetable substrates in two domestic compost heaps in Canberra (Australia). Large and consistent differences are found, both among the Drosophila and among the microbial classes, in their distributions across traps baited with different substrates. Moreover the distribution of each Drosophila species shows a unique set of strong associations with the microbial distributions. Thus the distributions of both D. simulans and D. melanogaster are found to be strongly negatively correlated with the abundance of bacteria while D. simulans is also strongly positively correlated with the titre of fermenter yeasts. D. immigrans is strongly positively correlated both with bacteria and with non-fermenter yeasts. D. hydei is positively correlated with nonfermentery yeasts and D. busckii is negatively correlated with fermenter yeasts. Moulds are the only microbial class not consistently associated with the distribution of any of the Drosophila species. The correlations with the other microbial classes are sufficient to explain the majority of the abundance differences of the Drosophila species among the trap types. It is therefore proposed that the clear partitioning of the fruit resources by the Drosophila is due to their differing primary interactions with the microflora.