Y‐chromosomal variation confirms independent domestications of swamp and river buffalo

Y‐chromosomal variation in the water buffalo was analysed by sequencing of DBY, ZFY and SRY gene segments. A clear separation of the paternal lineages of the river and swamp types parallels the differences between their maternal lineages and nuclear DNA. Sequence divergence was found to be comparable to the divergence of taurine cattle and zebu, and this divergence predated domestication, confirming that river and swamp buffalo originated from different wild populations. Within a sample of 23 Thai swamp buffaloes, we identified four haplotypes with different geographical distributions, two of which were shared by Thai wild buffaloes.     

Mitochondrial DNA and Y‐chromosome diversity in East Adriatic sheep

Variation in mitochondrial DNA (mtDNA) and Y‐chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon ( Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt‐3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt‐3 respectively. In contrast, Y‐chromosome analyses revealed extraordinary paternal allelic richness: HY‐6, 89.3%; HY‐8, 5.0%; HY‐18, 3.1%; HY‐7, 1.3%; and HY‐5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY‐18 (A‐oY1/135‐SRYM18), identified here for the first time, provides a link between the haplotype HY‐12 (A‐oY1/139‐SRYM18) found in a few rams in Turkey and haplotype HY‐9 (A‐oY1/131‐SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt‐55), and all were paternally monomorphic for haplotype HY‐6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y‐chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to ‘upgrade’ local populations using rams with different paternal origins.     

Y‐chromosome evidence supports widespread signatures of three‐species Canis hybridization in eastern North America

There has been considerable discussion on the origin of the red wolf and eastern wolf and their evolution independent of the gray wolf. We analyzed mitochondrial DNA (mtDNA) and a Y‐chromosome intron sequence in combination with Y‐chromosome microsatellites from wolves and coyotes within the range of extensive wolf–coyote hybridization, that is, eastern North America. The detection of divergent Y‐chromosome haplotypes in the historic range of the eastern wolf is concordant with earlier mtDNA findings, and the absence of these haplotypes in western coyotes supports the existence of the North American evolved eastern wolf (Canis lycaon). Having haplotypes observed exclusively in eastern North America as a result of insufficient sampling in the historic range of the coyote or that these lineages subsequently went extinct in western geographies is unlikely given that eastern‐specific mtDNA and Y‐chromosome haplotypes represent lineages divergent from those observed in extant western coyotes. By combining Y‐chromosome and mtDNA distributional patterns, we identified hybrid genomes of eastern wolf, coyote, gray wolf, and potentially dog origin in Canis populations of central and eastern North America. The natural contemporary eastern Canis populations represent an important example of widespread introgression resulting in hybrid genomes across the original C. lycaon range that appears to be facilitated by the eastern wolf acting as a conduit for hybridization. Applying conventional taxonomic nomenclature and species‐based conservation initiatives, particularly in human‐modified landscapes, may be counterproductive to the effective management of these hybrids and fails to consider their evolutionary potential.     

Novel Y‐chromosome polymorphisms in Chinese domestic yak

Y‐chromosome‐specific haplotypes (Y‐haplotypes) constructed using single nucleotide polymorphisms (Y‐SNPs) in the MSY (male‐specific region of the Y‐chromosome) are valuable in population genetic studies. But sequence variants in the yak MSY region have been poorly characterized so far. In this study, we screened a total of 16 Y‐chromosome‐specific gene segments from the ZFY, SRY, UTY, USP9Y, AMELY and OFD1Y genes to identify Y‐SNPs in domestic yaks. Six novel Y‐SNPs distributed in the USP9Y (g.223C>T), UTY19 (g.158A>C and g.169C>T), AMELY2 (g.261C>T), OFD1Y9 (g.165A>G) and SRY4 (g.104G>A) loci, which can define three Y‐haplotypes (YH1, YH2 and YH3) in yaks, were discovered. YH1 was the dominant and presumably most ancient haplotype based on the comparison of UTY19 locus with other bovid species. Interestingly, we found informative UTY19 markers (g.158A>C and g.169C>T) that can effectively distinguish the three yak Y‐haplotypes. The nucleotide diversity was 1.7 × 10^?4 ± 0.3 × 10^?4, indicating rich Y‐chromosome diversity in yaks. We identified two highly divergent lineages (YH1 and YH2 vs. YH3) that share similar frequencies (YH1 +  YH2: 0.82–0.89, YH3: 0.11–0.18) among all three populations. In agreement with previous mtDNA studies, we supported the hypothesis that the two highly divergent lineages (YH1 and YH2 vs. YH3) derived from a single gene pool, which can be explained by the reunion of at least two paternal populations with the divergent lineages already accumulated before domestication. We estimated a divergence time of 408 110 years between the two divergent lineages, which is consistent with the data from mitochondrial DNA in yaks.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Genome‐wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome‐wide sequence information, independent data were obtained from genotyping‐by‐sequencing and a target‐enrichment experiment that returned 244 low‐copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent‐based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat‐group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat‐group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.     

Biogeography,habitat transitions and hybridization in a radiation of South American silverside fishes revealed by mitochondrial and genomic RAD data

Rivers and lake systems in the southern cone of South America have been widely influenced by historical glaciations, carrying important implications for the evolution of aquatic organisms, including prompting transitions between marine and freshwater habitats and by triggering hybridization among incipient species via waterway connectivity and stream capture events. Silverside fishes (Odontesthes) in the region comprise a radiation of 19 marine and freshwater species that have been hypothesized on the basis of morphological or mitochondrial DNA data to have either transitioned repeatedly into continental waters from the sea or colonized marine habitats following freshwater diversification. New double digest restriction‐site associated DNA data presented here provide a robust framework to investigate the biogeographical history of and habitat transitions in Odontesthes. We show that Odontesthes silversides originally diversified in the Pacific but independently colonized the Atlantic three times, producing three independent marine‐to‐freshwater transitions. Our results also indicate recent introgression of marine mitochondrial haplotypes into two freshwater clades, with more recurring instances of hybridization among Atlantic‐ versus Pacific‐slope species. In Pacific freshwater drainages, hybridization with a marine species appears to be geographically isolated and may be related to glaciation events. Substantial structural differences of estuarine gradients between these two geographical areas may have influenced the frequency, intensity and evolutionary effects of hybridization events.     

Extensive hybridization and past introgression between divergent lineages in a quasi‐clonal hermaphroditic fish: Ramifications for species concepts and taxonomy

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self‐fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a ‘Central clade’ closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self‐fertilization for 1–3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.     

Genome‐wide association study confirms that the chromosome Z harbours a region responsible for rumplessness in Hongshan chickens

Rumplessness in Hongshan chickens has been reported as a novel sex‐linked characteristic. Re‐sequencing data suggest that a pseudogene on the Z chromosome, LOC431648, is affiliated with this phenotype. In this study, we chose 23 rumpless and 25 normal Hongshan chickens to localize the potential variation by means of a genome‐wide association study using a high density microarray. A region on the Z chromosome was found to be closely associated with rumplessness in Hongshan chickens. The region, located in gene LINGO2, was approximately 3 Mb away from pseudogene LOC431648. The function of this gene has not yet been studied in birds. Differential expression of the candidate genes in the tail feather follicles was not detected by q‐PCR, which suggests that the rumplessness trait could be attributed to other genetic mechanisms.     

Genotype × environment interaction analysis of North American shrub willow yield trials confirms superior performance of triploid hybrids

Development of dedicated bioenergy crop production systems will require accurate yield estimates, which will be important for determining many of the associated environmental and economic impacts of their production. Shrub willow (Salix spp) is being promoted in areas of the USA and Canada due to its adaption to cool climates and wide genetic diversity available for breeding improvement. Willow breeding in North America is in an early stage, and selection of elite genotypes for commercialization will require testing across broad geographic regions to gain an understanding of how shrub willow interacts with the environment. We analyzed a dataset of first‐rotation shrub willow yields of 16 genotypes across 10 trial environments in the USA and Canada for genotype‐by‐environment interactions using the additive main effects and multiplicative interactions (AMMI) model. Mean genotype yields ranged from 5.22 to 8.58 oven‐dry Mg ha^?1 yr^?1. Analysis of the main effect of genotype showed that one round of breeding improved yields by as much as 20% over check cultivars and that triploid hybrids, most notably Salix viminalis × S. miyabeana, exhibited superior yields. We also found important variability in genotypic response to environments, which suggests specific adaptability could be exploited among 16 genotypes for yield gains. Strong positive correlations were found between environment main effects and AMMI parameters and growing environment temperatures. These findings demonstrate yield improvements are possible in one generation and will be important for developing cultivar recommendations and for future breeding efforts.     

Tumour cell‐intrinsic CTLA4 regulates PD‐L1 expression in non‐small cell lung cancer

Cytotoxic T lymphocyte antigen 4 (CTLA4) and programmed cell death protein 1 (PD‐1) are immune checkpoint proteins expressed in T cells. Although CTLA4 expression was found in multiple tumours including non‐small cell lung cancer (NSCLC) tissues and cells, its function in tumour cells is unknown. Recently, PD‐1 was found to be expressed in melanoma cells and to promote tumorigenesis. We found that CTLA4 was expressed in a subset of NSCLC cell lines and in a subgroup of cancer cells within the lung cancer tissues. We further found that in NSCLC cells, anti‐CTLA4 antibody can induce PD‐L1 expression, which is mediated by CTLA4 and the EGFR pathway involving phosphorylation of MEK and ERK. In CTLA4 knockout cells, EGFR knockout cells or in the presence of an EGFR tyrosine kinase inhibitor, anti‐CTLA4 antibody was not able to induce PD‐L1 expression in NSCLC cells. Moreover, anti‐CTLA4 antibody promoted NSCLC cell proliferation in vitro and tumour growth in vivo in the absence of adaptive immunity. These results suggest that tumour cell‐intrinsic CTLA4 can regulate PD‐L1 expression and cell proliferation, and that anti‐CTLA4 antibody, by binding to the tumour cell‐intrinsic CTLA4, may result in the activation of the EGFR pathway in cancer cells.     

ICAM‐1 and MKL‐1 polymorphisms impose considerable impacts on coronary heart disease occurrence

This study was aimed to explore the correlation of intercellular adhesion molecule‐1 (ICAM‐1) K469E and megakaryoblastic leukaemia factor‐1 (MKL‐1) ?184C/T polymorphisms with the susceptibility to coronary heart disease (CHD) in the Chinese Han population. 100 CHD patients and 91 healthy people that had no blood connection with each other were enrolled in this case‐control study. ICAM‐1 and MKL‐1 polymorphisms were genotyped by polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) approach. Multiple logistic regression was used to analyse the correlation between polymorphisms of ICAM‐1 and MKL‐1 and CHD susceptibility. Differences of genotype and allele frequencies of the two SNPs between case and control groups were analysed by chi‐square test. Odds ratios (ORs) and 95% confidence intervals (CIs) were indicated relative susceptibility of CHD. The distributions of ICAM‐1 and MKL‐1 polymorphisms in each group conformed to Hardy‐Weinberg equilibrium (HWE). After adjusting for traditional risk factors, the TT genotype frequency of MKL‐1 ?184C/T polymorphism was found significantly higher in case group than in control group (P < .05). Meanwhile, T allele frequency increased in case group compared with control group, and the differences had statistical significance (P = .04, OR = 2.34, 95% CI = 1.34‐5.26). Logistic regression analysis in this study proved that smoking, hypertension, diabetes and triglyceride (TG) were all risk factors for CHD ICAM‐1 K469E polymorphism has no association with the onset of CHD. But MKL‐1 ?184C/T polymorphism is associated with the risk of CHD and T allele might be a susceptibility factor for CHD.     

CD1d‐mediated activation of group 3 innate lymphoid cells drives IL‐22 production

Innate lymphoid cells (ILCs) are a heterogeneous family of immune cells that play a critical role in a variety of immune processes including host defence against infection, wound healing and tissue repair. Whether these cells are involved in lipid‐dependent immunity remains unexplored. Here we show that murine ILCs from a variety of tissues express the lipid‐presenting molecule CD1d, with group 3 ILCs (ILC3s) showing the highest level of expression. Within the ILC3 family, natural cytotoxicity triggering receptor (NCR)^?CCR6⁺ cells displayed the highest levels of CD1d. Expression of CD1d on ILCs is functionally relevant as ILC3s can acquire lipids in vitro and in vivo and load lipids on CD1d to mediate presentation to the T‐cell receptor of invariant natural killer T (iNKT) cells. Conversely, engagement of CD1d in vitro and administration of lipid antigen in vivo induce ILC3 activation and production of IL‐22. Taken together, our data expose a previously unappreciated role for ILCs in CD1d‐mediated immunity, which can modulate tissue homeostasis and inflammatory responses.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.     

A single‐nucleotide polymorphism of miR‐196a‐2 and vitiligo: an association study and functional analysis in a Han Chinese population

Recent evidence indicates that oxidative stress and genetic factors play an important role in the pathogenesis of vitiligo. SNPs in miRNAs involved in oxidative stress could potentially influence the development of vitiligo. In this case–control study, we investigated the association of a functional SNP of rs11614913 in miR‐196a‐2 with risk of vitiligo. A significantly lower risk of vitiligo was associated with the rs11614913 miR‐196a‐2 CC genotype (adjusted OR, 0.77; CI, 0.60–0.98). In addition, TYRP1 gene expression was considerably down‐regulated by the rs11614913 C allele in miR‐196a‐2, which lowered the levels of intracellular reactive oxygen species (ROS) and reduced the proportion of early apoptosis in human melanocytes in response to H₂O₂ treatment. Our data suggest that the rs11614913 C allele in miR‐196a‐2 confers potential protection against oxidative effects on human melanocytes through the modulation of the target gene, TYRP1, which may account for the decreased risk of vitiligo in this study population.     

Valsartan independent of AT1 receptor inhibits tissue factor,TLR‐2 and ‐4 expression by regulation of Egr‐1 through activation of AMPK in diabetic conditions

Patients suffering from diabetes mellitus (DM) are at a severe risk of atherothrombosis. Early growth response (Egr)‐1 is well characterized as a central mediator in vascular pathophysiology. We tested whether valsartan independent of Ang II type 1 receptor (AT₁R) can reduce tissue factor (TF) and toll‐like receptor (TLR)‐2 and ‐4 by regulating Egr‐1 in THP‐1 cells and aorta in streptozotocin‐induced diabetic mice. High glucose (HG, 15 mM) increased expressions of Egr‐1, TF, TLR‐2 and ‐4 which were significantly reduced by valsartan. HG increased Egr‐1 expression by activation of PKC and ERK1/2 in THP‐1 cells. Valsartan increased AMPK phosphorylation in a concentration and time‐dependent manner via activation of LKB1. Valsartan inhibited Egr‐1 without activation of PKC or ERK1/2. The reduced expression of Egr‐1 by valsartan was reversed by either silencing Egr‐1, or compound C, or DN‐AMPK‐transfected cells. Valsartan inhibited binding of NF‐κB and Egr‐1 to TF promoter in HG condition. Furthermore, valsartan reduced inflammatory cytokine (TNF‐α, IL‐6 and IL‐1β) production and NF‐κB activity in HG‐activated THP‐1 cells. Interestingly, these effects of valsartan were not affected by either silencing AT₁R in THP‐1 cells or CHO cells, which were devoid of AT₁R. Importantly, administration of valsartan (20 mg/kg, i.p) for 8 weeks significantly reduced plasma TF activity, expression of Egr‐1, TLR‐2, ‐4 and TF in thoracic aorta and improved glucose tolerance of streptozotocin‐induced diabetic mice. Taken together, we concluded that valsartan may reduce atherothrombosis in diabetic conditions through AMPK/Egr‐1 regulation.     

Analysing variation in Drosophila aging across independent experimental studies: a meta‐analysis of survival data

Survival records of longevity experiments are a key component in research on aging. However, surprisingly there have been very few cross‐study analyses, besides comparisons of median lifespans or similar summary information. Here, we use a large set of full survival data from various studies to address questions in aging, which are beyond the scope of individual studies. We characterize survival differences between female and male flies of different genetic Drosophila strains, showing significant differences between strains. We further analyse the variation in survival of control cohorts recorded under highly similar conditions within different Drosophila strains. We found that overall transgenic constructs of the UAS/GAL4 expression system which should have no effect (e.g. a GAL4 construct alone) extend lifespan significantly in the w1118 strain. Using a large data set comprised of various studies, we found no evidence for larger lifespan extensions being associated with shorter lifespans of the control in Drosophila. This demonstrates that lifespan extending treatments are not purely rescuing weak backgrounds.     

lac‐1 and lag‐1 with ras‐1 affect aging and the biological clock in Neurospora crassa

Using an automated cell counting technique developed previously (Case et al., Ecology and Evolution 2014; 4: 3494), we explore the lifespan effects of lac‐1, a ceramide synthase gene paralogous to lag‐1 in Neurospora crassa in conjunction with the band bd (ras‐1) gene. We find that the replicative lifespan of a lac‐1^KO bd double mutants is short, about one race tube cycle, and this double mutant lacks a strong ~21‐hr clock cycle as shown by race tube and fluorometer analysis of fluorescent strains including lac‐1^KO. This short replicative lifespan phenotype is contrasted with a very long estimated chronological lifespan for lac‐1^KO bd double mutants from 247 to 462 days based on our regression analyses on log viability, and for the single mutant lac‐1^KO, 161 days. Both of these estimated lifespans are much higher than that of previously studied WT and bd single mutant strains. In a lac‐1 rescue and induction experiment, the expression of lac‐1⁺ as driven by a quinic acid‐dependent promoter actually decreases the median chronological lifespan of cells down to only 7 days, much lower than the 34‐day median lifespan found in control bd conidia also grown on quinic acid media, which we interpret as an effect of balancing selection acting on ceramide levels based on previous findings from the literature. Prior work has shown phytoceramides can act as a signal for apoptosis in stressed N. crassa cells. To test this hypothesis of balancing selection on phytoceramide levels, we examine the viability of WT, lag‐1^KO bd, and lac‐1^KO bd strains following the dual stresses of heat and glycolysis inhibition, along with phytoceramide treatments of different dosages. We find that the phytoceramide dosage–response curve is altered in the lag‐1^KO bd mutant, but not in the lac‐1^KO bd mutant. We conclude that phytoceramide production is responsible for the previously reported longevity effects in the lag‐1^KO bd mutant, but a different ceramide may be responsible for the longevity effect observed in the lac‐1^KO bd mutant.