Search

Genetic screens have been used extensively to probe interactions between nuclear genes and their impact on phenotypes. Probing interactions between mitochondrial genes and their phenotypic outcome, however, has not been possible due to a lack of tools to map the responsible polymorphisms. Here, using a toolkit we previously established in Drosophila, we isolate over 300 recombinant mitochondrial genomes and map a naturally occurring polymorphism at the cytochrome c oxidase III residue 109 (CoIII109) that fully rescues the lethality and other defects associated with a point mutation in cytochrome c oxidase I (CoIT300I). Through lipidomics profiling, biochemical assays and phenotypic analyses, we show that the CoIII109 polymorphism modulates cardiolipin binding to prevent complex IV instability caused by the CoIT300I mutation. This study demonstrates the feasibility of genetic interaction screens in animal mitochondrial DNA. It unwraps the complex intra-genomic interplays underlying disorders linked to mitochondrial DNA and how they influence disease expression. Chiang et al. map a genetic interaction in animal mitochondrial DNA by recombination. This reveals how polymorphisms in two complex IV subunits jointly affect cardiolipin binding to impact complex stability, organismal fitness and disease expression. [ABSTRACT FROM AUTHOR]

We and many others have noted the advantages of using heterogeneous (HS) animals to map genes and gene networks associated with both behavioral and non-behavioral phenotypes. Importantly, genetically complex Mus musculus crosses provide substantially increased resolution to examine old and new relationships between gene expression and behavior. Here we report on data obtained from two HS populations: the HS/NPT derived from eight inbred laboratory mouse strains and the HS-CC derived from the eight collaborative cross inbred mouse strains that includes three wild-derived strains. Our work has focused on the genes and gene networks associated with risk for excessive ethanol consumption, individual variation in ethanol consumption and the consequences, including escalation, of long-term ethanol consumption. Background data on the development of HS mice is provided, including advantages for the detection of expression quantitative trait loci. Examples are also provided of using HS animals to probe the genes associated with ethanol preference and binge ethanol consumption. [ABSTRACT FROM AUTHOR]

Animal models for cystic fibrosis (CF) have enhanced our understanding of the pathology and contributed to the development of new treatments. In the field of CF, many animal models have been developed and described. To our knowledge, thus far, none of the reviews of CF animal models has used a systematic methodology. A systematic approach to creating model overviews can lead to an objective, evidence-based choice of an animal model for new research questions. We searched Pubmed and Embase for the currently available animal models for CF. Two independent reviewers screened the results. We included all primary studies describing an animal model for CF. After duplicate removal, 12,304 publications were left. Because of the large number of models, in the current paper, only the genetic models are presented. A total of 636 publications were identified describing genetic animal models for CF in mice, pigs, ferrets, rats and zebrafish. Most of these models have an altered Cftr gene. An overview of basic model characteristics and outcome measures for these genetic models is provided, together with advice on using these data. As far as the authors are aware, this is one of the largest systematic mapping reviews on genetic animal models for CF. It can aid in selecting a suitable model and outcome measures. In general, the reporting quality of the included publications was poor. Further systematic reviews are warranted to determine the quality and translational value of these models further. [ABSTRACT FROM AUTHOR]

An unusual supernumerary chromosome has been reported for two related avian species, the zebra and Bengalese finches. This large, germline-restricted chromosome (GRC) is eliminated from somatic cells and spermatids and transmitted via oocytes only. Its origin, distribution among avian lineages, and function were mostly unknown so far. Using immunolocalization of key meiotic proteins, we found that GRCs of varying size and genetic content are present in all 16 songbird species investigated and absent from germline genomes of all eight examined bird species from other avian orders. Results of fluorescent in situ hybridization of microdissected GRC probes and their sequencing indicate that GRCs show little homology between songbird species and contain a variety of repetitive elements and unique sequences with paralogs in the somatic genome. Our data suggest that the GRC evolved in the common ancestor of all songbirds and underwent significant changes in the extant descendant lineages. [ABSTRACT FROM AUTHOR]

Since William Bateson's 1906 coinage of the term "genetics," the rise of mice as a model organism for biomedical research has gone hand in hand with genomic developments. In today's research environment, mice and rats make up the vast majority of all research subjects. While the advent of gene-editing tools such as CRISPR has made genetic manipulation of mice easier, these tools also signal a new trend toward an increased use of largeanimal models such as dogs, pigs, and nonhuman primates. Especially for neurological impairments, CRISPR gene editing offers the potential to generate large-animal models that better mimic human diseases. What are the ethical and regulatory implications of this trend? The professional and ethical framework for responsible conduct of animal research is widely recognized as the "three Rs": Reduction, Refinement, and Replacement. This Article points to the tension between reduction (decreasing the overall numbers of animals used) and relative replacement (the use of mice and rats instead of species with more "complex" capacities) that is implied by such a tradeoff. The Article argues, however, that a comparative analysis of regulatory frameworks in the United States and in the European Union shows that neither offers any substantial guidance to moderate a trend toward greater use of large-animal models. Nevertheless, we raise several ethical questions associated with the trend toward using relatively fewer animals but replacing less cognitively developed animals with those with potentially greater morally relevant capacities and moral standing. [ABSTRACT FROM AUTHOR]

Accuracy of genome-wide association studies, and the successful implementation of genomic selection depends on the level of linkage disequilibrium (LD) across the genome and also the persistence of LD phase between populations. In the present study LD between adjacent SNPs and LD decay between SNPs was calculated in three Iranian water buffalo populations. Persistence of LD phase was evaluated across these populations and effective population size (Ne) was estimated from corrected r2 information. A set of 404 individuals from three Iranian buffalo populations were genotyped with the Axiom Buffalo Genotyping 90K Array. Average r2 and |D'| between adjacent SNP pairs across all chromosomes was 0.27 and 0.66 for AZI, 0.29 and 0.68 for KHU, and 0.32 and 0.72 for MAZ. The LD between the SNPs decreased with increasing physical distance from 100Kb to 1Mb between markers, from 0.234 to 0.018 for AZI, 0.254 to 0.034 for KHU, and 0.297 to 0.119 for MAZ, respectively. These results indicate that a density of 90K SNP is sufficient for genomic analyses relying on long range LD (e.g. GWAS and genomic selection). The persistence of LD phase decreased with increasing marker distances across all the populations, but remained above 0.8 for AZI and KHU for marker distances up to 100Kb. For multi-breed genomic evaluation, the 90K SNP panel is suitable for AZI and KHU buffalo breeds. Estimated effective population sizes for AZI, KHU and MAZ were 477, 212 and 32, respectively, for recent generations. The estimated effective population sizes indicate that the MAZ is at risk and requires careful management. [ABSTRACT FROM AUTHOR]

Objective: Feed consumption contributes a large percentage for total production costs in the poultry industry. Detecting genes associated with feeding traits will be of benefit to improve our understanding of the molecular determinants for feed efficiency. The objective of this study was to identify candidate genes associated with feed conversion ratio (FCR) via genomewide association study (GWAS) using sequence data imputed from single nucleotide polymorphism (SNP) panel in a Chinese indigenous chicken population. Methods: A total of 435 Chinese indigenous chickens were phenotyped for FCR and were genotyped using a 600K SNP genotyping array. Twenty-four birds were selected for sequencing, and the 600K SNP panel data were imputed to whole sequence data with the 24 birds as the reference. The GWAS were performed with GEMMA software. Results: After quality control, 8,626,020 SNPs were used for sequence based GWAS, in which ten significant genomic regions were detected to be associated with FCR. Ten candidate genes, ubiquitin specific peptidase 44, leukotriene A4 hydrolase, ETS transcription factor, R-spondin 2, inhibitor of apoptosis protein 3, sosondowah ankyrin repeat domain family member D, calmodulin regulated spectrin associated protein family member 2, zinc finger and BTB domain containing 41, potassium sodium-activated channel subfamily T member 2, and member of RAS oncogene family were annotated. Several of them were within or near the reported FCR quantitative trait loci, and others were newly reported. Conclusion: Results from this study provide valuable prior information on chicken genomic breeding programs, and potentially improve our understanding of the molecular mechanism for feeding traits. [ABSTRACT FROM AUTHOR]

Objective: Internal organs indirectly affect economic performance and well-being of animals. Study of internal organs during later layer period will allow full utilization of layer hens. Hence, we conducted a genome-wide association study (GWAS) to identify potential quantitative trait loci or genes that potentially contribute to internal organ weight. Methods: A total of 1,512 chickens originating from White Leghorn and Dongxiang Blue- Shelled chickens were genotyped using high-density Affymetrix 600 K single nucleotide polymorphism (SNP) array. We conducted a GWAS, linkage disequilibrium analysis, and heritability estimated based on SNP information by using GEMMA, Haploview and GCTA software. Results: Our results displayed that internal organ weights show moderate to high (0.283 to 0.640) heritability. Variance partitioned across chromosomes and chromosome lengths had a linear relationship for liver weight and gizzard weight (R² = 0.493, 0.753). A total of 23 highly significant SNPs that associated with all internal organ weights were mainly located on Gallus gallus autosome (GGA) 1 and GGA4. Six SNPs on GGA2 affected heart weight. After the final analysis, five top SNPs were in or near genes 5-Hydroxytryptamine receptor 2A, general transcription factor IIF polypeptide 2, WD repeat and FYVE domain containing 2, non-SMC condensin I complex subunit G, and sonic hedgehog, which were considered as candidate genes having a pervasive role in internal organ weights. Conclusion: Our findings provide an understanding of the underlying genetic architecture of internal organs and are beneficial in the selection of chickens. [ABSTRACT FROM AUTHOR]

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock. [ABSTRACT FROM AUTHOR]

Abstract Copy number variation (CNV) related to complex traits, such as disease and quantitative phenotype, is considered an important and wealthy source of genetic and phenotypic diversity. It suggests that the copy number variation of function gene maybe leads to the phenotypic changes. Kupple like factor 3 (KLF3) gene is a vital transcription factor associated with the growth and development of muscle and adipose tissue. It has been mapped in a CNV region by animal genome re-sequencing. In this study, we detected the distribution diversity of KLF3 gene copy numbers in six Chinese cattle breeds (QC, NY, XN, PN, QDM and JX) and associated the phenotypic traits with it. Then, we analyzed the KLF3 gene transcription expression level in different tissues of Jiaxian (JX) cattle. Furthermore, we detected mRNA expression level of muscle and fat tissues of Jiaxian cattle (JX), Angus × Jiaxian (AJ). The results showed that the copy number in CNV loss was more frequent in QC than others. And we revealed a positive effect of KLF3 CNV on growth traits, such as body mass and heart girth (P < 0.05). In a word, we ascertained the significance between CNVs of KLF3 gene and growth traits in different cattle breeds, and our data indicates that the CNVs of KLF3 gene may as a marker for the future molecular breeding of Chinese beef cattle. Highlights • The CNV region was identified in Chinese cattle of KLF3 gene. • The KLF3 gene was higher in fat and liver than other tissues. • The KLF3 gene expression level was discrepant in same tissue of different cattle breeds. • There were significant differences in individual growth traits of different CNV types. [ABSTRACT FROM AUTHOR]

Genetic approaches have proved valuable to the study and conservation of endangered populations, especially for monitoring programs, and there is potential for further developments in this direction by extending analyses to the genomic level. We assembled the genome of the wolverine (Gulo gulo), a mustelid that in Scandinavia has recently recovered from a significant population decline, and obtained a 2.42 Gb draft sequence representing >85% of the genome and including >21,000 protein‐coding genes. We then performed whole‐genome resequencing of 10 Scandinavian wolverines for population genomic and demographic analyses. Genetic diversity was among the lowest detected in a red‐listed population (mean genome‐wide nucleotide diversity of 0.05%). Results of the demographic analyses indicated a long‐term decline of the effective population size (Ne) from 10,000 well before the last glaciation to <500 after this period. Current Ne appeared even lower. The genome‐wide FIS level was 0.089 (possibly signaling inbreeding), but this effect was not observed when analyzing a set of highly variable SNP markers, illustrating that such markers can give a biased picture of the overall character of genetic diversity. We found significant population structure, which has implications for population connectivity and conservation. We used an integrated microfluidic circuit chip technology to develop an SNP‐array consisting of 96 highly informative markers that, together with a multiplex pre‐amplification step, was successfully applied to low‐quality DNA from scat samples. Our findings will inform management, conservation, and genetic monitoring of wolverines and serve as a genomic roadmap that can be applied to other endangered species. The approach used here can be generally utilized in other systems, but we acknowledge the trade‐off between investing in genomic resources and direct conservation actions. Article impact statement: There is a direct link between genomics and management of a red‐listed population of wolverines. [ABSTRACT FROM AUTHOR]

Whiteflies are among the world’s most significant agricultural pests and chemical insecticides are extensively used to reduce crop damage to acceptable levels. However, nearly all insecticides pose a threat to the environment and alternative control methods, such as breeding of crop varieties that are inherently insect-resistant, are needed. Previously, a strong source of plant-age dependent resistance to the cabbage whitefly (Aleyrodes proletella) has been identified in the modern white cabbage (Brassica oleracea var. capitata) variety Rivera. However, nothing is known about the molecular mechanisms or the genes involved in this resistance. In the present study, a multidisciplinary approach combining transcriptome and metabolome profiling with genetic mapping was used to identify the molecular players of whitefly resistance in cabbage. Transcriptome profiles of young (susceptible) and older (resistant) Rivera plants were analyzed using RNA sequencing. While many genes involved in general processes were differentially expressed between both ages, several defense-related processes were overrepresented in the transcriptome profile of older plants. Hormone measurements revealed that jasmonic acid (JA) levels decreased upon whitefly infestation at both plant ages. Interestingly, abscisic acid (ABA) levels showed contrasting effects in response to whitefly infestation: ABA levels were reduced in young plants but induced in older plants upon whitefly feeding. Auxin levels were significantly lower in older plants compared with young plants, independent of whitefly presence, while glucosinolate levels were higher. Additionally, whitefly performance was monitored in an F2 population derived from a cross between Rivera and the susceptible white cabbage variety Christmas Drumhead. Significant QTL intervals were mapped on chromosome 2 and 9 for oviposition rate and whitefly adult survival, respectively. Several genes that were higher expressed in older plants and located in the identified QTL intervals were orthologous to Arabidopsis genes that have been related to ABA signaling, suggesting a role for ABA in the regulation of resistance towards whiteflies. Our results show that combining different omics approaches is a useful strategy to identify candidate genes underlying insect resistance. [ABSTRACT FROM AUTHOR]

Single-step genomic prediction models utilizing both genotyped and nongenotyped animals are likely to become the prevailing tool in genetic evaluations of livestock. Various single-step prediction models have been proposed, based either on estimation of individual marker effects or on direct prediction via a genomic relationship matrix. In this study, a classical pedigree-based animal model, a regular single-step genomic BLUP (ssGBLUP) model, algorithm for proven and young (APY) with 2 strategies for choosing core animals, and a single-step Bayesian regression (ssBR) model were compared for 305-d production traits (milk, fat, protein) in the Finnish red dairy cattle population. A residual polygenic effect with 10% of total genetic variance was included in the single-step models to reduce inflation of genomic predictions. Validation reliability was calculated as the squared Pearson correlation coefficient between genomically enhanced breeding value (GEBV) and yield deviation for masked records for 2,056 validation cows from the last year in the data set investigated. The results showed that gains of 0.02 to 0.04 on validation reliability were achieved by using single-step methods compared with the classical animal model. The regular ssGBLUP model and ssBR model with an extra polygenic effect yielded the same results. The APY methods yielded similar reliabilities as the regular ssGBLUP and ssBR. Exact prediction error variance of GEBV could be obtained by ssBR to avoid any approximation methods used for ssGBLUP when inversion left-hand side of mixed model equations is computationally infeasible for large data sets. [ABSTRACT FROM AUTHOR]

Abstract This study aims to improve the accuracy and consistency of T 1 and T 2 measurements using cardiac MR Fingerprinting (cMRF) by investigating and accounting for the effects of confounding factors including slice profile, inversion and T 2 preparation pulse efficiency, and B 1 +. The goal is to understand how measurements with different pulse sequences are affected by these factors. This can be used to determine which factors must be taken into account for accurate measurements, and which may be mitigated by the selection of an appropriate pulse sequence. Simulations were performed using a numerical cardiac phantom to assess the accuracy of over 600 cMRF sequences with different flip angles, TRs, and preparation pulses. A subset of sequences, including one with the lowest errors in T 1 and T 2 maps, was used in subsequent analyses. Errors due to non-ideal slice profile, preparation pulse efficiency, and B 1 + were quantified in Bloch simulations. Corrections for these effects were included in the dictionary generation and demonstrated in phantom and in vivo cardiac imaging at 3 T. Neglecting to model slice profile and preparation pulse efficiency led to underestimated T 1 and overestimated T 2 for most cMRF sequences. Sequences with smaller maximum flip angles were less affected by slice profile and B 1 +. Simulating all corrections in the dictionary improved the accuracy of T 1 and T 2 phantom measurements, regardless of acquisition pattern. More consistent myocardial T 1 and T 2 values were measured using different sequences after corrections. Based on these results, a pulse sequence which is minimally affected by confounding factors can be selected, and the appropriate residual corrections included for robust T 1 and T 2 mapping. [ABSTRACT FROM AUTHOR]

Population growth models are integral to ecological studies by providing estimates of population performance across space and time. Several models have been developed that estimate population growth through correlates of demographic traits, as measuring each parameter of the model can be prohibitive in experimental studies. Since differences in female size can accurately reflect changes in fecundity for many taxa, Livdahl and Sugihara developed a population growth index that incorporates size-fecundity relationships as a proxy for fecundity. To investigate the extent to which this model is robust to variation of this proxy, we tested if genetic (source population), temperature and resource treatments affect the size-fecundity relationship in Aedes albopictus (Skuse), the Asian tiger mosquito. We then determined if variation in the size-fecundity relationship alters the population growth estimates, lambda (λ’), when applied to Livdahl and Sugihara’s model. We performed 2 laboratory experiments in which we reared cohorts of four different geographic populations of A. albopictus across 5 temperature treatments (18, 21, 25, 18, 31°C) and three resource treatments (low, medium, high larval resources). We determined if the slope of the size-fecundity relationship varied by source population, temperature, or resource; and if variation in this relationship affects lambda (λ’) estimates in a competition study between A. albopictus and Culex pipiens (Linnaeus), the northern house mosquito. Temperature treatments significantly affected the size-fecundity relationship, resource level marginally affected the relationship, while source population had no effect. We found positive relationships between size and fecundity when mosquito larvae were reared at high temperatures and low resource levels but the relationship disappeared when mosquitoes were reared at a low temperature or with high levels of resources. The variation in the size-fecundity relationship produced from different temperatures resulted in statistically different lambda (λ’) estimates. However, these changes in lambda (λ’) did not alter the trends in the population performance across treatments or conclusions of the competition study. This study provides evidence that the population growth model is sensitive to variation in size-fecundity relationships and we recommend biologists apply the most compatible size-fecundity relationship to the models to obtain the most accurate estimates of population performance. [ABSTRACT FROM AUTHOR]

Landscape genetics investigations examine how the availability and configuration of habitat influence genetic structure of plants and animals. We used landscape genetics to evaluate the role that forest connectivity plays in determining genetic structure of the federally-threatened Northern Spotted Owl (Strix occidentalis caurina) using genotypes of 339 Northern Spotted Owls obtained for 10 microsatellite loci. Spatial clustering analyses identified a distinct genetic cluster at the southern extent of the region examined. This cluster could not be linked to landscape connectivity patterns and suggested that post-Pleistocene processes were involved with its development rather than contemporary landscape configuration. We also compared matrices of pairwise inter-individual genetic distances with resistance distances derived from a circuit-theory based framework. Resistance distances were obtained for an idealized raster map that reflected continuous unimpeded dispersal habitat across the landscape along with five empirically-derived raster maps reflecting the 1870’s, 1940’s, 1986, 1994, and 2012. Resistance distances from the idealized map served as surrogates for linear geographic distances. Relative to idealized conditions, resistance distances were ~250% higher in the 1940’s and ~200% higher from 1986 onward. Resistance distances from the 1870’s were ~40% higher than idealized conditions. Inter-individual genetic distances were most highly correlated with resistance distances from the idealized map rather than any of the empirical maps. Two hypotheses explain our results. First, our results may reflect temporal lags between the onset of large-scale habitat alterations and their novel effects on genetic structure in long-lived species such as Northern Spotted Owls. Second, because Northern Spotted Owls disperse over long distances, our results may indicate that forest habitat has never been sufficiently fragmented to the point where connectivity was disrupted. The second hypothesis could indicate that forest management practices mandated by the Northwest Forest Plan succeeded with one of its primary goals. However, our results do not represent a complete portrayal of the status of Northern Spotted Owls given detection of significant population declines and bottlenecks in other studies. Future investigations based on computer simulations may help distinguish between hypotheses. [ABSTRACT FROM AUTHOR]

Abstract: Dobzhansky–Muller (DM) incompatibilities involving sex chromosomes have been proposed to account for Haldane's rule (lowered fitness among hybrid offspring of the heterogametic sex) as well as Darwin's corollary (asymmetric fitness costs with respect to the direction of the cross). We performed simulation studies of a hybrid zone to investigate the effects of different types of DM incompatibilities on cline widths and positions of sex‐linked markers. From our simulations, X‐Y incompatibilities generate steep clines for both X‐linked and Y‐linked markers; random effects may produce strong noise in cline center positions when migration is high relative to fitness costs, but X‐ and Y‐centers always coincide strictly. X‐autosome and Y‐autosome incompatibilities also generate steep clines, but systematic shifts in cline centers occur when migration is high relative to selection, as a result of a dominance drive linked to Darwin's corollary. Interestingly, sex‐linked genes always show farther introgression than the associated autosomal genes. We discuss ways of disentangling the potentially confounding effects of sex biases in migration, we compare our results to those of a few documented contact zones, and we stress the need to study independent replicates of the same contact zone. [ABSTRACT FROM AUTHOR]

Background: Identification of genomic regions that have been targets of selection may shed light on the genetic history of livestock populations and help to identify variation controlling commercially important phenotypes. The Azeri and Kuzestani buffalos are the most common indigenous Iranian breeds which have been subjected to divergent selection and are well adapted to completely different regions. Examining the genetic structure of these populations may identify genomic regions associated with adaptation to the different environments and production goals. Results: A set of 385 water buffalo samples from Azeri (N = 262) and Khuzestani (N = 123) breeds were genotyped using the Axiom® Buffalo Genotyping 90 K Array. The unbiased fixation index method (FST) was used to detect signatures of selection. In total, 13 regions with outlier FST values (0.1%) were identified. Annotation of these regions using the UMD3.1 Bos taurus Genome Assembly was performed to find putative candidate genes and QTLs within the selected regions. Putative candidate genes identified include FBXO9, NDFIP1, ACTR3, ARHGAP26, SERPINF2, BOLA-DRB3, BOLA-DQB, CLN8, and MYOM2. Conclusions: Candidate genes identified in regions potentially under selection were associated with physiological pathways including milk production, cytoskeleton organization, growth, metabolic function, apoptosis and domestication-related changes include immune and nervous system development. The QTL identified are involved in economically important traits in buffalo related to milk composition, udder structure, somatic cell count, meat quality, and carcass and body weight. [ABSTRACT FROM AUTHOR]

In this study, data genotyping by sequence (GBS) was used to perform single step GWAS (ssGWAS) to identify SNPs associated with the litter traits in domestic pigs and search for candidate genes in the region of significant SNPs. After quality control, 167,355 high-quality SNPs from 532 pigs were obtained. Phenotypic traits on 2112 gilt litters from 532 pigs were recorded including total number born (TNB), number born alive (NBA), and litter weight born alive (LWB). A single-step genomic BLUP approach (ssGBLUP) was used to implement the genome-wide association analysis at a 5% genome-wide significance level. A total of 8, 23 and 20 significant SNPs were associated with TNB, NBA, and LWB, respectively, and these significant SNPs accounted for 62.78%, 79.75%, and 58.79% of genetic variance. Furthermore, 1 (SSC14: 16314857), 4 (SSC1: 81986236, SSC1: 66599775, SSC1: 161999013, and SSC1: 267883107), and 5 (SSC9: 29030061, SSC2: 32368561, SSC5: 110375350, SSC13: 45619882 and SSC13: 45647829) significant SNPs for TNB, NBA, and LWB were inferred to be novel loci. At SSC1, the AIM1 and FOXO3 genes were found to be associated with NBA; these genes increase ovarian reproductive capacity and follicle number and decrease gonadotropin levels. The genes SLC36A4 and INTU are involved in cell growth, cytogenesis and development were found to be associated with LWB. These significant SNPs can be used as an indication for regions in the Sus scrofa genome for variability in litter traits, but further studies are expected to confirm causative mutations. [ABSTRACT FROM AUTHOR]

Mycoplasma bovirhinis is a significant etiology in bovine pneumonia and mastitis, but our knowledge about the genetic and pathogenic mechanisms of M. bovirhinis is very limited. In this study, we sequenced the complete genome of M. bovirhinis strain GS01 isolated from the nasal swab of pneumonic calves in Gansu, China, and we found that its genome forms a 847,985 bp single circular chromosome with a GC content of 27.57% and with 707 protein-coding genes. The putative virulence determinants of M. bovirhinis were then analyzed. Results showed that three genomic islands and 16 putative virulence genes, including one adhesion gene enolase, seven surface lipoproteins, proteins involved in glycerol metabolism, and cation transporters, might be potential virulence factors. Glycerol and pyruvate metabolic pathways were defective. Comparative analysis revealed remarkable genome variations between GS01 and a recently reported HAZ141_2 strain, and extremely low homology with others mycoplasma species. Phylogenetic analysis demonstrated that M. bovirhinis was most genetically close to M. canis, distant from other bovine Mycoplasma species. Genomic dissection may provide useful information on the pathogenic mechanisms and genetics of M. bovirhinis. [ABSTRACT FROM AUTHOR]

Background: Malaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes. Results: In this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes. Conclusions: Our study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations. [ABSTRACT FROM AUTHOR]

In this study, we report evidence of a novel duplication of Melanocortin receptor 1 ( Mc1r) in the cavefish genome. This locus was discovered following the observation of excessive allelic diversity in a ∼820 bp fragment of Mc1r amplified via degenerate PCR from a natural population of Astyanax aeneus fish from Guerrero, Mexico. The cavefish genome reveals the presence of two closely related Mc1r open reading frames separated by a 1.46 kb intergenic region. One open reading frame corresponds to the previously reported Mc1r receptor, and the other open reading frame (duplicate copy) is 975 bp in length, encoding a receptor of 325 amino acids. Sequence similarity analyses position both copies in the syntenic region of the single Mc1r locus in 16 representative craniate genomes spanning bony fish (including Astyanax) to mammals, suggesting we discovered tandem duplicates of this important gene. The two Mc1r copies share ∼89% sequence similarity and, within Astyanax, are more similar to one another compared to other melanocortin family members. Future studies will inform the precise functional significance of the duplicated Mc1r locus and if this novel copy number variant may have adaptive significance for the Astyanax lineage. [ABSTRACT FROM AUTHOR]

The mangrove rivulus, Kryptolebias marmoratus, is one of only two self-fertilizing hermaphroditic fish species and inhabits mangrove forests. While selfing can be advantageous, it reduces heterozygosity and decreases genetic diversity. Studies using microsatellites found that there are variable levels of selfing among populations of K. marmoratus, but overall, there is a low rate of outcrossing and, therefore, low heterozygosity. In this study, we used whole-genome data to assess the levels of heterozygosity in different lineages of the mangrove rivulus and infer the phylogenetic relationships among those lineages. We sequenced whole genomes from 15 lineages that were completely homozygous at microsatellite loci and used single nucleotide polymorphisms (SNPs) to determine heterozygosity levels. More variation was uncovered than in studies using microsatellite data because of the resolution of full genome sequencing data. Moreover, missense polymorphisms were found most often in genes associated with immune function and reproduction. Inferred phylogenetic relationships suggest that lineages largely group by their geographic distribution. The use of whole-genome data provided further insight into genetic diversity in this unique species. Although this study was limited by the number of lineages that were available, these data suggest that there is previously undescribed variation within lineages of K. marmoratus that could have functional consequences and (or) inform us about the limits to selfing (e.g., genetic load, accumulation of deleterious mutations) and selection that might favor the maintenance of heterozygosity. These results highlight the need to sequence additional individuals within and among lineages. [ABSTRACT FROM AUTHOR]

A revision of mapping of main and alternative banding sequences in chromosome III (EF) has been made for 14 species of the Chironomus plumosus group. In total, new versions of mapping are presented for 18 banding sequences of arm E and 18 banding sequences of arm F. A new way of tracing the origins of banding sequences in chromosome III of the Ch. plumosus group in comparison with basic banding sequences of the genus Chironomus is suggested. The presented data indicate that h'pluE2 in arm E and p'borF2 in arm F are the closest to banding sequences of Ch. piger Strenzke, 1959 and thus should be considered the most ancient among banding sequences of chromosome III in the Ch. plumosus group. Phylogenetic relationships of banding sequences of chromosome III are discussed. [ABSTRACT FROM AUTHOR]

Background: The forest musk deer, Moschus berezovskii, is one of seven musk deer (Moschus spp.) and is distributed in Southwest China. Akin to other musk deer, the forest musk deer has been traditionally and is currently hunted for its musk (i.e., global perfume industry). Considerable hunting pressure and habitat loss have caused significant population declines. Consequently, the Chinese government commenced captive breeding programs for musk harvesting in the 1950s. However, the prevalence of fatal diseases is considerably restricting population increases. Disease severity and extent are exacerbated by inbreeding and genetic diversity declines in captive musk deer populations. It is essential that knowledge of captive and wild forest musk deer populations' immune system and genome be gained in order to improve their physical and genetic health. We have thus sequenced the whole genome of the forest musk deer, completed the genomic assembly and annotation, and performed preliminary bioinformatic analyses. Findings: A total of 407 Gb raw reads from whole-genome sequencing were generated using the Illumina HiSeq 4000 platform. The final genome assembly is around 2.72 Gb, with a contig N50 length of 22.6 kb and a scaffold N50 length of 2.85 Mb. We identified 24,352 genes and found that 42.05% of the genome is composed of repetitive elements. We also detected 1,236 olfactory receptor genes. The genome-wide phylogenetic tree indicated that the forest musk deer was within the order Artiodactyla, and it appeared as the sister clade of four members of Bovidae. In total, 576 genes were under positive selection in the forest musk deer lineage. Conclusions: We provide the first genome sequence and gene annotation for the forest musk deer. The availability of these resources will be very useful for the conservation and captive breeding of this endangered and economically important species and for reconstructing the evolutionary history of the order Artiodactyla. [ABSTRACT FROM AUTHOR]

Background: The Peruvian scallop, Argopecten purpuratus, is mainly cultured in southern Chile and Peru was introduced into China in the last century. Unlike other Argopecten scallops, the Peruvian scallop normally has a long life span of up to 7 to 10 years. Therefore, researchers have been using it to develop hybrid vigor. Here, we performed whole genome sequencing, assembly, and gene annotation of the Peruvian scallop, with an important aim to develop genomic resources for genetic breeding in scallops. Findings: A total of 463.19-Gb raw DNA reads were sequenced. A draft genome assembly of 724.78 Mb was generated (accounting for 81.87% of the estimated genome size of 885.29 Mb), with a contig N50 size of 80.11 kb and a scaffold N50 size of 1.02 Mb. Repeat sequences were calculated to reach 33.74% of the whole genome, and 26,256 protein-coding genes and 3,057 noncoding RNAs were predicted from the assembly. Conclusions: We generated a high-quality draft genome assembly of the Peruvian scallop, which will provide a solid resource for further genetic breeding and for the analysis of the evolutionary history of this economically important scallop. [ABSTRACT FROM AUTHOR]

Identification of genetic markers that affect economically important traits is of high value from a biological point of view, enabling the targeting of candidate genes and providing practical benefits for the industry such as wide-scale genomic selection. This study is one of the first to investigate the genetic background of economically important traits in dairy goats using the caprine 50K single nucleotide polymorphism (SNP) chip. The aim of the project was to perform a genomewide association study for milk yield and conformation of udder, teat, and feet and legs. A total of 137,235 milk yield records on 4,563 goats each scored for 10 conformation traits were available. Out of these, 2,381 goats were genotyped with the Illumina Caprine 50K BeadChip (Illumina Inc., San Diego, CA). A range of pseudo-phenotypes were used including deregressed breeding values and pseudo-estimated breeding values. Genome-wide association studies were performed using the multi-locus mixed model (MLMM) algorithm implemented in SNP & Variation Suite v7.7.8 (Golden Helix Inc., Bozeman, MT). A genome-wise significant [-log10(P-value) > 5.95] SNP for milk yield was identified on chromosome 19, with additional chromosomewise significant (-log10(P-value) > 4.46] SNP on chromosomes 4, 8, 14, and 29. Three genome-wise significant SNP for conformation of udder attachment, udder depth, and front legs were identified on chromosome 19, and chromosome-wise SNP were found on chromosomes 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 21, 23, and 27. The proportion of variance explained by the significant SNP was between 0.4 and 7.0% for milk yield and between 0.1 and 13.8% for conformation traits. This study is the first attempt to identify SNP associated with milk yield and conformation in dairy goats. Two genome-wise significant SNP for milk yield and 3 SNP for conformation of udder attachment, udder depth, and front legs were found. Our results suggest that conformation traits have a polygenic background because, for most of them, we did not identify any quantitative trait loci with major effect. [ABSTRACT FROM AUTHOR]

Individual wavenumbers of the infrared (IR) spectra of bovine milk have been shown to be moderately to highly heritable. The objective of this study was to identify genomic regions associated with individual milk IR wavenumbers. This is expected to provide information about the genetic background of milk composition and give insight in the relation between IR wavenumbers and milk components. For this purpose, a genome-wide association study was performed for a selected set of 50 individual IR wavenumbers measured on 1,748 Dutch Holstein cows. Significant associations were detected for 28 of the 50 wavenumbers. In total, 24 genomic regions distributed over 16 bovine chromosomes were identified. Major genomic regions associated with milk IR wavenumbers were identified on chromosomes 1, 5, 6, 14, 19, and 20. Most of these regions also showed significant associations with fat, protein, or lactose percentage. However, we also identified some new regions that were not associated with any one of these routinely collected milk composition traits. On chromosome 1, we identified 2 new genomic regions and hypothesized that they are related to variation in milk phosphorus content and orotic acid, respectively. On chromosome 20, we identified a new genomic region that seems to be related to citric acid. Identification of genomic regions associated with milk phosphorus content, orotic acid, and citric acid suggest that the milk IR spectra contain direct information on these milk components. Consequently milk IR analyses probably can be used to predict these milk components, which have low concentrations in milk; this can lead to novel applications of milk IR spectroscopy for dairy cattle breeding and herd management. [ABSTRACT FROM AUTHOR]

Genome-wide association testing facilitates the identification of genetic variants associated with complex traits. Mapping genes that promote genetic resistance to mastitis could reduce the cost of antibiotic use and enhance animal welfare and milk production by improving outcomes of breeding for udder health. Using imputed whole-genome sequence variants, we carried out association studies for 2 traits related to udder health, udder index, and milking speed in Nordic Holstein cattle. A total of 4,921 bulls genotyped with the BovineSNP50 BeadChip array were imputed to high-density genotypes (Illumina BovineHD Bead- Chip, Illumina, San Diego, CA) and, subsequently, to whole-genome sequence variants. An association analysis was carried out using a linear mixed model. Phenotypes used in the association analyses were deregressed breeding values. Multitrait meta-analysis was carried out for these 2 traits. We identified 10 and 8 chromosomes harboring markers that were significantly associated with udder index and milking speed, respectively. Strongest association signals were observed on chromosome 20 for udder index and chromosome 19 for milking speed. Multitrait meta-analysis identified 13 chromosomes harboring associated markers for the combination of udder index and milking speed. The associated region on chromosome 20 overlapped with earlier reported quantitative trait loci for similar traits in other cattle populations. Moreover, this region was located close to the FYB gene, which is involved in platelet activation and controls IL-2 expression; FYB is a strong candidate gene for udder health and worthy of further investigation. [ABSTRACT FROM AUTHOR]

The most common approach in dairy cattle to prove genotype by environment interactions is a multiple-trait model application, and considering the same traits in different environments as different traits. We enhanced such concepts by defining continuous phenotypic, genetic, and genomic herd descriptors, and applying random regression sire models. Traits of interest were test-day traits for milk yield, fat percentage, protein percentage, and somatic cell score, considering 267,393 records from 32,707 first-lactation Holstein cows. Cows were born in the years 2010 to 2013, and kept in 52 large-scale herds from 2 federal states of north-east Germany. The average number of genotyped cows per herd (45,613 single nucleotide polymorphism markers per cow) was 133.5 (range: 45 to 415 genotyped cows). Genomic herd descriptors were (1) the level of linkage disequilibrium (r²) within specific chromosome segments, and (2) the average allele frequency for single nucleotide polymorphisms in close distance to a functional mutation. Genetic herd descriptors were the (1) intra-herd inbreeding coefficient, and (2) the percentage of daughters from foreign sires. Phenotypic herd descriptors were (1) herd size, and (2) the herd mean for nonreturn rate. Most correlations among herd descriptors were close to 0, indicating independence of genomic, genetic, and phenotypic characteristics. Heritabilities for milk yield increased with increasing intraherd linkage disequilibrium, inbreeding, and herd size. Genetic correlations in same traits between adjacent levels of herd descriptors were close to 1, but declined for descriptor levels in greater distance. Genetic correlation declines were more obvious for somatic cell score, compared with test-day traits with larger heritabilities (fat percentage and protein percentage). Also, for milk yield, alterations of herd descriptor levels had an obvious effect on heritabilities and genetic correlations. By trend, multiple trait model results (based on created discrete herd classes) confirmed the random regression estimates. Identified alterations of breeding values in dependency of herd descriptors suggest utilization of specific sires for specific herd structures, offering new possibilities to improve sire selection strategies. Regarding genomic selection designs and genetic gain transfer into commercial herds, cow herds for the utilization in cow training sets should reflect the genomic, genetic, and phenotypic pattern of the broad population. [ABSTRACT FROM AUTHOR]

Equine-like G3P[8] rotavirus A strains with DS-1-like backbone genes have emerged since 2013. An equine-like RVA/Human-wt/JPN/15R429/2015/G3P[8] strain possessing I2-R2-C2-M2-A2-N2-T2-E2-H2 was detected in Japan in 2015. Its VP7 gene was ≥ 99.3% identical to those of equine-like G3P[4] strains detected in Japan, and the remaining 10 genes were 98.6-99.8% identical to G1P[8] double-gene reassortants detected in Japan, Thailand and the Philippines. Thus, 15R429 was likely generated through reassortment between the equine-like G3P[4] and G1P[8] reassortant strains. Notably, 15R429 was 98.5-99.8% identical across all 11 genes of the equine-like G3P[8] strains detected in Spain and Hungary in 2015. [ABSTRACT FROM AUTHOR]

Using combined genetic mapping, Illumina sequencing, bioinformatics analyses, and experimental validation, we identified 60 essential genes from 104 lethal mutations in two genomic regions of Caenorhabditis elegans totaling ~ 14 Mb on chromosome III(mid) and chromosome V(left). Five of the 60 genes had not previously been shown to have lethal phenotypes by RNA interference depletion. By analyzing the regions around the lethal missense mutations, we identified four putative new protein functional domains. Furthermore, functional characterization of the identified essential genes shows that most are enzymes, including helicases, tRNA synthetases, and kinases in addition to ribosomal proteins. Gene Ontology analysis indicated that essential genes often encode for enzymes that conduct nucleic acid binding activities during fundamental processes, such as intracellular DNA replication, transcription, and translation. Analysis of essential gene shows that they have fewer paralogs, encode proteins that are in protein interaction hubs, and are highly expressed relative to nonessential genes. All these essential gene traits in C. elegans are consistent with those of human disease genes. Most human orthologs (90%) of the essential genes in this study are related to human diseases. Therefore, functional characterization of essential genes underlines their importance as proxies for understanding the biological functions of human disease genes. [ABSTRACT FROM AUTHOR]

We investigated the genetic relatedness, antibiotic resistance and biofilm-producing ability of 114 strains of Salmonella , belonged to three serotypes (Corvallis, Brancaster and Albany), isolated from naturally contaminated poultry and their environment in wet markets and smale-scale processing plant from northern Malaysia. Pulsed-field gel electrophoresis revealed that Salmonella strains isolated from various wet markets were clonally related, suggesting the widespread dissemination of these three serotypes in northern Malaysia. All except one strain of Salmonella were resistant to more than two classes of antibiotics, hence regarded as multidrug resistant (MDR). Resistance to sulphonamide (96.5%), ampicillin (89.5%), tetracycline (85.1%), chloramphenicol (75.4%), trimethoprim (68.4%), trimethoprim-sulfamethoxazole (67.5%), streptomycin (58.8%) and nalidixic acid (44.4%) were observed. Resistance determinants, floR , cmlA , tetA , tetB , tetG , temB , bla PSE-1 , sul1 , sul2 , qnrA , qnrS , strA and aadA were detected by PCR among MDR Salmonella strains. Seventy-six strains (66.7%) harboured class-I integrons. The gene cassettes identified were dfrA1 , dfrA12 , aadA2 and an open reading frame orfC with unknown function. All Salmonella strains produced biofilm and 69.3% of them were strong biofilm-producers . Our findings suggested that most likely, persistent Salmonella colonises various sites in the processing environment by producing biofilm, which leads to their widespread dissemination in wet markets located in northern Malaysia. [ABSTRACT FROM AUTHOR]

Background: Milu, also known as P`ere David's deer (Elaphurus davidianus), was widely distributed in East Asia but recently experienced a severe bottleneck. Only 18 survived by the end of the 19th century, and the current population of 4500 individuals was propagated from just 11 kept by the 11th British Duke of Bedford. This species is known for its distinguishable appearance, the driving force behind which is still a mystery. To aid efforts to explore these phenomena, we constructed a draft genome of the species. Findings: In total, we generated 321.86 gigabases (Gb) of raw DNA sequence from whole-genome sequencing of a male milu deer using an Illumina HiSeq 2000 platform. Assembly yielded a final genome with a scaffold N50 size of 3.03 megabases (Mb) and a total length of 2.52 Gb. Moreover, we identified 20 125 protein-coding genes and 988.1 Mb of repetitive sequences. In addition, homology-based searches detected 280 rRNA, 1335 miRNA, 1441 snRNA, and 893 tRNA sequences in the milu genome. The divergence time between E. davidianus and Bos taurus was estimated to be about 28.20 million years ago (Mya). We identified 167 species-specific genes and 293 expanded gene families in the milu lineage. Conclusions: We report the first reference genome of milu, which will provide a valuable resource for studying the species' demographic history of severe bottleneck and the genetic mechanism(s) of special phenotypic evolution. [ABSTRACT FROM AUTHOR]

We present evidence that the dog hookworm (Ancylostoma caninum) is underutilised in the study of host-parasite interactions, particularly as a proxy for the human-hookworm relationship. The inability to passage hookworms through all life stages in vitro means that adult stage hookworms have to be harvested from the gut of their definitive hosts for ex vivo research. This makes study of the human-hookworm interface difficult for technical and ethical reasons. The historical association of humans, dogs and hookworms presents a unique triad of positive evolutionary pressure to drive the A. caninum-canine interaction to reflect that of the human-hookworm relationship. Here we discuss A. caninum as a proxy for human hookworm infection and situate this hookworm model within the current research agenda, including the various 'omics' applications and the search for next generation biologics to treat a plethora of human diseases. Historically, the dog hookworm has been well described on a physiological and biochemical level, with an increasing understanding of its role as a human zoonosis. With its similarity to human hookworm, the recent publications of hookworm genomes and other omics databases, as well as the ready availability of these parasites for ex vivo culture, the dog hookworm presents itself as a valuable tool for discovery and translational research. [ABSTRACT FROM AUTHOR]

Haploid cells are increasingly used for screening of complex pathways in animal genomes. Hemizygous mutations introduced through viral insertional mutagenesis can be directly selected for phenotypic changes. Here we present HaSAPPy a tool for analysing sequencing datasets of screens using insertional mutations in large pools of haploid cells. Candidate gene prediction is implemented through identification of enrichment of insertional mutations after selection by simultaneously evaluating several parameters. We have developed HaSAPPy for analysis of genetic screens for silencing factors of X chromosome inactivation in haploid mouse embryonic stem cells. To benchmark the performance, we further analyse several datasets of genetic screens in human haploid cells for which candidates have been validated. Our results support the effective candidate prediction strategy of HaSAPPy. HaSAPPy is implemented in Python, licensed under the MIT license, and is available from . [ABSTRACT FROM AUTHOR]

In the silkworm Bombyx mori, more than ten mutations exhibit an embryonic lethal phenotype. The lethal 19 (1-19) mutation is one of these, and embryos exhibiting this mutation die at blastokinesis. The 1-19 mutation is controlled by a single recessive gene mapped to chromosome 19. However, the responsible gene has not yet been identified. In this study, we performed morphological observations of the eggs and embryos using microscopy and mapped 1-19. Observations revealed that eggs with the 1-19 mutation die before or after the body pigmentation stage, but only embryos in eggs before the body pigmentation stage were dead at blastokinesis. This result suggests that the phenotype of the 1-19 mutation can be identified by sampling dead eggs at this stage. We next performed mapping of 1-19 using primer sets designed for chromosome 19. However, we could not identify the 1-19-linked region, because there were no regions where the genotypes of all F2 individuals were homozygous for the 1-19 mutant, suggesting that 1-19 does not map to chromosome 19. Therefore, we performed a linkage analysis using primer sets designed for each chromosome. Our analysis showed that 1-19 does not map to chromosome 19, but instead maps to chromosome 12. Further mapping indicated that the 1-19-linked region spanned approximately 10 Mb on chromosome 12. [ABSTRACT FROM AUTHOR]

Identifying the genetic bases for colour patterns has provided important insights into the control and expression of pigmentation and howthese characteristics influence fitness. However, much more is known about the genetic bases for traits based on melanin pigments than for traits based on another major class of pigments, carotenoids. Here, we use natural admixture in a hybrid zone between Audubon's and myrtle warblers (Setophaga coronata auduboni/S. c. coronata) to identify genomic regions associated with both types of pigmentation. Warblers are known for rapid speciation and dramatic differences in plumage. For each of five plumage coloration traits, we found highly significant associations with multiple single-nucleotide polymorphisms (SNPs) across the genome and these were clustered in discrete regions. Regions near significantly associated SNPs were enriched for genes associated with keratin filaments, fibrils that make up feathers. A carotenoid-based trait that differs between the taxa--throat colour--had more than a dozen genomic regions of association. One cluster of SNPs for this trait overlaps the Scavenger Receptor Class F Member 2 (SCARF2) gene. Other scavenger receptors are presumed to be expressed at target tissues and involved in the selective movement of carotenoids into the target cells, making SCARF2 a plausible new candidate for carotenoid processing. In addition, two melanin-based plumage traits--colours of the eye line and eye spot--show very strong associations with a single genomic region mapping to chromosome 20 in the zebra finch. These findings indicate that only a subset of the genomic regions differentiated between these two warblers are associated with the plumage differences between them and demonstrate the utility of reduced-representation genomic scans in hybrid zones. [ABSTRACT FROM AUTHOR]

Gene copy number variation ( CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays ( aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution. [ABSTRACT FROM AUTHOR]

Evolutionary ecologists aim to explain and predict evolutionary change under different selective regimes. Theory suggests that such evolutionary prediction should be more difficult for biomechanical systems in which different trait combinations generate the same functional output: 'many-to-one mapping.' Many-to-one mapping of phenotype to function enables multiple morphological solutions to meet the same adaptive challenges. Therefore, many-to-one mapping should undermine parallel morphological evolution, and hence evolutionary predictability, even when selection pressures are shared among populations. Studying 16 replicate pairs of lake- and stream-adapted threespine stickleback ( Gasterosteus aculeatus), we quantified three parts of the teleost feeding apparatus and used biomechanical models to calculate their expected functional outputs. The three feeding structures differed in their form-to-function relationship from one-to-one (lower jaw lever ratio) to increasingly many-to-one (buccal suction index, opercular 4-bar linkage). We tested for (1) weaker linear correlations between phenotype and calculated function, and (2) less parallel evolution across lake-stream pairs, in the many-to-one systems relative to the one-to-one system. We confirm both predictions, thus supporting the theoretical expectation that increasing many-to-one mapping undermines parallel evolution. Therefore, sole consideration of morphological variation within and among populations might not serve as a proxy for functional variation when multiple adaptive trait combinations exist. [ABSTRACT FROM AUTHOR]

MicroRNA (miRNA)-mediated gene silencing is a key mechanism regulating numerous biological processes such as development of organs and tumorigenesis. The expression of miRNA machinery genes linked to miRNA biogenesis and processing is finely regulated. Despite accumulating evidence for chicken miRNA in the female reproduction system, precise regulatory mechanisms are largely unknown. Therefore, the objective of this study was to determine changes in expression levels of miRNA machinery genes in developmental stages of the oviduct and ovarian carcinogenesis of laying hens. In the present study, differential expression of miRNA machinery genes during ovarian carcinogenesis was determined using cancerous and normal ovaries collected from normal laying hens and hens with cancer. Our results showed that 3 miRNA machinery genes (Dicer1, AGO3, and AGO4) were differentially expressed as laying hens' reproductive organs developed. These genes were simultaneously upregulated in cancerous ovaries compared with those in normal ovaries. Their transcripts were abundantly localized in glandular epithelial cells of cancerous ovaries. Our results indicate that Dicer1, AGO3, and AGO4 play critical roles in the development of reproductive organs and ovarian carcinogenesis in laying hens, suggesting that simultaneous overexpression of these genes might serve as a prognostic factor for ovarian cancer. [ABSTRACT FROM AUTHOR]

The biogeographical history of Southeast Asia is complicated due to the continuous emergences and disappearances of land bridges throughout the Pleistocene. Here, we use long-tailed macaques ( Macaca fascicularis ), which are widely distributed throughout the mainland and islands of Southeast Asia, as a model for better understanding the biogeographical patterns of diversification in this geographically complex region. A reliable intraspecific phylogeny including individuals from localities on oceanic islands, continental islands, and the mainland is needed to trace relatedness along with the pattern and timing of colonization in this region. We used high-throughput sequencing techniques to sequence mitochondrial genomes (mitogenomes) from 95 Southeast Asian M. fascicularis specimens housed at natural history museums around the world. To achieve a comprehensive picture, we more than tripled the mitogenome sample size for M. fascicularis from previous studies, and for the first time included documented samples from the Philippines and several small Indonesian islands. Confirming the result from a previous, recent intraspecific phylogeny for M. fascicularis , the newly reconstructed phylogeny of 135 specimens divides the samples into two major clades: Clade A includes haplotypes from the mainland and some from northern Sumatra, while Clade B includes all insular haplotypes along with lineages from southern Sumatra. This study resolves a previous disparity by revealing a disjunction in the origin of Sumatran macaques, with separate lineages originating within the two major clades, suggesting that at least two major migrations to Sumatra occurred. However, our dated phylogeny reveals that the two major clades split ∼1.88 Ma, which is earlier than in previously published phylogenies. Our new data reveal that most Philippine macaque lineages diverged from the Borneo stock within the last ∼0.06–0.43 Ma. Finally, our study provides insight into successful sequencing of DNA across museums and shotgun sequencing of DNA specimens as a method to sequence the mitogenome. [ABSTRACT FROM AUTHOR]

The chimpanzee is arguably the most important species for the study of human origins. A key resource for these studies is a high-quality reference genome assembly; however, as with most mammalian genomes, the current iteration of the chimpanzee reference genome assembly is highly fragmented. In the current iteration of the chimpanzee reference genome assembly (Pan tro 2.1.4), the sequence is scattered across more then 183 000 contigs, incorporating more than 159 000 gaps, with a genome-wide contig N50 of 51 Kbp. In this work, we produce an extensive and diverse array of sequencing datasets to rapidly assemble a new chimpanzee reference that surpasses previous iterations in bases represented and organized in large scaffolds. To this end, we show substantial improvements over the current release of the chimpanzee genome (Pan tro 2.1.4) by several metrics, such as increased contiguity by >750% and 300% on contigs and scaffolds, respectively, and closure of 77% of gaps in the Pan tro 2.1.4 assembly gaps spanning >850 Kbp of the novel coding sequence based on RNASeq data. We further report more than 2700 genes that had putatively erroneous frame-shift predictions to human in Pan tro 2.1.4 and show a substantial increase in the annotation of repetitive elements. We apply a simple 3-way hybrid approach to considerably improve the reference genome assembly for the chimpanzee, providing a valuable resource for the study of human origins. Furthermore, we produce extensive sequencing datasets that are all derived from the same cell line, generating a broad non-human benchmark dataset. [ABSTRACT FROM AUTHOR]